Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.

Mutations in Drosophila Myosin Rod Cause Defects in Myofibril Assembly

The roles of myosin during muscle contraction are well studied, but how different domains of this protein are involved in myofibril assembly in vivo is far less understood. The indirect flight muscles (IFMs) of Drosophila melanogaster provide a good model for understanding muscle development and function in vivo. We show that two missense mutations in the rod region of the myosin heavy-chain gene, Mhc, give rise to IFM defects and abnormal myofibrils. These defects likely result from thick filament abnormalities that manifest during early sarcomere development or later by hypercontraction. The thick filament defects are accompanied by marked reduction in accumulation of flightin, a myosin binding protein, and its phosphorylated forms, which are required to stabilise thick filaments. We investigated with purified rod fragments whether the mutations affect the coiled-coil structure, rod aggregate size or rod stability. No significant changes in these parameters were detected, except for rod thermodynamic stability in one mutation. Molecular dynamics simulations suggest that these mutations may produce localised rod instabilities. We conclude that the aberrant myofibrils are a result of thick filament defects, but that these in vivo effects cannot be detected in vitro using the biophysical techniques employed. The in vivo investigation of these mutant phenotypes in IFM development and function provides a useful platform for studying myosin rod and thick filament formation generically, with application to the aetiology of human myosin rod myopathies. (C) 2012 Elsevier Ltd. All rights reserved.

Transport in nanoporous zeolites: Relationships between sorbate size, entropy, and diffusivity

Molecular dynamics simulations have been performed on monatomic sorbates confined within zeolite NaY to obtain the dependence of entropy and self-diffusivity on the sorbate diameter. Previously, molecular dynamics simulations by Santikary and Yashonath J. Phys. Chem. 98, 6368 (1994)], theoretical analysis by Derouane J. Catal. 110, 58 (1988)] as well as experiments by Kemball Adv. Catal. 2, 233 (1950)] found that certain sorbates in certain adsorbents exhibit unusually high self-diffusivity. Experiments showed that the loss of entropy for certain sorbates in specific adsorbents was minimum. Kemball suggested that such sorbates will have high self-diffusivity in these adsorbents. Entropy of the adsorbed phase has been evaluated from the trajectory information by two alternative methods: two-phase and multiparticle expansion. The results show that anomalous maximum in entropy is also seen as a function of the sorbate diameter. Further, the experimental observation of Kemball that minimum loss of entropy is associated with maximum in self-diffusivity is found to be true for the system studied here. A suitably scaled dimensionless self-diffusivity shows an exponential dependence on the excess entropy of the adsorbed phase, analogous to excess entropy scaling rules seen in many bulk and confined fluids. The two trajectory-based estimators for the entropy show good semiquantitative agreement and provide some interesting microscopic insights into entropy changes associated with confinement.

Mechanical properties of ZnS nanowires and thin films: Microscopic origin of the dependence on size and growth direction

Mechanical properties of ZnS nanowires and thin films are studied as a function of size and growth direction using all-atom molecular dynamics simulations. Using the stress-strain relationship we extract Young's moduli of nanowires and thin films at room temperature. Our results show that Young's modulus of 0001] nanowires has strong size dependence. On the other hand, 01 (1) over bar0] nanowires do not exhibit a strong size dependence of Young's modulus in the size range we have investigated. We provide a microscopic understanding of this behavior on the basis of bond stretching and contraction due to the rearrangement of atoms in the surface layers. The ultimate tensile strengths of the nanowires do not show much size dependence. To investigate the mechanical behavior of ZnS in two dimensions, we calculate Young's modulus of thin films under tensile strain along the 0001] direction. Young's modulus of thin films converges to the bulk value more rapidly than that of the 0001] nanowire.

Accommodating Unwelcome Guests in Inorganic Layered Hosts: Inclusion of Chloranil in a Layered Double Hydroxide

The host-guest chemistry of most inorganic layered solids is limited to ion-exchange reactions. The guest species are either cations or anions to compensate for the charge deficit, either positive or negative, of the inorganic layers. Here, we outline a strategy to include neutral molecules like ortho- and para-chloranil, that are known to be good acceptors in donor-acceptor or charge-transfer complexes, within the galleries of a layered solid. We have succeeded in including neutral ortho- and para-chloranil molecules within the galleries of an Mg-Al layered double hydroxide (LDH) by using charge-transfer interactions with preintercalated p-aminobenzoate ions as the driving force. The p-aminobenzoate ions are introduced in the Mg-Al LDH via ion exchange. The intercalated LDH can adsorb ortho- and para-chloranil from chloroform solutions by forming charge-transfer complexes with the p-aminobenzoate anions present in the galleries. We use X-ray diffraction, spectroscopy, and molecular dynamics simulations to establish the nature of interactions and arrangement of the charge-transfer complex within the galleries of the layered double hydroxide.

Unraveling siRNA unzipping kinetics with graphene

