Solvation dynamics of tryptophan in water-dimethyl sulfoxide binary mixture: In search of molecular origin of composition dependent multiple anomalies

Experimental and simulation studies have uncovered at least two anomalous concentration regimes in water-dimethyl sulfoxide (DMSO) binary mixture whose precise origin has remained a subject of debate. In order to facilitate time domain experimental investigation of the dynamics of such binary mixtures, we explore strength or extent of influence of these anomalies in dipolar solvation dynamics by carrying out long molecular dynamics simulations over a wide range of DMSO concentration. The solvation time correlation function so calculated indeed displays strong composition dependent anomalies, reflected in pronounced non-exponential kinetics and non-monotonous composition dependence of the average solvation time constant. In particular, we find remarkable slow-down in the solvation dynamics around 10%-20% and 35%-50% mole percentage. We investigate microscopic origin of these two anomalies. The population distribution analyses of different structural morphology elucidate that these two slowing down are reflections of intriguing structural transformations in water-DMSO mixture. The structural transformations themselves can be explained in terms of a change in the relative coordination number of DMSO and water molecules, from 1DMSO:2H(2)O to 1H(2)O:1DMSO and 1H(2)O:2DMSO complex formation. Thus, while the emergence of first slow down (at 15% DMSO mole percentage) is due to the percolation among DMSO molecules supported by the water molecules (whose percolating network remains largely unaffected), the 2nd anomaly (centered on 40%-50%) is due to the formation of the network structure where the unit of 1DMSO:1H(2)O and 2DMSO:1H(2)O dominates to give rise to rich dynamical features. Through an analysis of partial solvation dynamics an interesting negative cross-correlation between water and DMSO is observed that makes an important contribution to relaxation at intermediate to longer times.

Simulations reveal that the HIV-1 gp120-CD4 complex dissociates via complex pathways and is a potential target of the polyamidoamine (PAMAM) dendrimer

The polyamidoamine (PAMAM) dendrimer prevents HIV-1 entry into target cells in vitro. Its mechanism of action, however, remains unclear and precludes the design of potent dendrimers targeting HIV-1 entry. We employed steered molecular dynamics simulations to examine whether the HIV-1 gp120-CD4 complex is a target of PAMAM. Our simulations mimicked single molecule force spectroscopy studies of the unbinding of the gp120-CD4 complex under the influence of a controlled external force. We found that the complex dissociates via complex pathways and defies the standard classification of adhesion molecules as catch and slip bonds. When the force loading rate was large, the complex behaved as a slip bond, weakening gradually. When the loading rate was small, the complex initially strengthened, akin to a catch bond, but eventually dissociated over shorter separations than with large loading rates. PAMAM docked to gp120 and destabilized the gp120-CD4 complex. The rupture force of the complex was lowered by PAMAM. PAMAM disrupted salt bridges and hydrogen bonds across the gp120-CD4 interface and altered the hydration pattern of the hydrophobic cavity in the interface. In addition, intriguingly, PAMAM suppressed the distinction in the dissociation pathways of the complex between the small and large loading rate regimes. Taken together, our simulations reveal that PAMAM targets the gp120-CD4 complex at two levels: it weakens the complex and also alters its dissociation pathway, potentially inhibiting HIV-1 entry.

Ozone and carbon trioxide synthesis by low energy ion implantation onto solid carbon dioxide and implications to astrochemistry

Ion implantation experiments were carried out on amorphous (30 K) and crystalline (80 K) solid CO2 using both reactive (D+, H+) and non-reactive (He+) ions, simulating different irradiation environments on satellite and dust grain surfaces. Such ion irradiation synthesized several new species in the ice including ozone (O-3), carbon trioxide (CO3), and carbon monoxide (CO) the main dissociation product of carbon dioxide. The yield of these products was found to be strongly dependent upon the ion used for irradiation and the sample temperature. Ion implantation changes the chemical composition of the ice with recorded infrared spectra clearly showing the coexistence of D-3h and C-2v isomers of CO3, for the first time, in ion irradiated CO2 ice. (C) 2013 AIP Publishing LLC.

Note: Effect of fluid phase compositions on the formation of substitutionally ordered solid phases from a binary mixture of oppositely charged colloidal suspensions

A binary mixture of oppositely charged colloidal particles can self-assemble into either a substitutionally ordered or substitutionally disordered crystalline phase depending on the nature and strength of interactions among the particles. An earlier study had mapped out favorable inter-particle interactions for the formation of substitutionally ordered crystalline phases from a fluid phase using Monte Carlo molecular simulations along with the Gibbs-Duhem integration technique. In this paper, those studies are extended to determine the effect of fluid phase composition on formation of substitutionally ordered solid phases.

