Self-Organization of n-Alkane Chains in Water: Length Dependent Crossover from Helix and Toroid to Molten Globule

We demonstrate a chain length dependent crossover in the structural properties of linear hydrocarbon (n-alkane) chains using detailed atomistic simulations in explicit water. We identify a number of exotic structures of the polymer chain through energy minimization of representative snapshots collected from molecular dynamics trajectory. While the collapsed state is ring-like (circular) for small chains (CnH2n+2; n <= 20) and spherical for very long ones (n = 100), we find the emergence of ordered helical structures at intermediate lengths (n similar to 40). We find different types of disordered helices and toroid-like structures at n = 60. We also report a sharp transition in the stability of the collapsed state as a function of the chain length through relevant free energy calculations. While the collapsed state is only marginally metastable for C20H42, a clear bistable free energy surface emerges only when the chain is about 30 monomers long. For n = 30, the polymer exhibits an intermittent oscillation between the collapsed and the coil structures, characteristic of two stable states separated by a small barrier.

Ions in water: Role of attractive interactions in size dependent diffusivity maximum

A molecular dynamics study of model ions in water is reported. The van der Waals diameter of both the cations and anions is varied. We have carried out two sets of simulations-with and without dispersion interaction-between the ion and water. Self-diffusivity of the ions exhibits an anomalous maximum as a function of the van der Waals diameter for both these sets. This existence of a maximum in self-diffusivity when there is no dispersion interaction between the ion and the water is attributed to the attractive term from electrostatic interactions. Detailed analysis of this effect shows that the solvent shell is more strongly defined in the presence of dispersion interactions. A smaller ion exhibits biexponential decay while a single exponential decay is seen for the ion with maximum diffusivity in the self-part of the intermediate scattering function. The solvent structure around the ion appears to determine much of the dynamics of the ion. Interesting trends are seen in the activation energies and these can be understood in terms of the levitation effect. (C) 2010 American Institute of Physics. doi:10.1063/1.3481656]

Dynamics of I-asterisk(P-2(1/2)) production inthe ultraviolet photodissociation of alkyn iodides

The relative quantum yields, phi*, for the production of I*(P-2(1/2)) at 266, 280, and similar to 305 nm are reported for a series of primary alkyl iodides using the technique of two-photon laser-induced fluorescence for the detection of I(P-2(3/2)) and I*(P-2(1/2)) atoms. Results are analyzed by invoking the impulsive energy disposal model, which summarizes the dynamics of dissociation as a single parameter. Comparison of our data with those calculated by a more sophisticated time-dependent quantum mechanical model is also made. Near the red edge of the alkyl iodide A band, absorption contribution from the (3)Q(1) state is important and the dynamics near the (3)Q(0)-(1)Q(1) curve-crossing region seem to be influenced by the kinematics of the dissociation process

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.

Temperature dependent transport studies in InN quantum dots grown by droplet epitaxy on silicon nitride/Si substrate

InN quantum dots (QDs) were fabricated on silicon nitride/Si (111) substrate by droplet epitaxy. Single-crystalline structure of InN QDs was verified by transmission electron microscopy, and the chemical bonding configurations of InN QDs were examined by x-ray photoelectron spectroscopy. Photoluminescence measurement shows a slight blue shift compared to the bulk InN, arising from size dependent quantum confinement effect. The interdigitated electrode pattern was created and current-voltage (I-V) characteristics of InN QDs were studied in a metal-semiconductor-metal configuration in the temperature range of 80-300K. The I-V characteristics of lateral grown InN QDs were explained by using the trap model. (C) 2011 American Institute of Physics. [doi:10.1063/1.3651762]

Chirality dependent elastic properties of single-walled boron nitride nanotubes under uniaxial and torsional loading

The elastic behavior of single-walled boron nitride nanotubes is studied under axial and torsional loading. Molecular dynamics simulation is carried out with a tersoff potential for modeling the interatomic interactions. Different chiral configurations with similar diameter are considered to study the effect of chirality on the elastic and shear moduli. Furthermore, the effects of tube length on elastic modulus are also studied by considering different aspects ratios. It is observed that both elastic and shear moduli depend upon the chirality of a nanotube. For aspect ratios less than 15, the elastic modulus reduces monotonically with an increase in the chiral angle. For chiral nanotubes, the torsional response shows a dependence on the direction of loading. The difference between the shear moduli against and along the chiral twist directions is maximum for chiral angle of 15 degrees, and zero for zigzag (0 degrees) and armchair (30 degrees) configurations. (C) 2014 AIP Publishing LLC.

Anomalous power law decay in solvation dynamics of DNA: a mode coupling theory analysis of ion contribution

