Structure and transport properties of a plastic crystal ion conductor : Diethyl(methyl)(isobutyl)phosphonium hexafluorophosphate

Understanding the ion transport behavior of organic ionic plastic crystals (OIPCs) is crucial for their potential application as solid electrolytes in various electrochemical devices such as lithium batteries. In the present work, the ion transport mechanism is elucidated by analyzing experimental data (single-crystal XRD, multinuclear solid-state NMR, DSC, ionic conductivity, and SEM) as well as the theoretical simulations (second moment-based solid static NMR line width simulations) for the OIPC diethyl(methyl)(isobutyl)phosphonium hexafluorophosphate ([P1,2,2,4][PF6]). This material displays rich phase behavior and advantageous ionic conductivities, with three solid–solid phase transitions and a highly “plastic” and conductive final solid phase in which the conductivity reaches 10–3 S cm–1. The crystal structure shows unique channel-like packing of the cations, which may allow the anions to diffuse more easily than the cations at lower temperatures. The strongly phase-dependent static NMR line widths of the 1H, 19F, and 31P nuclei in this material have been well simulated by different levels of molecular motions in different phases. Thus, drawing together of the analytical and computational techniques has allowed the construction of a transport mechanism for [P1,2,2,4][PF6]. It is also anticipated that utilization of these techniques will allow a more detailed understanding of the transport mechanisms of other plastic crystal electrolyte materials.

Motional heterogeneities in siloxane/poly(ethylene glycol) ormolyte nanocomposites studied by C-13 solid-state exchange NMR

C-13 exchange solid-state NMR methods were used to study two families of siloxane/poly-(ethylene glycol) hybrid materials: Types I and II, where the polymer chains interact with the inorganic phase through physical (hydrogen bonds or van der Waals forces) or chemical (covalent bonds) interactions, respectively. These methods were employed to analyze the effects of the interactions between the organic and inorganic phases on the polymer dynamics in the milliseconds to seconds time scale, which occurs at temperatures below the motional narrowing of the NMR line width and around the polymer glass transition. Motional heterogeneities associated with these interactions and evidence of both small and large amplitude motions were directly observed for both types of hybrids. The results revealed that the hindrance to the slow molecular motions of the polymer chains due to the siloxane structures depends on the chain length and the nature of the interaction between the organic and inorganic phases.

Untersuchung der Strukturdynamik arbeitender F1-ATPase und ATP-Synthase aus Micrococcus luteus in Einzelschußexperimenten mit Synchrotronstrahlung

ZusammenfassungDie ATP-Synthase koppelt im Energiestoffwechsel der Zellen den Protonentransport über die biologische Membran mit der Synthese des energiespeichernden Moleküls ATP aus ADP und Phosphat. ATP-Synthasen bestehen aus 2 Subkomplexen, wobei der katalytische F1-Teil von der membranständigen Domäne abgelöst werden kann und nur zur ATP-Hydrolyse fähig ist. Der hochkooperative Reaktionsmechanismus der dreizentrigen ATP-Synthasen ist weitgehend unklar.Im Rahmen dieser Arbeit wurde der ATP-Synthasekomplex und ihr wasserlösliches katalytisches F1-Fragment aus Micrococcus luteus in präparativem Maßstab mittels chromatographischer Trennmethoden isoliert. Die Überprüfung der Funktionalität beider Enzyme erfolgte mit enzymatischen Methoden. Durch zeitaufgelöste Röntgenkleinwinkelstreuung wurde die Strukturdynamik der arbeitenden ATP-Synthase und ihres F1-Fragmentes aus Micrococcus luteus im Laufe des ATP-Hydrolysezyklus untersucht. Diese Methode diente zum Nachweis weiträumiger Konformationsänderungen innerhalb der arbeitenden Enzyme unter nativen physiologischen Bedingungen. Die zeitaufgelösten Streuexperimente fanden an der ESRF (Europäische Synchrotronstrahlungsquelle) in Grenoble (F) statt. Dort wurden für beide Enzyme im Laufe des ATP-Hydrolysezykus molekulare Bewegungen nachgewiesen. Als Referenz zu den zeitaufgelösten Messungen dienten statische Messungen zur Strukturuntersuchung der Proteine am schwächeren DESY. Anhand dieser Strukturdaten wurden Molekülmodelle der F1-ATPase und ATP-Synthase aus Micrococcus luteus konstruiert. Das Molekülmodell der F1-ATPase war die Grundlage zur Modellierung einzelner Teilschritte des ATP-Hydrolysezyklus bei 20°C. Die experimentellen Daten wurden mit einer Kippbewegung der membranseitigen Domäne der katalytischen b-Untereinheiten der F1-ATPase während des ATP-Hydrolysezyklus interpretiert.

