Biomedical subjects

James T Kindt

Publications and source records attributed to James T Kindt.

8 recordsLinked to original sources

PubMed2006

Equilibrium distributions of dipalmitoyl phosphatidylcholine and dilauroyl phosphatidylcholine in a mixed lipid bilayer: atomistic semigrand canonical ensemble simulations.

Conventional molecular dynamics (MD) simulations are seriously limited by the slow rate of diffusive mixing in their ability to predict lateral distributions of different lipid types within mixed-lipid bilayers using atomistic models. A method to overcome this limitation, using configuration-bias Monte Carlo (MC) "mutation" moves to transform lipids from one type to another in dynamic equilibrium, is demonstrated in binary fluid-phase mixtures of lipids whose tails differ in length by four carbons. The hybrid MC-MD method operates within a semigrand canonical ensemble, so that an equilibrium composition of the mixture is determined by a constant difference in chemical potential (Delta(mu)) chosen for the components. Within several nanoseconds, bilayer structures initiated as pure dipalmitoyl phosphatidylcholine (DPPC) or pure dilauroyl phosphatidylcholine (DLPC) converge to a common composition and structure in independent simulations conducted at the same Delta(mu). Trends in bilayer thickness, area per lipid, density distributions across the bilayer, and order parameters have been investigated at three mixture compositions and compared with results from the pure bilayers at 323 K. The mixtures showed a moderate increase in DPPC acyl tail sites crossing the bilayer midplane relative to pure DPPC. Correlations between lateral positions of the two lipid types within or across the bilayer were found to be weak or absent. While the lateral distribution is consistent with nearly ideal mixing, the dependence of composition on Delta(mu) indicates a positive excess free energy of mixing.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2006

Theoretical analysis of polydispersity in the nematic phase of self-assembled semiflexible chains.

In recent simulations of semiflexible equilibrium polymers [X. Lu and J. T. Kindt, J. Chem. Phys. 120, 10328 (2004)] a roughly biexponential distribution of chain lengths was observed in the nematic phase. In this study, we show that a theory representing the nematic phase as an equilibrium mixture of randomly oriented chains below a critical length Lc and chains above Lc obeying a Gaussian orientational distribution reproduces this distribution qualitatively. The agreement between predicted and simulated phase boundaries is improved, with a narrower isotropic-nematic coexistence region, compared with a monoexponential chain distribution. We find, furthermore, that the critical length Lc scales as the inverse of monomer concentration, irrespective of the bond strength of the chains, but that at the phase boundary, the fraction of monomers found in disordered chains peaks at a certain bond strength.

Journal Article↗

PubMed2006

Coarse-grained model simulations of mixed-lipid systems: composition and line tension of a stabilized bilayer edge.

Bilayer disks and ribbons composed of a mixture of short- and long-tail phospholipids have been studied by molecular dynamics with a coarse-grained model. The effects of system composition on the edge structure, composition, and line tension were analyzed. Increases in the fraction of short-tail lipids tend to decrease the line tension (i.e., stabilize the edge) but not eliminate it. The short-tail lipid is generally enriched at the curved rim forming the bilayer edge, with an excess of 3 to 4 molecules per nanometer (relative to the bulk), but complete segregation was not observed. In all mixtures, a region depleted in the short-tail component occurs just before the edge, corresponding to a bulge in the bilayer thickness. The bulge and depletion are more prominent as the bilayer composition shifts toward a majority of short-tail lipids. In one case, a net excess of long-tail lipids at the edge was demonstrated, suggesting that certain circumstances give rise to a "segregation inversion" in which the long-tail lipid behaves as an edge stabilizer.

Journal Article↗

PubMed2006

Complementarity and clustering in a simple model mixed bilayer.

A bilayer of uniform thickness containing a mixture of long and short lipids is simulated using a parallel hard-rod model to illustrate the effect of transbilayer repulsions between the tails of the long component. Monte Carlo simulations show considerable entropy-driven clustering within each layer. Demixing reaches a maximum at the highest packing fraction of the liquid state and decreases as the system orders. The formation of complementary clusters of long and short rods on opposite sides of the bilayer increases translational freedom within each cluster by reducing constraints imposed by the opposing leaflet, an effect that becomes less important as rods lock into facing hexagonally ordered arrays.

Computer Simulation↗

PubMed2005

Simulation and theory of self-assembly and network formation in reversibly cross-linked equilibrium polymers.

