Monte Carlo simulations of patch models with applications to soft matter.

Patch models have been employed to describe anisotropic interactions in disparate soft matter systems. In this work, we present a unified study of a patch model to explore particle self-assembly in both monodisperse and polydisperse systems, with applications to both proteins and colloids. In the first case, we obtained a temperature-density phase diagram for a model of the protein polyglutamine from Monte Carlo simulations. These simulations evinced clusters in the gas phase and, via a comparison with the corresponding coarse-grained PLUM model for this system, we verified that dense clustering in the gas phase falls into the supersaturation region of the saturation curve. In the second case, we have investigated the effect of size polydispersity on the phase behavior of binary colloidal mixtures. It was found that the width of gas-liquid phase coexistence increases with increasing polydispersity and that, in addition to the aforementioned particle clustering, small particles decorate patches on large particles, thereby creating large-particle bridges. Our aim is to compare and contrast self-assembly in these two prototypical systems.

Early stage aggregation of a coarse-grained model of polyglutamine.

In this paper, we study the early stages of aggregation of a model of polyglutamine (polyQ) for different repeat lengths (number of glutamine amino acid groups in the chain). In particular, we use the Large-scale Atomic/Molecular Massively Parallel Simulator to study a generic coarse-grained model proposed by Bereau and Deserno. We focus on the primary nucleation mechanism involved and find that our results for the initial self-assembly process are consistent with the two-dimensional classical nucleation theory of Kashchiev and Auer. More specifically, we find that with decreasing supersaturation, the oligomer fibril (protofibril) transforms from a one-dimensional ß sheet to two-, three-, and higher layer ß sheets as the critical nucleus size increases. We also show that the results are consistent with several predictions of their theory, including the dependence of the critical nucleus size on the supersaturation. Our results for the time dependence of the mass aggregation are in reasonable agreement with an approximate analytical solution of the filament theory by Knowles and collaborators that corresponds to an additional secondary nucleation arising from filament fragmentation. Finally, we study the dependence of the critical nucleus size on the repeat length of polyQ. We find that for the larger length polyglutamine chain that we study, the critical nucleus is a monomer, in agreement with experiment and in contrast to the case for the smaller chain, for which the smallest critical nucleus size is four.

Phase diagram of a model of the protein amelogenin.

There has been considerable recent interest in the self-assembly and phase behavior of models of colloidal and protein particles with anisotropic interactions. One example of particular interest is amelogenin, an important protein involved in the formation of dental enamel. Amelogenin is primarily hydrophobic with a 25-residue charged C-terminus tail. This protein undergoes a hierarchical assembly process that is crucial to mineral deposition, and experimental work has demonstrated that the deletion of the C-terminus tail prevents this self-assembly. A simplified model of amelogenin has been proposed in which the protein is treated as a hydrophobic sphere, interacting via the Asakura-Oosawa (AO) potential, with a tethered point charge on its surface. In this paper, we examine the effect of the Coulomb interaction between the point charges in altering the phase diagram of the AO model. For the parameter case specific to amelogenin, we find that the previous in vitro experimental and model conditions correspond to the system being near the low-density edge of the metastable region of the phase diagram. Our study illustrates more generally the importance of understanding the phase diagram for proteins, in that the kinetic pathway for self-assembly and the resulting aggregate morphology depends on the location of the initial state in the phase diagram.

Phase behavior of patchy spheroidal fluids.

We employ Gibbs-ensemble Monte Carlo computer simulation to assess the impact of shape anisotropy and particle interaction anisotropy on the phase behavior of a colloidal (or, by extension, protein) fluid comprising patchy ellipsoidal particles, with an emphasis on critical behavior. More specifically, we obtain the fluid-fluid equilibrium phase diagram of hard prolate ellipsoids having Kern-Frenkel surface patches under a variety of conditions and study the critical behavior of these fluids as a function of particle shape parameters. It is found that the dependence of the critical temperature on aspect ratio for particles having the same volume can be described approximately in terms of patch solid angles. In addition, ordering in the fluid that is associated with particle elongation is also found to be an important factor in dictating phase behavior.