Microwave, infrared-microwave double resonance, and theoretical studies of C2H4 center dot center dot center dot H2S complex

In this manuscript, rotational spectra of four new isotopologues of the S-H center dot center dot center dot pi bonded C2H4 center dot center dot center dot H2S complex, i.e., C2D4 center dot center dot center dot H2S, C2D4 center dot center dot center dot D2S, C2D4 center dot center dot center dot HDS, and (CCH4)-C-13 center dot center dot center dot H2S have been reported and analyzed. All isotopologues except C2D4 center dot center dot center dot HDS show a four line pattern whereas a doubling of the transition frequencies was observed for C2D4 center dot center dot center dot HDS. These results together with our previous report on the title complex M. Goswami, P. K. Mandal, D. J. Ramdass, and E. Arunan, Chem. Phys. Lett. 393(1-3), 22-27 (2004)] confirm that both subunits (C2H4 and H2S) are involved in large amplitude motions leading to a splitting of each rotational transition to a quartet. Further, the results also confirm that the motions which are responsible for the observed splittings involve both monomers. Molecular symmetry group analysis, considering the interchange of equivalent H atoms in H2S and C2H4 could explain the observed four line pattern and their intensities in the microwave spectrum. In addition, hydride stretching fundamentals of the complex were measured using coherence-converted population transfer Fourier Transform Microwave-infrared (IR-MW double resonance) experiments in the S-H and C-H stretch regions. Changes in the tunneling splittings upon vibrational excitation are consistent with the isotopic dependence of pure rotational transitions. A complexation shift of 2.7-6.5 cm(-1) has been observed in the two fundamental S-H stretching modes of the H2S monomer in the complex. Vibrational pre-dissociation in the bound S-H stretch has been detected whereas the instrument-limited line-shapes in other S-H and C-H stretches indicate slower pre-dissociation rate. Some local perturbations in the vibrational spectra have been observed. Two combination bands have been observed corresponding to both the S-H stretching fundamentals and what appears to be the intermolecular stretching mode at 55 cm(-1). The tunneling splitting involved in the rotation of C2H4 unit has been deduced to be 1.5 GHz from the IR-MW results. In addition, ab initio barrier heights derived for different motions of the monomers support the experimental results and provide further insight into the motions causing the splitting. (C) 2013 AIP Publishing LLC.

System Reliability Analysis of Flexible Pavements

The uncertainty in material properties and traffic characterization in the design of flexible pavements has led to significant efforts in recent years to incorporate reliability methods and probabilistic design procedures for the design, rehabilitation, and maintenance of pavements. In the mechanistic-empirical (ME) design of pavements, despite the fact that there are multiple failure modes, the design criteria applied in the majority of analytical pavement design methods guard only against fatigue cracking and subgrade rutting, which are usually considered as independent failure events. This study carries out the reliability analysis for a flexible pavement section for these failure criteria based on the first-order reliability method (FORM) and the second-order reliability method (SORM) techniques and the crude Monte Carlo simulation. Through a sensitivity analysis, the most critical parameter affecting the design reliability for both fatigue and rutting failure criteria was identified as the surface layer thickness. However, reliability analysis in pavement design is most useful if it can be efficiently and accurately applied to components of pavement design and the combination of these components in an overall system analysis. The study shows that for the pavement section considered, there is a high degree of dependence between the two failure modes, and demonstrates that the probability of simultaneous occurrence of failures can be almost as high as the probability of component failures. Thus, the need to consider the system reliability in the pavement analysis is highlighted, and the study indicates that the improvement of pavement performance should be tackled in the light of reducing this undesirable event of simultaneous failure and not merely the consideration of the more critical failure mode. Furthermore, this probability of simultaneous occurrence of failures is seen to increase considerably with small increments in the mean traffic loads, which also results in wider system reliability bounds. The study also advocates the use of narrow bounds to the probability of failure, which provides a better estimate of the probability of failure, as validated from the results obtained from Monte Carlo simulation (MCS).

Influence of Spatial Variation of Hydraulic Conductivity of Municipal Solid Waste on Performance of Bioreactor Landfill

