Unzipping and binding of small interfering RNA with single walled carbon nanotube: A platform for small interfering RNA delivery

In an effort to design efficient platform for siRNA delivery, we combine all atom classical and quantum simulations to study the binding of small interfering RNA (siRNA) by pristine single wall carbon nanotube (SWCNT). Our results show that siRNA strongly binds to SWCNT surface via unzipping its base-pairs and the propensity of unzipping increases with the increase in the diameter of the SWCNTs. The unzipping and subsequent wrapping events are initiated and driven by van der Waals interactions between the aromatic rings of siRNA nucleobases and the SWCNT surface. However, molecular dynamics (MD) simulations of double strand DNA (dsDNA) of the same sequence show that the dsDNA undergoes much less unzipping and wrapping on the SWCNT in the simulation time scale of 70 ns. This interesting difference is due to smaller interaction energy of thymidine of dsDNA with the SWCNT compared to that of uridine of siRNA, as calculated by dispersion corrected density functional theory (DFT) methods. After the optimal binding of siRNA to SWCNT, the complex is very stable which serves as one of the major mechanisms of siRNA delivery for biomedical applications. Since siRNA has to undergo unwinding process with the effect of RNA-induced silencing complex, our proposed delivery mechanism by SWCNT possesses potential advantages in achieving RNA interference. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3682780]

Theoretical analysis of photoinduced H-atom elimination in thiophenol

The photoinduced hydrogen elimination reaction in thiophenol via the conical intersections of the dissociative (1)pi sigma* excited state with the bound (1)pi pi* excited state and the electronic ground state has been investigated with ab initio electronic-structure calculations and time-dependent quantum wave-packet calculations. A screening of the coupling constants of the symmetry-allowed coupling modes at the (1)pi pi*-(1)pi sigma* and (1)pi sigma*-S-0 conical intersection shows that the SH torsional mode is by far the most important coupling mode at both conical intersections. A model including three intersecting potential-energy surfaces (S-0, (1)pi pi*, (1)pi sigma*) and two nuclear degrees of freedom (SH stretch and SH torsion) has been constructed on the basis of ab initio complete-active-space self-consistent field and multireference second-order perturbation theory calculations. The nonadiabatic quantum wave-packet dynamics initiated by optical excitation of the (1)pi pi* and (1)pi sigma* states has been explored for this three-state two-coordinate model. The photodissociation dynamics is characterized in terms of snapshots of time-dependent wave packets, time-dependent electronic population probabilities, and the branching ratio of the (2)sigma/(2)pi electronic states of the thiophenoxyl radical. The dependence of the timescale of the photodissociation process and the branching ratio on the initial excitation of the SH stretching and SH torsional vibrations has been analyzed. It is shown that the node structure, which is imposed on the nuclear wave packets by the initial vibrational preparation as well as by the transitions through the conical intersections, has a profound effect on the photodissociation dynamics. The effect of additional weak coupling modes of CC twist (nu(16a)) and ring-distortion (nu(16b)) character has been investigated with three-dimensional and four-dimensional time-dependent wave-packet calculations, and has been found to be minor. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4709608]

Relaxation and jump dynamics of water at the mica interface

The orientational relaxation dynamics of water confined between mica surfaces is investigated using molecular dynamics simulations. The study illustrates the wide heterogeneity that exists in the dynamics of water adjacent to a strongly hydrophilic surface such as mica. Analysis of the survival probabilities in different layers is carried out by normalizing the corresponding relaxation times with bulk water layers of similar thickness. A 10-fold increase in the survival times is observed for water directly in contact with the mica surface and a non-monotonic variation in the survival times is observed moving away from the mica surface to the bulk-like interior. The orientational relaxation time is highest for water in the contact layer, decreasing monotonically away from the surface. In all cases the ratio of the relaxation times of the 1st and 2nd rank Legendre polynomials of the HH bond vector is found to lie between 1.5 and 1.9 indicating that the reorientational relaxation in the different water layers is governed by jump dynamics. The orientational dynamics of water in the contact layer is particularly novel and is found to undergo distinct two-dimensional hydrogen bond jump reorientational dynamics with an average waiting time of 4.97 ps. The waiting time distribution is found to possess a long tail extending beyond 15 ps. Unlike previously observed jump dynamics in bulk water and other surfaces, jump events in the mica contact layer occur between hydrogen bonds formed by the water molecule and acceptor oxygens on the mica surface. Despite slowing down of the water orientational relaxation near the surface, life-times of water in the hydration shell of the K ion are comparable to that observed in bulk salt solutions. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4717710]