Attenuated T2 relaxation by mutual cancellation of dipole–dipole coupling and chemical shift anisotropy indicates an avenue to NMR structures of very large biological macromolecules in solution

Fast transverse relaxation of 1H, 15N, and 13C by dipole-dipole coupling (DD) and chemical shift anisotropy (CSA) modulated by rotational molecular motions has a dominant impact on the size limit for biomacromolecular structures that can be studied by NMR spectroscopy in solution. Transverse relaxation-optimized spectroscopy (TROSY) is an approach for suppression of transverse relaxation in multidimensional NMR experiments, which is based on constructive use of interference between DD coupling and CSA. For example, a TROSY-type two-dimensional 1H,15N-correlation experiment with a uniformly 15N-labeled protein in a DNA complex of molecular mass 17 kDa at a 1H frequency of 750 MHz showed that 15N relaxation during 15N chemical shift evolution and 1HN relaxation during signal acquisition both are significantly reduced by mutual compensation of the DD and CSA interactions. The reduction of the linewidths when compared with a conventional two-dimensional 1H,15N-correlation experiment was 60% and 40%, respectively, and the residual linewidths were 5 Hz for 15N and 15 Hz for 1HN at 4°C. Because the ratio of the DD and CSA relaxation rates is nearly independent of the molecular size, a similar percentagewise reduction of the overall transverse relaxation rates is expected for larger proteins. For a 15N-labeled protein of 150 kDa at 750 MHz and 20°C one predicts residual linewidths of 10 Hz for 15N and 45 Hz for 1HN, and for the corresponding uniformly 15N,2H-labeled protein the residual linewidths are predicted to be smaller than 5 Hz and 15 Hz, respectively. The TROSY principle should benefit a variety of multidimensional solution NMR experiments, especially with future use of yet somewhat higher polarizing magnetic fields than are presently available, and thus largely eliminate one of the key factors that limit work with larger molecules.

Splitting the fast and slow motions in molecular dynamics simulations based on the change of cold potential well bottom

Abstract. The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.

Splitting the Fast and Slow Motions in Molecular Dynamics Simulations Based on the Change of Cold Potential Well Bottom

The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.

Continuum mechanics modelling of microfilament networks with different architectures based on molecular investigation of single F-actin

The actin microfilament plays a critical role in many cellular processes including embryonic development, wound healing, immune response, and tissue development. It is commonly organized in the form of networks whose mechanical properties change with changes in their architecture due to cell evolution processes. This paper presents a new nonlinear continuum mechanics model of single filamentous actin (F-actin) that is based on nanoscale molecular simulations. Following this continuum model of the single F-actin, mechanical properties of differently architected lamellipodia are studied. The results provide insight that can contribute to the understanding of the cell edge motions of living cells.

NMR studies of molecular dynamics in tetramethylammonium tetrabromo cadmate