The current study analyzes the leachate distribution in the Orchard Hills Landfill, Davis Junction, Illinois, using a two-phase flow model to assess the influence of variability in hydraulic conductivity on the effectiveness of the existing leachate recirculation system and its operations through reliability analysis. Numerical modeling, using finite-difference code, is performed with due consideration to the spatial variation of hydraulic conductivity of the municipal solid waste (MSW). The inhomogeneous and anisotropic waste condition is assumed because it is a more realistic representation of the MSW. For the reliability analysis, the landfill is divided into 10 MSW layers with different mean values of vertical and horizontal hydraulic conductivities (decreasing from top to bottom), and the parametric study is performed by taking the coefficients of variation (COVs) as 50, 100, 150, and 200%. Monte Carlo simulations are performed to obtain statistical information (mean and COV) of output parameters of the (1) wetted area of the MSW, (2) maximum induced pore pressure, and (3) leachate outflow. The results of the reliability analysis are used to determine the influence of hydraulic conductivity on the effectiveness of the leachate recirculation and are discussed in the light of a deterministic approach. The study is useful in understanding the efficiency of the leachate recirculation system. (C) 2013 American Society of Civil Engineers.

Seismic Bearing Capacity of Foundations on Cohesionless Slopes

By applying the lower bound theorem of limit analysis in conjunction with finite elements and nonlinear optimization, the bearing capacity factor N has been computed for a rough strip footing by incorporating pseudostatic horizontal seismic body forces. As compared with different existing approaches, the present analysis is more rigorous, because it does not require an assumption of either the failure mechanism or the variation of the ratio of the shear to the normal stress along the footing-soil interface. The magnitude of N decreases considerably with an increase in the horizontal seismic acceleration coefficient (kh). With an increase in kh, a continuous spread in the extent of the plastic zone toward the direction of the horizontal seismic body force is noted. The results obtained from this paper have been found to compare well with the solutions reported in the literature. (C) 2013 American Society of Civil Engineers.

Dynamics of edge dislocations in a sheared lamellar mesophase

The dynamics and interactions of edge dislocations in a nearly aligned sheared lamellar mesophase is analysed to provide insights into the relationship between disorder and rheology. First, the mesoscale permeation and momentum equations for the displacement field in the presence of external forces are derived from the model H equations for the concentration and momentum field. The secondary flow generated due to the mean shear around an isolated defect is calculated, and the excess viscosity due to the presence of the defect is determined from the excess energy dissipation due to the secondary flow. The excess viscosity for an isolated defect is found to increase with system size in the cross-stream direction as L-3/2 for an isolated defect, though this divergence is cut-off due to interactions in a defect suspension. As the defects are sheared past each other due to the mean flow, the Peach-Koehler force due to elastic interaction between pairs of defects is found to cause no net displacement relative to each other as they approach from large separation to the distance of closest approach. The equivalent force due to viscous interactions is found to increase the separation for defects of opposite sign, and decrease the separation for defects of same sign. During defect interactions, we find that there is no buckling instability due to dilation of layers for systems of realistic size. However, there is another mechanism, which is the velocity difference generated across a slightly deformed bilayer due to the mean shear, which could result in the creation of new defects. (C) 2013 AIP Publishing LLC.

Layerwise decomposition of water dynamics in reverse micelles: A simulation study of two-dimensional infrared spectrum

We present computer simulation study of two-dimensional infrared spectroscopy (2D-IR) of water confined in reverse micelles (RMs) of various sizes. The present study is motivated by the need to understand the altered dynamics of confined water by performing layerwise decomposition of water, with an aim to quantify the relative contributions of different layers water molecules to the calculated 2D-IR spectrum. The 0-1 transition spectra clearly show substantial elongation, due to in-homogeneous broadening and incomplete spectral diffusion, along the diagonal in the surface water layer of different sized RMs. Fitting of the frequency fluctuation correlation functions reveal that the motion of the surface water molecules is sub-diffusive and indicate the constrained nature of their dynamics. This is further supported by two peak nature of the angular analogue of van Hove correlation function. With increasing system size, the water molecules become more diffusive in nature and spectral diffusion almost completes in the central layer of the larger size RMs. Comparisons between experiments and simulations establish the correspondence between the spectral decomposition available in experiments with the spatial decomposition available in simulations. Simulations also allow a quantitative exploration of the relative role of water, sodium ions, and sulfonate head groups in vibrational dephasing. Interestingly, the negative cross correlation between force on oxygen and hydrogen of O-H bond in bulk water significantly decreases in the surface layer of each RM. This negative cross correlation gradually increases in the central water pool with increasing RMs size and this is found to be partly responsible for the faster relaxation rate of water in the central pool. (C) 2013 AIP Publishing LLC.

Stability of fluctuating and transient aggregates of amphiphilic solutes in aqueous binary mixtures: Studies of dimethylsulfoxide, ethanol, and tert-butyl alcohol