Crowding, molecular volume and plasticity: An assessment involving crystallography, NMR and simulations

The discrepancy between the X-ray and NMR structures of Mycobacterium tuberculosis peptidyl-tRNA hydrolase in relation to the functionally important plasticity of the molecule led to molecular dynamics simulations. The X-ray and the NMR studies along with the simulations indicated an inverse correlation between crowding and molecular volume. A detailed comparison of proteins for which X-ray and the NMR structures appears to confirm this correlation. In consonance with the reported results of the investigations in cellular compartments and aqueous solution, the comparison indicates that the crowding results in compaction of the molecule as well as change in its shape, which could specifically involve regions of the molecule important in function. Crowding could thus influence the action of proteins through modulation of the functionally important plasticity of the molecule. Selvaraj M, Ahmad R, Varshney U and Vijayan M 2012 Crowding, molecular volume and plasticity: An assessment involving crystallography, NMR and simulations. J. Biosci. 37 953-963] DOI 10.1007/s12038-012-9276-5

Simulation of Influence of Bilayer Melting on Dynamics and Thermodynamics of Interfacial Water

We investigate the effect of bilayer melting transition on thermodynamics and dynamics of interfacial water using molecular dynamics simulation with the two-phase thermodynamic model. We show that the diffusivity of interface water depicts a dynamic crossover at the chain melting transition following an Arrhenius behavior until the transition temperature. The corresponding change in the diffusion coefficient from the bulk to the interface water is comparable with experimental observations found recently for water near 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) vesicles Phys. Chem. Chem. Phys. 13, 7732 (2011)]. The entropy and potential energy of interfacial water show distinct changes at the bilayer melting transition, indicating a strong correlation in the thermodynamic state of water and the accompanying first-order phase transition of the bilayer membrane. DOI: 10.1103/PhysRevLett.110.018303

Charge-transfer-driven inclusion of neutral TCNQ molecules in the galleries of a layered double hydroxide: an experimental and computational study

The layered double hydroxides (LDH) or anionic clays are an important class of ion-exchange materials. They consist of positively charged brucite-like inorganic sheets with charge-compensating exchangeable anions in the interlamellar space. Here we show how neutral TCNQ (7,7,8,8-tetracyanoquinodimethane) molecules can be included within the galleries of an LDH. To do so, we exploit the fact that TCNQ is a good electron acceptor that forms donor acceptor complexes with a variety of donors. The electron donor aniline was intercalated into a Mg-Al LDH as p-aminobenzoate (AB) ions by a conventional ion-exchange reaction. We show here that neutral TCNQ molecules may be driven into the galleries of the layered solid by charge-transfer complex formation with the intercalated p-aminobenzoate anions. We use diffraction and spectroscopic measurements in combination with molecular dynamics simulations and quantum chemical calculations to establish the nature of interactions and arrangement of the charge-transfer complex within the galleries of the layered double hydroxide. Electrostatic interactions between the TCNQ molecules and the anchored AB ions, subsequent to charge transfer, are the driving force for the inclusion of TCNQ molecules in the galleries of the LDH.

Engineering Gold Nanoparticle Interaction by PAMAM Dendrimer

Bare faceted gold nanoparticles (AuNPs) have a tendency to aggregate through a preferred attachment of the 111] surfaces. We have used fully atomistic classical molecular dynamics simulations to obtain a quantitative estimate of this surface interaction using umbrella sampling (US) at various temperatures. To tune this surface interaction, we use polyamidoamine (PAMAM) dendrimer to coat the gold surface under various conditions. We observe a spontaneous adsorption of the protonated as well as nonprotonated PAMAM dendrimer on the AuNP surface. The adsorbed dendrimer on the nanoparticle surface strongly alters the interaction between the nanoparticles. We calculate the interaction between dendrimercoated AuNPs using US and show how the interaction between two bare faceted AuNPs can be tuned as a function of dendrimer concentration and charge (pH dependent) With appropriate choice of the dendrimer concentration and charge, two strongly interacting AuNPs can be made effectively noninteracting. Our simulation results demonstrate a strategy to tune the nanoparticle interaction, which can help in engineering self-assembly of such nanoparticles.