Internal motions in a A2BX4 compound (tetramethylammonium tetrabromo cadmate) have been investigated using proton spin—lattice relaxation time (T1) and second moment (M2) measurements in the temperature range 77 to 400 K. T1 measurements at three Larmor frequencies (10, 20 and 30 MHz) show isotropic tumbling of the tetramethylammonium group, random reorientation of methyl groups and spin—rotation interaction, and the corresponding parameters have been computed. The cw spectrum is narrow throughout the temperature range and shows side bands at the lowest temperature. This observation, along with the free-induction-decay behavior at these temperatures, is interpreted as the onset of a coherent motion, e.g. methyl group quantum tunnelling.

The hydrogen bond: a molecular beam microwave spectroscopist's view with a universal appeal

In this manuscript, we propose a criterion for a weakly bound complex formed in a supersonic beam to be characterized as a `hydrogen bonded complex'. For a `hydrogen bonded complex', the zero point energy along any large amplitude vibrational coordinate that destroys the orientational preference for the hydrogen bond should be significantly below the barrier along that coordinate so that there is at least one bound level. These are vibrational modes that do not lead to the breakdown of the complex as a whole. If the zero point level is higher than the barrier, the `hydrogen bond' would not be able to stabilize the orientation which favors it and it is no longer sensible to characterize a complex as hydrogen bonded. Four complexes, Ar-2-H2O, Ar-2-H2S, C2H4-H2O and C2H4-H2S, were chosen for investigations. Zero point energies and barriers for large amplitude motions were calculated at a reasonable level of calculation, MP2(full)/aug-cc-pVTZ, for all these complexes. Atoms in molecules (AIM) theoretical analyses of these complexes were carried out as well. All these complexes would be considered hydrogen bonded according to the AIM theoretical criteria suggested by Koch and Popelier for C-H center dot center dot center dot O hydrogen bonds (U. Koch and P. L. A. Popelier, J. Phys. Chem., 1995, 99, 9747), which has been widely and, at times, incorrectly used for all types of contacts involving H. It is shown that, according to the criterion proposed here, the Ar-2-H2O/H2S complexes are not hydrogen bonded even at zero kelvin and C2H4-H2O/H2S complexes are. This analysis can naturally be extended to all temperatures. It can explain the recent experimental observations on crystal structures of H2S at various conditions and the crossed beam scattering studies on rare gases with H2O and H2S.

Study of molecular reorientation and quantum rotational tunneling in tetramethylammonium selenate by H-1 NMR

H-1 NMR spin-lattice relaxation time measurements have been carried out in [(CH3)(4)N](2)SeO4 in the temperature range 389-6.6K to understand the possible phase transitions, internal motions and quantum rotational tunneling. A broad T, minimum observed around 280K is attributed to the simultaneous motions of CH3 and (CH3)(4)N groups. Magnetization recovery is found to be stretched exponential below 72 K with varying stretched exponent. Low-temperature T-1 behavior is interpreted in terms of methyl groups undergoing quantum rotational tunneling. (c) 2007 Elsevier Inc. All rights reserved.

Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Crystal Structure of Staphylococcus aureus Metallopeptidase (Sapep) Reveals Large Domain Motions between the Manganese-bound and Apo-states

Proteases belonging to the M20 family are characterized by diverse substrate specificity and participate in several metabolic pathways. The Staphylococcus aureus metallopeptidase, Sapep, is a member of the aminoacylase-I/M20 protein family. This protein is a Mn2+-dependent dipeptidase. The crystal structure of this protein in the Mn2+-bound form and in the open, metal-free state suggests that large interdomain movements could potentially regulate the activity of this enzyme. We note that the extended inactive conformation is stabilized by a disulfide bond in the vicinity of the active site. Although these cysteines, Cys(155) and Cys(178), are not active site residues, the reduced form of this enzyme is substantially more active as a dipeptidase. These findings acquire further relevance given a recent observation that this enzyme is only active in methicillin-resistant S. aureus. The structural and biochemical features of this enzyme provide a template for the design of novel methicillin-resistant S. aureus-specific therapeutics.