Using all atom molecular dynamics simulations, we report spontaneous unzipping and strong binding of small interfering RNA (siRNA) on graphene. Our dispersion corrected density functional theory based calculations suggest that nucleosides of RNA have stronger attractive interactions with graphene as compared to DNA residues. These stronger interactions force the double stranded siRNA to spontaneously unzip and bind to the graphene surface. Unzipping always nucleates at one end of the siRNA and propagates to the other end after few base-pairs get unzipped. While both the ends get unzipped, the middle part remains in double stranded form because of torsional constraint. Unzipping probability distributions fitted to single exponential function give unzipping time (tau) of the order of few nanoseconds which decrease exponentially with temperature. From the temperature variation of unzipping time we estimate the energy barrier to unzipping. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4742189]

Translocation and encapsulation of siRNA inside carbon nanotubes

We report spontaneous translocation of small interfering RNA (siRNA) inside carbon nanotubes (CNTs) of various diameters and chirality using all atom molecular dynamics simulations with explicit solvent. We use umbrella sampling method to calculate the free energy landscape of the siRNA entry and translocation event. Free energy profiles show that siRNA gains free energy while translocating inside CNT, and barrier for siRNA exit from CNT ranges from 40 to 110 kcal/mol depending on CNT chirality and salt concentration. The translocation time tau decreases with the increase of CNT diameter with a critical diameter of 24 angstrom for the translocation. In contrast, double strand DNA of the same sequence does not translocate inside CNT due to large free energy barrier for the translocation. This study helps in understanding the nucleic acid transport through nanopores at microscopic level and may help designing carbon nanotube based sensor for siRNA. (C) 2013 American Institute of Physics. http://dx.doi.org/10.1063/1.4773302]

Application of a modified jogged-screw model for creep of titanium aluminides: evaluation of the key substructural parameters

A modification of the jogged-screw model has been adopted recently by the authors to explain observations of 1/2[110]-type jogged-screw dislocations in equiaxed Ti-48Al under creep conditions. The aim of this study has been to verify and validate the parameters and functional dependencies that have been assumed in this previous work. The original solution has been reformulated to take into account the finite length of the moving jog. This is a better approximation of the tall jog. The substructural model parameters have been further investigated in light of the Finite Length Moving Line (FLML) source approximation. The original model assumes that the critical jog height (beyond which the jog is not dragged) is inversely proportional to the applied stress. By accounting for the fact that there are three competing mechanisms (jog dragging, dipole dragging, dipole bypass) possible, we can arrive at a modified critical jog height. The critical jog height was found to be more strongly stress dependent than assumed previously. The original model assumes the jog spacing to be invariant over the stress range. However, dynamic simulation using a line tension model has shown that the jog spacing is inversely proportional to the applied stress. This has also been confirmed by TEM measurements of jog spacings over a range of stresses. Taylor's expression assumed previously to provide the dependence of dislocation density on the applied stress, has now been confirmed by actual dislocation density measurements. Combining all of these parameters and dependencies, derived both from experiment and theory, leads to an excellent prediction of creep rates and stress exponents. The further application of this model to other materials, and the important role of atomistic and dislocation dynamics simulations in its continued development is also discussed.

Chemical unfolding of chicken villin headpiece in aqueous dimethyl sulfoxide solution: cosolvent concentration dependence, pathway, and microscopic mechanism

Unfolding of a protein often proceeds through partial unfolded intermediate states (PUIS). PUIS have been detected in several experimental and simulation studies. However, complete analyses of transitions between different PUIS and the unfolding trajectory are sparse. To understand such dynamical processes, we study chemical unfolding of a small protein, chicken villin head piece (HP-36), in aqueous dimethyl sulfoxide (DMSO) solution. We carry out molecular dynamics simulations at various solution compositions under ambient conditions. In each concentration, the initial step of unfolding involves separation of two adjacent native contacts, between phenyl alanine residues (11-18 and 7-18). This first step induces, under appropriate conditions, subsequent separation among other hydrophobic contacts, signifying a high degree of cooperativity in the unfolding process. The observed sequence of structural changes in HP-36 on increasing DMSO concentration and the observed sequence of PUIS, are in approximate agreement with earlier simulation results (in pure water) and experimental observations on unfolding of HP-36. Peculiar to water-DMSO mixture, an intervening structural transformation (around 15% of DMSO) in the binary mixture solvent retards the progression of unfolding as composition is increased. This is reflected in a remarkable nonmonotonic composition dependence of RMSD, radius of gyration and the fraction of native contacts. At 30% mole fraction of DMSO, we find the extended randomly coiled structure of the unfolded protein. The molecular mechanism of DMSO induced unfolding process is attributed to the initial preferential solvation of the hydrophobic side chain atoms through the methyl groups of DMSO, followed by the hydrogen bonding of the oxygen atom of DMSO to the exposed backbone NH groups of HP-36.

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.

Ionomer based blend as water vapor barrier material for organic device encapsulation

