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.

Dimethyl sulfoxide induced structural transformations and non-monotonic concentration dependence of conformational fluctuation around active site of lysozyme

Experimental studies have observed significant changes in both structure and function of lysozyme (and other proteins) on addition of a small amount of dimethyl sulfoxide (DMSO) in aqueous solution. Our atomistic molecular dynamic simulations of lysozyme in water-DMSO reveal the following sequence of changes on increasing DMSO concentration. (i) At the initial stage (around 5% DMSO concentration) protein's conformational flexibility gets markedly suppressed. From study of radial distribution functions, we attribute this to the preferential solvation of exposed protein hydrophobic residues by the methyl groups of DMSO. (ii) In the next stage (10-15% DMSO concentration range), lysozome partially unfolds accompanied by an increase both in fluctuation and in exposed protein surface area. (iii) Between 15-20% concentration ranges, both conformational fluctuation and solvent accessible protein surface area suddenly decrease again indicating the formation of an intermediate collapse state. These results are in good agreement with near-UV circular dichroism (CD) and fluorescence studies. We explain this apparently surprising behavior in terms of a structural transformation which involves clustering among the methyl groups of DMSO. (iv) Beyond 20% concentration of DMSO, the protein starts its final sojourn towards the unfolding state with further increase in conformational fluctuation and loss in native contacts. Most importantly, analysis of contact map and fluctuation near the active site reveal that both partial unfolding and conformational fluctuations are centered mostly on the hydrophobic core of active site of lysozyme. Our results could offer a general explanation and universal picture of the anomalous behavior of protein structure-function observed in the presence of cosolvents (DMSO, ethanol, tertiary butyl alcohol, dioxane) at their low concentrations. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3694268]

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.

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. © 2012 Wiley Periodicals, Inc.

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.

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. (c) 2012 Wiley Periodicals, Inc.

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]

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.

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.