Dependence of structure of dehydrated PVA on the crystallinity of the original PVA

Poly(vinyl alcohol) (PVA) films having different crystallinities were prepared by elongating PVA films to different degrees and heating the films with and without elongation treatment at several temperatures between 60 and 200°C. Then, they were dehydrated by heating from 80°C to 230 or 330°C in hydrogen chloride gas. Infrared spectral measurements were made on the dehydrated PVA films obtained. Absorbances of the absorption bands due to several groups seen in the infrared spectra depended only on the density of original PVA films, regardless of its degree of elongation. From these dependences, it was found that the dehydration reaction progressed more easily in the crystalline region than in the amorphous region, that the aromatic rings with four or five adjacent hydrogen atoms were formed mainly in the amorphous region and those with two adjacent hydrogen atoms in the crystalline region, and that the aromatic rings with two isolated hydrogen atoms were formed in both the amorphous and crystalline regions. Relative concentrations of the groups of which the dehydrated PVA were composed were estimated. Comparison of the numbers of carbon atoms among the composing groups indicated that the main groups were the methylene group and the aromatic ring.     

Understanding the recognition mechanism of protein-RNA complexes using energy based approach

Protein-RNA interactions perform diverse functions within the cell. Understanding the recognition mechanism of protein-RNA complexes is a challenging task in molecular and computational biology. In this work, we have developed an energy based approach for identifying the binding sites and important residues for binding in protein-RNA complexes. The new approach considers the repulsive interactions as well as the effect of distance between the atoms in protein and RNA in terms of interaction energy, which are not considered in traditional distance based methods to identify the binding sites. We found that the positively charged, polar and aromatic residues are important for binding. These residues influence to form electrostatic, hydrogen bonding and stacking interactions. Our observation has been verified with the experimental binding specificity of protein-RNA complexes and found good agreement with experiments. Further, the propensities of residues/nucleotides in the binding sites of proteins/RNA and their atomic contributions have been derived. Based on these results we have proposed a novel mechanism for the recognition of protein-RNA complexes: the charged and polar residues in proteins initiate recognition with RNA by making electrostatic and hydrogen bonding interactions between them; the aromatic side chains tend to form aromatic-aromatic interactions and the hydrophobic residues aid to stabilize the complex.     

A first principles molecular dynamics study of lithium atom solvation in binary liquid mixture of water and ammonia: structural, electronic, and dynamical properties

The preferential solvation of solutes in mixed solvent systems is an interesting phenomenon that plays important roles in solubility and kinetics. In the present study, solvation of a lithium atom in aqueous ammonia solution has been investigated from first principles molecular dynamics simulations. Solvation of alkali metal atoms, like lithium, in aqueous and ammonia media is particularly interesting because the alkali metal atoms release their valence electrons in these media so as to produce solvated electrons and metal counterions. In the present work, first principles simulations are performed employing the Car-Parrinello molecular dynamics method. Spontaneous ionization of the Li atom is found to occur in the mixed solvent system. From the radial distribution functions, it is found that the Li(+) ion is preferentially solvated by water and the coordination number is mostly four in its first solvation shell and exchange of water molecules between the first and second solvation shells is essentially negligible in the time scale of our simulations. The Li(+) ion and the unbound electron are well separated and screened by the polar solvent molecules. Also the unbound electron is primarily captured by the hydrogens of water molecules. The diffusion rates of Li(+) ion and water molecules in its first solvation shell are found to be rather slow. In the bulk phase, the diffusion of water is found to be slower than that of ammonia molecules because of strong ammonia-water hydrogen bonds that participate in solvating ammonia molecules in the mixture. The ratio of first and second rank orientational correlation functions deviate from 3, which suggests a deviation from the ideal Debye-type orientational diffusion. It is found that the hydrogen bond lifetimes of ammonia-ammonia pairs is very short. However, ammonia-water H-bonds are found to be quite strong when ammonia acts as an acceptor and these hydrogen bonds are found to live longer than even water-water hydrogen bonds.     

Mechanism of the formation of a dehydrated ion by an unusual loss of oxygen at the 4'-carbonyl group of abscisic acid methyl ester in electron ionization mass spectrometry

