Hongxia Guo

Simulation Study of Protein-Mediated Vesicle Fusion

A dissipative particle dynamics simulation method is used to probe the mechanism of protein-mediated membrane fusion. The coarse-grained models for proteins are designed based on the function of fusion proteins. Attractive forces have been introduced to produce a protein complex. The formation of protein complexes provides mechanical forces to bring membranes in proximity and trigger their merging. The whole fusion process is in good agreement with the scaffold hypothesis. Additionally, if self-defined interactions are also imposed on transmembrane segments (TMSs), their association will yield an unstable fusion pore, which can incorporate lipid and water to accomplish membrane fusion. It indicates the formation of a protein-lined pore will promote the stalk-pore transition and accelerate the fusion process.

Shear-induced parallel-to-perpendicular orientation transition in the amphiphilic lamellar phase: A nonequilibrium molecular-dynamics simulation study

The present work is devoted to a study of the shear-induced parallel-to-perpendicular orientation transition in the lamellar system by the large-scale nonequilibrium molecular-dynamics (NEMD) simulation. An effective generic model-A2B2 tetramer for amphiphilies is used. The NEMD simulation produces unambiguous evidence that undulation instability along the vorticity direction sets in well above a critical shear rate and grows in magnitude as the shear rate is further increased. At a certain high shear rate, the coherent undulation instability grows so large that defects are nucleated and the global lamellar monodomain breaks into several aligned lamellar domains. Subsequently layers in these domains rotate into the perpendicular orientation with the rotation of chains towards the y direction, merge into a global perpendicular-aligned lamellar monodomain, and organize into a perfect well-aligned perpendicular lamellar phase by the migration and annihilation of edge dislocations and disclinations. The macroscopic observable viscosity as a function of time or shear rate is correlated with the structural response such as the mesoscopic domain morphology and the microscopic chain conformation. The onset of undulation instability concurs with the start-up of shear-thinning behavior. During the orientation transformation at the high shear rate, the complex time-dependent thixotropic behavior is observed. The smaller viscosity in the perpendicular lamellar phase gives an energetic reason for the shear-induced orientation transition.

Molecular dynamics simulation of the phase behavior of lamellar amphiphilic model systems

Using efficient simplified continuum amphiphile models-AB dimers, A2B2 tetramers, and A4B4 octamers-as models for diblock copolymers, surfactants and liquid crystals, the phase behavior of the corresponding lamellar systems is investigated by large-scale parallelized dissipative particle dynamics type of molecular dynamics simulation. We not only observe the first order nature of the order–disorder transition, but also find that this first order phase transition becomes weaker as the block length of the models increases. The theoretically predicted “chain stretching” effect is reproduced as the like-monomer attractive potential well depth is increased. At the order–disorder transition point, a singular discontinuous jump in the stretching amplitude is accompanied with a distinctive orientation-induced stretching effect in the layer normal direction and an orientation induced-compressing effect in the layer plane directions. The configuration of the lamellar phase is very sensitive to the commensurability ...

Amphiphilic lamellar model systems under dilation and compression: Molecular dynamics study

Using an effective simplified continuum model for diblock copolymers, surfactants, and liquid crystals, the response of lamellar systems upon dilation and compression along the layer normal is investigated by large-scale parallelized dissipative particle dynamics like stabilized molecular dynamics simulation. We demonstrate dilation induced long-lived undulations and give clear evidence for introduction or elimination of a layer by the creation and movement of dislocation under higher dilation or compression strain.

Nonequilibrium molecular dynamics simulation study on the orientation transition in the amphiphilic lamellar phase under shear flow

