Statistically interacting vacancy particles

Free-energy landscapes and insertion pathways for peptides in membrane environment

Free-energy landscapes for short peptides-specifically for variants of the pH low insertion peptide (pHLIP)-in the heterogeneous environment of a lipid bilayer or cell membrane are constructed, taking into account a set of dominant interactions and the conformational preferences of the peptide backbone. Our methodology interprets broken internal H-bonds along the backbone of a polypeptide as statistically interacting quasiparticles, activated from the helix reference state. The favored conformation depends on the local environment (ranging from polar to nonpolar), specifically on the availability of external H-bonds (with H_{2}O molecules or lipid headgroups) to replace internal H-bonds. The dominant side-chain contribution is accounted for by residue-specific transfer free energies between polar and nonpolar environments. The free-energy landscape is sensitive to the level of pH in the aqueous environment surrounding the membrane. For high pH, we identify pathways of descending free energy that suggest a coexistence of membrane-adsorbed peptides with peptides in solution. A drop in pH raises the degree of protonation of negatively charged residues and thus increases the hydrophobicity of peptide segments near the C terminus. For low pH, we identify insertion pathways between the membrane-adsorbed state and a stable trans-membrane state with the C terminus having crossed the membrane.

Mechanical response to tension and torque of molecular chains via statistically interacting particles associated with extension, contraction, twist, and supercoiling

A methodology for the statistical mechanical analysis of polymeric chains under tension introduced previously is extended to include torque. The response of individual bonds between monomers or of entire groups of monomers to a combination of tension and torque involves, in the framework of this method of analysis, the (thermal or mechanical) activation of a specific mix of statistically interacting particles carrying quanta of extension or contraction and quanta of twist or supercoiling. The methodology, which is elucidated in applications of increasing complexity, is capable of describing the conversion between twist chirality and plectonemic chirality in quasistatic processes. The control variables are force or extension and torque or linkage (a combination of twist and writhe). The versatility of this approach is demonstrated in two applications relevant and promising for double-stranded DNA under controlled tension and torque. One application describes conformational transformations between (native) B-DNA, (underwound) S-DNA, and (overwound) P-DNA in accord with experimental data. The other application describes how the conversion between a twisted chain and a supercoiled chain accommodates variations of linkage and excess length in a buckling transition.

Molecular chains under tension: Thermal and mechanical activation of statistically interacting extension and contraction particles

This work introduces a methodology for the statistical mechanical analysis of polymeric chains under tension controlled by optical or magnetic tweezers at thermal equilibrium with an embedding fluid medium. The response of single bonds between monomers or of entire groups of monomers to tension is governed by the activation of statistically interacting particles representing quanta of extension or contraction. This method of analysis is capable of describing thermal unbending of the freely jointed or wormlike chain kind, linear or nonlinear contour elasticity, and structural transformations including effects of cooperativity. The versatility of this approach is demonstrated in an application to double-stranded DNA undergoing torsionally unconstrained stretching across three regimes of mechanical response including an overstretching transition. The three-regime force-extension characteristic, derived from a single free-energy expression, accurately matches empirical evidence.

Density profiles of a self-gravitating lattice gas in one, two, and three dimensions

We consider a lattice gas in spaces of dimensionality D=1,2,3. The particles are subject to a hardcore exclusion interaction and an attractive pair interaction that satisfies Gauss' law as do Newtonian gravity in D=3, a logarithmic potential in D=2, and a distance-independent force in D=1. Under mild additional assumptions regarding symmetry and fluctuations we investigate equilibrium states of self-gravitating material clusters, in particular radial density profiles for closed and open systems. We present exact analytic results in several instances and high-precision numerical data in others. The density profile of a cluster with finite mass is found to exhibit exponential decay in D=1 and power-law decay in D=2 with temperature-dependent exponents in both cases. In D=2 the gas evaporates in a continuous transition at a nonzero critical temperature. We describe clusters of infinite mass in D=3 with a density profile consisting of three layers (core, shell, halo) and an algebraic large-distance asymptotic decay. In D=3 a cluster of finite mass can be stabilized at T>0 via confinement to a sphere of finite radius. In some parameter regime, the gas thus enclosed undergoes a discontinuous transition between distinct density profiles. For the free energy needed to identify the equilibrium state we introduce a construction of gravitational self-energy that works in all D for the lattice gas. The decay rate of the density profile of an open cluster is shown to transform via a stretched exponential for 1<D<2, whereas it crosses over from one power-law at intermediate distances to a different power-law at larger distances for 2<D<3.

Calculation of displacement correlation tensor indicating vortical cooperative motion in two-dimensional colloidal liquids

As an indicator of cooperative motion in a system of Brownian particles that models two-dimensional colloidal liquids, a displacement correlation tensor is calculated analytically and compared with numerical results. The key idea for the analytical calculation is to relate the displacement correlation tensor, which is a kind of four-point space-time correlation, to the Lagrangian two-time correlation of the deformation gradient tensor. Tensorial treatment of the statistical quantities, including the displacement correlation itself, allows capturing the vortical structure of the cooperative motion. The calculated displacement correlation also implies a negative long-time tail in the velocity autocorrelation, which is a manifestation of the cage effect. Both the longitudinal and transverse components of the displacement correlation are found to be expressible in terms of a similarity variable, suggesting that the cages are nested to form a self-similar structure in the space-time.

Monodisperse hard rods in external potentials

We consider linear arrays of cells of volume V(c) populated by monodisperse rods of size σV(c),σ=1,2,..., subject to hardcore exclusion interaction. Each rod experiences a position-dependent external potential. In one application we also examine effects of contact forces between rods. We employ two distinct methods of exact analysis with complementary strengths and different limits of spatial resolution to calculate profiles of pressure and density on mesoscopic and microscopic length scales at thermal equilibrium. One method uses density functionals and the other statistically interacting vacancy particles. The applications worked out include gravity, power-law traps, and hard walls. We identify oscillations in the profiles on a microscopic length scale and show how they are systematically averaged out on a well-defined mesoscopic length scale to establish full consistency between the two approaches. The continuum limit, realized as V(c)→0,σ→∞ at nonzero and finite σV(c), connects our highest-resolution results with known exact results for monodisperse rods in a continuum. We also compare the pressure profiles obtained from density functionals with the average microscopic pressure profiles derived from the pair distribution function.

Coil-helix transition of polypeptide at water-lipid interface

We present the exact solution of a microscopic statistical mechanical model for the transformation of a long polypeptide between an unstructured coil conformation and an α-helix conformation. The polypeptide is assumed to be adsorbed to the interface between a polar and a non-polar environment such as realized by water and the lipid bilayer of a membrane. The interfacial coil-helix transformation is the first stage in the folding process of helical membrane proteins. Depending on the values of model parameters, the conformation changes as a crossover, a discontinuous transition, or a continuous transition with helicity in the role of order parameter. Our model is constructed as a system of statistically interacting quasiparticles that are activated from the helix pseudo-vacuum. The particles represent links between adjacent residues in coil conformation that form a self-avoiding random walk in two dimensions. Explicit results are presented for helicity, entropy, heat capacity, and the average numbers and sizes of sboth coil and helix segments.