Complexation of siRNA with Dendrimer: A Molecular Modeling Approach

This paper reports the structural behavior and thermodynamics of the complexation of siRNA with poly(amidoamine) (PAMAM) dendrimers of generation 3 (G3) and 4 (G4) through fully atomistic molecular dynamics (MD) simulations accompanied by free energy calculations and inherent structure determination. We have also done simulation with one siRNA and two dendrimers (2 x G3 or 2xG4) to get the microscopic picture of various binding modes. Our simulation results reveal the formation of stable siRNA-dendrimer complex over nanosecond time scale. With the increase in dendrimcr generation, the charge ratio increases and hence the binding energy between siRNA and dendrimer also increases in accordance with available experimental measurements. Calculated radial distribution functions of amines groups of various subgenerations in a given generation of dendrimer and phosphate in backbone of siRNA reveals that one dendrimer of generation 4 shows better binding with siRNA almost wrapping the dendrimer when compared to the binding with lower generation dendrimer like G3. In contrast, two dendrimers of generation 4 show binding without siRNA wrapping the den-rimer because of repulsion between two dendrimers. The counterion distribution around the complex and the water molecules in the hydration shell of siRNA give microscopic picture of the binding dynamics. We see a clear correlation between water. counterions motions and the complexation i.e. the water molecules and counterions which condensed around siRNA are moved away from the siRNA backbone when dendrimer start binding to the siRNA back hone. As siRNA wraps/bind to the dendrimer counterions originally condensed onto siRNA (Na-1) and dendrimer (Cl-) get released. We give a quantitative estimate of the entropy of counterions and show that there is gain in entropy due to counterions release during the complexation. Furthermore, the free energy of complexation of IG3 and IG4 at two different salt concentrations shows that increase in salt concentration leads to the weakening of the binding affinity of siRNA and dendrimer.

Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase

The catalytic conversion of adenosine triphosphate (ATP) and adenosine monophosphate (AMP) to adenosine diphosphate (ADP) by adenylate kinase (ADK) involves large amplitude, ligand induced domain motions, involving the opening and the closing of ATP binding domain (LID) and AMP binding domain (NMP) domains, during the repeated catalytic cycle. We discover and analyze an interesting dynamical coupling between the motion of the two domains during the opening, using large scale atomistic molecular dynamics trajectory analysis, covariance analysis, and multidimensional free energy calculations with explicit water. Initially, the LID domain must open by a certain amount before the NMP domain can begin to open. Dynamical correlation map shows interesting cross-peak between LID and NMP domain which suggests the presence of correlated motion between them. This is also reflected in our calculated two-dimensional free energy surface contour diagram which has an interesting elliptic shape, revealing a strong correlation between the opening of the LID domain and that of the NMP domain. Our free energy surface of the LID domain motion is rugged due to interaction with water and the signature of ruggedness is evident in the observed root mean square deviation variation and its fluctuation time correlation functions. We develop a correlated dynamical disorder-type theoretical model to explain the observed dynamic coupling between the motion of the two domains in ADK. Our model correctly reproduces several features of the cross-correlation observed in simulations. (C) 2011 American Institute of Physics. doi:10.1063/1.3516588]

Structure and order of the discotic compound 2,3,6,7,10,11-hexakispentyloxytriphenylene as revealed by diffraction and molecular simulation studies

The aggregate structure of the discotic compound 2,3,6,7,10,11-hexakispentyloxytriphenylene (HPT) was studied both for the crystalline state and the liquid crystalline state by using electron crystallography and a molecular simulation approach. In the crystalline state, HPT was found to adopt an orthorhombic P-2212 space group with cell parameters a = 36.73 Angstrom, b = 27.99 Angstrom and c = 4.91 Angstrom. Molecular packing calculations were conducted to elucidate the molecular conformation and mutual orientational characteristics in the different states. Phase transitions and relationships are discussed from a structural point of view.