Kinetics of aggregation of an anisotropic model of self-assembling molecules.

We study the kinetics of aggregation of a two site model of interacting spherical molecules. A given site on one molecule can interact with one or more sites on other neighboring molecules. The sites represent the result of a simple coarse graining of putative amino acid residues or two specifically designed sites on a colloidal particle. We study the kinetics and equilibrium morphology for a fixed angle between the two sites, for several angles between 30° and 150°. In the model, the sites interact via an attractive Asakura-Oosawa potential and the molecules have the usual hard sphere repulsion interaction. We find a transition from a micelle-like morphology at small angles to a rod-like morphology at intermediate angles and to a gel-like structure at values of the angle greater than about ninety degrees. However, at 150 degrees, after a long induction time during which there is no aggregation, we observe a nucleation and growth process that leads to a final spherical-like aggregate. Our results show that this angle is a control parameter for the kinetics and equilibrium properties of the system.

The impact of anisotropy and interaction range on the self-assembly of Janus ellipsoids.

We assess the roles of anisotropy and interaction range on the self-assembly of Janus colloidal particles. In particular, Monte Carlo simulation is employed to investigate the propensity for the formation of aggregates in a spheroidal model of a colloid having a relatively short-ranged interaction that is consistent with experimentally realizable systems. By monitoring the equilibrium distribution of aggregates as a function of temperature and density, we identify a "micelle" transition temperature and discuss its dependence on particle shape. We find that, unlike systems with longer ranged interactions, this system does not form micelles below a transition temperature at low density. Rather, larger clusters comprising 20-40 particles characterize the transition. We then examine the dependence of the second virial coefficient on particle shape and well width to determine how these important system parameters affect aggregation. Finally, we discuss possible strategies suggested by this work to promote self-assembly for the encapsulation of particles.

Multicenter Australian trial of islet transplantation: improving accessibility and outcomes.

Whilst initial rates of insulin independence following islet transplantation are encouraging, long-term function using the Edmonton Protocol remains a concern. The aim of this single-arm, multicenter study was to evaluate an immunosuppressive protocol of initial antithymocyte globulin (ATG), tacrolimus and mycophenolate mofetil (MMF) followed by switching to sirolimus and MMF. Islets were cultured for 24 h prior to transplantation. The primary end-point was an HbA1c of <7% and cessation of severe hypoglycemia. Seventeen recipients were followed for ≥ 12 months. Nine islet preparations were transported interstate for transplantation. Similar outcomes were achieved at all three centers. Fourteen of the 17 (82%) recipients achieved the primary end-point. Nine (53%) recipients achieved insulin independence for a median of 26 months (range 7-39 months) and 6 (35%) remain insulin independent. All recipients were C-peptide positive for at least 3 months. All subjects with unstimulated C-peptide >0.2 nmol/L had cessation of severe hypoglycemia. Nine of the 17 recipients tolerated switching from tacrolimus to sirolimus with similar graft outcomes. There was a small but significant reduction in renal function in the first 12 months. The combination of islet culture, ATG, tacrolimus and MMF is a viable alternative for islet transplantation.

Kinetics of nanochain formation in a simplified model of amelogenin biomacromolecules.

We show that the kinetics of nanochain formation of amelogenin molecules is well described by a combination of translational and rotational diffusion of a simplified anisotropic bipolar model consisting of hydrophobic spherical colloid particles and a point charge located on each particle surface. The colloid particles interact via a standard depletion attraction whereas the point charges interact through a screened Coulomb repulsion. We study the kinetics via a Brownian dynamics simulation of both translational and rotational motions and show that the anisotropy brought in by the charge dramatically affects the kinetic pathway of cluster formation and our simple model captures the main features of the experimental observations.

Synergistic effects of genetic beta cell dysfunction and maternal glucose intolerance on offspring metabolic phenotype in mice.