Several time dependent fluorescence Stokes shift (TDFSS) experiments have reported a slow power law decay in the hydration dynamics of a DNA molecule. Such a power law has neither been observed in computer simulations nor in some other TDFSS experiments. Here we observe that a slow decay may originate from collective ion contribution because in experiments DNA is immersed in a buffer solution, and also from groove bound water and lastly from DNA dynamics itself. In this work we first express the solvation time correlation function in terms of dynamic structure factors of the solution. We use mode coupling theory to calculate analytically the time dependence of collective ionic contribution. A power law decay in seen to originate from an interplay between long-range probe-ion direct correlation function and ion-ion dynamic structure factor. Although the power law decay is reminiscent of Debye-Falkenhagen effect, yet solvation dynamics is dominated by ion atmosphere relaxation times at longer length scales (small wave number) than in electrolyte friction. We further discuss why this power law may not originate from water motions which have been computed by molecular dynamics simulations. Finally, we propose several experiments to check the prediction of the present theoretical work.

Composition dependent multiple structural transformations of myoglobin in aqueous ethanol solution: A combined experimental and theoretical study

Experimental studies (circular dichroism and ultra-violet (UV) absorption spectra) and large scale atomistic molecular dynamics simulations (accompanied by order parameter analyses) are combined to establish a number of remarkable (and unforeseen) structural transformations of protein myoglobin in aqueous ethanol mixture at various ethanol concentrations. The following results are particularly striking. (1) Two well-defined structural regimes, one at x(EtOH) similar to 0.05 and the other at x(EtOH) similar to 0.25, characterized by formation of distinct partially folded conformations and separated by a unique partially unfolded intermediate state at x(EtOH) similar to 0.15, are identified. (2) Existence of non-monotonic composition dependence of (i) radius of gyration, (ii) long range contact order, (iii) residue specific solvent accessible surface area of tryptophan, and (iv) circular dichroism spectra and UV-absorption peaks are observed. Interestingly at x(EtOH) similar to 0.15, time averaged value of the contact order parameter of the protein reaches a minimum, implying that this conformational state can be identified as a molten globule state. Multiple structural transformations well known in water-ethanol binary mixture appear to have considerably stronger effects on conformation and dynamics of the protein. We compare the present results with studies in water-dimethyl sulfoxide mixture where also distinct structural transformations are observed along with variation of co-solvent composition. (C) 2015 AIP Publishing LLC.

How does contrasting dependence of impurity-atom diffusivity on the density of host disordered medium arise?

We report results of molecular dynamics investigations into neutral impurity diffusing within an amorphous solid as a function of the size of the diffusant and density of the host amorphous matrix. We find that self diffusivity exhibits an anomalous maximum as a function of the size of the impurity species. An analysis of properties of the impurity atom with maximum diffusivity shows that it is associated with lower mean square force, reduced backscattering of velocity autocorrelation function, near-exponential decay of the intermediate scattering function (as compared to stretched-exponential decay for other sizes of the impurity species) and lower activation energy. These results demonstrate the existence of size-dependent diffusivity maximum in disordered solids. Further, we show that the diffusivity maximum is observed at lower impurity diameters with increase in density. This is explained in terms of the Levitation parameter and the void structure of the amorphous solid. We demonstrate that these results imply contrasting dependence of self diffusivity (D) on the density of the amorphous matrix, p. D increases with p for small sizes of the impurity but shows an increase followed by a decrease for intermediate sizes of the impurity atom. For large sizes of the impurity atom, D decreases with increase in p. These contrasting dependence arises naturally from the existence of Levitation Effect.

Dynamical Transition of Water in the Grooves of DNA Duplex at Low Temperature

At low temperature (below its freezing/melting temperature), liquid water under confinement is known to exhibit anomalous dynamical features. Here we study structure and dynamics of water in the grooves of a long DNA duplex using molecular dynamics simulations with TIP5P potential at low temperature. We find signatures of a dynamical transition in both translational and orientational dynamics of water molecules in both the major and the minor grooves of a DNA duplex. The transition occurs at a slightly higher temperature (TGL ≈ 255 K) than the temperature at which the bulk water is found to undergo a dynamical transition, which for the TIP5P potential is at 247 K. Groove water, however, exhibits markedly different temperature dependence of its properties from the bulk. Entropy calculations reveal that the minor groove water is ordered even at room temperature, and the transition at T ≈ 255 K can be characterized as a strong-to-strong dynamical transition. Confinement of water in the grooves of DNA favors the formation of a low density four-coordinated state (as a consequence of enthalpy−entropy balance) that makes the liquid−liquid transition stronger. The low temperature water is characterized by pronounced tetrahedral order, as manifested in the sharp rise near 109° in the O−O−O angle distribution. We find that the Adams−Gibbs relation between configurational entropy and translational diffusion holds quite well when the two quantities are plotted together in a master plot for different region of aqueous DNA duplex (bulk, major, and minor grooves) at different temperatures. The activation energy for the transfer of water molecules between different regions of DNA is found to be weakly dependent on temperature.