In aqueous binary mixtures, amphiphilic solutes such as dimethylsulfoxide (DMSO), ethanol, tertbutyl alcohol (TBA), etc., are known to form aggregates (or large clusters) at small to intermediate solute concentrations. These aggregates are transient in nature. Although the system remains homogeneous on macroscopic length and time scales, the microheterogeneous aggregation may profoundly affect the properties of the mixture in several distinct ways, particularly if the survival times of the aggregates are longer than density relaxation times of the binary liquid. Here we propose a theoretical scheme to quantify the lifetime and thus the stability of these microheterogeneous clusters, and apply the scheme to calculate the same for water-ethanol, water-DMSO, and water-TBA mixtures. We show that the lifetime of these clusters can range from less than a picosecond (ps) for ethanol clusters to few tens of ps for DMSO and TBA clusters. This helps explaining the absence of a strong composition dependent anomaly in water-ethanol mixtures but the presence of the same in water-DMSO and water-TBA mixtures. (C) 2013 AIP Publishing LLC.

New insights into designing metallacarborane based room temperature hydrogen storage media

Metallacarboranes are promising towards realizing room temperature hydrogen storage media because of the presence of both transition metal and carbon atoms. In metallacarborane clusters, the transition metal adsorbs hydrogen molecules and carbon can link these clusters to form metal organic framework, which can serve as a complete storage medium. Using first principles density functional calculations, we chalk out the underlying principles of designing an efficient metallacarborane based hydrogen storage media. The storage capacity of hydrogen depends upon the number of available transition metal d-orbitals, number of carbons, and dopant atoms in the cluster. These factors control the amount of charge transfer from metal to the cluster, thereby affecting the number of adsorbed hydrogen molecules. This correlation between the charge transfer and storage capacity is general in nature, and can be applied to designing efficient hydrogen storage systems. Following this strategy, a search for the best metallacarborane was carried out in which Sc based monocarborane was found to be the most promising H-2 sorbent material with a 9 wt.% of reversible storage at ambient pressure and temperature. (C) 2013 AIP Publishing LLC.

Communication: Benzene dimer-The free energy landscape

Establishing the relative orientation of the two benzene molecules in the dimer has remained an enigmatic challenge. Consensus has narrowed the choice of structures to either a T-shape, that may be tilted, or a parallel displaced arrangement, but the relatively small energy differences makes identifying the global minimum difficult. Here we report an ab initio Car-Parrinello Molecular Dynamics based metadynamics computation of the free-energy landscape of the benzene dimer. Our calculations show that although competing structures may be isoenergetic, free energy always favors a tilted T-shape geometry at all temperatures where the bound benzene dimer exist. (C) 2013 AIP Publishing LLC.

Solid-solid collapse transition in a two dimensional model molecular system

Solid-solid collapse transition in open framework structures is ubiquitous in nature. The real difficulty in understanding detailed microscopic aspects of such transitions in molecular systems arises from the interplay between different energy and length scales involved in molecular systems, often mediated through a solvent. In this work we employ Monte-Carlo simulation to study the collapse transition in a model molecular system interacting via both isotropic as well as anisotropic interactions having different length and energy scales. The model we use is known as Mercedes-Benz (MB), which, for a specific set of parameters, sustains two solid phases: honeycomb and oblique. In order to study the temperature induced collapse transition, we start with a metastable honeycomb solid and induce transition by increasing temperature. High density oblique solid so formed has two characteristic length scales corresponding to isotropic and anisotropic parts of interaction potential. Contrary to the common belief and classical nucleation theory, interestingly, we find linear strip-like nucleating clusters having significantly different order and average coordination number than the bulk stable phase. In the early stage of growth, the cluster grows as a linear strip, followed by branched and ring-like strips. The geometry of growing cluster is a consequence of the delicate balance between two types of interactions, which enables the dominance of stabilizing energy over destabilizing surface energy. The nucleus of stable oblique phase is wetted by intermediate order particles, which minimizes the surface free energy. In the case of pressure induced transition at low temperature the collapsed state is a disordered solid. The disordered solid phase has diverse local quasi-stable structures along with oblique-solid like domains. (C) 2013 AIP Publishing LLC.

IRES mediated translational regulation of p53 isoforms

p53 is a well known tumor suppressor protein that plays a critical role in cell cycle arrest and apoptosis. It has several isoforms which are produced by transcriptional and posttranscriptional regulatory mechanisms. p53 mRNA has been demonstrated to be translated into two isoforms, full-length p53 (FL-p53) and a truncated isoform N-p53 by the use of alternative translation initiation sites. The mechanism of translation regulation of these two isoforms was further elucidated by the discovery of IRES elements in the p53 mRNA. These two IRESs were shown to regulate the translation of p53 and N-p53 in a distinct cell-cycle phase-dependent manner. This review focuses on the current understanding of the regulation of p53 IRES mediated translation and the role of cis and trans acting factors that influence expression of p53 isoforms. (C) 2013 John Wiley & Sons, Ltd.