Influence of a Counterion on the Ion Atmosphere of an Anion: A Molecular Dynamics Study of LiX and CsX (X = F-, Cl-, I-) in Methanol

We report molecular dynamics (MD) simulations to explore the influence of a counterion on the structure and dynamics of cationic and anionic solvation shells for various ions in methanol at 298 K. We show that the variation in ionic size of either the cation or the anion in an ion pair influences the solvation structure of the other ion as well as the diffusivity in an electrolyte solution of methanol. The extent of ionic association between the cation and its counteranion of different ionic sizes has been investigated by analyzing the radial distribution functions (RDFs) and the orientation of methanol molecules in the first solvation shell (FSS) of ions. It is shown that the methanol in the FSS of the anion as well the cation exhibit quite different radial and orientational structures as compared to methanol which lie in the FSS of either the anion or the cation but not both. We find that the coordination number (CN) of F-, Cr-, and I- ions decreases with increasing size of the anion which is contrary to the trend reported for the anions in H2O. The mean residence time (MRT) of methanol molecules in the FSS of ions has been calculated using the stable states picture (SSP) approach. It is seen that the ion-counterion interaction has a considerable influence on the MRT of methanol molecules in the FSS of ions. We also discuss the stability order of the ion-counterion using the potentials of mean force (PMFs) for ion pairs with ions of different sizes. The PMF plots reveal that the Li+-F- pair (small-small) is highly stable and the Li+-I- pair is least stable (small-large) in electrolyte solutions.

Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pK(a) Values from Metadynamics Sampling

Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation deprotonation reaction of the 20 canonical alpha amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metad-ynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pK(a) values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pK(a) values with a mean relative error, with respect to experimental results, of 0.2 pK(a) units.

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

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

Molecular structure of the discotic liquid crystalline phase of hexa-peri-hexabenzocoronene/oligothiophene hybrid and their charge transport properties

Using atomistic molecular dynamics simulation, we study the discotic columnar liquid crystalline (LC) phases formed by a new organic compound having hexa-peri-Hexabenzocoronene (HBC) core with six pendant oligothiophene units recently synthesized by Nan Hu et al. Adv. Mater. 26, 2066 (2014)]. This HBC core based LC phase was shown to have electric field responsive behavior and has important applications in organic electronics. Our simulation results confirm the hexagonal arrangement of columnar LC phase with a lattice spacing consistent with that obtained from small angle X-ray diffraction data. We have also calculated various positional and orientational correlation functions to characterize the ordering of the molecules in the columnar arrangement. The molecules in a column are arranged with an average twist of 25 degrees having an average inter-molecular separation of similar to 5 angstrom. Interestingly, we find an overall tilt angle of 43 degrees between the columnar axis and HBC core. We also simulate the charge transport through this columnar phase and report the numerical value of charge carrier mobility for this liquid crystal phase. The charge carrier mobility is strongly influenced by the twist angle and average spacing of the molecules in the column. (C) 2015 AIP Publishing LLC.

The putative role of some conserved water molecules in the structure and function of human transthyretin

Human transthyretin (hTTR) is a multifunctional protein that is involved in several neurodegenerative diseases. Besides the transportation of thyroxin and vitamin A, it is also involved in the proteolysis of apolipoprotein A1 and A beta peptide. Extensive analyses of 32 high-resolution X-ray and neutron diffraction structures of hTTR followed by molecular-dynamics simulation studies using a set of 15 selected structures affirmed the presence of 44 conserved water molecules in its dimeric structure. They are found to play several important roles in the structure and function of the protein. Eight water molecules stabilize the dimeric structure through an extensive hydrogen-bonding network. The absence of some of these water molecules in highly acidic conditions (pH <= 4.0) severely affects the interfacial hydrogen-bond network, which may destabilize the native tetrameric structure, leading to its dissociation. Three pairs of conserved water molecules contribute to maintaining the geometry of the ligand-binding cavities. Some other water molecules control the orientation and dynamics of different structural elements of hTTR. This systematic study of the location, absence, networking and interactions of the conserved water molecules may shed some light on various structural and functional aspects of the protein. The present study may also provide some rational clues about the conserved water-mediated architecture and stability of hTTR.

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

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

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (L-alpha) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the L-alpha phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the L-alpha phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.

