Aggregation of an Amyloidogenic Peptide on Gold Surfaces

Solid surfaces have been shown to affect the aggregation and assembly of many biomolecular systems. One important example is the formation of protein fibrils, which can occur on a range of biological and synthetic surfaces. The rate of fibrillation depends on both the protein structure and the surface chemistry, with the different molecular and oligomer structures adopted by proteins on surfaces likely to be crucial. In this paper, the aggregation of the model amyloidogenic peptide, Aβ(16-22), corresponding to a hydrophobic segment of the amyloid beta protein on a gold surface is studied using molecular dynamics simulation. Previous simulations of this peptide on gold surfaces have shown that it adopts conformations on surfaces that are quite different from those in bulk solution. These simulations show that this then leads to significant differences in the oligomer structures formed in solution and on gold surfaces. In particular, oligomers formed on the surface are low in beta-strands so are unlike the structures formed in bulk solution. When oligomers formed in solution adsorb onto gold surfaces they can then restructure themselves. This can then help explain the inhibition of Aβ(16-22) fibrillation by gold surfaces and nanoparticles seen experimentally.

Effect of Surface Hydrophobicity on the Adsorption of a Pilus-Derived Adhesin-like Peptide

Bacterial colonization of abiotic surfaces such as those of medical implants, membrane filters, and everyday household items is a process of tremendous importance for public health. Bacteria use adhesive cell surface structures called adhesins to establish contact with abiotic surfaces. Among them, protein filaments called type IV pili are particularly important and found in many Gram-negative pathogens such as Pseudomonas aeruginosa. Understanding the interaction of such adhesin proteins with different abiotic surfaces at the molecular level thus represents a fundamental prerequisite for impeding bacterial colonization and preventing the spread of infectious diseases. In this work, we investigate the interaction of a synthetic adhesin-like peptide, PAK128-144ox, derived from the type IV pilus of P. aeruginosa with hydrophilic and hydrophobic self-assembled monolayers (SAMs). Using a combination of molecular dynamics (MD) simulations, quartz crystal microbalance with dissipation monitoring (QCM-D), and spectroscopic investigations, we find that PAK128-144ox has a higher affinity for hydrophobic than for hydrophilic surfaces. Additionally, PAK128-144ox adsorption on the hydrophobic SAM is furthermore accompanied by a strong increase in α-helix content. Our results show a clear influence of surface hydrophobicity and further indicate that PAK128-144ox adsorption on the hydrophobic surface is enthalpically favored, while on the hydrophilic surface, entropic contributions are more significant. However, our spectroscopic investigations also suggest aggregation of the peptide under the employed experimental conditions, which is not considered in the MD simulations and should be addressed in more detail in future studies.

No-show rates to a sleep clinic: drivers and determinants

STUDY OBJECTIVES

Attendance to sleep clinic appointments is imperative to diagnose sleep-related disorders and to offer appropriate treatment. As part of our quality assurance program, we assessed predictors of no-show rates at our sleep clinic. We hypothesize that no-show rates can be predicted by demographics, appointment type (new vs established) and timing, and insurance status.

METHODS

We performed a 10-month, retrospective chart review of patients scheduled at Saint Louis University's SLUCare Sleep Disorders Center. Multivariable logistic regression was used to determine which factors were independently associated with no-show.

RESULTS

A total of 2,532 clinical visits were reviewed, and the overall no-show rate was 21.2%. Factors associated with a higher incidence of no-show rates included younger age (17-40 years: 21.5%; 41-64 years: 23.5%; ≥65 years: 14.0%; P < .0001), appointment type (new: 30.5% vs established: 18.3%; P < .0001), and insurance status (no insurance: 24.6% vs public: 22.6% vs private: 15.9%; P < .0001). Multivariable logistic regression confirmed the independent association between no-show and age ≤ 40 years (adjusted odds ratio = 1.72; 95% confidence interval: 1.44, 2.20), new patient status (adjusted odds ratio = 1.78; 95% confidence interval: 1.44, 2.20), and absence of health insurance (adjusted odds ratio = 1.62; 95% confidence interval: 1.24, 2.11). Sex, appointment time, day of the week, and season did not significantly influence no-show rates.