Effect of Stress and Interface defects on Photo Luminescence of Si nano-crystals embedded in SiO2

Effect of stress and interface defects on photo luminescence property of a silicon nano-crystal (Si-nc) embedded in amorphous silicon dioxide (a-SiO2) are studied in this paper using a self-consistent quantum-continuum based modeling framework. Si-ncs or quantum dots show photoluminescence at room temperature. Whether its origin is due to Si-nc/a-SiO2 interface defects or quantum confinement of carriers in Si-nc is still an outstanding question. Earlier reports have shown that stresses greater than 12 GPa change the indirect energy band gap structure of bulk Si to a direct energy band gap structure. Such stresses are observed very often in nanostructures and these stresses influence the carrier confinement energy significantly. Hence, it is important to determine the effect of stress in addition to the structure of interface defects on photoluminescence property of Si-nc. In the present work, first a Si-nc embedded in a-SiO2 is constructed using molecular dynamics simulation framework considering the actual conditions they are grown so that the interface and residual stress in the structure evolves naturally during formation. We observe that the structure thus created has an interface of about 1 nm thick consisting of 41.95% of defective states mostly Sin+ (n = 0 to 3) coordination states. Further, both the Si-nc core and the embedding matrix are observed to be under a compressive strain. This residual strain field is applied in an effective mass k.p Hamiltonian formulation to determine the energy states of the carriers. The photo luminescence property computed based on the carrier confinement energy and interface energy states associated with defects will be analysed in details in the paper.

Molecular dynamic simulations of elastomer structure and its influence on anisotropic stress under time-varying strain

We investigate the evolution of polymer structure and its influence on uniaxial anisotropic stress under time-varying uniaxial strain, and the role of external control variables such as temperature, strain rate, chain length, and density, using molecular dynamics simulation. At temperatures higher than glass transition, stress anisotropy in the system is reduced even though the bond stretch is greater at higher temperatures. There is a significant increase in the stress level with increasing density. At higher densities, the uncoiling of the chains is suppressed and the major contribution to the deformation is by internal deformation of the chains. At faster rates of loading stress anisotropy increases. The deformation mechanism is mostly due to bond stretch and bond bending rather than overall shape and size. Stress levels increase with longer chain length. There is a critical value of the functionality of the cross-linkers beyond which the uniaxial stress developed increases caused primarily by bond stretching due to increased constraint on the motion of the monomers. Stacking of the chains in the system also plays a dominant role in the behaviour in terms of excluded volume interactions. Low density, high temperature, low values of functionality of cross-linkers, and short chain length facilitate chain uncoiling and chain slipping in cross-linked polymers.

Solvent Sensitivity of Protein Unfolding: Dynamical Study of Chicken Villin Headpiece Subdomain in Water Ethanol Binary Mixture

We carry out a series of long atomistic molecular dynamics simulations to study the unfolding of a small protein, chicken villin headpiece (HP-36), in water-ethanol (EtOH) binary mixture. The prime objective of this work is to explore the sensitivity of protein unfolding dynamics toward increasing concentration of the cosolvent and unravel essential features of intermediates formed in search of a dynamical pathway toward unfolding. In water ethanol binary mixtures, HP-36 is found to unfold partially, under ambient conditions, that otherwise requires temperature as high as similar to 600 K to denature in pure aqueous solvent. However, an interesting course of pathway is observed to be followed in the process, guided by the formation of unique intermediates. The first step of unfolding is essentially the separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely, Phe-7, Phe-11, and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. The initial steps are similar to temperature-induced unfolding as well as chemical unfolding using DMSO as cosolvent. Subsequent unfolding steps follow a unique path. As water-ethanol shows composition-dependent anomalies, so do the details of unfolding dynamics. With an increase in cosolvent concentration, different partially unfolded intermediates are found to be formed. This is reflected in a remarkable nonmonotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states can be attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol. We further quantify the local dynamics of unfolding by using a Marcus-type theory.

Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

In addition to the biologically active monomer of the protein insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The insulin monomer contains two distinct surfaces, namely, the dimer forming surface (DFS) and the hexamer forming surface (HFS), that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behavior of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. In particular, we find that more structured water molecules, with higher residence times (similar to 300-500 ps), are present near HFS than those near DFS. A significant slowing down is observed in the decay of associated rotational auto time correlation functions of O-H bond vector of water in the vicinity of HFS. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly an outward H-atoms directed, clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inward H-atoms directed orientation owing to the flat curvature of hydrophobic surface and the interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces that reveal the significant differences between the two hydration layers.

