Monte Carlo Simulation of the Assembly of bis‐Biotinylated DNA and Streptavidin

We present Monte Carlo simulations of the self‐assembly of bivalent bis‐biotinylated DNA molecules with the tetravalent biotin‐binding protein streptavidin (STV). By fitting the STV binding probabilities for the four possible valencies, the modelling correctly reproduces the dependencies of various network parameters experimentally observed in an earlier study. The combined results from the experimental and theoretical studies suggest that the binding probability for divalent STV formation is about 50 times larger than for the formation of trivalent and about 200 times larger than for tetravalent STV. In accordance with the experimental results, the modelling also indicates that the mixture of an equimolar ratio of DNA and STV leads to a maximum in size of the oligomeric DNA–STV clusters formed. Furthermore, we found a percolation transition in which the DNA cluster size increases rapidly with increasing DNA concentration resulting in the formation of a single supercluster at elevated concentrations. This behaviour coincides with the occurrence of an immobile band previously observed in electrophoretic experiments, indicating the formation of extremely large DNA–STV aggregate networks.     

Anion mediated polytype selectivity among the basic salts of Co(II)

Basic salts of Co(II) crystallize in the rhombohedral structure. Two different polytypes, 3R₁ and 3R₂, with distinct stacking sequences of the metal hydroxide slabs, are possible within the rhombohedral structure. These polytypes are generated by simple translation of successive layers by (2/3, 1/3, z) or (1/3, 2/3, z). The symmetry of the anion and the mode of coordination influences polytype selection. Cobalt hydroxynitrate crystallizes in the structure of the 3R₂ polytype while the hydroxytartarate, hydroxychloride and α-cobalt hydroxide crystallize in the structure of the 3R₁ polytype. Cobalt hydroxysulfate is turbostratically disordered. The turbostratic disorder is a direct consequence of the mismatch between the crystallographically defined interlayer sites generated within the crystal and the tetrahedral symmetry of the SO₄²⁻ ions.     

Estimation of the absolute internal-rotation entropy of molecules with two torsional degrees of freedom from stochastic simulations

A method of statistical estimation is applied to the problem of evaluating the absolute entropy of internal rotation in a molecule with two torsional degrees of freedom. The configurational part of the entropy is obtained as that of the joint probability density of an arbitrary form represented by a two-dimensional Fourier series, the coefficients of which are statistically estimated using a sample of the torsional angles of the molecule obtained by a stochastic simulation. The internal rotors in the molecule are assumed to be attached to a common frame, and their reduced moments of inertia are initially calculated as functions of the two torsional angles, but averaged over all the remaining internal degrees of freedom using the stochastic-simulation sample of the atomic configurations of the molecule. The torsional-angle dependence of the reduced moments of inertia can be also averaged out, and the absolute internal-rotation entropy of the molecule is obtained in a good approximation as the sum of the configurational entropy and a kinetic contribution fully determined by the averaged reduced moments of inertia. The method is illustrated using Monte Carlo simulations of isomers of stilbene and halogenated derivatives of propane. The two torsional angles in cis-stilbene are found to be much more strongly correlated than those in trans-stilbene, while the degree of the angular correlation in propane increases strongly on substitution of hydrogen atoms with chlorine.     

Simple algorithm for isothermal-isobaric molecular dynamics

We review principles of non-Hamiltonian statistical mechanics and present a new set of equations and integration algorithm for isothermal-isobaric dynamics. The chief advantage of the present scheme is that it is somewhat simpler than previous methods. We perform numerical simulations to test the accuracy of the algorithm and compare its stability to that of a "gold standard," a symplectic integrator for Hamiltonian dynamics of the same system. The stability of the isothermal-isobaric algorithm is comparable to the stability of the symplectic integrator.     

Fine structures of zeolite-Linde-L (LTL): surface structures, growth unit and defects

