Levitation Effect: Role of Symmetry and Dependence of Diffusivity on the Bond Length of Homonuclear and Heteronuclear Diatomic Species

Molecular dynamics investigation of model diatomic species confined to the alpha-cages of zeolite NaY is reported. The dependence of self-diffusivity on the bond length of the diatomic species has been investigated. Three different sets of runs have been carried out. In the first set, the two atoms of the diatomic molecule interact with the zeolite atoms with equal strength (example, O-2, the symmetric case). In the second and third sets which correspond to asymmetric cases, the two atoms of the diatomic molecule interact with unequal strengths (example, CO). The result for the symmetric case exhibits a well-defined maximum in self-diffusivity for an intermediate bond length. In contrast to this, the intermediate asymmetry leads to a less pronounced maximum. For the large asymmetric case, the maximum is completely absent. These findings are analyzed by computing a number of related properties. These results provide a direct confirmation at the microscopic level of the suggestion by Derouane that the supermobility observed experimentally by Kemball has its origin in the mutual cancellation of forces. The maximum in diffusivity from molecular dynamics is seen at the value predicted by the levitation effect. Further, these findings suggest a role for symmetry in the existence of a diffusivity maximum as a function of diameter of the diffusant often referred to as the levitation effect. The nature of the required symmetry for the existence of anomalous diffusivity is interaction symmetry which is different from that normally encountered in crystallography.

Interaction of single-walled carbon nanotubes with poly(propyl ether imine) dendrimers

We study the complexation of nontoxic, native poly(propyl ether imine) dendrimers with single-walled carbon nanotubes (SWNTs). The interaction was monitored by measuring the quenching of inherent fluorescence of the dendrimer. The dendrimer-nanotube binding also resulted in the increased electrical resistance of the hole doped SWNT, due to charge-transfer interaction between dendrimer and nanotube. This charge-transfer interaction was further corroborated by observing a shift in frequency of the tangential Raman modes of SWNT. We also report the effect of acidic and neutral pH conditions on the binding affinities. Experimental studies were supplemented by all atom molecular dynamics simulations to provide a microscopic picture of the dendrimer-nanotube complex. The complexation was achieved through charge transfer and hydrophobic interactions, aided by multitude of oxygen, nitrogen, and n-propyl moieties of the dendrimer. (C) 2011 American Institute of Physics. doi:10.1063/1.3561308]

Orientational relaxation in dipolar systems: How much do we understand the role of correlations

Dipolar systems, both liquids and solids, constitute a class of naturally abundant systems that are important in all branches of natural science. The study of orientational relaxation provides a powerful method to understand the microscopic properties of these systems and, fortunately, there are many experimental tools to study orientational relaxation in the condensed phases. However, even after many years of intense research, our understanding of orientational relaxation in dipolar systems has remained largely imperfect. A major hurdle towards achieving a comprehensive understanding is the long range and complex nature of dipolar interactions which also made reliable theoretical study extremely difficult. These difficulties have led to the development of continuum model based theories, which although they provide simple, elegant expressions for quantities of interest, are mostly unsatisfactory as they totally neglect the molecularity of inter-molecular interactions. The situation has improved in recent years because of renewed studies, led by computer simulations. In this review, we shall address some of the recent advances, with emphasis on the work done in our laboratory at Bangalore. The reasons for the failure of the continuum model, as revealed by the recent Brownian dynamics simulations of the dipolar lattice, are discussed. The main reason is that the continuum model predicts too fast a decay of the torque-torque correlation function. On the other hand, a perturbative calculation, based on Zwanzig's projection operator technique, provides a fairly satisfactory description of the single particle orientational dynamics for not too strongly polar dipolar systems. A recently developed molecular hydrodynamic theory that properly includes the effects of intermolecular orientational pair correlations provides an even better description of the single-particle orientational dynamics. We also discuss the rank dependence of the dielectric friction. The other topics reviewed here includes dielectric relaxation and solvation dynamics, as they are intimately connected with orientational relaxation. Recent molecular dynamics simulations of the dipolar lattice are also discussed. The main theme of the present review is to understand the effects of intermolecular interactions on orientational relaxation. The presence of strong orientational pair correlation leads to a strong coupling between the single particle and the collective dynamics. This coupling can lead to rich dynamical properties, some of which are detailed here, while a major part remains yet unexplored.