CONCLUSIONS

Independent predictors of no-show appointments included younger age, new patient status, and lack of health insurance. Our findings will aid future efforts to identify patients with high predictors of nonadherence. Further studies are needed to develop methods to decrease no-show rates once high-risk appointments have been identified.

Measuring the Impact of PEGylation on a Protein-Polysaccharide Interaction

PEGylation is the most widely used half-life extension strategy for protein therapeutics. While it imparts a range of attractive attributes PEGylation can impede protein binding and reduce efficacy. A model system to probe the effects of PEGylation on protein binding has practical applications. Here, we present a system based on complex formation between a hexavalent lectin (RSL) and the globular polysaccharide Ficoll PM70 (a type of glycocluster). Mutants of the lectin were used to generate conjugates with 3, 6, or 12 PEG (1 kDa) chains. Using NMR spectroscopy we monitored how the degree of PEGylation impacted the lectin-Ficoll interaction. The binding propensity was observed to decrease with increasing polymer density. Apparently, the extended PEG chains sterically impede the lectin-Ficoll binding. This deduction was supported by molecular dynamics simulations of the protein-polymer conjugates. The implications for protein-surface interactions are discussed.

Atomistic Study of Zwitterionic Peptoid Antifouling Brushes

Using molecular dynamics (MD) simulations, we study the molecular behavior and hydration properties of a set of zwitterionic "peptoid" brushes, grafted on a rutile surface, that has been previously reported to exhibit excellent resistance against protein adsorption and cell attachment. Peptoids are novel poly( N-substituted glycine) peptide mimics with the side chains attached to amide nitrogens. They constitute a unique model polymer system because hundreds of side chains have been demonstrated, and the exact chain length and sequence order of the residues/monomers may be specified in experiments. In this report, we vary the brush grafting density as well as the side chain/polymer molecular volume. We include in our study polysarcosine as an uncharged comparison with a small polymer chain cross-section. Sarcosine is the simplest peptoid residue with only a nominally hydrophobic methyl group as side chain, but is also reported to exhibit high antifouling performance. Overall, we show in detail how molecular volume and hydration effects are intertwined in a zwitterionic polymer brush. For example, the zwitterionic design significantly promotes extended chain conformations and could actually lower the overall electrostatic potential. Some properties promoted by the balanced charges, such as chain flexibility and hydration, increase more prominently at "low" to "intermediate" chain densities. These and other observations should provide insight on the molecular behavior of peptoids and inform the design of zwitterionic antifouling polymer brushes.

Adsorption and conformations of lysozyme and α-lactalbumin at a water-octane interface

As proteins contain both hydrophobic and hydrophilic amino acids, they will readily adsorb onto interfaces between water and hydrophobic fluids such as oil. This adsorption normally causes changes in the protein structure, which can result in loss of protein function and irreversible adsorption, leading to the formation of protein interfacial films. While this can be advantageous in some applications (e.g., food technology), in most cases it limits our ability to exploit protein functionality at interfaces. To understand and control protein interfacial adsorption and function, it is necessary to understand the microscopic conformation of proteins at liquid interfaces. In this paper, molecular dynamics simulations are used to investigate the adsorption and conformation of two similar proteins, lysozyme and α-lactalbumin, at a water-octane interface. While they both adsorb onto the interface, α-lactalbumin does so in a specific orientation, mediated by two amphipathic helices, while lysozyme adsorbs in a non-specific manner. Using replica exchange simulations, both proteins are found to possess a number of distinct interfacial conformations, with compact states similar to the solution conformation being most common for both proteins. Decomposing the different contributions to the protein energy at oil-water interfaces suggests that conformational change for α-lactalbumin, unlike lysozyme, is driven by favourable protein-oil interactions. Revealing these differences between the factors that govern the conformational change at interfaces in otherwise similar proteins can give insight into the control of protein interfacial adsorption, aggregation, and function.