Methyl ester of abscisic acid (ABA), a plant hormone, gives a dehydrated ion at m/z 260 in electron ionization mass spectrometry (EI-MS). This dehydrated ion had been considered to be derived only from the elimination of the tertiary hydroxyl group at C-1'. We found that 34% of the dehydrated ion was formed by elimination of the oxygen atom at the 4'-carbonyl group, and the remaining 66% by elimination of the 1'-hydroxyl group. This unusual elimination of the carbonyl oxygen was shown with [4'-(18)O]ABA methyl ester. Involvement of the 4'-carbonyl oxygen in dehydration was observed in methyl ester of phaseic acid (PA), a natural metabolite of ABA, but not in 1'-deoxy-ABA methyl ester or isophorone. This suggested that the 1'-hydroxyl group was necessary for the elimination of the 4'-carbonyl oxygen. ABA methyl esters labeled with stable isotopes showed that hydrogen atoms at the 1'-hydroxyl group and at C-4 or -5 or -3' or - 5' or -7' were eliminated with the 4'-carbonyl oxygen. These results allow us to propose a formation mechanism of the dehydrated ion derived from the elimination of 4'-carbonyl oxygen and hydrogen atoms at C-4 and 1'-oxygen in ABA methyl ester as follows: first, ionization at the 1'-hydroxyl group occurs to give an ion radical, and the proton at the 1'-oxygen migrates to the 4'-carbonyl oxygen after the bond fission between C-1'-C-6'; second, migration of the proton at C-4 to the 1'-oxygen is followed by migration of the protons at C-5 and C-7' to C-4 and C-5, respectively; finally, the proton at the 1'-oxygen migrates to the 4'-hydroxyl group, and H(2)O at C-4' is eliminated to give the dehydrated ion. Our findings point out that a dehydrated ion is not always derived from the elimination of a hydroxyl group.     

Langmuir-Blodgett film characteristics and phospholipid membrane ion conduction : Part 2. Ethylenic Acyl Chain Oxidation

The influence of polar moieties located in the non-polar hydrocarbon zone in bilayer lipid membranes on ion conduction is described. Natural egg-derived phosphatidyl choline combined with cholesterol produces membranes containing several ethylenic residues in the hydrocarbon interior. Catalytic oxidation of these residues by ultraviolet radiation provides a significant density of permanently-bound polar species within the membrane. The extent of such oxidation is correlated with the Arrhenius energy barrier to ion conduction for bilayer membranes and molecular packing characteristics obtained from Langmuir-Blodgett monolayer compression experiments. The results confirm that membrane ion current is almost entirely controlled by molecular packing.     

Is the mobility of the pore walls and water molecules in the selectivity filter of KcsA channel functionally important?

We performed in-depth analysis of the forces which act on the K(+) ions in the selectivity filter of the KcsA channel in order to estimate the relative importance of static and dynamic influence of the filter wall and water molecules on ion permeation and selectivity. The forces were computed using the trajectories of all-atom molecular dynamics simulations. It is shown that the dynamics of the selectivity filter contributes about 3% to the net force acting on the ions and can be neglected in the studies focused on the macroscopic properties of the channel, such as the current. Among the filter atoms, only the pore-forming carbonyl groups can be considered as dynamic in the studies of microscopic events of conduction, while the dynamic effects from all other atoms are negligible. We also show that the dynamics of the water molecules in the filter can not be neglected. The fluctuating forces from the water molecules can be as strong as net forces from the pore walls and can effectively drive the ions through the local energy barriers in the filter.     

A molecular dynamics study of chloride binding by the cryptand SC24

The capture of chloride from water by the tetraprotonated form of the spherical macrotricyclic molecule SC24 was studied using molecular dynamics simulation methods. This model ionophore represents a broad class of molecules which remove ions from water. Two binding sites for the chloride were found, one inside and one outside the ligand. These sites are separated by a potential energy barrier of approximately 20 kcal mol⁻¹. The major contribution to this barrier comes from dehydration of the chloride. The large, unfavorable dehydration effect is compensated for by an increase in electrostatic attraction between the oppositely charged chloride and cryptand, and by energetically favorable rearrangements of water structure. Additional assistance in crossing the barrier and completing the dehydration of the ion is provided by the shift of three positively charged hydrogen atoms of the cryptand towards the chloride. This structural flexibility of the cryptand, leads to a decrease in the energy barrier, whereas, its structural rigidity is partially responsible for its selectivity.     

An extraction system for low-energy hydrogen ions formed by electron impact