Modeling of angiogenin-3 '-NMP complex

Angiogenin belongs to the Ribonuclease superfamily and has a weak enzymatic activity that is crucial for its biological function of stimulating blood vessel growth. Structural studies on ligand bound Angiogenin will go a long way in understanding the mechanism of the protein as well as help in designing drugs against it. In this study we present the first available structure of nucleotide ligand bound Angiogenin obtained by computer modeling. The importance of this study in itself notwithstanding, is a precursor to modeling a full dinucleotide substrate onto Angiogenin. Bovine Angiogenin, the structure of which has been solved at a high resolution, was earlier subjected to Molecular Dynamics simulations for a nanosecond. The MD structures offer better starting points for docking as they offer lesser obstruction than the crystal structure to ligand binding. The MD structure with the least serious short contacts was modeled to obtain a steric free Angiogenin - 3' mononucleotide complex structure. The structures were energetically minimized and subjected to a brief spell of Molecular Dynamics. The results of the simulation show that all the li,ligand-Angiogenin interactions and hydrogen bonds are retained, redeeming the structure and docking procedure. Further, following ligand - protein interactions in the case of the ligands 3'-CMP and 3'-UMP we were able to speculate on how Angiogenin, a predominantly prymidine specific ribonuclease prefers Cytosine to Uracil in the first base position.

Neutral Nanosheets that Gel: Exfoliated Layered Double Hydroxides in Toluene

A simple strategy to exfoliate inorganic layered double hydroxide (LDH) solids to their ultimate constituent, intact single layers of nanometer thickness and micrometer size, is presented. The procedure involves intercalation of an ionic surfactant that forms a hydrophobic anchored surfactant bilayer in the galleries of the solid followed by simply stirring the intercalated solid in toluene. The method is rapid but at the same time gentle enough to produce exfoliated nanosheets of regular morphology that are electrically neutral and form stable gels at higher concentrations. In this Letter, we describe the phenomena and use molecular dynamics simulations to show that exfoliation of the LDH in toluene is a consequence of the modification of the cohesive dispersive interactions between surfactant chains anchored on opposing inorganic sheets by the toluene molecules. The toluene molecules function as a molecular glue, holding the surfactant-anchored LDH sheets together, leading to gel formation.

Inverse surface melting in confined clusters: Ar-13 in zeolite L

Monte Carlo and molecular dynamics simulations on an Ar-13 cluster in zeolite L have been carried out at a series of temperatures to understand the rigid-nonrigid transition corresponding to the solid-liquid transition exhibited by the free Ar-13 cluster. The icosahedral geometry of the free cluster is no longer preferred when the cluster is confined in the zeolite. The root-mean-squared pair distance fluctuation, delta, exhibits a sharp, well-defined rigid-nonrigid transition at 17 K as compared to 27 K for the free cluster. Multiple peaks in the distribution of short-time averages of the guest-host interaction energy indicate coexistence of two phases.; It is shown that this transition is associated with the inner atoms becoming mobile at 17 K even while the outer layer atoms, which are in close proximity to the zeolitic wall, continue to be comparatively immobile. This may be contrasted with the melting of large free clusters of 40 or more atoms which exhibit surface melting. Guest-host interactions seem to play a predominant role in determining the properties of confined clusters. We demonstrate that the volume of the cluster increases rather sharply at 17 and 27 K respectively for the confined and the free cluster. Power spectra suggest that the motion of the inner atoms is generally parallel to the atoms which form the cage wall.

Levitation effect and its dependence on sorbate concentration

Earlier work has reported the existence of a diffusion anomaly in porous solids at dilute sorbate concentrations. In this work we have carried out molecular dynamics simulations at higher sorbate concetrations. Results indicate the persistence of a diffusion anomaly even at significantly higher sorbate concentrations, which means that this anomaly can be used for separation of mixtures under conditions prevailing in industries.

Temperature dependence of the levitation effect - Implications for separation of multicomponent mixtures