Macroscopic chiral symmetry breaking in monolayers of achiral nonconvex platelets

The fabrication of chiral structures using achiral building blocks is a fundamental problem that remains a challenge in materials science. In this work we present a molecular dynamics simulation study of nonconvex polygonal platelets, interacting via soft-repulsive interactions, that are confined in two-dimensional space. These particle models are designed to promote, even at moderate densities, a natural offset displacement between the edges of neighbouring particles. In particular we demonstrate that nonconvex platelets exhibit macroscopic chiral symmetry breaking when the symmetry of the particles equals (or is multiple of) the number of nearest neighbours in the condensed crystalline phase, corresponding to the situation of platelets with 4-, 6-, and 12-fold symmetries.

Adsorption of the natural protein surfactant Rsn-2 onto liquid interfaces

To stabilize foams, droplets and films at liquid interfaces a range of protein biosurfactants have evolved in nature. Compared to synthetic surfactants, these combine surface activity with biocompatibility and low solution aggregation. One recently studied example is Rsn-2, a component of the foam nest of the frog Engystomops pustulosus, which has been predicted to undergo a clamshell-like opening transition at the air-water interface. Using atomistic molecular dynamics simulations and surface tension measurements we study the adsorption of Rsn-2 onto air-water and cyclohexane-water interfaces. The protein adsorbs readily at both interfaces, with adsorption mediated by the hydrophobic N-terminus. At the cyclohexane-water interface the clamshell opens, due to the favourable interaction between hydrophobic residues and cyclohexane molecules and the penetration of cyclohexane molecules into the protein core. Simulations of deletion mutants showed that removal of the N-terminus inhibits interfacial adsorption, which is consistent with the surface tension measurements. Deletion of the hydrophilic C-terminus also affects adsorption, suggesting that this plays a role in orienting the protein at the interface. The characterisation of the interfacial behaviour gives insight into the factors that control the interfacial adsorption of proteins, which may inform new applications of this and similar proteins in areas including drug delivery and food technology and may also be used in the design of synthetic molecules showing similar changes in conformation at interfaces.

Conformations of Myoglobin-Derived Peptides at the Air-Water Interface

The conformational change exhibited by proteins at liquid interfaces, such as the air-water and oil-water interfaces, has long been of interest both for understanding protein structure outside of native environments and for applications in areas including food technology and pharmaceuticals. Using molecular simulation, this article studies the conformations of two peptides derived from myoglobin, for which the emulsification behavior has been studied. Both peptides were found to readily adsorb onto the air-water interface, with one of these (experimentally, the more effective stabilizer) adopting a flat, extended conformation and the other peptide remaining close to its solution conformation.

Phase behaviour of self-assembled monolayers controlled by tuning physisorbed and chemisorbed states: A lattice-model view

The self-assembly of molecules on surfaces into 2D structures is important for the bottom-up fabrication of functional nanomaterials, and the self-assembled structure depends on the interplay between molecule-molecule interactions and molecule-surface interactions. Halogenated benzene derivatives on platinum have been shown to have two distinct adsorption states: a physisorbed state and a chemisorbed state, and the interplay between the two can be expected to have a profound effect on the self-assembly and phase behaviour of these systems. We developed a lattice model that explicitly includes both adsorption states, with representative interactions parameterised using density functional theory calculations. This model was used in Monte Carlo simulations to investigate pattern formation of hexahalogenated benzene molecules on the platinum surface. Molecules that prefer the physisorbed state were found to self-assemble with ease, depending on the interactions between physisorbed molecules. In contrast, molecules that preferentially chemisorb tend to get arrested in disordered phases. However, changing the interactions between chemisorbed and physisorbed molecules affects the phase behaviour. We propose functionalising molecules in order to tune their adsorption states, as an innovative way to control monolayer structure, leading to a promising avenue for directed assembly of novel 2D structures.

Scaling Behavior of Polymers at Liquid/Liquid Interfaces