A system has been developed for extracting near-zero kinetic energy H⁻ and D⁻ ions formed by dissociative electron attachment. It is the essential part of a new set-up for vibrational spectroscopy of hydrogen molecules. A magnetic field is used to collimate the probing electron beam. Ions produced by electron collision with the target molecules are collected by the combined action of this field and an electrostatic field penetrating into the interaction region. Highly effective extraction is achieved by taking into account the correct out-of plane displacement of ion trajectories which is usually neglected in similar arrangements. The extraction conditions are mass dependent so that by proper tuning, mass selection of detected ions is achieved. The new system is also used for detecting positive ions created by electron collisions with hydrogen atoms and molecules.     

Molecular-dynamics simulation of the effect of ions on a liquid-liquid interface for a partially miscible mixture

Molecular-dynamics simulations were performed to model the effect of added salt ions on the liquid-liquid interface in a partially miscible system. Simulations of the interface between saturated phases of a model 1-hexanol+water system show a bilayer structure of 1-hexanol molecules at the interface with -OH heads of the first layer directed into the water phase and the opposite orientation for the second layer. The alignment of the polar -OH groups at the interface stabilizes a charge separation of sodium and chloride ions when salt is introduced into the aqueous phase, producing an electrical double layer. Chloride ions aggregate nearer the interface and sodium ions move toward the bulk water phase, consistent with the explanation that the -OH alignment presents a region of partial positive charges to which the hydrated chloride atoms are attracted. Ions near the interface were found to be less solvated than those in the bulk phase. An electric field was also applied to drive ions through the interface. Ions crossing the interface tended to shed water molecules as they entered the hexanol bilayer, leaving a trail of water molecules. Stabilization and facilitated transport of the ion by interactions with the second layer of hexanol molecules appeared to be an important step in the mechanism of sodium ion transport.     

Detecting a sodium chloride ion pair in ice using terahertz time-domain spectroscopy.

The hydrogen bond resonance of a sodium chloride (NaCl) ion pair trapped in aqueous ice has been observed by transmission terahertz time-domain spectroscopy. The absorption peak of a sodium chloride ion pair in ice is 1.65 THz at 83 K. By investigating the interaction of the cation and anion with other chemical compounds, we deduce that the absorption peak originates from the hydrogen bond resonance of sodium chloride and water molecules. The charge redistribution that occurs when other ion pairs are added to aqueous salt solution changes the absorption spectrum. Furthermore, the results also indicate that simple molecules such as sodium halides have fingerprints in the terahertz region when the ions are trapped in ice. NaCl ion pairs in seawater and in Ringer's solution were examined.     

Web interface for Brownian dynamics simulation of ion transport and its applications to beta-barrel pores

Brownian dynamics (BD) based on accurate potential of mean force is an efficient and accurate method for simulating ion transport through wide ion channels. Here, a web-based graphical user interface (GUI) is presented for carrying out grand canonical Monte Carlo (GCMC) BD simulations of channel proteins: http://www.charmm-gui.org/input/gcmcbd. The webserver is designed to help users avoid most of the technical difficulties and issues encountered in setting up and simulating complex pore systems. GCMC/BD simulation results for three proteins, the voltage dependent anion channel (VDAC), α-Hemolysin (α-HL), and the protective antigen pore of the anthrax toxin (PA), are presented to illustrate the system setup, input preparation, and typical output (conductance, ion density profile, ion selectivity, and ion asymmetry). Two models for the input diffusion constants for potassium and chloride ions in the pore are compared: scaling of the bulk diffusion constants by 0.5, as deduced from previous all-atom molecular dynamics simulations of VDAC, and a hydrodynamics based model (HD) of diffusion through a tube. The HD model yields excellent agreement with experimental conductances for VDAC and α-HL, while scaling bulk diffusion constants by 0.5 leads to underestimates of 10-20%. For PA, simulated ion conduction values overestimate experimental values by a factor of 1.5-7 (depending on His protonation state and the transmembrane potential), implying that the currently available computational model of this protein requires further structural refinement.     

Analysis of ligand-bound water molecules in high-resolution crystal structures of protein-ligand complexes