By the extensive large-scale nonequilibrium molecular dynamics simulation on an effective generic model-A2B2 tetramer for amphiphiles, we investigate the shear-induced parallel to perpendicular orientation transition in the lamellar phase as a function of segregation degree and shear rate. Under low rate shear flow the evolution of parallel lamellar configurations at different segregation strengths shows a similar kinetic pathway independent of the segregation degree. While under high rate shear flow in which the lifetime of undulation instability exceeds the characteristic time of the applied shear flow, the kinetic pathway of the shear-induced parallel-to-perpendicular orientation transition in lamellar systems is the segregation degree dependent. Comparing the temporal mesoscopic domain morphology, the microscopic chain conformation, and macroscopic observable-viscosity changes with the experimentally proposed mechanisms, we find that the undulation instability, partial breakup of monodomain, grain rotation, and recombination combined with defect migration and annihilation are the kinetic pathway for the parallel-to-perpendicular orientation transition in the lamellar phase in or near the intermediate segregation limit, and that the undulation instability, domain dissolution, and reformation along the preferred direction combined with defect migration and annihilation are the kinetic pathway for the parallel-to-perpendicular orientation transition in the lamellar phase close to the order-to-disorder phase transition point. A detailed underlying microscopic picture of the alignment process illustrates that the orientation transition is driven by the alignment of molecules with shear flow. The orientation diagram that characterizes the steady-state orientations as a function of shear rate and attractive potential depth is built, in which the attractive potential depth takes the role of an inverse temperature, somewhat like the Flory-Huggins interaction parameter. The microscopic mechanism of the critical orientation transition condition is discussed.

Kinetics of the shear-induced isotropic-to-lamellar transition of an amphiphilic model system: A nonequilibrium molecular dynamics simulation study

The shear-induced isotropic-to-lamellar phase transition in the amphiphilic systems in the vicinity of the quiescent order-to-disorder transition point is investigated by the large-scale parallel nonequilibrium molecular dynamics simulations of simple amphiphilic model systems. There is a shear-induced upward shift of the ordering temperature. The initial isotropic phase orders into a lamellar phase perpendicular to the shear vorticity. The phase diagram as a function of temperature and shear rate is established. The dependency of the ordering transition on interaction strength and shear rate is rationalized by the competition between shear rate and chain relaxation. The time evolution of morphology reveals that the shear-induced ordering proceeds via nucleation and growth, a signature of a first-order phase transition. At low shear rate, a single ordered domain grows after an incubation period. With increasing shear rate ordering speeds up, but eventually develops in a lamellar system with disordered shear bands. The time dependence of the order parameter follows that of the mean-squared end-to-end distance, shear viscosity, and bulk pressure, and follows an Avrami scheme with an Avrami exponent between 2 and 4.

A computer simulation study of the segregation of amphiphiles in binary immiscible matrices: Short asymmetric copolymers in short homopolymers

The bulk and interfacial properties of ternary mixtures with asymmetric amphiphiles (A2B8) in A2 and B2 matrices and in A2 and B10 matrices are investigated by the dissipative particle dynamics type of molecular-dynamics simulations. The monomer concentrations of A2B8(phiA2B8) studied are below the critical micelle concentration (phiA2B8(cmc)) for the formation of micelles in the presence of an adsorbed amphiphilic monolayer at the interface. Macrophase separation from the mixed phase to the segregated state with A-rich and B-rich coexisting phases and the segregation of A2B8 at the interface are thermodynamically gradual but are accompanied by a pronounced stretching and orientation of the constituent chains. The segregation of A2B8 at the interface broadens the interfacial region and reduces the interfacial tension. The chain conformation of the asymmetric amphiphilic molecules and the interfacial properties are dominated by the majority block in the amphiphilic chain and dependent on the composition of the matrix in contact with the majority block. In the A2 and B2 matrices, the B8 blocks in A2B8 chains at the interface resemble a wet brush swollen by short B2 chains. Swelling is responsible for the pronounced stretching and orienting of the amphiphilic chains and the reduced interfacial amphiphile enrichment. At the same interfacial amphiphile excess, however, swollen amphiphiles are more efficient in reducing the interfacial tension than nonswollen amphiphiles.