The dynamics of a polymer chain confined in a soft 2D slit formed by two immiscible liquids is studied by means of molecular dynamics simulations. We show that the scaling behavior of a polymer confined between two liquids does not follow that predicted for polymers adsorbed on solid or soft surfaces such as lipid bilayers. Indeed, our results show that in the diffusive regime the polymer behaves like in bulk solution, following the Zimm model, and with the hydrodynamic interactions dominating its dynamics. Although the presence of the interface does not affect the long-time diffusion properties, it has an influence on the dynamics at short time scale, where for low molecular weight polymers the subdiffusive regime almost disappears. Simulations carried out when the liquid interface is sandwiched between two solid walls show that, when the confinement is a few times larger than the blob size, the Rouse dynamics is recovered.

Aggregation of nanoparticles on one and two-component bilayer membranes

Using dissipative particle dynamics simulations the aggregation of nanoparticles on single and two-component bilayers is investigated. For a uniform bilayer the aggregation of nanoparticles depends strongly on the location of the particles in the bilayer; particles residing on the bilayer exterior cluster strongly under the influence of bilayer-mediated interactions, whereas the interaction between the particles in the bilayer interior is significantly weaker leading to more loosely bound, dynamic aggregates. The aggregation of nanoparticles on two-component bilayers composed of immiscible components changes due to competition between nanoparticle clustering and their adsorption on the boundary between the bilayer components. This reduces the size of the nanoparticle clusters formed on the bilayer exterior, with the clusters adhering onto the boundary between the bilayer components. Due to their weaker attraction nanoparticles in the interior of a mixed bilayer no longer aggregate and instead form strings along the boundary between the two bilayer components. Nanoparticles with an affinity to one bilayer component nucleate small domains of their favoured component around themselves. For asymmetric mixtures this leads to a notable change in the aggregation behaviour of the nanoparticles.

Mechanochemical Tuning of the Binaphthyl Conformation at the Air-Water Interface

Gradual and reversible tuning of the torsion angle of an amphiphilic chiral binaphthyl, from -90° to -80°, was achieved by application of a mechanical force to its molecular monolayer at the air-water interface. This 2D interface was an ideal location for mechanochemistry for molecular tuning and its experimental and theoretical analysis, since this lowered dimension enables high orientation of molecules and large variation in the area. A small mechanical energy (<1 kcal mol(-1) ) was applied to the monolayer, causing a large variation (>50 %) in the area of the monolayer and modification of binaphthyl conformation. Single-molecule simulations revealed that mechanical energy was converted proportionally to torsional energy. Molecular dynamics simulations of the monolayer indicated that the global average torsion angle of a monolayer was gradually shifted.

Thermodynamics of linear and star polymers at fluid interfaces

Performing molecular dynamics simulations on model systems we study the structural changes and thermodynamic stability of polymers of varying topology (linear and star-shaped) at interface between two liquids. We find that homopolymers are attracted to the interface in both good and poor solvent conditions showing that they are surface active molecules even though not amphiphilic. In most cases changing polymer topology had only a minor effect on the desorption free energy. A noticeable dependence on polymer topology is only seen for relatively high molecular weight polymers at interface between two good solvents. Examining separately the enthalpic and entropic components of the desorption free energy suggests that its largest contribution is the decrease in the enthalpic part of interfacial free energy caused by the adsorption of the polymer at the interface. Finally we propose a simple method to qualitatively predict the trend of the interfacial free energy as a function of the polymer molecular weight.

Conformations and effective interactions of polymer-coated nanoparticles at liquid interfaces

