Tractable molecular theory of transport of Lennard-Jones fluids in nanopores

We present here a tractable theory of transport of simple fluids in cylindrical nanopores, which is applicable over a wide range of densities and pore sizes. In the Henry law low-density region the theory considers the trajectories of molecules oscillating between diffuse wall collisions, while at higher densities beyond this region the contribution from viscous flow becomes significant and is included through our recent approach utilizing a local average density model. The model is validated by means of equilibrium as well nonequilibrium molecular dynamics simulations of supercritical methane transport in cylindrical silica pores over a wide range of temperature, density, and pore size. The model for the Henry law region is exact and found to yield an excellent match with simulations at all conditions, including the single-file region of very small pore size where it is shown to provide the density-independent collective transport coefficient. It is also shown that in the absence of dispersive interactions the model reduces to the classical Knudsen result, but in the presence of such interactions the latter model drastically overpredicts the transport coefficient. For larger micropores beyond the single-file region the transport coefficient is reduced at high density because of intermolecular interactions and hindrance to particle crossings leading to a large decrease in surface slip that is not well represented by the model. However, for mesopores the transport coefficient increases monotonically with density, over the range studied, and is very well predicted by the theory, though at very high density the contribution from surface slip is slightly overpredicted. It is also seen that the concept of activated diffusion, commonly associated with diffusion in small pores, is fundamentally invalid for smooth pores, and the apparent activation energy is not simply related to the minimum pore potential or the adsorption energy as generally assumed. (C) 2004 American Institute of Physics.

Modeling molecular transport in slit pores

We examine the transport of methane in microporous carbon by performing equilibrium and nonequilibrium molecular dynamics simulations over a range of pore sizes, densities, and temperatures. We interpret these simulation results using two models of the transport process. At low densities, we consider a molecular flow model, in which intermolecular interactions are neglected, and find excellent agreement between transport diffusion coefficients determined from simulation, and those predicted by the model. Simulation results indicate that the model can be applied up to fluid densities of the order to 0.1-1 nm(-3). Above these densities, we consider a slip flow model, combining hydrodynamic theory with a slip condition at the solid-fluid interface. As the diffusion coefficient at low densities can be accurately determined by the molecular flow model, we also consider a model where the slip condition is supplied by the molecular flow model. We find that both density-dependent models provide a useful means of estimating the transport coefficient that compares well with simulation. (C) 2004 American Institute of Physics.

Molecular recognition of an RNA trafficking element by heterogeneous nuclear ribonucleoprotein A2

Heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a multitasking protein involved in RNA packaging, alternative splicing of pre-mRNA. telomere maintenance, cytoplasmic RNA trafficking, and translation. It binds short segments of single-stranded nucleic acids, including the A2RE11 RNA element that is necessary and sufficient for cytoplasmic transport of a subset of rnRNAs in oligodendrocytes and neurons. We have explored the structures of hnRNP A2, its RNA recognition motifs (RRMs) and Gly-rich module, and the RRM complexes with A2RE11. Circular dichroism spectroscopy showed that the secondary structure of the first 189 residues of hnRNP A2 parallels that of the tandem beta alpha beta beta alpha beta RRMs of its paralogue, hnRNP A1, previously deduced from X-ray diffraction studies. The unusual GRD was shown to have substantial beta-sheet and beta-turn structure. Sedimentation equilibrium and circular dichroism results were consistent with the tandem RRM region being monomeric and supported earlier evidence for the binding of two A2RE11 oligoribonucleotides to this domain, in contrast to the protein dimer formed by the complex of hnRNP A1 with the telomeric ssDNA repeat. A three-dimensional structure for the N-terminal, two-RRM-containing segment of hnRNP A2 was derived by homology modeling. This structure was used to derive a model for the complex with A2RE11 using the previously described interaction of pairs of stacked nucleotides with aromatic residues on the RRM beta-sheet platforms, conserved in other RRM-RNA complexes, together with biochemical data and molecular dynamics-based observations of inter-RRM mobility.

Quantum effects on adsorption and diffusion of hydrogen and deuterium in microporous materials

Monte Carlo and molecular dynamics simulations and neutron scattering experiments are used to study the adsorption and diffusion of hydrogen and deuterium in zeolite Rho in the temperature range of 30-150 K. In the molecular simulations, quantum effects are incorporated via the Feynman-Hibbs variational approach. We suggest a new set of potential parameters for hydrogen, which can be used when Feynman-Hibbs variational approach is used for quantum corrections. The dynamic properties obtained from molecular dynamics simulations are in excellent agreement with the experimental results and show significant quantum effects on the transport at very low temperature. The molecular dynamics simulation results show that the quantum effect is very sensitive to pore dimensions and under suitable conditions can lead to a reverse kinetic molecular sieving with deuterium diffusing faster than hydrogen.