We have performed a comprehensive analysis of water molecules at the protein-ligand interfaces observed in 392 high-resolution crystal structures. There are a total of 1829 ligand-bound water molecules in these 392 complexes; 18% are surface water molecules, and 72% are interfacial water molecules. The number of ligand-bound water molecules in each complex structure ranges from 0 to 21 and has an average of 4.6. Of these interfacial water molecules, 76% are considered to be bridging water molecules, characterized by having polar interactions with both ligand and protein atoms. Among a number of factors that may influence the number of ligand-bound water molecules, the polar van der Waals (vdw) surface area of ligands has the highest Pearson linear correlation coefficient of 0.63. Our regression analysis predicted that one more ligand-bound water molecule is expected for every additional 24 A2 in the polar vdw surface area of the ligand. In contrast to the observation that the resolution is the primary factor influencing the number of water molecules in crystallographic models of proteins, we found that there is only a weak relationship between the number of ligand-bound water molecules and the resolution of the crystal structures. An analysis of the isotropic B factors of buried ligand-bound water molecules suggested that, when water molecules have fewer than two polar interactions with the protein-ligand complex, they are more mobile than protein atoms in the crystal structures; when they have more than three polar interactions, they are significantly less mobile than protein atoms.     

A Halogen-Bond-Driven Artificial Chloride-Selective Channel Constructed from 5-Iodoisophthalamide-based Molecules

Despite considerable emphasis on advancing artificial ion channels, progress is constrained by the limited availability of small molecules with the necessary attributes of self-assembly and ion selectivity. In this study, a library of small molecules based on 5-haloisophthalamide and a non-halogenated isophthalamide were examined for their ion transport properties across the lipid bilayer membranes, and the finding demonstrates that the di-hexyl-substituted 5-iodoisophthalamide derivative exhibits the highest level of activity. Furthermore, it was established that the highest active compound facilitates the selective chloride transport that occurs via an antiport-mediated mechanism. The crystal structure of the compound unveils a distinctive self-assembly of molecules, forming a zig-zag channel pore that is well-suited for the permeation of anions. Planar bilayer conductance measurements proved the formation of chloride selective channels. A molecular dynamics simulation study, relying on the self-assembled component derived from the crystal structure, affirmed the paramount significance of intermolecular hydrogen bonding in the formation of supramolecular barrel-rosette structures that span the bilayer. Furthermore, it was demonstrated that the transport of chloride across the lipid bilayer membrane is facilitated by the synergistic effects of halogen bonding and hydrogen bonding within the channel.     

Computer Simulation of Cation-Exchange Membrane Structure: An Elementary Act of Hydrated Ion Transport

A fragment of the structure of a sulfo cation exchanger, which is the basis of most cation-exchange membranes, is calculated by an ab initio method. An analysis of interatomic bonds in the structure shows that, to detach a mobile ion from a fixed ion, it is necessary to break the hydrogen bond between hydration water molecules of the counterion in addition to overcoming the electrostatic attraction. As the hydrogen bond cleavage work for simple hydrated ions is ten times the electrostatic attraction energy, an elementary act of ion transport in a cation-exchange membrane is considered mostly as the hydrogen bond transfer reaction.     

Water inside a hydrophobic cavitand molecule

The structure and dynamics of water inside a water-soluble, bowl-shaped cavitand molecule with a hydrophobic interior are studied using molecular dynamics computer simulations. The simulations find that the number of inside water molecules is about 4.5, but it fluctuates from being completely empty to full on a time scale of tens of nanoseconds. The transition from empty to full is energetically favorable and entropically unfavorable. The water molecules inside have fewer hydrogen bonds than the bulk and in general weaker interactions; the lower energy results from the nearest-neighbor interactions with the cavitand atoms and the water molecules at the entrance of the cavitand, interactions that are lost upon dewetting. An analysis of translational and rotational motion suggests that the lower entropy of the inside water molecules is due to decreased translational entropy, which outweighs an increased orientational entropy. The cavitand molecule acts as a host binding hydrophobic guests, and dewetting can be induced by the presence of a hydrophobic guest molecule about 3 A above the entrance. At this position, the guest displaces the water molecules which stabilize the inside water molecules and the empty cavitand becomes more stable than the full.     

The Role of Aromatic Residues in Stabilizing the Secondary and Tertiary Structure of Avian Pancreatic Polypeptide