A simulation model of hard spheres capable of reversible assembly into chains, which then may reversibly cross-link into networks, has been studied through grand canonical Monte Carlo simulation. Effects of varying intra- and interchain bond strengths, chain flexibilities, and restrictions on cross-linking angle were investigated. Observations including chain-length distributions and phase separation could be captured in most cases using a simple model theory. The coupling of chain growth to cross-linking was shown to be highly sensitive to the treatment of cross-linking by chain ends. In some systems, ladderlike domains of several cross-links joining two chains were common, resulting from cooperativity in the cross-linking. Extended to account for this phenomenon, the model theory predicts that such cooperativity will suppress phase separation in weakly polymerizing chains and at high cross-link concentration. In the present model, cross-linking stabilizes the isotropic phase with respect to the nematic phase, causing a shift in the isotropic-nematic transition to higher monomer concentration than in simple equilibrium polymers.

Computer Simulation↗

PubMed2004

Monte Carlo simulation of the self-assembly and phase behavior of semiflexible equilibrium polymers.

Grand canonical Monte Carlo simulations of a simple model semiflexible equilibrium polymer system, consisting of hard sphere monomers reversibly self-assembling into chains of arbitrary length, have been performed using a novel sampling method to add or remove multiple monomers during a single MC move. Systems with two different persistence lengths and a range of bond association constants have been studied. We find first-order lyotropic phase transitions between isotropic and nematic phases near the concentrations predicted by a statistical thermodynamic theory, but with significantly narrower coexistence regions. A possible contribution to the discrepancy between theory and simulation is that the length distribution of chains in the nematic phase is bi-exponential, differing from the simple exponential distribution found in the isotropic phase and predicted from a mean-field treatment of the nematic. The additional short length-scale characterizing the distribution appears to arise from the lower orientational order of short chains. The dependence of this length-scale on chemical potential, bond association constant, and total monomer concentration has been examined.

Journal Article↗

PubMed2004

Molecular dynamics simulations of the lipid bilayer edge.

Phospholipid bilayers have been intensively studied by molecular dynamics (MD) simulation in recent years. The properties of bilayer edges are important in determining the structure and stability of pores formed in vesicles and biomembranes. In this work, we use molecular dynamics simulation to investigate the structure, dynamics, and line tension of the edges of bilayer ribbons composed of pure dimyristoylphosphatidylcholine (DMPC) or palmitoyl-oleoylphosphatidylethanolamine (POPE). As expected, we observe a significant reorganization of lipids at and near the edges. The treatment of electrostatic effects is shown to have a qualitative impact on the structure and stability of the edge, and significant differences are observed in the dynamics and structure of edges formed by DMPC and palmitoyl-oleoylphosphatidylethanolamine. From the pressure anisotropy in the simulation box, we calculate a line tension of approximately 10-30 pN for the DMPC edge, in qualitative agreement with experimental estimates for similar lipids.

Anisotropy↗

PubMed2003

Forces and pressures in DNA packaging and release from viral capsids.

In a previous communication (Kindt et al., 2001) we reported preliminary results of Brownian dynamics simulation and analytical theory which address the packaging and ejection forces involving DNA in bacteriophage capsids. In the present work we provide a systematic formulation of the underlying theory, featuring the energetic and structural aspects of the strongly confined DNA. The free energy of the DNA chain is expressed as a sum of contributions from its encapsidated and released portions, each expressed as a sum of bending and interstrand energies but subjected to different boundary conditions. The equilibrium structure and energy of the capsid-confined and free chain portions are determined, for each ejected length, by variational minimization of the free energy with respect to their shape profiles and interaxial spacings. Numerical results are derived for a model system mimicking the lambda-phage. We find that the fully encapsidated genome is highly compressed and strongly bent, forming a spool-like condensate, storing enormous elastic energy. The elastic stress is rapidly released during the first stage of DNA injection, indicating the large force (tens of pico Newtons) needed to complete the (inverse) loading process. The second injection stage sets in when approximately 1/3 of the genome has been released, and the interaxial distance has nearly reached its equilibrium value (corresponding to that of a relaxed torus in solution); concomitantly the encapsidated genome begins a gradual morphological transformation from a spool to a torus. We also calculate the loading force, the average pressure on the capsid's walls, and the anisotropic pressure profile within the capsid. The results are interpreted in terms of the (competing) bending and interaction components of the packing energy, and are shown to be in good agreement with available experimental data.

Bacteriophage lambda↗