Fractal nematic colloids.

Fractals are remarkable examples of self-similarity where a structure or dynamic pattern is repeated over multiple spatial or time scales. However, little is known about how fractal stimuli such as fractal surfaces interact with their local environment if it exhibits order. Here we show geometry-induced formation of fractal defect states in Koch nematic colloids, exhibiting fractal self-similarity better than 90% over three orders of magnitude in the length scales, from micrometers to nanometres. We produce polymer Koch-shaped hollow colloidal prisms of three successive fractal iterations by direct laser writing, and characterize their coupling with the nematic by polarization microscopy and numerical modelling. Explicit generation of topological defect pairs is found, with the number of defects following exponential-law dependence and reaching few 100 already at fractal iteration four. This work demonstrates a route for generation of fractal topological defect states in responsive soft matter.

Understanding the nanoparticle-protein corona complexes using computational and experimental methods.

Nanoparticles (NP) have capability to adsorb proteins from biological fluids and form protein layer, which is called protein corona. As the cell sees corona coated NPs, the protein corona can dictate biological response to NPs. The composition of protein corona is varied by physicochemical properties of NPs including size, shape, surface chemistry. Processing of protein adsorption is dynamic phenomena; to that end, a protein may desorb or leave a surface vacancy that is rapidly filled by another protein and cause changes in the corona composition mainly by the Vroman effect. In this review, we discuss the interaction between NP and proteins and the available techniques for identification of NP-bound proteins. Also we review current developed computational methods for understanding the NP-protein complex interactions.

Stiffer double-stranded DNA in two-dimensional confinement due to bending anisotropy.

Using analytical approach and Monte Carlo (MC) simulations, we study the elastic behavior of the intrinsically twisted elastic ribbons with bending anisotropy, such as double-stranded DNA (dsDNA), in two-dimensional (2D) confinement. We show that, due to the bending anisotropy, the persistence length of dsDNA in 2D conformations is always greater than three-dimensional (3D) conformations. This result is in consistence with the measured values for DNA persistence length in 2D and 3D in equal biological conditions. We also show that in two dimensions, an anisotropic, intrinsically twisted polymer exhibits an implicit twist-bend coupling, which leads to the transient curvature increasing with a half helical turn periodicity along the bent polymer.

Extreme bendability of DNA double helix due to bending asymmetry.

Experimental data of the DNA cyclization (J-factor) at short length scales exceed the theoretical expectation based on the wormlike chain (WLC) model by several orders of magnitude. Here, we propose that asymmetric bending rigidity of the double helix in the groove direction can be responsible for extreme bendability of DNA at short length scales and it also facilitates DNA loop formation at these lengths. To account for the bending asymmetry, we consider the asymmetric elastic rod (AER) model which has been introduced and parametrized in an earlier study [B. Eslami-Mossallam and M. R. Ejtehadi, Phys. Rev. E 80, 011919 (2009)]. Exploiting a coarse grained representation of the DNA molecule at base pair (bp) level and using the Monte Carlo simulation method in combination with the umbrella sampling technique, we calculate the loop formation probability of DNA in the AER model. We show that the DNA molecule has a larger J-factor compared to the WLC model which is in excellent agreement with recent experimental data.

Particle selection through topographic templates in nematic colloids.

Liquid crystal colloids have been proposed as suitable candidates for responsive photonic crystals. Large scale growth of such colloidal systems is, however, a challenge and recently template-assisted assembly has been proposed to guide the growth of colloidal crystals, with controlled symmetries, in nematic liquid crystals. Known for their long-range anisotropic interactions, these colloidal systems are stabilized typically at the center of the cells due to strong particle-wall repulsion from the confining substrates. This behaviour is dramatically changed in the presence of topographic patterning. Here we propose the use of topographic modulation of surfaces to select and localize particles in nematic colloids. By considering convex and concave deformations of one of the confining surfaces we show that the colloid-flat surface repulsion may be enhanced or switched into an attraction. In particular, we find that when the colloidal particles have the same anchoring conditions as the patterned surfaces, they are strongly attracted to concave dimples, while if they exhibit different anchoring conditions they are pinned at the top of convex protrusions. Although dominated by elastic interactions the first mechanism is reminiscent of the depletion induced attraction or of the key-lock mechanism, while the second is specific to liquid crystal colloids. These long-ranged, highly tunable, surface-colloid interactions contribute to the development of template-assisted assembly of large colloidal crystals, with well defined symmetries, as required for applications.

