Structure-Property Relationships in Liquid Crystals: Can Modelling do Better than Empiricism?

Abstract

This paper reviews the relationships observed experimentally between the physical properties of liquid crystals and the molecular structures of the constituent molecules, and reports molecular mechanics calculations designed to provide a predictive and interpretative basis for structure/property relationships in liquid crystals. The calculations are of the minimum energy configurations of dimers of interacting liquid crystal molecules, and the geometry of the dimers, relative orientation of molecular dipoles and the extent of parallel correlation are related to liquid crystal properties. A large number of structural types are examined and the results discussed in terms of dipole correlation, apolar angular correlation, chiral twist sense and transition temperatures of liquid crystals.     

Molecular dynamics simulation of proteins under high pressure: Structure,function and thermodynamics

BackgroundMolecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.Scope of reviewFirst, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.Major conclusionsRecent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.General significanceMD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.     

Computer Simulation of Liquid Methanol II. System Size Effects

Abstract

A series of molecular dynamics simulations of liquid methanol has been carried out on a supernode transputer array. Four system sizes from 125 to 512 molecules have been considered, in order to study the effect of system size on the calculated structural, orientational and dynamic properties. The dielectric constant and the dielectric relaxation time are compared with experimental data.     

On the Implementation of Friedman Boundary Conditions in Liquid Water Simulations

Abstract

We have applied the image approximation to the reaction field as suggested by H.L. Friedman [Mol. Phys., 29, 1533 (1975)] by investigating appropriate cavity sizes and system parameters for use in molecular simulations. The energy of and the structure around a central simple point charge (SPC) water molecule in a dielectric cavity was found to be in good agreement with the properties of a liquid sample. To confine the water molecules within the cavity, we introduced a short-range repulsion between a real charge and its image as the Lennard-Jones repulsive potential between oxygen atoms of the SPC potential. For a system of 65 water molecules a cavity radius of 10.45 Å is appropriate; this radius is altered to 12.00 Å for a cavity surrounding 113 molecules. The effect of the boundary is restricted to the outer-most water layer which is in contact with the dielectric continuum.     

Influence of Cooperativity on Hydrogen Bond Networks

Abstract

Cooperative effects are known to strongly affect the geometrical, energetic and vibrational properties of hydrogen bonded systems. In particular, such effects strongly favor molecular arrangements where each molecule is simultaneously a donor and an acceptor of hydrogen bonds (HBs), regardless of the chemical nature of the monomer subunits. In the particular case of water systems, it has been shown that the more a molecule is a proton donor in HBs, the more the HBs where it is a proton acceptor are reinforced. Such a property could be at the origin of the equilibrium between the two species of hydrogen bonded water molecules in liquid water (one with a strong hydrogen bonding character, and one with a weaker one), as experimentally evidenced and as a molecular dynamic study of the small (H₂O)₂₄ cluster clearly suggests.     

A Strong to Fragile Transition in a Model of Liquid Silica

Abstract

