James T. Kindt

DNA packaging and ejection forces in bacteriophage

We calculate the forces required to package (or, equivalently, acting to eject) DNA into (from) a bacteriophage capsid, as a function of the loaded (ejected) length, under conditions for which the DNA is either self-repelling or self-attracting. Through computer simulation and analytical theory, we find the loading force to increase more than 10-fold (to tens of piconewtons) during the final third of the loading process; correspondingly, the internal pressure drops 10-fold to a few atmospheres (matching the osmotic pressure in the cell) upon ejection of just a small fraction of the phage genome. We also determine an evolution of the arrangement of packaged DNA from toroidal to spool-like structures.

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 capsids 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.

Far-infrared absorption spectra of water, ammonia, and chloroform calculated from instantaneous normal mode theory

Instantaneous normal mode (INM) theory was used to calculate absolute far-infrared absorption spectra of water, ammonia, and chloroform. Three procedures for weighting the INM density of states to yield absorption intensities were tested against spectra based on dipole time correlation functions generated from molecular dynamics (MD) simulations. Weighting method I, which utilizes only the rotational character of a mode in determining its contribution to absorption, performed slightly better than method II, a more exact treatment which incorporates the extent to which a mode is IR-active. Method III, which includes the contributions of induced dipoles, was successful in describing the influence of induced dipoles on the far-infrared spectra of the model liquids. The contribution to absorption of unstable modes with imaginary frequencies was found to be significant at low frequencies, and was treated by a simple approximation. Agreement between INM theory and the MD analysis was quite good for chloroform a...

Electron-hole pair contributions to scattering, sticking, and surface diffusion: CO on Cu(100)

To assess the importance of coupling to electron-hole pair (ehp) excitations for molecular sticking, scattering, and diffusion dynamics at metal surfaces, simulations of the CO/Cu(100) system were performed using the “molecular dynamics with electronic frictions” method. Over a range of incident translational energies, energy losses to ehp excitations produce a moderate increase in sticking probability and account for 5%–10% of initial translational energy in scattered molecules, significantly less than phonon losses. Vibrational excitation and deexcitation of scattered molecules, while remaining a minor pathway for energy flow, is strongly affected by the inclusion of ehp excitations. Finally, although equilibrium diffusion constants are unaffected by the inclusion of coupling to ehp, it causes a significant quenching of transient mobility following adsorption of translationally hot molecules.

Theory for determination of the low-frequency time-dependent response function in liquids using time-resolved terahertz pulse spectroscopy

With the current rapid improvement in techniques for generating and detecting femtosecond terahertz (fs-THz) pulses, it is now possible to detect time-dependent perturbations to low-frequency intermolecular far-infrared modes during solvation events. The interpretation of these time-resolved terahertz spectroscopy experiments has generally relied on a collapse of the information available in the fs-THz pulse to give a one-dimensional time-dependent function, whose physical significance is not well defined. We describe a method to exploit the full electric field profile of the pulse as a function of pump/probe delay time to obtain a two-dimensional function that describes the system’s time-dependent dielectric response. We discuss the physical significance of this response function in terms of the time evolution of solute/solvent modes following photoexcitation, and give examples to show how its qualitative features relate to the measurable signal.

Atomistic Simulation of Cholesterol Effects on Miscibility of Saturated and Unsaturated Phospholipids: Implications for Liquid-Ordered/Liquid-Disordered Phase Coexistence

Mixed MD/MC simulation at fixed difference in chemical potential (Δμ) between two lipid types provides a computational indicator of the relative affinities of the two lipids for different environments. Applying this technique to ternary DPPC/DOPC/cholesterol bilayers yields a DPPC/DOPC ratio that increases with increasing cholesterol content at fixed Δμ, consistent with the known enrichment of DPPC and cholesterol-rich in liquid-ordered phase domains in the fluid-fluid coexistence region of the ternary phase diagram. Comparison of the cholesterol-dependence of PC compositions at constant Δμ with experimentally measured coexistence tie line end point compositions affords a direct test of the faithfulness of the atomistic model to experimental phase behavior. DPPC/DOPC ratios show little or no dependence on cholesterol content at or below 16% cholesterol in the DOPC-rich region of the composition diagram, indicating cooperativity in the favorable interaction between DPPC and cholesterol. The relative affinity of DPPC and DOPC for high cholesterol bilayer environments in simulations is explicitly shown to depend on the degree of cholesterol alignment with the bilayer normal, suggesting that a source of the cooperativity is the composition dependence of cholesterol tilt angle distributions.