Potency and pharmacokinetics of broad spectrum and isoform-specific p110 gamma and delta inhibitors in cancers

Emerging data on cancer suggesting that target-based therapy is promising strategy in cancer treatment. PI3K-AKT pathway is extensively studied in many cancers; several inhibitors target this pathway in different levels. Recent finding on this pathway uncovered the therapeutic applications of PI3K-specific inhibitors; PI3K, AKT, and mTORC broad spectrum inhibitors. Noticeably, class I PI3K isoforms, p110 and p110 catalytic subunits have rational therapeutic application than other isoforms. Therefore, three classes of inhibitors: isoform-specific, dual-specific and broad spectrum were selected for molecular docking and dynamics. First, p110 structure was modelled; active site was analyzed. Then, molecular docking of each class of inhibitors were studied; the docked complexes were further used in 1.2ns molecular dynamics simulation to report the potency of each class of inhibitor. Remarkably, both the studies retained the similar kind of protein ligand interactions. GDC-0941, XL-147 (broad spectrum); TG100-115 (dual-specific); and AS-252424, PIK-294 (isoform-specific) were found to be potential inhibitors of p110 and p110, respectively. In addition to that pharmacokinetic properties are within recommended ranges. Finally, molecular phylogeny revealed that p110 and p110 are evolutionarily divergent; they probably need separate strategies for drug development.

Molecular dynamics and simulations study on the vibrational and electronic solvatochromism of benzophenone

Solvent plays a key role in diverse physico-chemical and biological processes. Therefore, understanding solute-solvent interactions at the molecular level of detail is of utmost importance. A comprehensive solvatochromic analysis of benzophenone (Bzp) was carried out in various solvents using Raman and electronic spectroscopy, in conjunction with Density Functional Theory (DFT) calculations of supramolecular solute-solvent clusters generated using classical Molecular Dynamics Simulations (c-MDSs). The >C=O stretching frequency undergoes a bathochromic shift with solvent polarity. Interestingly, in protic solvents this peak appears as a doublet: c-MDS and ad hoc explicit solvent ab initio calculations suggest that the lower and higher frequency peaks are associated with the hydrogen bonded and dangling carbonyl group of Bzp, respectively. Additionally, the dangling carbonyl in methanol (MeOH) solvent is 4 cm(-1) blue-shifted relative to acetonitrile solvent, despite their similar dipolarity/polarizability. This suggests that the cybotactic region of the dangling carbonyl group in MeOH is very different from its bulk solvent structure. Therefore, we propose that this blue-shift of the dangling carbonyl originates in the hydrophobic solvation shell around it resulting from extended hydrogen bonding network of the protic solvents. Furthermore, the 1(1)n pi* (band I) and 1(1)pi pi* (band II) electronic transitions show a hypsochromic and bathochromic shift, respectively. In particular, these shifts in protic solvents are due to differences in their excited state-hydrogen bonding mechanisms. Additionally, a linear relationship is obtained for band I and the >C=O stretching frequency (cm(-1)), which suggests that the different excitation wavelengths in band I correspond to different solvation states. Therefore, we hypothesize that the variation in excitation wavelengths in band I could arise from different solvation states leading to varying solvation dynamics. This will have implications for ultrafast processes associated with electron-transfer, charge transfer, and also the photophysical aspects of excited states. (C) 2016 AIP Publishing LLC.

Diffusion of flexible, charged, nanoscopic molecules in solution: Size and pH dependence for PAMAM dendrimer

In order to understand self-diffusion (D) of a charged, flexible, and porous nanoscopic molecule in water, we carry out very long, fully atomistic molecular dynamics simulation of PAMAM dendrimer up to eight generations in explicit salt water under varying pH. We find that while the radius of gyration (R-g) varies as N-1/3, the self-diffusion constant (D) scales, surprisingly, as N-alpha, with alpha=0.39 at high pH and 0.5 at neutral pH, indicating a dramatic breakdown of Stokes-Einstein relation for diffusion of charged nanoscopic molecules. The variation in D as a function of radius of gyration demonstrates the importance of treating water and ions explicitly in the diffusion process of a flexible nanoscopic molecule. In agreement with recent experiments, the self-diffusion constant increases with pH, revealing the importance of dielectric friction in the diffusion process. The shape of a dendrimer is found to fluctuate on a nanosecond time scale. We argue that this flexibility (and also the porosity) of the dendrimer may play an important role in determining the mean square displacement of the dendrimer and the breakdown of the Stokes-Einstein relation between diffusion constant and the radius.

