Sharp interface limits of phase-field models

The use of continuum phase-field models to describe the motion of well-defined interfaces is discussed for a class of phenomena, that includes order/disorder transitions, spinodal decomposition and Ostwald ripening, dendritic growth, and the solidification of eutectic alloys. The projection operator method is used to extract the ``sharp interface limit'' from phase field models which have interfaces that are diffuse on a length scale $\xi$. In particular,phase-field equations are mapped onto sharp interface equations in the limits $\xi \kappa \ll 1$ and $\xi v/D \ll 1$, where $\kappa$ and $v$ are respectively the interface curvature and velocity and $D$ is the diffusion constant in the bulk. The calculations provide one general set of sharp interface equations that incorporate the Gibbs-Thomson condition, the Allen-Cahn equation and the Kardar-Parisi-Zhang equation.Comment: 17 pages, 9 figure

Entire solutions of hydrodynamical equations with exponential dissipation

We consider a modification of the three-dimensional Navier--Stokes equations and other hydrodynamical evolution equations with space-periodic initial conditions in which the usual Laplacian of the dissipation operator is replaced by an operator whose Fourier symbol grows exponentially as \ue ^{|k|/\kd} at high wavenumbers $|k|$. Using estimates in suitable classes of analytic functions, we show that the solutions with initially finite energy become immediately entire in the space variables and that the Fourier coefficients decay faster than \ue ^{-C(k/\kd) \ln (|k|/\kd)} for any $C<1/(2\ln 2)$. The same result holds for the one-dimensional Burgers equation with exponential dissipation but can be improved: heuristic arguments and very precise simulations, analyzed by the method of asymptotic extrapolation of van der Hoeven, indicate that the leading-order asymptotics is precisely of the above form with $C= C_\star =1/\ln2$. The same behavior with a universal constant $C_\star$ is conjectured for the Navier--Stokes equations with exponential dissipation in any space dimension. This universality prevents the strong growth of intermittency in the far dissipation range which is obtained for ordinary Navier--Stokes turbulence. Possible applications to improved spectral simulations are briefly discussed.Comment: 29 pages, 3 figures, Comm. Math. Phys., in pres

Nature of complex singularities for the 2D Euler equation

A detailed study of complex-space singularities of the two-dimensional incompressible Euler equation is performed in the short-time asymptotic r\'egime when such singularities are very far from the real domain; this allows an exact recursive determination of arbitrarily many spatial Fourier coefficients. Using high-precision arithmetic we find that the Fourier coefficients of the stream function are given over more than two decades of wavenumbers by \hat F(\k) = C(\theta) k^{-\alpha} \ue ^ {-k \delta(\theta)}, where \k = k(\cos \theta, \sin \theta). The prefactor exponent $\alpha$, typically between 5/2 and 8/3, is determined with an accuracy better than 0.01. It depends on the initial condition but not on $\theta$. The vorticity diverges as $s^{-\beta}$, where $\alpha+\beta= 7/2$ and $s$ is the distance to the (complex) singular manifold. This new type of non-universal singularity is permitted by the strong reduction of nonlinearity (depletion) which is associated to incompressibility. Spectral calculations show that the scaling reported above persists well beyond the time of validity of the short-time asymptotics. A simple model in which the vorticity is treated as a passive scalar is shown analytically to have universal singularities with exponent $\alpha =5/2$.Comment: 22 pages, 24 figures, published version; a version of the paper with higher-quality figures is available at http://www.obs-nice.fr/etc7/euler.pd

Analytic solutions and Singularity formation for the Peakon b--Family equations

Using the Abstract Cauchy-Kowalewski Theorem we prove that the $b$-family equation admits, locally in time, a unique analytic solution. Moreover, if the initial data is real analytic and it belongs to $H^s$ with $s > 3/2$, and the momentum density $u_0 - u_{0,{xx}}$ does not change sign, we prove that the solution stays analytic globally in time, for $b\geq 1$. Using pseudospectral numerical methods, we study, also, the singularity formation for the $b$-family equations with the singularity tracking method. This method allows us to follow the process of the singularity formation in the complex plane as the singularity approaches the real axis, estimating the rate of decay of the Fourier spectrum

Lattice Boltzmann simulations of soft matter systems

This article concerns numerical simulations of the dynamics of particles immersed in a continuum solvent. As prototypical systems, we consider colloidal dispersions of spherical particles and solutions of uncharged polymers. After a brief explanation of the concept of hydrodynamic interactions, we give a general overview over the various simulation methods that have been developed to cope with the resulting computational problems. We then focus on the approach we have developed, which couples a system of particles to a lattice Boltzmann model representing the solvent degrees of freedom. The standard D3Q19 lattice Boltzmann model is derived and explained in depth, followed by a detailed discussion of complementary methods for the coupling of solvent and solute. Colloidal dispersions are best described in terms of extended particles with appropriate boundary conditions at the surfaces, while particles with internal degrees of freedom are easier to simulate as an arrangement of mass points with frictional coupling to the solvent. In both cases, particular care has been taken to simulate thermal fluctuations in a consistent way. The usefulness of this methodology is illustrated by studies from our own research, where the dynamics of colloidal and polymeric systems has been investigated in both equilibrium and nonequilibrium situations.Comment: Review article, submitted to Advances in Polymer Science. 16 figures, 76 page

Molecular Dynamics Simulation of Phosphorylated KID Post-Translational Modification