High-resolution electron microscopy (HREM) has been used to image the surface structure of nano- and micrometer-sized synthetic crystals of zeolite-Linde-L (LTL). Columnar holes and rotational, nano-sized, wheel-like defects were observed within the crystals, where the hole has a minimum size equal to that of the rotational defect. Predictions of surface structure from atomistic computer simulation concur with the observations from HREM and provide insight into the crystal growth mechanism of perfect and defective LTL. Analysis of the energetics of the formation of rotational defect structures reveals that the driving force for defect creation is thermodynamic and furthermore, the rotational defects could be created in high concentrations. Formation of a columnar hole is found to be slightly energetically unfavourable and therefore we speculate that the incidence of both rotational and nano-sized vacancy defects is strongly dependent on kinetic factors and reaction conditions. The morphology of nano- and microcrystalline LTL is contradistinct and we use insights from simulation to propose an explanation of the disparity in crystal shape.     

Subtleties in data analysis related to the size of critical region

We comment on the analysis of the critical behavior of a layered driven diffusive system recently done by Achahbar and Marro. We discuss why we believe their method of taking the thermodynamic limit and determining the order-parameter exponent leads to unreliable estimates.     

Formation and morphology of latex monolayers. Computer simulation studies

The results of computer simulations of monolayers created from monodisperse latex particles are presented and discussed. Layers are characterized by the normalized coverage,NC (the actual coverage of the surface related to its maximum possible coverage with particles), and by the average number of neighbors,ANN, calculated as the number of particles being in contact with a given one and averaged over all the particles on the surface. Variable parameters used in simulations include: the rate of particles deposition, the probability of lateral movements, the probability of desorption of particles adsorbed on the surface, the probability of covalent immobilization of adsorbed particles, and the on-sphere slip parameter, OSS (characterizing the scattering of a falling particle on the particles being already attached to the surface). Morphology of monolayers is qualitatively characterized by relations betweenANN andNC. It is shown that for all monolayers formed without adhesion (and without repulsion) between the particles adsorbed on the surface the dependence ofANN vs.NC is described by a characteristic master curve (regardless of the values of probabilities of desorption and lateral movements of particles). For the monolayers created including adhesive forces between the adsorbed particles the plots ofANN vs.NC lie above the master curve, while similar plots obtained for the layers made of particles showing various types of repulsive interactions are always placed below it. Thus, the dependencies ofANN vs.NC, derived from computer simulations, can be used for the determination of the character of the interparticle interactions in the real systems.     

Renormalization group for Markov chains and application to metastability

In this paper we introduce a new renormalization group method for the study of the long-time behavior of Markov chains with finite state space and with transition probabilities exponentially small in an external parameter. A general approach of metastability problems emerges from this analysis and is discussed in detail in the case of a two-dimensional Ising system at low temperature,     

Atomistic Simulation Study of Calcium,Phosphate and Fluoride Ion Association to the Teleopeptide‐Tails of Collagen – Initial Steps to Biomineral Formation

The attachment of single ions to putative adsorption sites in the tails of collagen fibers is investigated by means of molecular dynamics simulations and discussed with respect to the very early steps of apatite/collagen biomineral formation. Our studies clearly demonstrate an increased flexibility of the tails of the triple‐helical collagen protein. Apart from the termini of the backbone, several side chains were also observed to be freely accessible to ion attachment from aqueous solution. The teleopeptide was systematically scanned for suitable adsorption sites for calcium, phosphate and fluoride ions. Association of these ions was then explored from potential of mean force calculations. The resulting energy profiles reveal a variety of favorable protein‐ion bonds and hint at the suitability of the collagen tails to promote apatite aggregation.     

Computer experiments on crystalline nylons: structural analysis of nylons with large aliphatic segments

A theoretical study based on force-field calculations has been performed to investigate the structural preferences of crystalline even nylons n with large and very large aliphatic segments. Atomistic energy calculations and Monte Carlo simulations were carried out considering the conventional and forms of nylons 10, 12, 18, 24, and 32. Results indicated that the form is the most favored for nylons 10, 12, 18, and 24. However, the structure was unstable for nylon 32, a polymer in which the density of hydrogen bonds is almost negligible. In this case, the arrangement is energetically more favored than the one.