AIMS/HYPOTHESIS: Diabetes in pregnancy is linked to development of obesity in the offspring, but the mechanisms are not fully understood. Gestational diabetes mellitus (GDM) occurs when beta cells are unable to compensate for the normal insulin resistance of late pregnancy. In this study, we used a murine model of beta cell dysfunction to examine the effects of maternal GDM on phenotype in male offspring with and without an inherited predisposition for beta cell dysfunction. METHODS: Beta cell-specific aryl-hydrocarbon receptor nuclear translocator-null (ßArnt) mice develop GDM from beta cell dysfunction. ßArnt and control female mice were used to induce GDM and non-diabetic pregnancies, respectively. RESULTS: Offspring from GDM pregnancies became spontaneously obese on a normal-chow diet. They were heavier than offspring from non-diabetic pregnancies, with increased body fat. Respiratory exchange ratio (RER) was higher, indicating decreased capacity to switch to lipid oxidation. Metabolic rate in GDM offspring was decreased prior to onset of obesity. The phenotype was more pronounced in ßArnt GDM offspring than in GDM offspring of control genotype, demonstrating an interaction between genotype and pregnancy exposure. ßArnt GDM offspring had increased hypothalamic neuropeptide Y (Npy) and decreased pro-opiomelanocortin (Pomc) expression. Weight, body fat, insulin sensitivity and RER in all mice, and hypothalamic Npy in ßArnt mice were significantly correlated with AUC of maternal late pregnancy glucose tolerance tests (p < 0.01), but not with litter size, maternal weight, triacylglycerol or pre-pregnancy glycaemia. CONCLUSIONS/INTERPRETATION: In ßArnt mice, exposure to GDM and inheritance of genetic beta cell dysfunction had additive effects on male offspring obesity; severity of the offspring phenotype correlated with maternal glycaemia.

Pathways of cluster growth and kinetic slowing down in a model of short-range attractive colloids.

We present the results of an extensive 3D Brownian dynamics simulation of the self-assembly of colloidal particles for a short-range attractive model that is quenched below its metastable critical point. In particular, results are obtained in the small-volume-fraction, low-temperature region in which we find so-called sticky beads that diffuse around the system, without reaching a final large cluster on the timescale of our simulation. For larger volume fractions in this low-temperature regime, a gel forms as the result of kinetically slowed down spinodal decomposition, as shown earlier for other short-range attractive models (Foffi, G.; De Michele, C.; Sciortino, F.; Tartaglia, P. Phys. Rev. Lett. 2005, 94, 078301. Zaccarelli, E. J. Phys.: Condens. Matter 2007, 19, 323101). We also show that for quenches below the critical point but above the intersection of the binodal with the glass line, two-step crystallization takes place. For sufficiently small volume fractions, the first step is the nucleation of dense fluid drops, followed by the second step of crystallization within these drops, as first proposed for a model of protein crystallization for quenches just above the metastable critical point (ten Wolde, P. R.; Frenkel, D. Science 1997, 277, 1975). For larger values of the volume fraction, the initial step is spinodal decomposition that leads to the formation of an interconnected network of low- and high-density fluids. The second step is crystallization that takes place within the dense fluid phase.

Statistical mechanics of helical wormlike chain model.

We investigate the statistical mechanics of polymers with bending and torsional elasticity described by the helical wormlike model. Noticing that the energy function is factorizable, we provide a numerical method to solve the model using a transfer matrix formulation. The tangent-tangent and binormal-binormal correlation functions have been calculated and displayed rich profiles which are sensitive to the combination of the temperature and the equilibrium torsion. Their behaviors indicate that there is no finite temperature Lifshitz point between the disordered and helical phases. The asymptotic behavior at low temperature has been investigated theoretically and the predictions fit the numerical results very well. Our analysis could be used to understand the statics of dsDNA and other chiral polymers.

Brownian dynamics simulation of insulin microsphere formation from break-up of a fractal network.