We investigate conformations and effective interactions of polymer-coated nanoparticles adsorbed at a model liquid-liquid interface via molecular dynamics simulations. The polymer shells strongly deform at the interface, with the shape governed by a balance between maximizing the decrease in interfacial area between the two solvent components, minimizing unfavorable contact between polymer and solvent, and maximizing the conformational entropy of the polymers. Using potential of mean force calculations, we compute the effective interaction between the nanoparticles at the liquid-liquid interface. We find that it differs quantitatively from the bulk and is significantly affected by the length of the polymer chains and by the solvent quality. Under good solvent conditions, the effective interactions are always repulsive and soft for long chains. The repulsion range decreases as the solvent quality decreases. In particular, under poor solvent conditions, short chains may fail to induce steric repulsion, leading to a net attraction between the nanoparticles, whereas with long-enough chains the effective interaction potential may feature an additional repulsive shoulder at intermediate distances.

Molecular simulation of hydrophobin adsorption at an oil-water interface

Hydrophobins are small, amphiphilic proteins expressed by strains of filamentous fungi. They fulfill a number of biological functions, often related to adsorption at hydrophobic interfaces, and have been investigated for a number of applications in materials science and biotechnology. In order to understand the biological function and applications of these proteins, a microscopic picture of the adsorption of these proteins at interfaces is needed. Using molecular dynamics simulations with a chemically detailed coarse-grained potential, the behavior of typical hydrophobins at the water-octane interface is studied. Calculation of the interfacial adsorption strengths indicates that the adsorption is essentially irreversible, with adsorption strengths of the order of 100 k(B)T (comparable to values determined for synthetic nanoparticles but significantly larger than small molecule surfactants and biomolecules). The protein structure at the interface is unchanged at the interface, which is consistent with the biological function of these proteins. Comparison of native proteins with pseudoproteins that consist of uniform particles shows that the surface structure of these proteins has a large effect on the interfacial adsorption strengths, as does the flexibility of the protein.

Structural variability and dynamics of the P3HT/PCBM interface and its effects on the electronic structure and the charge-transfer rates in solar cells

Using a range of realistic interface geometries obtained from a molecular dynamics simulation we study the effects of different microscopic atomic arrangements on the electronic structure and charge transfer rates of the prototypical photovoltaic interface between P3HT (poly(3-hexylthiophene)) and PCBM ([6,6]-phenyl-C(61)-butyric acid methyl ester). The electronic structures of charge-transfer (CT) states belong to two groups that can be denoted as "charge-separated" and "charge-bridging" states. For the former group of structures, which may lead to fully separated charges, the ranges and the average values of internal reorganization energy, the electronic coupling and the charge separated states energy are evaluated. A range and distribution of absolute charge separation (CS) and recombination (CR) rates are computed using the Marcus-Levich-Jortner rate equation. Due to the variety of P3HT/PCBM interface structures, a very broad range of CS (7.7 × 10(9)-1.8 × 10(12) s(-1)) and CR (2.5 × 10(5)-1.1 × 10(10) s(-1)) "instantaneous" rates are computed. However, the energetic parameters affecting the rate evolve in time due to the dynamic nature of the interface with a characteristic timescale of about 10 ns. For this reason the slowest CR instantaneous rates are not observed and the minimum CR rate observed is determined by the rate of conformational rearrangement at the interface. The combination of these observations provides a more general framework for the interpretation of experimental spectroscopic data, suggesting that the analysis based on simple first order rates may be insufficient to describe charge transfer in organic solar cell interfaces.

Charge transport in semiconductors with multiscale conformational dynamics

In partially ordered organic semiconductors, the characteristic times of nuclear motion are comparable to those of charge carrier dynamics. It is impossible to describe charge transport using either static disorder models or temperature averaged electronic Hamiltonians. We build a model Hamiltonian which allows the study of charge transport whenever carrier and nuclear dynamics are not easily separable. Performing nanoseconds long molecular dynamics of a columnar mesophase of a discotic liquid crystal and evaluating electronic couplings, we identify realistic parameters of the Hamiltonian. All modes which are coupled to the electron dynamics can be described in the model Hamiltonian by a limited number of Langevin oscillators. This method can be applied to systems with both slow (nanoseconds) and fast (hundreds of femtoseconds) nuclear motions, i.e., with both dynamic and static disorder.

Interaction of nanoparticles with ideal liquid-liquid interfaces