Blends of poly (ethylene-co-methacrylic acid) (PEMA) and poly (vinyl alcohol-co-ethylene) (EVOH) were studied for encapsulating Schottky structured organic devices. A calcium degradation test was used to determine water vapor transmission rates and to determine the moisture barrier performance of neat and blend films. Moisture barrier analysis for the neat and blend compositions was discussed concerning the interactions in the blend, diffusivity of water molecules through the unit cell systems, and the occupiable free volumes available in the unit cells using molecular dynamics simulations. The experimental results of water vapor permeation were correlated with diffusion behavior predicted from molecular dynamics simulations results. The effectiveness of the blend as a suitable barrier material in increasing the lifetime of an encapsulated Schottky structured organic device was determined.

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.

Mechanistic features of Salmonella typhimurium propionate kinase (TdcD): Insights from kinetic and crystallographic studies

Short-chain fatty acids (SCFAs) play a major role in carbon cycle and can be utilized as a source of carbon and energy by bacteria. Salmonella typhimurium propionate kinase (StTdcD) catalyzes reversible transfer of the gamma-phosphate of ATP to propionate during L-threonine degradation to propionate. Kinetic analysis revealed that StTdcD possesses broad ligand specificity and could be activated by various SCFAs (propionate > acetate approximate to butyrate), nucleotides (ATP approximate to GTP > CTP approximate to TTP; dATP > dGTP > dCTP) and metal ions (Mg2+ approximate to Mn2+ > Co2+). Inhibition of StTdcD by tricarboxylic acid (TCA) cycle intermediates such as citrate, succinate, alpha-ketoglutarate and malate suggests that the enzyme could be under plausible feedback regulation. Crystal structures of StTdcD bound to PO4 (phosphate), AMP, ATP, Ap4 (adenosine tetraphosphate), GMP, GDP, GTP, CMP and CTP revealed that binding of nucleotide mainly involves hydrophobic interactions with the base moiety and could account for the broad biochemical specificity observed between the enzyme and nucleotides. Modeling and site-directed mutagenesis studies suggest Ala88 to be an important residue involved in determining the rate of catalysis with SCFA substrates. Molecular dynamics simulations on monomeric and dimeric forms of StTdcD revealed plausible open and closed states, and also suggested role for dimerization in stabilizing segment 235-290 involved in interfacial interactions and ligand binding. Observation of an ethylene glycol molecule bound sufficiently close to the gamma-phosphate in StTdcD complexes with triphosphate nucleotides supports direct in-line phosphoryl transfer. (C) 2013 Elsevier B.V. All rights reserved.

Conformational Diversity in Contryphans from Conus Venom: cis-trans Isomerisation and Aromatic/Proline Interactions in the 23-Membered Ring of a 7-Residue Peptide Disulfide Loop

Conformational diversity or shapeshifting in cyclic peptide natural products can, in principle, confer a single molecular entity with the property of binding to multiple receptors. Conformational equilibria have been probed in the contryphans, which are peptides derived from Conus venom possessing a 23-membered cyclic disulfide moiety. The natural sequences derived from Conus inscriptus, GCV(D)LYPWC* (In936) and Conus loroisii, GCP(D)WDPWC* (Lo959) differ in the number of proline residues within the macrocyclic ring. Structural characterisation of distinct conformational states arising from cis-trans equilibria about Xxx-Pro bonds is reported. Isomerisation about the C2-P3 bond is observed in the case of Lo959 and about the Y5-P6 bond in In936. Evidence is presented for as many as four distinct species in the case of the synthetic analogue V3P In936. The Tyr-Pro-Trp segment in In936 is characterised by distinct sidechain orientations as a consequence of aromatic/proline interactions as evidenced by specific sidechain-sidechain nuclear Overhauser effects and ring current shifted proton chemical shifts. Molecular dynamics simulations suggest that Tyr5 and Trp7 sidechain conformations are correlated and depend on the geometry of the Xxx-Pro bond. Thermodynamic parameters are derived for the cis trans equilibrium for In936. Studies on synthetic analogues provide insights into the role of sequence effects in modulating isomerisation about Xxx-Pro bonds.

Increasing cleavage specificity and activity of restriction endonuclease KpnI

Restriction enzyme KpnI is a HNH superfamily endonuclease requiring divalent metal ions for DNA cleavage but not for binding. The active site of KpnI can accommodate metal ions of different atomic radii for DNA cleavage. Although Mg2+ ion higher than 500 mu M mediates promiscuous activity, Ca2+ suppresses the promiscuity and induces high cleavage fidelity. Here, we report that a conservative mutation of the metal-coordinating residue D148 to Glu results in the elimination of the Ca2+-mediated cleavage but imparting high cleavage fidelity with Mg2+. High cleavage fidelity of the mutant D148E is achieved through better discrimination of the target site at the binding and cleavage steps. Biochemical experiments and molecular dynamics simulations suggest that the mutation inhibits Ca2+-mediated cleavage activity by altering the geometry of the Ca2+-bound HNH active site. Although the D148E mutant reduces the specific activity of the enzyme, we identified a suppressor mutation that increases the turnover rate to restore the specific activity of the high fidelity mutant to the wild-type level. Our results show that active site plasticity in coordinating different metal ions is related to KpnI promiscuous activity, and tinkering the metal ion coordination is a plausible way to reduce promiscuous activity of metalloenzymes.