Motivated by a recent experiment on insulin microsphere formation where polyethylene glycol (PEG) is used as the precipitating agent, we have developed a simple theoretical model that can predict the formation of a fractal network of insulin monomers and the subsequent break-up of the fractal network into microsphere aggregates. In our approach the effect of PEG on insulin is modeled via a standard depletion attraction mechanism via the Asakura-Oosawa model. We show that even in the context of this simple model, it is possible to mimic important aspects of the insulin experiment in a brownian dynamics simulation. We simulate the effect of changing temperature in our model by changing the well depth of the Asakura-Oosawa potential. A fractal network is observed in a "deep quench" of the system, followed by a "heating" that results in a break-up of the network and subsequent formation of microspheres.

A simple model of directional interactions for proteins.

We study a simple two patch model of globular protein solutions. The model consists of a hard sphere interaction, together with a weak isotropic attraction, decorated with stronger attractive patch-patch interactions. The isotropic and anisotropic attractions are modeled with square well potentials with an interaction range of 1.2sigma, where sigma is the hard sphere diameter. We have calculated its phase diagram and found a metastable fluid-fluid phase separation curve. We have also obtained the different stable crystal structures at various temperatures. The orientationally ordered, body-centered crystal (bcc-o) lattice is the crystal structure that has the lowest energy; it is stable at lower temperatures and moderate pressures. The orientationally ordered face-centered crystal (fcc) lattice has a higher energy and becomes stable with respect to bcc-o at high pressures and lower temperatures. Finally, at high temperatures when the entropy effect becomes important, the orientational order is lost and an orientationally disordered fcc lattice becomes stable.

Ion specific effects on phase transitions in protein solutions.

A recent Monte Carlo simulation determined the potential of mean force between two lysozyme molecules in various aqueous solutions [M. Lund, Phys. Rev. Lett. 100, 258105 (2008)]. The study involved a combination of explicit solvent and continuum model simulations and showed that there are significant ion-specific protein-protein interactions due to hydrophobic patches on the protein surfaces. In this paper we use the results of their study to determine the phase diagram for lysozyme for aqueous solutions of NaCl and NaI. Two of the three phase diagrams have a stable fluid-fluid critical point, while the third has a slightly metastable critical point. This results from a secondary extremum in the potential associated with a repulsive interaction. This repulsive interaction reduces the effective range of the attractive interaction and produces a metastable critical point. We compare the results of one of these phase diagrams with that for a model that includes ion-dispersion forces, but does not contain solvent structural effects.

Kinetics and microstructure associated with nonisothermal nucleation and growth processes.

The impact of nonisothermal conditions on the kinetics and microstructure associated with a first-order phase transformation is investigated using simulation and a correlation function formalism. More specifically, equal-time correlation functions and microstructural descriptors of grain area are calculated using an N-fold Monte Carlo simulation and compared directly with their theoretical counterparts for a range of heating rates. Finally, connections between these results and those obtained experimentally via calorimetry and microscopy are discussed.

Phase diagram of a model of nanoparticles in electrolyte solutions.

We obtain accurate fluid-fluid coexistence curves for a recent simple model of interacting nanoparticles that includes the effects of ion-dispersion forces. It has been proposed that these ion-dispersion forces provide at least a partial explanation for the Hofmeister effect [M. Bostrom et al. Phys. Rev. Lett. 87, 168103 (2001)]. We study a model of aluminum oxide nanoparticle [Deniz et al., Colloids Surf. A 319, 98 (2008)] for three different electrolyte solutions with added salt type being sodium chloride, sodium iodide, and a nonpolarizable salt. We observe that the fluid-fluid coexistence curves depend substantially on the identity of added salt; this provides an efficient way of tuning the phase behavior of nanoparticles. The methods we employ include finite-size scaling (FSS), multicanonical histogram reweighting, and Gibbs ensemble methods. We show that, as expected, all three cases belong to the Ising universality class. The scaling fields and critical point parameters are obtained in the thermodynamic limit of infinite system size by extrapolation of our FSS results.

Hofmeister effect and the phase diagram of lysozyme.