Using molecular simulations the interaction between a noncharged nanoparticle and an ideal liquid-liquid interface is studied. The free energy profile as function of nanoparticle-interface separation is determined using Wang-Landau sampling. Comparison between the simulation results and macroscopic theories shows that the latter give a poor description of the free energy profile. In particular, they underestimate both the range of interaction between the particle and the interface and its strength, with the discrepancy lessening as the particle radius increases. On increasing the solvent chemical potential the interaction strength increases and the interaction range decreases due to the increase in interfacial tension and consequent decrease in interfacial width.

Elastic constants of hard thin platelets by Monte Carlo simulation and virial expansion

In this paper we present an investigation into the calculation of the Frank elastic constants of hard platelets via molecular simulation and virial expansion beyond second order. Monte Carlo simulations were carried out and director fluctuations measured as a function of wave vector k, giving the elastic constants through a fit in the low-k limit. Additionally, the virial expansion coefficients of the elastic constants up to sixth order were calculated, and the validity of the theory determined by comparison with the simulation results. The simulation results are also compared with experimental measurements on colloidal suspensions of platelike particles.

Structure and stability of isotropic states of hard platelet fluids

We study the thermodynamics and the pair structure of hard, infinitely thin, circular platelets in the isotropic phase. Monte Carlo simulation results indicate a rich spatial structure of the spherical expansion components of the direct correlation function, including nonmonotonical variation of some of the components with density. Integral equation theory is shown to reproduce the main features observed in simulations. The hypernetted chain closure, as well as its extended versions that include the bridge function up to second and third order in density, perform better than both the Percus-Yevick closure and Verlet bridge function approximation. Using a recent fundamental measure density functional theory, an analytic expression for the direct correlation function is obtained as the sum of the Mayer bond and a term proportional to the density and the intersection length of two platelets. This is shown to give a reasonable estimate of the structure found in simulations, but to fail to capture the nonmonotonic variation with density. We also carry out a density functional stability analysis of the isotropic phase with respect to nematic ordering and show that the limiting density is consistent with that where the Kerr coefficient vanishes. As a reference system, we compare to simulation results for hard oblate spheroids with small, but nonzero elongations, demonstrating that the case of vanishingly thin platelets is approached smoothly.

Modelling charge transport in organic semiconductors: from quantum dynamics to soft matter

The charge carrier dynamics in organic semiconductors has been traditionally discussed with the models used in inorganic crystalline and amorphous solids but this analogy has severe limitations because of the more complicated role of nuclear motions in organic materials. In this perspective, we discuss how a new approach to the modelling of charge transport is emerging from the alliance between the conventional quantum chemical methods and the methods more traditionally used in soft-matter modelling. After describing the conventional limit cases of charge transport we discuss the problems arising from the comparison of the theory with the experimental and computational results. Several recent applications of numerical methods based on the propagation of the wavefunction or kinetic Monte Carlo methods on soft semiconducting materials are reviewed.

Forces between cylindrical nanoparticles in a liquid crystal

Using classical density functional theory, the forces between two cylindrical nanoparticles in a liquid crystal solvent are calculated. Both the nematic and isotropic phases of the solvent are considered. In the nematic phase, the interaction is highly anisotropic. At short range, changes in the defect structure around the cylinders leads to a complex interaction between them. In the isotropic phase, an attractive interaction arises due to overlap between halos of ordered fluid adsorbed on the surfaces of the cylinders.

Structure of molecular liquids: hard rod-disk mixtures

The structure of hard rod-disk mixtures is studied using Monte Carlo simulations and integral equation theory, for a range of densities in the isotropic phase. By extension of methods used in single component fluids, the pair correlation functions of the molecules are calculated and comparisons between simulation and integral equation theory, using a number of different closure relations, are made. Comparison is also made for thermodynamic data and phase behavior.

Liquid-crystal-mediated force between a cylindrical nanoparticle and substrate