Avian Pancreatic Polypeptide is a 36 residue protein that exhibits a tertiary fold. Results of previous experimental and computational studies indicate that the structure of aPP is stabilized more by non-bonded interactions than by the hydrophobic effect. Aromatic residues are known to participate in a variety of long range non-bonded interactions, with both backbone atoms and the atoms of other side-chains, which could be responsible, in part, for the stability of both the local secondary structure and the tertiary fold. The effect of these aromatic interactions on the stability of aPP was calculated using BHandHLYP/cc-pVTZ. Aromatic residues were shown to participate in multiple hydrogen bonded and weakly polar interactions in the secondary structure. The energies of the weakly polar interactions are comparable with those of hydrogen bonds. Aromatic residues were also shown to participate in multiple weakly polar interactions across the tertiary fold, again with energies similar to those of hydrogen bonds.     

Ab Initio quantum mechanical charge field study of hydrated bicarbonate ion: Structural and dynamical properties

The ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism was applied to simulate the bicarbonate ion, HCO₃^?, in aqueous solution. The difference in coordination numbers obtained by summation over atoms (6.6) and for the solvent‐accessible surface (5.4) indicates the sharing of some water molecules between the individual atomic hydration shells. It also proved the importance to consider the hydration of the chemically different atoms individually for the evaluation of structural and dynamical properties of the ion. The orientation of water molecules in the hydration shell was visualized by the θ–tilt surface plot. The mean residence time in the surroundings of the HCO₃^? ion classify it generally as a structure‐breaking ion, but the analysis of the individual ion‐water hydrogen bonds revealed a more complex behavior of the different coordination sites. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Molecular dynamics simulations for water and ions in protein crystals

The spatial and temporal properties of water and ions in bionanoporous materials-protein crystals-have been investigated using molecular dynamics simulations. Three protein crystals are considered systematically with different morphologies and chemical topologies: tetragonal lysozyme, orthorhombic lysozyme, and tetragonal thermolysin. It is found that the thermal fluctuations of C(alpha) atoms in the secondary structures of protein molecules are relatively weak due to hydrogen bonding. The solvent-accessible surface area per residue is nearly identical in the three protein crystals; the hydrophobic and hydrophilic residues in each crystal possess approximately the same solvent-accessible surface area. Water distributes heterogeneously and has different local structures within the biological nanopores of the three protein crystals. The mobility of water and ions in the crystals is enhanced as the porosity increases and also by the fluctuations of protein atoms particularly in the two lysozyme crystals. Anisotropic diffusion is found preferentially along the pore axis, as experimentally observed. The anisotropy of the three crystals increases in the order: tetragonal thermolysin < tetragonal lysozyme < orthorhombic lysozyme.     

Interactions between hydrophobic and ionic solutes in aqueous guanidinium chloride and urea solutions: lessons for protein denaturation mechanism

In order to clarify the mechanism of denaturant-induced unfolding of proteins we have calculated the interactions between hydrophobic and ionic species in aqueous guanidinium chloride and urea solutions using molecular dynamics simulations. Hydrophobic association is not significantly changed in urea or guanidinium chloride solutions. The strength of interaction between ion pairs is greatly diminished by the guanidinium ion. Although the changes in electrostatic interactions in urea are small, examination of structures, using appropriate pair functions, of urea and water around the solutes show strong hydrogen bonding between urea's carbonyl oxygen and the positively charged solute. Our results strongly suggest protein denaturation occurs by the direct interaction model according to which the most commonly used denaturants unfold proteins by altering electrostatic interactions either by solvating the charged residues or by engaging in hydrogen bonds with the protein backbone. To further validate the direct interaction model we show that, in urea and guanidinium chloride solutions, unfolding of an unusually stable helix (H1) from mouse PrPC (residues 144-153) occurs by hydrogen bonding of denaturants to charged side chains and backbone carbonyl groups.     

Physical quantity of residue electrostatic energy in flavin mononucleotide binding protein dimer

The electrostatic (ES) energy of each residue was for the first time quantitatively evaluated in a flavin mononucleotide binding protein (FBP). A residue electrostatic energy (RES) was obtained as the sum of the ES energies between atoms in each residue and all other atoms in the FBP dimer using atomic coordinates obtained by a molecular dynamics (MD) simulation. ES is one of the most important energies among the interaction energies in a protein. It is determined from the RES, the residues which mainly contribute to stabilize the structure of each subunit, and the binding energy between two subunits can be estimated. The RES of all residues in subunit A (Sub A) and subunit B (Sub B) were attractive forces, even though the residues contain net negative or positive charges. This reveals that the ES energies of any of the residues can contribute to stabilize the protein structure. The total binding ES energy over all residues among the subunits was distributed between −0.2 to −1.2 eV (mean = −0.67 eV) from the MD simulation time.