De novo inference of protein function from coarse-grained dynamics

Inference of molecular function of proteins is the fundamental task in the quest for understanding cellular processes. The task is getting increasingly difficult with thousands of new proteins discovered each day. The difficulty arises primarily due to lack of high-throughput experimental technique for assessing protein molecular function, a lacunae that computational approaches are trying hard to fill. The latter too faces a major bottleneck in absence of clear evidence based on evolutionary information. Here we propose a de novo approach to annotate protein molecular function through structural dynamics match for a pair of segments from two dissimilar proteins, which may share even <10% sequence identity. To screen these matches, corresponding 1 mu s coarse-grained (CG) molecular dynamics trajectories were used to compute normalized root-mean-square-fluctuation graphs and select mobile segments, which were, thereafter, matched for all pairs using unweighted three-dimensional autocorrelation vectors. Our in-house custom-built forcefield (FF), extensively validated against dynamics information obtained from experimental nuclear magnetic resonance data, was specifically used to generate the CG dynamics trajectories. The test for correspondence of dynamics-signature of protein segments and function revealed 87% true positive rate and 93.5% true negative rate, on a dataset of 60 experimentally validated proteins, including moonlighting proteins and those with novel functional motifs. A random test against 315 unique fold/function proteins for a negative test gave >99% true recall. A blind prediction on a novel protein appears consistent with additional evidences retrieved therein. This is the first proof-of-principle of generalized use of structural dynamics for inferring protein molecular function leveraging our custom-made CG FF, useful to all. (C) 2014 Wiley Periodicals, Inc.

DNA Assisted Self-Assembly of PAMAM Dendrimers

We report DNA assisted self-assembly of polyamidoamine (PAMAM) dendrimers using all atom Molecular Dynamics (MD) simulations and present a molecular level picture of a DNA-linked PAMAM dendrimer nanocluster, which was first experimentally reported by Choi et al. (Nano Lett., 2004, 4, 391-397). We have used single stranded DNA (ssDNA) to direct the self-assembly process. To explore the effect of pH on this mechanism, we have used both the protonated (low pH) and nonprotonated (high pH) dendrimers. In all cases studied here, we observe that the DNA strand on one dendrimer unit drives self-assembly as it binds to the complementary DNA strand present on the other dendrimer unit, leading to the formation of a DNA-linked dendrimer dimeric complex. However, this binding process strongly depends on the charge of the dendrimer and length of the ssDNA. We observe that the complex with a nonprotonated dendrimer can maintain a DNA length dependent inter-dendrimer distance. In contrast, for complexes with a protonated dendrimer, the inter-dendrimer distance is independent of the DNA length. We attribute this observation to the electrostatic complexation of a negatively charged DNA strand with the positively charged protonated dendrimer.

Hydrophobic hydration driven self-assembly of curcumin in water: Similarities to nucleation and growth under large metastability, and an analysis of water dynamics at heterogeneous surfaces

As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of curcumin monomers from 4 to 16 monomers. We find that the curcumin monomers form clusters in a very definite pattern where they tend to aggregate both in parallel and anti-parallel orientation of the phenyl rings, often seen in the formation of beta-sheet in proteins. A considerable enhancement in the population of parallel alignments is observed with increasing the system size from 12 to 16 curcumin monomers. Due to the prevalence of such parallel alignment for large system size, a more closely packed cluster is formed with maximum number of hydrophobic contacts. We also follow the pathway of cluster growth, in particular the transition from the initial segregated to the final aggregated state. We find the existence of a metastable structural intermediate involving a number of intermediate-sized clusters dispersed in the solution. We have constructed a free energy landscape of aggregation where the metatsable state has been identified. The course of aggregation bears similarity to nucleation and growth in highly metastable state. The final aggregated form remains stable with the total exclusion of water from its sequestered hydrophobic core. We also investigate water structure near the cluster surface along with their orientation. We find that water molecules form a distorted tetrahedral geometry in the 1st solvation layer of the cluster, interacting rather strongly with the hydrophilic groups at the surface of the curcumin. The dynamics of such quasi-bound water molecules near the surface of curcumin cluster is considerably slower than the bulk signifying a restricted motion as often found in protein hydration layer. (C) 2014 AIP Publishing LLC.

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.