Using classical density functional theory, the structure of a molecular fluid around a cylindrical nanoparticle near a solid substrate is studied. The solvent-mediated force between the nanoparticle and the substrate is calculated in both the nematic and isotropic phases of the solvent. In the nematic phase, the force is short ranged and arises due to interaction between high-density regions near the substrate and nanoparticle. In the isotropic phase, the formation of a nematic bridge between the substrate and nanoparticle gives rise to an attractive force between them. The potential between the nanoparticle and substrate as a function of separation calculated numerically is compared to that calculated from the Derjaguin approximation. In the isotropic phase these are found to be in reasonable agreement at low separations, while the agreement is poorer in the nematic phase.

Structure of molecular liquids: closure relations for hard spheroids

We present the results of Monte Carlo simulations of hard spheroids of revolution of different elongations. Both prolate and oblate shapes are examined. A systematic study of the bridge function b(1,2), and direct comparison with the indirect correlation function gamma(1,2)=h(1,2)-c(1,2) at densities spanning the isotropic fluid range, allow us to evaluate the accuracy of various proposed closure relations for integral equations.

Structure of a liquid crystalline fluid around a macroparticle: Density functional theory

The structure of a molecular liquid, in both the nematic liquid crystalline and isotropic phases, around a cylindrical macroparticle, is studied using density functional theory. In the nematic phase the structure of the fluid is highly anisotropic with respect to the director, in agreement with results from simulation and phenomenological theories. On going into the isotropic phase the structure becomes rotationally invariant around the macroparticle with an oriented layer at the surface.

Structure of molecular liquids: cavity and bridge functions of the hard spheroid fluid

We present methodologies for calculating the direct correlation function c(1,2), the cavity function y(1,2), and the bridge function b(1,2), for molecular liquids, from Monte Carlo simulations. As an example we present results for the isotropic hard spheroid fluid with elongation e = 3. The simulation data are compared with the results from integral equation theory. In particular, we solve the Percus-Yevick and hypernetted chain equations. In addition, we calculate the first two terms in the virial expansion of the bridge function and incorporate this into the closure. At low densities, the bridge functions calculated by theory and from simulation are in good agreement, lending support to the correctness of our numerical procedures. At higher densities, the hypernetted chain results are brought into closer agreement with simulation by incorporating the approximate bridge function, but significant discrepancies remain.

Monte Carlo simulations of liquid crystals near rough walls

The effect of surface roughness on the structure of liquid crystalline fluids near solid substrates is studied by Monte Carlo simulations. The liquid crystal is modeled as a fluid of soft ellipsoidal molecules and the substrate is modeled as a hard wall that excludes the centers of mass of the fluid molecules. Surface roughness is introduced by embedding a number of molecules with random positions and orientations within the wall. It is found that the density and order near the wall are reduced as the wall becomes rougher, i.e., the number of embedded molecules is increased). Anchoring coefficients are determined from fluctuations in the reciprocal space order tensor. It is found that the anchoring strength decreases with increasing surface roughness.

Calculation of flexoelectric coefficients for a nematic liquid crystal by atomistic simulation

Equilibrium molecular dynamics calculations have been performed for the liquid crystal molecule n-4-(trans-4-n-pentylcyclohexyl)benzonitrile (PCH5) using a fully atomistic model. Simulation data have been obtained for a series of temperatures in the nematic phase. The simulation data have been used to calculate the flexoelectric coefficients e(s) and e(b) using the linear response formalism of Osipov and Nemtsov [M. A. Osipov and V. B. Nemtsov, Sov. Phys. Crstallogr. 31, 125 (1986)]. The temperature and order parameter dependence of e(s) and e(b) are examined, as are separate contributions from different intermolecular interactions. Values of e(s) and e(b) calculated from simulation are consistent with those found from experiment.

A density-functional theory study of the confined soft ellipsoid fluid

