A cellulosic liquid crystal pool for cellulose nanocrystals: structure and molecular dynamics at high shear rates

Cellulose and its derivatives, such as hydroxypropylcellulose (HPC) have been studied for a long time but they are still not well understood particularly in liquid crystalline solutions. These systems can be at the origin of networks with properties similar to liquid crystalline (LC) elastomers. The films produced from LC solutions can be manipulated by the action of moisture allowing for instance the development of a soft motor (Geng et al., 2013) driven by humidity. Cellulose nanocrystals (CNC), which combine cellulose properties with the specific characteristics of nanoscale materials, have been mainly studied for their potential as a reinforcing agent. Suspensions of CNC can also self-order originating a liquid-crystalline chiral nematic phases. Considering the liquid crystalline features that both LC-HPC and CNC can acquire, we prepared LC-HPC/CNC solutions with different CNC contents (1,2 and 5 wt.%). The effect of the CNC into the LC-HPC matrix was determined by coupling rheology and NMR spectroscopy - Rheo-NMR a technique tailored to analyse orientational order in sheared systems. (C) 2015 Elsevier Ltd. All rights reserved.

The effects of solvent composition on the affinity of a peptide towards hair keratin: experimental and molecular dynamics data

The study of the interaction between hair filaments and formulations or peptides is of utmost importance in fields like cosmetic research. Keratin intermediate filaments structure is not fully described, limiting the molecular dynamics (MD) studies in this field although its high potential to improve the area. We developed a computational model of a truncated protofibril, simulated its behavior in alcoholic based formulations and with one peptide. The simulations showed a strong interaction between the benzyl alcohol molecules of the formulations and the model, leading to the disorganization of the keratin chains, which regress with the removal of the alcohol molecules. This behavior can explain the increase of peptide uptake in hair shafts evidenced in fluorescence microscopy pictures. The model developed is valid to computationally reproduce the interaction between hair and alcoholic formulations and provide a robust base for new MD studies about hair properties. It is shown that the MD simulations can improve hair cosmetic research, improving the uptake of a compound of interest.

Molecular dynamics of the neuronal Ca 2+ -binding proteins Caldendrin and Calneurons

Magdeburg, Univ., Fak. für Naturwiss., Diss., 2009

Methods for analyzing the influence of molecular dynamics on neuronal activity

Magdeburg, Univ., Fak. für Informatik, Diss., 2015

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex.

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.

Assessing protein stability of the dimeric DNA-binding domain of E2 human papillomavirus 18 with molecular dynamics

The objective of this study is to understand the structural flexibility and curvature of the E2 protein of human papillomavirus type 18 using molecular dynamics (6 ns). E2 is required for viral DNA replication and its disruption could be an anti-viral strategy. E2 is a dimer, with each monomer folding into a stable open-faced β-sandwich. We calculated the mobility of the E2 dimer and found that it was asymmetric. These different mobilities of E2 monomers suggest that drugs or vaccines could be targeted to the interface between the two monomers.

Molecular determinants of desensitization in an ENaC/degenerin channel.

Epithelial Na(+) channel (ENaC)/degenerin family members are involved in mechanosensation, blood pressure control, pain sensation, and the expression of fear. Several of these channel types display a form of desensitization that allows the channel to limit Na(+) influx during prolonged stimulation. We used site-directed mutagenesis and chemical modification, functional analysis, and molecular dynamics simulations to investigate the role of the lower palm domain of the acid-sensing ion channel 1, a member of the ENaC/degenerin family. The lower palm domains of this trimeric channel are arranged around a central vestibule, at ∼20 Å above the plasma membrane and are covalently linked to the transmembrane channel parts. We show that the lower palm domains approach one another during desensitization. Residues in the palm co-determine the pH dependence of desensitization, its kinetics, and the stability of the desensitized state. Mutations of palm residues impair desensitization by preventing the closing movement of the palm. Overexpression of desensitization-impaired channel mutants in central neurons allowed--in contrast to overexpression of wild type--a sustained signaling response to rapid pH fluctuations. We identify and describe here the function of an important regulatory domain that most likely has a conserved role in ENaC/degenerin channels.

Influence of ionization state on the activation of temocapril by hCES1: a molecular-dynamics study.

Temocapril is a prodrug whose hydrolysis by carboxylesterase 1 (CES1) yields the active ACE inhibitor temocaprilat. This molecular-dynamics (MD) study uses a resolved structure of the human CES1 (hCES1) to investigate some mechanistic details of temocapril hydrolysis. The ionization constants of temocapril (pK1 and pK3) and temocaprilat (pK1, pK2, and pK3) were determined experimentally and computationally using commercial algorithms. The constants so obtained were in good agreement and revealed that temocapril exists mainly in three ionic forms (a cation, a zwitterion, and an anion), whereas temocaprilat exists in four major ionic forms (a cation, a zwitterion, an anion, and a dianion). All these ionic forms were used as ligands in 5-ns MS simulations. While the cationic and zwitterionic forms of temocapril were involved in an ion-pair bond with Glu255 suggestive of an inhibitor behavior, the anionic form remained in a productive interaction with the catalytic center. As for temocaprilat, its cation appeared trapped by Glu255, while its zwitterion and anion made a slow departure from the catalytic site and a partial egress from the protein. Only its dianion was effectively removed from the catalytic site and attracted to the protein surface by Lys residues. A detailed mechanism of product egress emerges from the simulations.