Structure and dynamics of ASC2, a pyrin domain-only protein that regulates inflammatory signaling

Pyrin domain (PYD)-containing proteins are key components of pathways that regulate inflammation, apoptosis, and cytokine processing. Their importance is further evidenced by the consequences of mutations in these proteins that give rise to autoimmune and hyperinflammatory syndromes. PYDs, like other members of the death domain ( DD) superfamily, are postulated to mediate homotypic interactions that assemble and regulate the activity of signaling complexes. However, PYDs are presently the least well characterized of all four DD subfamilies. Here we report the three-dimensional structure and dynamic properties of ASC2, a PYD-only protein that functions as a modulator of multidomain PYD-containing proteins involved in NF-KB and caspase-1 activation. ASC2 adopts a six-helix bundle structure with a prominent loop, comprising 13 amino acid residues, between helices two and three. This loop represents a divergent feature of PYDs from other domains with the DD fold. Detailed analysis of backbone N-15 NMR relaxation data using both the Lipari-Szabo model-free and reduced spectral density function formalisms revealed no evidence of contiguous stretches of polypeptide chain with dramatically increased internal motion, except at the extreme N and C termini. Some mobility in the fast, picosecond to nanosecond timescale, was seen in helix 3 and the preceding alpha 2-alpha 3 loop, in stark contrast to the complete disorder seen in the corresponding region of the NALP1 PYD. Our results suggest that extensive conformational flexibility in helix 3 and the alpha 2-alpha 3 loop is not a general feature of pyrin domains. Further, a transition from complete disorder to order of the alpha 2-alpha 3 loop upon binding, as suggested for NALP1, is unlikely to be a common attribute of pyrin domain interactions.

Implementation of particle-scale rotation in the 3-d Lattice Solid Model

In this study, 3-D Lattice Solid Model (LSMearth or LSM) was extended by introducing particle-scale rotation. In the new model, for each 3-D particle, we introduce six degrees of freedom: Three for translational motion, and three for orientation. Six kinds of relative motions are permitted between two neighboring particles, and six interactions are transferred, i.e., radial, two shearing forces, twisting and two bending torques. By using quaternion algebra, relative rotation between two particles is decomposed into two sequence-independent rotations such that all interactions due to the relative motions between interactive rigid bodies can be uniquely decided. After incorporating this mechanism and introducing bond breaking under torsion and bending into the LSM, several tests on 2-D and 3-D rock failure under uni-axial compression are carried out. Compared with the simulations without the single particle rotational mechanism, the new simulation results match more closely experimental results of rock fracture and hence, are encouraging. Since more parameters are introduced, an approach for choosing the new parameters is presented.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

Computational mechanics of molecular systems:quantifying high-dimensional dynamics by distribution of Poincaré recurrence times

A framework that connects computational mechanics and molecular dynamics has been developed and described. As the key parts of the framework, the problem of symbolising molecular trajectory and the associated interrelation between microscopic phase space variables and macroscopic observables of the molecular system are considered. Following Shalizi and Moore, it is shown that causal states, the constituent parts of the main construct of computational mechanics, the e-machine, define areas of the phase space that are optimal in the sense of transferring information from the micro-variables to the macro-observables. We have demonstrated that, based on the decay of their Poincare´ return times, these areas can be divided into two classes that characterise the separation of the phase space into resonant and chaotic areas. The first class is characterised by predominantly short time returns, typical to quasi-periodic or periodic trajectories. This class includes a countable number of areas corresponding to resonances. The second class includes trajectories with chaotic behaviour characterised by the exponential decay of return times in accordance with the Poincare´ theorem.

Identifying and correcting non-Markov states in peptide conformational dynamics

Conformational transitions in proteins define their biological activity and can be investigated in detail using the Markov state model. The fundamental assumption on the transitions between the states, their Markov property, is critical in this framework. We test this assumption by analyzing the transitions obtained directly from the dynamics of a molecular dynamics simulated peptide valine-proline-alanine-leucine and states defined phenomenologically using clustering in dihedral space. We find that the transitions are Markovian at the time scale of ˜ 50 ps and longer. However, at the time scale of 30–40 ps the dynamics loses its Markov property. Our methodology reveals the mechanism that leads to non-Markov behavior. It also provides a way of regrouping the conformations into new states that now possess the required Markov property of their dynamics.