Confinement dynamics of a semiflexible chain inside nano-spheres.

We study the conformations of a semiflexible chain, confined in nano-scaled spherical cavities, under two distinct processes of confinement. Radial contraction and packaging are employed as two confining procedures. The former method is performed by gradually decreasing the diameter of a spherical shell which envelopes a confined chain. The latter procedure is carried out by injecting the chain inside a spherical shell through a hole on the shell surface. The chain is modeled with a rigid body molecular dynamics simulation and its parameters are adjusted to DNA base-pair elasticity. Directional order parameter is employed to analyze and compare the confined chain and the conformations of the chain for two different sizes of the spheres are studied in both procedures. It is shown that for the confined chains in the sphere sizes of our study, they appear in spiral or tennis-ball structures, and the tennis-ball structure is more likely to be observed in more compact confinements. Our results also show that the dynamical procedure of confinement and the rate of the confinement are influential parameters of the structure of the chain inside spherical cavities.

Gating and conduction of nano-channel forming proteins: a computational approach.

Monitoring conformational changes in ion channels is essential to understand their gating mechanism. Here, we explore the structural dynamics of four outer membrane proteins with different structures and functions in the slowest nonzero modes of vibration. Normal mode analysis was performed on the modified elastic network model of channel in the membrane. According to our results, when membrane proteins were analyzed in the dominant mode, the composed pores, TolC and α-hemolysin showed large motions at the intramembrane ß-barrel region while, in other porins, OmpA and OmpF, largest motions observed in the region of external flexible loops. A criterion based on equipartition theorem was used to measure the possible amplitude of vibration in channel forming proteins. The current approach complements theoretical and experimental techniques including HOLE, Molecular Dynamics (MD), and voltage clamp used to address the channel's structure and dynamics and provides the means to conduct a theoretical simultaneous study of the structure and function of the channel. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:3.

Definition of the persistence length in the coarse-grained models of DNA elasticity.

By considering the detailed structure of DNA in the base pair level, two possible definitions of the persistence length are compared. One definition is related to the orientation of the terminal base pairs, and the other is based on the vectors which connect two adjacent base pairs at each end of the molecule. It is shown that although these definitions approach each other for long DNA molecules, they are dramatically different on short length scales. We show analytically that the difference mostly comes from the shear flexibility of the molecule and can be used to measure the shear modulus of DNA.

Directed motion of C60 on a graphene sheet subjected to a temperature gradient.

Nonequilibrium molecular dynamics simulations are used to study the motion of a C(60) molecule on a graphene sheet subjected to a temperature gradient. The C(60) molecule is actuated and moves along the system while it just randomly dances along the perpendicular direction. Increasing the temperature gradient increases the directed velocity of C(60). It is found that the free energy decreases as the C(60) molecule moves toward the cold end. The driving mechanism based on the temperature gradient suggests the construction of nanoscale graphene-based motors.

Contribution of nonlocal interactions to DNA elasticity.

A nonlocal harmonic elastic rod model is proposed to describe the elastic behavior of short DNA molecules. We show that the nonlocal interactions contribute to effective bending energy of the molecule and affect its apparent persistence length. It is also shown that the anomalous behavior which has been observed in all-atom molecular dynamic simulations [A. K. Mazur, Biophys. J. 134, 4507 (2006)] can be a consequence of both nonlocal interactions between DNA base pairs and the intrinsic curvature of DNA.

Membrane interactions control residue fluctuations of outer membrane porins.

Bacterial outer membrane porins (Omp) that have robust ß -barrel structures, show potential applications for nanomedicine devices in synthetic membranes and single molecule detection biosensors. Here, we explore the conformational dynamics of a set of 22 outer membrane porins, classified into five major groups: general porins, specific porins, transport Omps, poreless Omps and composed pores. Normal mode analysis, based on mechanical vibration theory and elastic network model, is performed to study the fluctuations of residues of aforementioned porins around their equilibrium positions. We find that a simple modification in this model considering weak interaction between protein and membrane, dramatically enhance the stability of results and improve the correlation coefficient between computational output and experimental results.

Diffusive motion of C60 on a graphene sheet.