Intra and Inter-Molecular Communications Through Protein Structure Network

Communication within and across proteins is crucial for the biological functioning of proteins. Experiments such as mutational studies on proteins provide important information on the amino acids, which are crucial for their function. However, the protein structures are complex and it is unlikely that the entire responsibility of the function rests on only a few amino acids. A large fraction of the protein is expected to participate in its function at some level or other. Thus, it is relevant to consider the protein structures as a completely connected network and then deduce the properties, which are related to the global network features. In this direction, our laboratory has been engaged in representing the protein structure as a network of non-covalent connections and we have investigated a variety of problems in structural biology, such as the identification of functional and folding clusters, determinants of quaternary association and characterization of the network properties of protein structures. We have also addressed a few important issues related to protein dynamics, such as the process of oligomerization in multimers, mechanism on protein folding, and ligand induced communications (allosteric effect). In this review we highlight some of the investigations which we have carried out in the recent past. A review on protein structure graphs was presented earlier, in which the focus was on the graphs and graph spectral properties and their implementation in the study of protein structure graphs/networks (PSN). In this article, we briefly summarize the relevant parts of the methodology and the focus is on the advancement brought out in the understanding of protein structure-function relationships through structure networks. The investigations of structural/biological problems are divided into two parts, in which the first part deals with the analysis of PSNs based on static structures obtained from x-ray crystallography. The second part highlights the changes in the network, associated with biological functions, which are deduced from the network analysis on the structures obtained from molecular dynamics simulations.

Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W I and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the Ob atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.

Stress-induced phase transformation and pseudo-elastic/pseudo-plastic recovery in intermetallic Ni-Al nanowires

Extensive molecular dynamics (MD) simulations have been performed in a B2-NiAl nanowire using an embedded atom method (EAM) potential. We show a stress induced B2 -> body-centered-tetragonal (BCT) phase transformation and a novel temperature and cross-section dependent pseudo-elastic/pseudo-plastic recovery from such an unstable BCT phase with a recoverable strain of similar to 30% as compared to 5-8% in polycrystalline materials. Such a temperature and cross-section dependent pseudo-elastic/pseudo-plastic strain recovery can be useful in various interesting applications of shape memory and strain sensing in nanoscale devices. Effects of size, temperature, and strain rate on the structural and mechanical properties have also been analyzed in detail. For a given size of the nanowire the yield stress of both the B2 and the BCT phases is found to decrease with increasing temperature, whereas for a given temperature and strain rate the yield stress of both the B2 and the BCT phase is found to increase with increase in the cross-sectional dimensions of the nanowire. A constant elastic modulus of similar to 80 GPa of the B2 phase is observed in the temperature range of 200-500 K for nanowires of cross-sectional dimensions in the range of 17.22-28.712 angstrom, whereas the elastic modulus of the BCT phase shows a decreasing trend with an increase in the temperature.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

Role of the triplet state in the green emission peak of polyfluorene films: A time evolution study

The blue emission of ethyl-hexyl substituted polyfluorene (PF2/6) films is accompanied by a low energy green emission peak around 500 nm in inert atmosphere. The intensity of this 500 nm peak is large in electroluminescence (EL) compared to photoluminescence (PL)measurements. Furthermore, the green emission intensity reduces dramatically in the presence of molecular oxygen. To understand this, we have modeled various nonradiative processes by time dependent quantum many body methods. These are (i) intersystem crossing to study conversion of excited singlets to triplets leading to a phosphorescence emission, (ii) electron-hole recombination (e-hR) process in the presence of a paramagnetic impurity to follow the yield of triplets in a polyene system doped with paramagnetic metal atom, and (iii) quenching of excited triplet states in the presence of oxygen molecules to understand the low intensity of EL emission in ambient atmosphere, when compared with that in nitrogen atmosphere. We have employed the Pariser-Parr-Pople Hamiltonian to model the molecules and have invoked electron-electron repulsions beyond zero differential approximation while treating interactions between the organic molecule and the rest of the system. Our time evolution methods show that there is a large cross section for triplet formation in the e-hR process in the presence of paramagnetic impurity with degenerate orbitals. The triplet yield through e-hR process far exceeds that in the intersystem crossing pathway, clearly pointing to the large intensity of the 500 nm peak in EL compared to PL measurements. We have also modeled the triplet quenching process by a paramagnetic oxygen molecule which shows a sizable quenching cross section especially for systems with large sizes. These studies show that the most probable origin of the experimentally observed low energy EL emission is the triplets.