Persistence in a random bond ising model of socio-econo dynamics

We study the persistence phenomenon in a socio-econo dynamics model using computer simulations at a nite temperature on hypercubic lattices in dimensions up to ve. The model includes a \social" local eld which contains the magnetization at time t. The nearest neighbour quenched interactions are drawn from a binary distribution which is a function of the bond concentration, p. The decay of the persistence probability in the model depends on both the spatial dimension and p. We nd no evidence of \blocking" in this model. We also discuss the implications of our results for possible applications in the social and economic elds. It is suggested that the absence, or otherwise, of blocking could be used as a criterion to decide on the validity of a given model in dierent scenarios.

Complexity of classical dynamics of molecular systems. II. Finite statistical complexity of a water-Na+ system

The computational mechanics approach has been applied to the orientational behavior of water molecules in a molecular dynamics simulated water–Na + system. The distinctively different statistical complexity of water molecules in the bulk and in the first solvation shell of the ion is demonstrated. It is shown that the molecules undergo more complex orientational motion when surrounded by other water molecules compared to those constrained by the electric field of the ion. However the spatial coordinates of the oxygen atom shows the opposite complexity behavior in that complexity is higher for the solvation shell molecules. New information about the dynamics of water molecules in the solvation shell is provided that is additional to that given by traditional methods of analysis.

Theoretical study of chloroperoxidase catalyzed chlorination of beta-cyclopentanedione and role of water in the chlorination mechanism

Chloroperoxidase (CPO) is a potential biocatalyst for use in asymmetric synthesis. The mechanisms of CPO catalysis are therefore of interest. The halogenation reaction, one of several chemical reactions that CPO catalyzes, is not fully understood and is the subject of this dissertation. The mechanism by which CPO catalyzes halogenation is disputed. It has been postulated that halogenation of substrates occurs at the active site. Alternatively, it has been proposed that hypochlorous acid, produced at the active site via oxidation of chloride, is released prior to reaction, so that halogenation occurs in solution. The free-solution mechanism is supported by the observation that halogenation of most substrates often occurs non-stereospecifically. On the other hand, the enzyme-bound mechanism is supported by the observation that some large substrates undergo halogenation stereospecifically. The major purpose of this research is to compare chlorination of the substrate β-cyclopentanedione in the two environments. One study was of the reaction with limited hydration because such a level of hydration is typical of the active site. For this work, a purely quantum mechanical approach was used. To model the aqueous environment, the limited hydration environment approach is not appropriate. Instead, reaction precursor conformations were obtained from a solvated molecular dynamics simulation, and reaction of potentially reactive molecular encounters was modeled with a hybrid quantum mechanical/molecular mechanical approach. Extensive work developing parameters for small molecules was pre-requisite for the molecular dynamics simulation. It is observed that a limited and optimized (active-site-like) hydration environment leads to a lower energetic barrier than the fully solvated model representative of the aqueous environment at room temperature, suggesting that the stable water network near the active site is likely to facilitate the chlorination mechanism. The influence of the solvent environment on the reaction barrier is critical. It is observed that stabilization of the catalytic water by other solvent molecules lowers the barrier for keto-enol tautomerization. Placement of water molecules is more important than the number of water molecules in such studies. The fully-solvated model demonstrates that reaction proceeds when the instantaneous dynamical water environment is close to optimal for stabilizing the transition state.

A study of electrochemical transport and diffuse charge dynamics using a Langevin equation

This thesis aims to develop new numerical and computational tools to study electrochemical transport and diffuse charge dynamics at small scales. Previous efforts at modeling electrokinetic phenomena at scales where the noncontinuum effects become significant have included continuum models based on the Poisson-Nernst-Planck equations and atomic simulations using molecular dynamics algorithms. Neither of them is easy to use or conducive to electrokinetic transport modeling in strong confinement or over long time scales. This work introduces a new approach based on a Langevin equation for diffuse charge dynamics in nanofluidic devices, which incorporates features from both continuum and atomistic methods. The model is then extended to include steric effects resulting from finite ion size, and applied to the phenomenon of double layer charging in a symmetric binary electrolyte between parallel-plate blocking electrodes, between which a voltage is applied. Finally, the results of this approach are compared to those of the continuum model based on the Poisson-Nernst-Planck equations.