The transport properties of an ionic model for liquid silica [1] at high temperatures and pressure are investigated using molecular dynamics simulations. With increasing pressure, a clear change from “strong” to “fragile” behaviour (according to Angell's classification of glass-forming liquids) is observed, albeit only on the small viscosity range that can be explored in MD simulations. This change is related to structural changes, from an almost perfect four-fold coordination to an imperfect five or six-fold coordination.     

Isobaric and Isothermal Molecular Dynamics Simulations of Diatomic Systems

Abstract

Isobaric molecular dynamics simulations were carried out for diatomic systems using different algorithms available in the literature. Two-centered Lennard-Jones potentials with and without quadrupolar interactions were used. Thermodynamic properties obtained from the isobaric algorithms compared very well with those of an equivalent simulation in the microcanonical ensemble; however, some differences were observed when similar comparisons were carried out for dynamic properties. More specifically, the constant pressure constraint affects the translational dynamics of the system because of the non-negligible differences between the momenta and the instantaneous velocities of the molecules.

Furthermore, the following studies were carried out using isobaric MD simulations: 1. Low temperature spontaneous FCC-orthorhombic (and vice versa) transition of a diatomic system with quadrupolar interactions as a function of the molecular bond length. 2. Effect of quadrupolar interaction on isobaric melting of a model diatomic system. 3. Effect of pressure on melting properties of a model diatomic system with quadrupolar interactions.     

An integrative method for evaluating the biological effects of nanoparticle-protein corona

BackgroundNanoplastics in the environment can enter the human body through gastrointestinal intake, dermal contact, and pulmonary inhalation, posing a threat to human health. Protein molecules in body fluids will quickly adsorb on the surfaces of the nanoplastics, forming a protein corona, which has implications for the interaction of the nanoplastics with cells and the metabolic pathways of the nanoplastic within cells. For years, practical tools such as dynamic light scattering, transmission electron microscopy, and liquid chromatography have been developed to understand the protein corona of nanoparticles (NPs), either in vitro or in cellular or molecular level. However, an integrated approach to understand the nanoparticles-protein corona is still lacking.MethodsUsing the most frequently observed environmental nanoplastics, polystyrene nanoplastics (PS), as a standard, we established an integrative structural characterization platform, a biophysical and biochemical evaluation method to investigate the effect of surface charge on protein corona composition. The cellular and molecular mechanisms were also explored through in vitro cellular experiments.ResultsThe first integrative method for characterizing biological properties of NPs-protein corona has been established. This method comprehensively covers the critical aspects to understand NPs-protein corona interactions, from structure to function.ConclusionsThe integrative method for nanoplastics microstructure characterization can be applied to the structural characterization of nanoparticles in nanoscale, which is of universal significance from in vitro characterization to cellular experiments and then to molecular mechanism studies.General significanceThis strategy has high reliability and repeatability and can be applied both in environment and nanomedicine safety assessment.     

2. A Theory for the Displacement of Proteins and Viruses with Polyethylene Glycol

ABSTRACT

The displacement action of polyethylene glycol of different molecular weights may be linked to the ability of the polymers to form coiled particles in solution. From conclusions drawn from their sedimentating properties in centrifugal fields the polyethylene glycols of low molecular weights, as expected, are less randomly coiled than those of higher molecular weight. It is suggested that protein molecules have the ability to diffuse into the coils of the polyethylene glycol from which they are excluded when the random coiling increases with increasing polymer concentration. From considerations based on the interaction of the polymer filament with the displaced particle the distribution of the substance between the coils and the intermolecular spaces may be predicted semi-quantitatively.     

From Molecular to Process Simulation: Novel Approaches to the Prediction of Phase Equilibria and PVT Behavior Based on Molecular/Quantum Mechanics and Molecular Dynamics Simulations

Abstract

Actually, in modern process simulators, more than 75% of the code implemented is dedicated to physical properties estimation, calculation and predictions. Data banks storing pure component parameters and binary interaction parameters for phase equilibrium calculations are extensively used and continuously implemented in actual process simulators. This gives an idea of the important role physical properties availability plays in process simulation.

In this paper we propose a new way for coupling molecular and process simulation. The basic machinery is to resort to molecular/quantum mechanics and molecular dynamics simulation techniques for generating the parameters of some equations of state that will subsequently be used for the prediction of phase equilibria and PVT behavior of small and polymeric molecules as well. This information, in turn, will be used as input in the process simulator, thus creating a final and well-defined bridge between molecular and process simulations in chemical engineering.     

Study on properties of liquid ammonia via molecular dynamics simulation based on ABEEMσπ polarisable force field

Abstract

The molecular dynamics simulation has been performed to investigate the charge distribution, structural and dynamical properties of liquid ammonia at 273 K using a polarisable force field of the atom-bond electronegativity equalisation method (ABEEMσπ). One ammonia molecule in this model has eight charge sites, one N atomic site, three H atomic sites, three N–H bond sites and one lone-pair electron site. ABEEMσπ model can present the quantitative site charges of molecular ammonias in liquid and their changing in response to their surroundings. The radial distribution functions and dynamical properties are in fair agreement with the available experimental data. The first peak of g_NN(r) appears at N–N distance of ~3.50 ± 0.05 Å where most hydrogen bonds are formed. The average coordination number of the first shell is 13.0 ± 0.1 among which a central ammonia molecule intimately connects 3 ~ 4 ammonia molecules by hydrogen bonds. The power spectrum shows the vibrations of hydrogen bonds. For a reference, a simple estimation of the average hydrogen bonding energy in liquid ammonia is 6.5 ± 0.1 kcal/mol larger than 3.8 ± 0.3 kcal/mol in dimer ammonia. Our simulation results provide more detailed information about liquid ammonia.     

The Calculation of 3D Solubility Parameters Using Molecular Models

Abstract

In this paper we describe the use of molecular mechanics models to examine detailed intermolecular interactions within the liquid state of a common nonionic surfactant system, nonyl phenol ethoxylate (NPE). Using constant energy molecular dynamics simulations we have studied the relative strengths of dispersive interactions versus polar interactions and have estimated three dimensional solubility parameters for NPE systems as a function of temperature and ethylene oxide content. The predictions at 300 K are in good agreement with three dimensional solubility parameters predicted using group contribution tables. Models of the amorphous liquid state were represented by single molecular structures of NPE in a periodic cell. The solubility parameters predicted with these models were in good agreement with those values derived from models having eight NPE molecules packed into a cell with the exception of the electrostatic interactions, which are the most sensitive to system size effects.     

The nanoscopic molecular pathway through human skin

BackgroundKnowledge regarding the barrier properties of human skin is important for understanding skin pathology, developing of transdermal drug delivery systems and computational skin absorption models; however, the molecular pathways through human skin remains to be fully investigated on a nanoscopic level. In particular the nanoscopic pathway of molecules passing the intercellular lipid bilayers separating the corneocytes in the stratum corneum (SC) is not fully elucidated.MethodsUsing stimulated emission depletion microscopy (STED) and Förster resonance energy transfer (FRET) the molecular pathways through the SC, the main barrier of the skin, are determined for lipophilic and water-soluble molecules at a nanoscopic resolution.ResultsUsing STED and confocal microscopy, water-soluble dyes, were observed to be present in both the corneocytes and in the intercellular lipid matrix, whereas the lipophilic dyes were predominately in the intercellular lipid bilayers. FRET was observed in the SC between the lipophilic and water-soluble dyes, the existence of a minimum possible distance between acceptor and donor molecules of 4.0 ± 0.1 nm was found.ConclusionsThe results indicate that lipophilic molecules penetrate the stratum corneum via the intercellular lipids bilayers separating the corneocytes in the SC, while the more water-soluble molecules penetrate the stratum corneum via the transcellular route through the corneocytes and intercellular lipid bilayers via the polar head groups of lipid molecules in the bilayers.General significanceKnowledge of the nanoscopic molecular pathways through human skin will help understand the skin barrier function and will be of use for computational skin absorption models and transdermal drug delivery strategies.     

The Interactions of Enzyme and Chemical Probes with Inverted Repeats in Supercoiled DNA

Abstract

In negatively supercoiled DNA molecules some inverted repeat sequences adopt a perturbed conformation which is characterised by the following properties. They are centrally hypersensitive to single-strand-specific nucleases such as SI, and to a much lower extent the flanking regions may also be sensitive. They are also hypersensitive to modification by bromoacetaldehyde, particularly in their flanking region. They may be resistant to endo- nucleolysis by restriction enzymes and are cleaved (resolved) by a T4 resolving enzyme. All these properties can only be consistently explained by a model in which the inverted repeat adopts a cruciform structure. This property has been shown to depend sharply on a superhelix density, and the transition to nuclease sensitivity is accompanied by a marked alteration in the overall molecular geometry as judged by frictional properties. The probable dynamics of these structures are discussed.     

Computational study about the derivatives of pyrrole as high-energy-density compounds

ABSTRACT

A series of nitro-derivatives of pyrrole were designed by replacing the hydrogen atoms on the pyrrole ring with nitro system. In order to investigate the thermodynamic and kinetic stability of these molecules, the enthalpy of formation, bond dissociation energy and bond order are calculated. The results show that most of the molecules we designed have sufficient thermodynamic and kinetic stability. Furthermore, the determination factors are confirmed in detail for the thermal stability of title molecules. The detonation velocity and detonation pressure of these molecules have been calculated by Kamlet–Jacobs equation. The calculated results show that these molecules have excellent detonation properties. Considering the stability and detonation properties, four nitro-derivatives of pyrrole (1,2,3-trinitro-1H-pyrrole, 2,3,4-trinitro-1H-pyrrole, 2,3,5-trinitro-1H-pyrrole and 2,3,4,5-tetranitro-1H-pyrrole) are screened out as potential high-energy-density molecules to further study.     

Phase Equilibria of Quadrupolar Fluids by Simulation in the Gibbs Ensemble

Abstract

Vapour-liquid phase diagrams for pure fluids and mixtures of molecules with Lennard-Jones plus quadrupole-quadrupole interaction potentials were determined by Monte Carlo simulation in the Gibbs ensemble [1]. This is the first reported application of the method to molecular fluids. We have demonstrated that the Gibbs method works reliably for strongly interacting molecular fluids at liquid densities. Pure fluid calculations were performed for reduced quadrupole strengths, Q* = Q/(εσ⁵)^1/2 equal to 1 and √2, typical of molecules like C₂H₂ and C₂H₄. It was found that the critical temperature of the quadrupolar fluid increased rapidly with increasing quadrupolar strength, in good agreement with previous computer simulation and theoretical results. A single mixture with components characterized by identical Lennard-Jones parameters and Q*₁ = + 1, Q*₂ = - 1 was studied at three temperatures. A negative azeotrope was observed at the lowest temperature studied, as seen experimentally in the CO₂/C₂H₂ mixture. The perturbation theory calculations are in good agreement with the simulation results for all properties except coexisting liquid densities. The results illustrate some of the strengths and limitations of perturbation theories based on the Padé approximant for the free energy of polar fluids.     

Evaporation and Condensation at a Liquid Surface of Methanol

Abstract

Evaporation and condensation processes at a liquid surface of methanol were investigated at room temperature with a microcanonical molecular dynamics computer simulation technique. The condensation coefficient (the number ratio of condensed molecules to incident ones) was estimated by comparing two types of autocorrelation functions, and found to be less than unity, which is in qualitative agreement with experiments. A variety of complex dynamic phenomena were observed at the surface.     

Systolic Loop Methods for Molecular Dynamics Simulation,Generalised for Macromolecules

Abstract

Systolic loop programs have been shown to be very efficient for molecular dynamics simulations of liquid systems on distributed memory parallel computers. The original methods address the case where the number of molecules simulated exceeds the number of processors used. Simulations of large flexible molecules often do not meet this condition, requiring the three- and four-body terms used to model chemical bonds within a molecule to be distributed over several processors. This paper discusses how the systolic loop methods may be generalised to accommodate such systems, and describes the implementation of a computer program for simulation of protein dynamics. Performance figures are given for this program running typical simulations on a Meiko Computing Surface using different number of processors.     

Application of ITC in foods: A powerful tool for understanding the gastrointestinal fate of lipophilic compounds

BackgroundIsothermal titration calorimetry (ITC) is a biophysical technique widely used to study molecular interactions in biological and non-biological systems. It can provide important information about molecular interactions (such as binding constant, number of binding sites, free energy, enthalpy, and entropy) simply by measuring the heat absorbed or released during an interaction between two liquid solutions.Scope of the reviewIn this review, we present an overview of ITC applications in food science, with particular focus on understanding the fate of lipids within the human gastrointestinal tract. In this area, ITC can be used to study micellization of bile salts, inclusion complex formation, the interaction of surface-active molecules with proteins, carbohydrates and lipids, and the interactions of lipid droplets.Major conclusionsITC is an extremely powerful tool for measuring molecular interactions in food systems, and can provide valuable information about many types of interactions involving food components such as proteins, carbohydrates, lipids, surfactants, and minerals. For systems at equilibrium, ITC can provide fundamental thermodynamic parameters that can be used to establish the physiochemical origin of molecular interactions.General significanceIt is expected that ITC will continue to be utilized as a means of providing fundamental information about complex materials such as those found in foods. This knowledge may be used to create functional foods designed to behave in the gastrointestinal tract in a manner that will improve human health and well-being. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.