BACKGROUND:Kinase-inducible domain (KID) as transcriptional activator can stimulate target gene expression in signal transduction by associating with KID interacting domain (KIX). NMR spectra suggest that apo-KID is an unstructured protein. After post-translational modification by phosphorylation, KID undergoes a transition from disordered to well folded protein upon binding to KIX. However, the mechanism of folding coupled to binding is poorly understood. METHODOLOGY:To get an insight into the mechanism, we have performed ten trajectories of explicit-solvent molecular dynamics (MD) for both bound and apo phosphorylated KID (pKID). Ten MD simulations are sufficient to capture the average properties in the protein folding and unfolding. CONCLUSIONS:Room-temperature MD simulations suggest that pKID becomes more rigid and stable upon the KIX-binding. Kinetic analysis of high-temperature MD simulations shows that bound pKID and apo-pKID unfold via a three-state and a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound pKID folds in the order of KIX access, initiation of pKID tertiary folding, folding of helix alpha(B), folding of helix alpha(A), completion of pKID tertiary folding, and finalization of pKID-KIX binding. Our data show that the folding pathways of apo-pKID are different from the bound state: the foldings of helices alpha(A) and alpha(B) are swapped. Here we also show that Asn139, Asp140 and Leu141 with large Phi-values are key residues in the folding of bound pKID. Our results are in good agreement with NMR experimental observations and provide significant insight into the general mechanisms of binding induced protein folding and other conformational adjustment in post-translational modification

Polar or Apolar—The Role of Polarity for Urea-Induced Protein Denaturation

Urea-induced protein denaturation is widely used to study protein folding and stability; however, the molecular mechanism and driving forces of this process are not yet fully understood. In particular, it is unclear whether either hydrophobic or polar interactions between urea molecules and residues at the protein surface drive denaturation. To address this question, here, many molecular dynamics simulations totalling ca. 7 µs of the CI2 protein in aqueous solution served to perform a computational thought experiment, in which we varied the polarity of urea. For apolar driving forces, hypopolar urea should show increased denaturation power; for polar driving forces, hyperpolar urea should be the stronger denaturant. Indeed, protein unfolding was observed in all simulations with decreased urea polarity. Hyperpolar urea, in contrast, turned out to stabilize the native state. Moreover, the differential interaction preferences between urea and the 20 amino acids turned out to be enhanced for hypopolar urea and suppressed (or even inverted) for hyperpolar urea. These results strongly suggest that apolar urea–protein interactions, and not polar interactions, are the dominant driving force for denaturation. Further, the observed interactions provide a detailed picture of the underlying molecular driving forces. Our simulations finally allowed characterization of CI2 unfolding pathways. Unfolding proceeds sequentially with alternating loss of secondary or tertiary structure. After the transition state, unfolding pathways show large structural heterogeneity

An Evolutionary Trade-Off between Protein Turnover Rate and Protein Aggregation Favors a Higher Aggregation Propensity in Fast Degrading Proteins

We previously showed the existence of selective pressure against protein aggregation by the enrichment of aggregation-opposing ‘gatekeeper’ residues at strategic places along the sequence of proteins. Here we analyzed the relationship between protein lifetime and protein aggregation by combining experimentally determined turnover rates, expression data, structural data and chaperone interaction data on a set of more than 500 proteins. We find that selective pressure on protein sequences against aggregation is not homogeneous but that short-living proteins on average have a higher aggregation propensity and fewer chaperone interactions than long-living proteins. We also find that short-living proteins are more often associated to deposition diseases. These findings suggest that the efficient degradation of high-turnover proteins is sufficient to preclude aggregation, but also that factors that inhibit proteasomal activity, such as physiological ageing, will primarily affect the aggregation of short-living proteins

Unfolding Simulations of Holomyoglobin from Four Mammals: Identification of Intermediates and β-Sheet Formation from Partially Unfolded States

Myoglobin (Mb) is a centrally important, widely studied mammalian protein. While much work has investigated multi-step unfolding of apoMb using acid or denaturant, holomyoglobin unfolding is poorly understood despite its biological relevance. We present here the first systematic unfolding simulations of holoMb and the first comparative study of unfolding of protein orthologs from different species (sperm whale, pig, horse, and harbor seal). We also provide new interpretations of experimental mean molecular ellipticities of myoglobin intermediates, notably correcting for random coil and number of helices in intermediates. The simulated holoproteins at 310 K displayed structures and dynamics in agreement with crystal structures (R g ~1.48-1.51 nm, helicity ~75%). At 400 K, heme was not lost, but some helix loss was observed in pig and horse, suggesting that these helices are less stable in terrestrial species. At 500 K, heme was lost within 1.0-3.7 ns. All four proteins displayed exponentially decaying helix structure within 20 ns. The C- and F-helices were lost quickly in all cases. Heme delayed helix loss, and sperm whale myoglobin exhibited highest retention of heme and D/E helices. Persistence of conformation (RMSD), secondary structure, and ellipticity between 2-11 ns was interpreted as intermediates of holoMb unfolding in all four species. The intermediates resemble those of apoMb notably in A and H helices, but differ substantially in the D-, E- and F-helices, which interact with heme. The identified mechanisms cast light on the role of metal/cofactor in poorly understood holoMb unfolding. We also observed β-sheet formation of several myoglobins at 500 K as seen experimentally, occurring after disruption of helices to a partially unfolded, globally disordered state; heme reduced this tendency and sperm-whale did not display any sheet propensity during the simulations