Bilayer edge and curvature effects on partitioning of lipids by tail length: atomistic simulations.

The partitioning of lipids among different microenvironments in a bilayer is of considerable relevance to characterization of composition variations in biomembranes. Atomistic simulation has been ill-suited to model equilibrated lipid mixtures because the time required for diffusive exchange of lipids among microenvironments exceeds typical submicrosecond molecular dynamics trajectories. A method to facilitate local composition fluctuations, using Monte Carlo mutations to change lipid structures within the semigrand-canonical ensemble (at a fixed difference in component chemical potentials, Deltamu), was recently implemented to address this challenge. This technique was applied here to mixtures of dimyristoylphosphatidylcholine and a shorter-tail lipid, either symmetric (didecanoylphosphatidylcholine (DDPC)) or asymmetric (hexanoyl-myristoylphosphatidylcholine), arranged in two types of structure: bilayer ribbons and buckled bilayers. In ribbons, the shorter-tail component showed a clear enrichment at the highly curved rim, more so for hexanoyl-myristoylphosphatidylcholine than for DDPC. Results on buckled bilayers were variable. Overall, the DDPC content of buckled bilayers tended to exceed by several percent the DDPC content of flat ones simulated at the same Deltamu, but only for mixtures with low overall DDPC content. Within the buckled bilayer structure, no correlation could be resolved between the sign or magnitude of the local curvature of a leaflet and the mean local lipid composition. Results are discussed in terms of packing constraints, surface area/volume ratios, and curvature elasticity.

Mechanism of structural transformations induced by antimicrobial peptides in lipid membranes

It has long been suggested that pore formation is responsible for the increase in membrane permeability by antimicrobial peptides (AMPs). To better understand the mechanism of AMP activity, the disruption of model membrane by protegrin-1 (PG-1), a cationic antimicrobial peptide, was studied using atomic force microscopy. We present here the direct visualization of the full range of structural transformations in supported lipid bilayer patches induced by PG-1 on zwitterionic 1,2-dimyristoyl-snglycero-phospho-choline (DMPC) membranes. When PG-1 is added to DMPC, the peptide first induces edge instability at low concentrations, then pore-like surface defects at intermediate concentrations, and finally wormlike structures with a specific length scale at high concentrations. The formation of these structures can be understood using a mesophase framework of a binary mixture of lipids and peptides, where PG-1 acts as a line-active agent. Atomistic molecular dynamics simulations on lipid bilayer ribbons with PG-1 molecules placed at the edge or interior positions are carried out to calculate the effect of PG-1 in reducing line tension. Further investigation of the placement of PG-1 and its association with defects in the bilayer is carried out using unbiased assembly of a PG-1 containing bilayer from a random mixture of PG-1, DMPC, and water. A generalized model of AMP induced structural transformations is also presented in this work. This article is part of a Special Issue entitled: Membrane protein structure and function.

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.

Quantum dots encapsulated within phospholipid membranes: phase-dependent structure, photostability, and site-selective functionalization.

Lipid vesicle encapsulation is an efficient approach to transfer quantum dots (QDs) into aqueous solutions, which is important for renewable energy applications and biological imaging. However, little is known about the molecular organization at the interface between a QD and lipid membrane. To address this issue, we investigated the properties of 3.0 nm CdSe QDs encapsulated within phospholipid membranes displaying a range of phase transition temperatures (Tm). Theoretical and experimental results indicate that the QD locally alters membrane structure, and in turn, the physical state (phase) of the membrane controls the optical and chemical properties of the QDs. Using photoluminescence, ICP-MS, optical microscopy, and ligand exchange studies, we found that the Tm of the membrane controls optical and chemical properties of lipid vesicle-embedded QDs. Importantly, QDs encapsulated within gel-phase membranes were ultrastable, providing the most photostable non-core/shell QDs in aqueous solution reported to date. Atomistic molecular dynamics simulations support these observations and indicate that membranes are locally disordered displaying greater disordered organization near the particle-solution interface. Using this asymmetry in membrane organization near the particle, we identify a new approach for site-selective modification of QDs by specifically functionalizing the QD surface facing the outer lipid leaflet to generate gold nanoparticle-QD assemblies programmed by Watson-Crick base-pairing.