Sufficiently long molecular dynamics simulations have been carried out on spherical monatomic sorbates in NaY zeolite, interacting via simple Lennard-Jones potentials, to investigate the dependence of the levitation effect on the temperature. Simulations carried out in the range 100-300 K suggest that the anomalous peak in the diffusion coefficient (observed when the levitation parameter, gamma, is near unity) decreases in intensity with increase in temperature. The rate of cage-to-cage migrations also exhibits a similar trend. The activation energy obtained from Arrhenius plots is found to exhibit a minimum when the diffusion coefficient is a maximum, corresponding to the gamma approximate to 1 sorbate diameter. In the linear or normal regime, the activation energy increases with increase in sorbate diameter until it shows a sharp decrease in the anomalous regime. Locations and energies of the adsorption sites and their dependence on the sorbate size gives interesting insight into the nature of the underlying potential-energy surface and further explain the observed trend in the activation energy with sorbate size. Cage residence times, tau(c), show little or no change with temperature for the sorbate with diameter corresponding to gamma approximate to 1, whereas there is a significant decrease in tau(c) with increase in temperature for sorbates in the linear regime. The implications of the present study for the separation of mixtures of sorbates are discussed.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.

Structural basis for poor uracil excision from hairpin DNA

Two-dimensional NMR and molecular dynamics simulations have been used to determine the three-dimensional structures of two hairpin DNA structures: d-CTAGAG GATCCUTTTGGATCCT (abbreviated as U1-hairpin) and d-CTAGAGGATCCTTUTGGATCCT (abbreviated as U3-hairpin). The (1) H resonances of both of these hairpin structures have been assigned almost completely. NMR restrained molecular dynamics and energy minimization procedures have been used to describe the three-dimensional structures of these hairpins. This study and concurrent NMR structural studies on two other d-CTAGAGGA TCCTUTTGGATCCT (abbreviated as U2-hairpin) and d-CTAGAGGATCCTTTUGGATCCT (abbreviated as U4-hairpin) have shed light upon various interactions reported between Echerichia coli uracil DNA glycosylase (UDG) and uracil-containing DNA. The backbone torsion angles, which partially influence the local conformation of U12 and U14 in U1 and U3-hairpins, respectively, are probably locked in the trans conformation as in the case of U-13 in the U2-hairpin. Such a stretched-out backbone conformation in the vicinity of U-12 and U-14 is thought to be the reason why the K-m value is poor for U1- and U3-hairpins as it is for the U2-hairpin. Furthermore, the bases U-12 and U-14 in both U1- and U3-hairpins adopt an anti conformation, in contrast with the base conformation of U-13 in the U2-hairpin, which adopts a syn conformation. The clear discrepancy observed in the U-base orientation with respect to the sugar moieties could explain why the V-max value is 10- to 20-fold higher for the U1- and U3-hairpins compared with the U2-hairpin. Taken together, these observations support our interpretation that the unfavourable backbone results in a poor K-m value, whereas the unfavourable nucleotide conformation results in a poor V-max value. These two parameters therefore make the U1- and U3-hairpins better substrates for UDG compared with the U2-hairpin, as reported earlier [Kumar, N. V. & Varshney, U. (1997) Nucleic Acids Res. 25, 2336-2343.].

Study of translational and rotational mobility and orientational preference of ethane in one-dimensional channels

Structural and dynamical properties of ethane in one-dimensional channels of AlPO4-5 and carbon nanotube have been investigated at dilute concentration with the help of molecular dynamics simulation. Density distributions and orientational structure of ethane have been analyzed. Repulsive interactions seem to play an important role when ethane is located in the narrow part of the AlPO4-5 channel. In AlPO4-5, parallel orientation is predominant over perpendicular orientation except when ethane is located in the broader part of the channel. Unlike in the case of single-file diffusion, our results in carbon nanotube show that at dilute concentrations the mean squared displacement, mu(2)(t) approximate to t(alpha), alpha = 1.8. The autocorrelation function for the z-component of angular velocity of ethane in space-fixed frame of reference shows a pronounced negative correlation. This is attributed to the restriction in the movement of ethane along the x- and y- directions. It is seen that the ratio of reorientational correlation times does not follow the Debye model for confined ethane but it is closer to the predictions of the Debye model for bulk ethane.

Probing the Allosteric Mechanism in Pyrrolysyl-tRNA Synthetase Using Energy-Weighted Network Formalism