The phase diagrams for lysozyme are calculated for two different precipitant salts, NaCl and NaSCN, using a potential of mean force that takes into account contributions from ion-dispersion forces [M. Boström, J. Phys. Chem. B 110, 24757 (2006)]. Our results are consistent with a recent perturbation theory calculation (referenced above) in that the phase diagram for lysozyme with NaCl is quite different than for lysozyme with NaSCN for the same molar concentration (0.2M) . However, in contrast to the perturbation theory calculation, we find that the lysozyme phase diagram with NaCl has a metastable fluid-fluid coexistence curve and that the metastability gap in the case of NaSCN is much larger than predicted by perturbation theory.

Ethnicity influences cardiovascular outcomes and complications in patients with type 2 diabetes.

AIM: To determine whether cardiovascular outcomes in type 2 diabetes (T2D) differ according to ethnicity, and whether ethnicity influences the effect of gender on these outcomes in Caucasians, East-Southeast-Asians, Middle-Easterners, South-Asians and Pacific-Islanders. METHODS: We compared demographics, HbA1c, lipid profile, renal function markers, and prevalence of macrovascular and microvascular complications between ethnic groups. Cross-sectional data was prospectively collected from 204 consecutive patients at Westmead Hospital's T2D clinic from April-October 2015. Univariate analysis was performed using chi-squared test for categorical data, and Mann-Whitney-U or Kruskal-Wallis test for continuous data. RESULTS: Compared to Caucasians, South-Asians were diagnosed younger, were currently younger, had lower body-mass-index (BMI) and better renal function but higher rates of non-ST-elevation myocardial infarction (STEMI, 21.7% versus 3.5%, p<0.05). East-Southeast-Asians had lower BMI but more nephropathy than Caucasians (59% versus 39%, p<0.05). East-Southeast-Asian males had fewer CVD than Caucasians, but this protection was absent in East-Southeast-Asian females. Middle-Easterners had more non-STEMI than Caucasians (5.3% vs 3.5%, p<0.05). Middle-Eastern females were not at lower CVD risk than males. Caucasians had most PVD (20% versus 6%, p<0.05). CONCLUSIONS: Ethnicity influences rates of diabetes-related complications. Female CVD protection is altered in some groups. Ethnicity should be considered in assessing CVD and complications risk.

1,25-Dihydroxycholecalciferol (calcitriol) modifies uptake and release of 25-hydroxycholecalciferol in skeletal muscle cells in culture.

The major circulating metabolite of vitamin D3, 25-hydroxycholecalciferol [25(OH)D], has a remarkably long half-life in blood for a (seco)steroid. Data from our studies and others are consistent with the hypothesis that there is a role for skeletal muscle in the maintenance of vitamin D status. Muscle cells internalise vitamin D-binding protein (DBP) from the circulation by means of a megalin/cubilin plasma membrane transport mechanism. The internalised DBP molecules then bind to actin and thus provide an intracellular array of high affinity binding sites for its specific ligand, 25(OH)D. There is evidence that the residence time for DBP in muscle cells is short and that it undergoes proteolytic degradation, releasing bound 25(OH)D. The processes of internalisation of DBP and its intracellular residence time, bound to actin, appear to be regulated. To explore whether 1,25-dihydroxycholecalciferol (calcitriol) has any effect on this process, cell cultures of myotubes and primary skeletal muscle fibers were incubated in a medium containing 10-10M calcitriol but with no added DBP. After 3h pre-incubation with calcitriol, the net uptake of 25(OH)D by these calcitriol-treated cells over a further 4h was significantly greater than that in vehicle-treated control cells. This was accompanied by a significant increase in intracellular DBP protein. However, after 16h of pre-incubation with calcitriol, the muscle cells showed a significantly depressed ability to accumulate 25(OH)D compared to control cells over a further 4 or 16hours. These effects of pre-incubation with calcitriol were abolished in fibers from VDR-knockout mice. The effect was also abolished by the addition of 4,4'-diisothiocyano-2,2'-stilbenedisulfonic acid (DIDS), which inhibits chloride channel opening. Incubation of C2 myotubes with calcitriol also significantly reduced retention of previously accumulated 25(OH)D after 4 or 8h. It is concluded from these in vitro studies that calcitriol can modify the DBP-dependent uptake and release of 25(OH)D by skeletal muscle cells in a manner that suggests some inducible change in the function of these cells.