Correlation between conductivity or diffusivity and activation energy in amorphous solids

There exist many investigations of ionic transport in a variety of glasses. These studies exhibit strong correlation between ionic conductivity and activation energy: Typically, it is found that higher conductivity is associated with lower activation energies and vice versa. Although there are explanations for this at a phenomenological level, there is no consistent physical picture to explain the correlation between conductivity and activation energy. We have carried out molecular dynamics simulation as a function of the size of the impurity atom or diffusant (both neutral and charged) in a host amorphous matrix. We find that there is a maximum in self-diffusivity as a function of the size of the impurity atom suggesting that there is an appropriate size for which the diffusivity is maximum. The activation energy is found to be the lowest for this size of the impurity. A similar maximum has been previously found in other condensed phases, such as confined fluids and dense liquids, and has its origin in the levitation effect. The implications of this result for understanding ionic conductivity in glasses are discussed. Our results suggest that there is a relation between microscopic structure of the amorphous solid, diffusivity or conductivity, and activation energy. The nature of this relationship is discussed in terms of the levitation parameter showing that diffusivity is maximum when the size of the neck or doorway radius is comparable with the size of the diffusant. Our computational results here are in excellent agreement with independent experimental results of Nascimento et al. [Braz. J. Phys. 35, 626 (2005)] that structural features of the glass are important in determining the ionic conductivity.

Structural Polymorphism in d(T).d(A)*d(T) Triple Helices

DNA triple helices containing two thymine strands and one adenine strand have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that the structure with a parallel orientation of the third strand, in which the third strand base forms Hoogsteen hydrogen bonds with the adenine base in the Watson-Crick duplex, is energetically the most favourable. An antiparallel orientation of the third strand is also possible, in which the third strand base hydrogen bonds to both the bases in the Watson-Crick duplex. This structure is energetically comparable to the parallel structure. For the parallel triplex a 200ps molecular dynamics simulation starting from two different starting structures indicates that at 300K significant structural heterogeneity exists in this tripler structure. The results are compared with existing structural data on this class of triplexes derived from theoretical and NMR techniques.

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.

Strain induced phase transition in CdSe nanowires: Effect of size and temperature

Using all-atom molecular dynamics simulation, we have studied the effect of size and temperature on the strain induced phase transition of wurtzite CdSe nanowires. The wurtzite structure transforms into a five-fold coordinated structure under uniaxial strain along the c axis. Our results show that lower temperature and smaller size of the nanowires stabilize the five-fold coordinated phase which is not a stable structure in bulk CdSe. High reversibility of this transformation with a very small heat loss will make these nanowires suitable for building efficient nanodevices. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4734990]

Effect of pressure on the ionic conductivity of Li+ and Cl- ions in water

A molecular dynamics simulation study of aqueous solution of LiCl is reported as a function of pressure. Experimental measurements of conductivity of Li+ ion as a function of pressure shows an increase in conductivity with pressure. Our simulations are able to reproduce the observed trend in conductivity. A number of relevant properties have been computed in order to understand the reasons for the increase in conductivity with pressure. These include radial distribution function, void and neck distributions, hydration or coordination numbers, diffusivity, velocity autocorrelation functions, angles between ion-oxygen and dipole of water as well as OH vector, mean residence time for water in the hydration shell, etc. These show that the increase in pressure acts as a structure breaker. The decay of the self part of the intermediate scattering function at small wave number k shows a bi-exponential decay at 1 bar which changes to single exponential decay at higher pressures. The k dependence of the ratio of the self part of the full width at half maximum of the dynamic structure factor to 2Dk(2) exhibits trends which suggest that the void structure of water is playing a role. These support the view that the changes in void and neck distributions in water can account for changes in conductivity or diffusivity of Li+ with pressure. These results can be understood in terms of the levitation effect. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4756909]