The motion of a C60 molecule over a graphene sheet at finite temperature is investigated both theoretically and computationally. We show that a graphene sheet generates a van der Waals laterally periodic potential, which directly influences the motion of external objects in its proximity. The translational motion of a C60 molecule near a graphene sheet is found to be diffusive in the lateral directions, while in the perpendicular direction, the motion may be described as diffusion in an effective harmonic potential which is determined from the distribution function of the position of the C60 molecule. We also examine the rotational diffusion of C60 and show that its motion over the graphene sheet is not a rolling motion.

Asymmetric elastic rod model for DNA.

In this paper we consider the anharmonic corrections to the anisotropic elastic rod model for DNA. Our model accounts for the difference between the bending energies of positive and negative rolls, which comes from the asymmetric structure of the DNA molecule. We will show that the model can explain the high flexibility of DNA at small length scales, as well as kink formation at high deformation limit.

Effective potential of longitudinal interactions between microtubule protofilaments.

An effective potential for longitudinal interactions between adjacent protofilaments in a microtubule is introduced. Our proposed interaction potential is a periodic and continuous function of the offset between two protofilaments, which also incorporates the bending energy of protofilaments. This potential produces the results of atomistic simulations. Further, using the potential, a Monte Carlo simulation gives results for the skew angles of observed structures that are in good agreement with experiments.

Stretching an anisotropic DNA.

We present a perturbation theory to find the response of an anisotropic DNA to the external tension. It is shown that the anisotropy has a nonzero but small contribution to the force-extension curve of the DNA. Thus an anisotropic DNA behaves like an isotropic one with an effective bending constant equal to the harmonic average of its soft and hard bending constants.

Coarse-grained interaction potentials for anisotropic molecules.

We have proposed an efficient parametrization method for a recent variant of the Gay Berne potential for dissimilar and biaxial particles [Phys. Rev. E 67, 041710 (2003)] and demonstrated it for a set of small organic molecules. Compared with the previously proposed coarse-grained models, the new potential exhibits a superior performance in close contact and large distant interactions. The repercussions of thermal vibrations and elasticity have been studied through a statistical method. The study justifies that the potential of mean force is representable with the same functional form, extending the application of this coarse-grained description to a broader range of molecules. Moreover, the advantage of employing coarse-grained models over truncated atomistic summations with large distance cutoffs has been briefly studied.

Conservation of statistical results under the reduction of pair-contact interactions to solvation interactions.

We show that the hydrophobicity of sequences is the leading term in Miyazawa-Jernigan interactions. Being the source of additive (solvation) terms in pair-contact interactions, they were used to reduce the energy parameters while resulting in a clear vector manipulation of energy. The reduced (additive) potential performs considerably successful in predicting the statistical properties of arbitrary structures. The evaluated designabilities of the structures by both models are highly correlated. Suggesting geometrically nondegenerate vectors (structures) as proteinlike structures, the additive model is a powerful tool for protein design. Moreover, a crossing point in the log-linear diagram of designability ranking shows that about 1/e of the structures have designabilities above the average, independent on the used model.

Protein folding rates correlate with heterogeneity of folding mechanism.

By observing trends in the folding kinetics of experimental 2-state proteins at their transition midpoints, and by observing trends in the barrier heights of numerous simulations of coarse-grained, C(alpha) model Go proteins, we show that folding rates correlate with the degree of heterogeneity in the formation of native contacts. Statistically significant correlations are observed between folding rates and measures of heterogeneity inherent in the native topology, as well as between rates and the variance in the distribution of either experimentally measured or simulated phi values.

Three-body interactions improve the prediction of rate and mechanism in protein folding models.

Here we study the effects of many-body interactions on rate and mechanism in protein folding by using the results of molecular dynamics simulations on numerous coarse-grained Calpha-model single-domain proteins. After adding three-body interactions explicitly as a perturbation to a Go-like Hamiltonian with native pairwise interactions only, we have found (i) a significantly increased correlation with experimental phi values and folding rates, (ii) a stronger correlation of folding rate with contact order, matching the experimental range in rates when the fraction of three-body energy in the native state is approximately 20%, and (iii) a considerably larger amount of three-body energy present in chymotripsin inhibitor than in the other proteins studied.

Interaction potentials for soft and hard ellipsoids.

Using results from colloid science we derive interaction potentials for computer simulations of mixtures of soft or hard ellipsoids of arbitrary shape and size. Our results are in many respects reminicent of potentials of the Gay-Berne type but have a well-defined microscopic interpretation and no adjustable parameters. Since our potentials require the calculation of similar variables, the modification of existing simulation codes for Gay-Berne potentials is straightforward. The computational performance should remain unaffected.