Comparison between molecular dynamics abd Monte Carlo simulations of an ordering process in a binary alloy

Ordering in a binary alloy is studied by means of a molecular-dynamics (MD) algorithm which allows to reach the domain growth regime. Results are compared with Monte Carlo simulations using a realistic vacancy-atom (MC-VA) mechanism. At low temperatures fast growth with a dynamical exponent x>1/2 is found for MD and MC-VA. The study of a nonequilibrium ordering process with the two methods shows the importance of the nonhomogeneity of the excitations in the system for determining its macroscopic kinetics.

Molecular dynamics study of orientational cooperativity in water

Recent experiments on liquid water show collective dipole orientation fluctuations dramatically slower than expected (with relaxation time >tation, the self-dipole randomization time tr, which is an upper limit on ta; we find that tr5ta. Third, to check if there are correlated domains of dipoles in water which have large relaxation times compared to the individual dipoles, we calculate the randomization time tbox of the site-dipole field, the net dipole moment formed by a set of molecules belonging to a box of edge Lbox. We find that the site-dipole randomization time tbox2.5ta for Lbox3 , i.e., it is shorter than the same quantity calculated for the self-dipole. Finally, we find that the orientational correlation length is short even at low T.

Computer simulation study of liquid lithium at 470 and 843 K

Both structural and dynamical properties of 7Li at 470 and 843 K are studied by molecular dynamics simulation and the results are comapred with the available experimental data. Two effective interatomic potentials are used, i.e., a potential derived from the Ashcroft pseudopotential [Phys. Lett. 23, 48 (1966)] and a recently proposed potential deduced from the neutral pseudoatom method [J. Phys.: Condens. Matter 5, 4283 (1993)]. Although the shape of the two potential functions is very different, the majority of the properties calculated from them are very similar. The differences among the results using the two interaction models are carefully discussed.

Molecular dynamics study of the longitudinal modes in disparate-mass binaryliquid mixtures.

A series of molecular dynamics simulations of simple liquid binary mixtures of soft spheres with disparate-mass particles were carried out to investigate the origin of the marked differences between the dynamic structure factors of some liquid binary mixtures such as the Li0.7Mg0.3 and Li0.8Pb0.2 alloys. It is shown that the facility for observing peaks associated with fast-propagating modes in the partial Li-Li dynamic structure factor of Li0.8Pb0.2 should be mainly attributed to the structure of this alloy, which is characterized by an incipient ABAB ordering as found in molten salts. The longitudinal dispersion relations at intermediate wave vectors obtained from the longitudinal current spectra are very similar for the two alloys and reflect the existence of both fast-and slow-propagating modes of kinetic character associated with light and heavy particles, respectively. The influence of the hardness of the repulsive potential cores as well as the composition of the mixture on the longitudinal collective modes is also discussed.

Understanding and modulating the competitive surface-adsorption of proteins through coarse-grained molecular dynamics simulations

It is now well accepted that cellular responses to materials in a biological medium reflect greatly the adsorbed biomolecular layer, rather than the material itself. Here, we study by molecular dynamics simulations the competitive protein adsorption on a surface (Vroman effect), i.e. the non-monotonic behavior of the amount of protein adsorbed on a surface in contact with plasma as functions of contact time and plasma concentration. We find a complex behavior, with regimes during which small and large proteins are not necessarily competing between them, but are both competing with others in solution ("cooperative" adsorption). We show how the Vroman effect can be understood, controlled and inverted.

Diffusive transport properties in monovalent and divalent metal-ion halide melts: A computer simulation study

Self- and cross-velocity correlation functions and related transport coefficients of molten salts are studied by molecular-dynamics simulation. Six representative systems are considered, i.e., NaCl and KCl alkali halides, CuCl and CuBr noble-metal halides, and SrCl2 and ZnCl2 divalent metal-ion halides. Computer simulation results are compared with experimental self-diffusion coefficients and electrical conductivities. Special attention is paid to dynamic cross correlations and their dependence on the Coulomb interactions as well as on the size and mass differences between anions and cations.

Fourier transform convolution integrals applied to generalized Born molecular volume.

Generalized Born methods are currently among the solvation models most commonly used for biological applications. We reformulate the generalized Born molecular volume method initially described by (Lee et al, 2003, J Phys Chem, 116, 10606; Lee et al, 2003, J Comp Chem, 24, 1348) using fast Fourier transform convolution integrals. Changes in the initial method are discussed and analyzed. Finally, the method is extensively checked with snapshots from common molecular modeling applications: binding free energy computations and docking. Biologically relevant test systems are chosen, including 855-36091 atoms. It is clearly demonstrated that, precision-wise, the proposed method performs as good as the original, and could better benefit from hardware accelerated boards.