Pyrrolysyl-tRNA synthetase (PyIRS) is an atypical enzyme responsible for charging tRNA(Pyl) with pyrrolysine, despite lacking precise tRNA anticodon recognition. This dimeric protein exhibits allosteric regulation of function, like any other tRNA synthetases. In this study we examine the paths of allosteric communication at the atomic level, through energy-weighted networks of Desulfitobacterium hafniense PyIRS (DhPyIRS) and its complexes with tRNA(Pyl) and activated pyrrolysine. We performed molecular dynamics simulations of the structures of these complexes to obtain an ensemble conformation-population perspective. Weighted graph parameters relevant to identifying key players and ties in the context of social networks such as edge/node betweenness, closeness index, and the concept of funneling are explored in identifying key residues and interactions leading to shortest paths of communication in the structure networks of DhPylRS. Further, the changes in the status of important residues and connections and the costs of communication due to ligand induced perturbations are evaluated. The optimal, suboptimal, and preexisting paths are also investigated. Many of these parameters have exhibited an enhanced asymmetry between the two subunits of the dimeric protein, especially in the pretransfer complex, leading us to conclude that encoding of function goes beyond the sequence/structure of proteins. The local and global perturbations mediated by appropriate ligands and their influence on the equilibrium ensemble of conformations also have a significant role to play in the functioning of proteins. Taking a comprehensive view of these observations, we propose that the origin of many functional aspects (allostery rand half-sites reactivity in the case of DhPyIRS) lies in subtle rearrangements of interactions and dynamics at a global level.

Conformational basis for substrate recognition and regulation of catalytic activity in Staphylococcus aureus nucleoside di-phosphate kinase

Nucleoside diphosphate kinases (NDK) are characterized by high catalytic turnover rates and diverse substrate specificity. These features make this enzyme an effective activator of a pro-drug an application that has been actively pursued for a variety of therapeutic strategies. The catalytic mechanism of this enzyme is governed by a conserved histidine that coordinates a magnesium ion at the active site. Despite substantial structural and biochemical information on NDK, the mechanistic feature of the phospho-transfer that leads to auto-phosphorylation remains unclear. While the role of the histidine residue is well documented, the other active site residues, in particular the conserved serine remains poorly characterized. Studies on some homologues suggest no role for the serine residue at the active site, while others suggest a crucial role for this serine in the regulation and quaternary association of this enzyme in some species. Here we report the biochemical features of the Staphylococcus aureus NDK and the mutant enzymes. We also describe the crystal structures of the apo-NDK, as a transition state mimic with vanadate and in complex with different nucleotide substrates. These structures formed the basis for molecular dynamics simulations to understand the broad substrate specificity of this enzyme and the role of active site residues in the phospho-transfer mechanism and oligomerization. Put together, these data suggest that concerted changes in the conformation of specific residues facilitate the stabilization of nucleotide complexes thereby enabling the steps involved in the ping-pong reaction mechanism without large changes to the overall structure of this enzyme. (C) 2011 Elsevier B.V. All rights reserved.

Structural Rigidity of Paranemic Crossover and Juxtapose DNA Nanostructures

Crossover motifs are integral components for designing DNA-based nanostructures and nanomechanical devices due to their enhanced rigidity compared to the normal B-DNA. Although the structural rigidity of the double helix B-DNA has been investigated extensively using both experimental and theoretical tools, to date there is no quantitative information about structural rigidity and the mechanical strength of parallel crossover DNA motifs. We have used fully atomistic molecular dynamics simulations in explicit solvent to get the force-extension curve of parallel DNA nanostructures to characterize their mechanical rigidity. In the presence of monovalent Na(+) ions, we find that the stretch modulus (gamma(1)) of the paranemic crossover and its topoisomer JX DNA structure is significantly higher (similar to 30%) compared to normal B-DNA of the same sequence and length. However, this is in contrast to the original expectation that these motifs are almost twice as rigid compared to the double-stranded B-DNA. When the DNA motif is surrounded by a solvent with Mg(2+) counterions, we find an enhanced rigidity compared to Na(+) environment due to the electrostatic screening effects arising from the divalent nature of Mg(2+) ions. To our knowledge, this is the first direct determination of the mechanical strength of these crossover motifs, which can be useful for the design of suitable DNA for DNA-based nanostructures and nanomechanical devices with improved structural rigidity.