A system of soft ellipsoid molecules confined between two planar walls is studied using classical density-functional theory. Both the isotropic and nematic phases are considered. The excess free energy is evaluated using two different Ansätze and the intermolecular interaction is incorporated using two different direct correlation functions (DCF's). The first is a numerical DCF obtained from simulations of bulk soft ellipsoid fluids and the second is taken from the Parsons-Lee theory. In both the isotropic and nematic phases the numerical DCF gives density and order parameter profiles in reasonable agreement with simulation. The Parsons-Lee DCF also gives reasonable agreement in the isotropic phase but poor agreement in the nematic phase.

Parametrization and validation of a force field for liquid-crystal forming molecules

First principles density functional calculations have been carried out to determine the structures and conformational energies of a series of liquid-crystal fragment molecules. The calculations have been used to derive a molecular mechanics force field that describes a subset of commonly occurring liquid-crystal molecules. The force field has been used to carry out molecular dynamics simulations of the bulk phase for these fragment molecules. Computed densities and heats of vaporization are in good agreement with experimental data. These results should be useful in future molecular dynamics simulations of liquid-crystal systems.

Coronary artery fistulas: long-term results of surgical correction

BACKGROUND

Opinions vary as to whether operation should be offered patients with coronary artery fistula, particularly to those who are asymptomatic. Published studies lacked long-term follow-up data.

METHODS

We studied 41 patients with coronary artery fistula operated in our unit in the past 30 years with restudies including coronary angiograms in those who agreed to the investigation.

RESULTS

There was no operative mortality and operative morbidity was low. The mean duration of follow-up was 9.1 years and 96.9% of the patients were asymptomatic. Twenty-one patients had a coronary angiogram. The native coronary artery either remained dilated and tortuous, or more frequently had thromboses with a short proximal stump. (None of these patients had evidence of myocardial ischemia.) Four patients had demonstrable recurrence fistula but without hemodynamic disturbance.

CONCLUSIONS

We advocate operation for all patients with coronary artery fistulas and demonstrable shunting in view of minimal operative risks. Small asymptomatic fistulas without demonstrable shunting should be left alone. The relatively high incidence of residual or recurrent fistula makes long-term follow-up mandatory.

Surgical management of blunt traumatic rupture of the descending thoracic aorta

Acute rupture of the descending thoracic aorta following blunt trauma is a life-threatening injury that requires emergent operative intervention. From February 1989 to January 1997, 4 patients with multiple injuries including traumatic rupture in the region of the aortic isthmus were surgically treated at our institution. Diagnosis was confirmed in all patients by aortogram prior to aortic repair. One patient had preoperative paraplegia secondary to an unstable spinal fracture. All patients underwent repair under cardiopulmonary bypass (3 partial, 1 total with hypothermic arrest). The site of rupture was resected and replaced with an interposition artificial graft. There was no perioperative mortality and no additional occurrence of paraplegia. Our experience and a review of the literature indicate that for survivors of traumatic aortic rupture, excellent outcomes can be achieved only if the diagnosis is made early and the surgical treatment is prompt.

A non-fatal impalement injury of the thorax

Impalement is an uncommon injury with only occasional reports in the literature. There are even fewer reports of impalement injuries limited to the thorax. We report herein the case of a 24-year-old man who survived impalement injury of the left side of the thorax with a steel rod while working at a construction site. The great vessels of the thorax were spared but the second thoracic vertebra was fractured resulting in complete paralysis of the left lower limb. The precise nature and extent of the injury were determined pre-operatively by computed tomography and aortography. The important principles of surgical management contributing to the successful outcome are described, these being minimal manipulation of the impalement object before and during transport, careful pre-operative planning and a multidisciplinary approach.

Smooth muscle tumors of the inferior vena cava and right heart

Four Chinese patients with extensive smooth muscle tumors of the inferior vena cava (IVC) and right heart are described. Two patients had primary IVC leiomyosarcomas and one had both IVC and cardiac metastases from uterine leiomyosarcoma. The remaining patient had extra-uterine intravenous leiomyomatosis, with extensive tumor calcification. Clinical presentation was variable, consisting of the Budd-Chiari syndrome, a new heart murmur, and signs of cerebral embolism. All the lesions were defined radiologically prior to surgical resection of cardiac tumors.