Path-integral and Coarse-graining Strategies for Complex Molecular Phenomena

Molecular simulation provides a powerful tool for connecting molecular-level processes to physical observables. However, the facility to make those connections relies upon the application and development of theoretical methods that permit appropriate descriptions of the systems or processes to be studied. In this thesis, we utilize molecular simulation to study and predict two phenomena with very different theoretical challenges, beginning with (1) lithium-ion transport behavior in polymers and following with (2) equilibrium isotope effects with relevance to position-specific and clumped isotope studies. In the case of ion transport in polymers, there is motivation to use molecular simulation to provide guidance in polymer electrolyte design, but the length and timescales relevant for ion diffusion in polymers preclude the use of direct molecular dynamics simulation to compute ion diffusivities in more than a handful of candidate systems. In the case of equilibrium isotope effects, the thermodynamic driving forces for isotopic fractionation are often fundamentally quantum mechanical in nature, and the high precision of experimental instruments demands correspondingly accurate theoretical approaches. Herein, we describe respectively coarse-graining and path-integral strategies to address outstanding questions in these two subject areas.

Kohärente transiente und effiziente Endzustandskontrolle in Atomen und Molekülen

Es ist ein lang gehegter Traum in der Chemie, den Ablauf einer chemischen Reaktion zu kontrollieren und das Aufbrechen und Bilden chemischer Bindungen zu steuern. Diesem Ziel verschreibt sich auch das Forschungsgebiet der Femtochemie. Hier werden Femtosekunden Laserpulse eingesetzt um auf dem Quantenlevel molekulare Dynamiken auf ihren intrinsischen Zeitskalen zu kontrollieren und das System selektiv und effizient von einem Anfangs- in einen Zielzustand zu überführen. Der Wunsch, mit geformten Femtosekunden Laserpulsen Kontrolle über transiente Dynamiken und finale Populationen auszuüben, zu beobachten und zu verstehen, bildet auch die Motivation für diese Arbeit. Hierzu wurden mit Hilfe der Photoelektronenspektroskopie Untersuchungen zur Wechselwirkung atomarer und molekularer Prototypsysteme mit intensiven, geformten Femtosekunden Laserpulsen durchgeführt. Die Verwendung von Modelsystemen ermöglicht es, grundlegende Mechanismen der kohärenten Kontrolle in intensiven Laserfeldern zu analysieren, ohne dass sie durch komplexe Wechselwirkungen verschleiert werden. Zunächst wurde die Wechselwirkung von Kaliumatomen mit gechirpten Femtosekunden Laserpulsen untersucht. In den Experimenten wurden sowohl transiente Dynamiken als auch die Endbesetzungen der elektronischen Zustände abgebildet. In den folgenden Experimenten wurde das Quantenkontrollszenario SPODS auf die gekoppelte Elektronen-Kern-Dynamik in Molekülen übertragen. Die Kontrolle basiert auf der Erzeugung und Manipulation von Ladungsoszillationen durch Pulssequenzen. Der letzte Teil widmet sich der Entwicklung adiabatischer Kontrollmechanismen in Molekülen. Bei den Experimenten wurden gechirpte Airypulse eingesetzt um robuste Starkfeldanregung in molekularen Systemen zu induzieren. In Zukunft wird die Erforschung immer komplexerer Moleküle im Rahmen der transienten Kontrolle im Fokus stehen. Dabei werden nicht nur die effiziente Besetzung gebundener Zustände von Interesse sein, sondern auch die gezielte Dissoziation in spezifische Fragmente, photoinduzierte Isomerisierungsreaktionen oder die Kontrolle über transiente Dynamiken, die Einfluss auf andere molekulare Eigenschaften haben. Vor dem Hintergrund dieses übergeordneten Wunsches, photochemische Reaktionen immer komplexerer Moleküle, bis hin zu großen, biologisch relevanten Molekülen, zu kontrollieren, ist es umso wichtiger, die zugrundeliegenden Anregungsmechanismen in einfachen Systemen nachzuvollziehen. In den hier präsentierten Experimenten wurde gezeigt, wie die simultane Beobachtung der bekleideten und der stationären Zustände in atomaren Systemen zu einem umfassenden Bild der lichtinduzierte Dynamiken führen kann. Die gewonnenen Erkenntnisse können auf die Steuerung gekoppelter Dynamiken übertragen werden, durch die Kontrolle auch in molekularen Systemen möglich wird.