Combining structure and dynamics: non-denaturing high-pressure effect on lysozyme in solution

Small-angle X-ray scattering (SAXS) and elastic and quasi-elastic neutron scattering techniques were used to investigate the high-pressure-induced changes on interactions, the low-resolution structure and the dynamics of lysozyme in solution. SAXS data, analysed using a global-fit procedure based on a new approach for hydrated protein form factor description, indicate that lysozyme completely maintains its globular structure up to 1500 bar, but significant modi. cations in the protein-protein interaction potential occur at approximately 600-1000 bar. Moreover, the mass density of the protein hydration water shows a clear discontinuity within this pressure range. Neutron scattering experiments indicate that the global and the local lysozyme dynamics change at a similar threshold pressure. A clear evolution of the internal protein dynamics from diffusing to more localized motions has also been probed. Protein structure and dynamics results have then been discussed in the context of protein-water interface and hydration water dynamics. According to SAXS results, the new configuration of water in the first hydration layer induced by pressure is suggested to be at the origin of the observed local mobility changes.

2-Aminopurine non-radiative decay and emission in aqueous solution: A theoretical study

The minimum energy path along the lowest-lying pi pi* excited state of 2-aminopurine was calculated to elucidate the mechanisms of radiationless decay and emission in water. The sequential Monte Carlo quantum mechanics approach with a multiconfigurational and perturbative description of the wave function was employed to compute the minimum, transition state, and conical intersection. It was found that the barrier in the potential energy surface to access the conical intersection funnel increases in aqueous environment, making the system prone to enlarge the emission yield. These results rationalize the observed enhancement of emission in 2-aminopurine upon increasing of the solvent polarity. (c) 2008 Elsevier B.V. All rights reserved.

Non-LTE Monte Carlo radiative transfer. II. Nonisothermal solutions for viscous Keplerian disks

We discuss the basic hydrodynamics that determines the density structure of the disks around hot stars. Observational evidence supports the idea that these disks are Keplerian (rotationally supported) gaseous disks. A popular scenario in the literature, which naturally leads to the formation of Keplerian disks, is the viscous decretion model. According to this scenario, the disks are hydrostatically supported in the vertical direction, while the radial structure is governed by the viscous transport. This suggests that the temperature is one primary factor that governs the disk density structure. In a previous study we demonstrated, using three-dimensional non-LTE Monte Carlo simulations, that viscous Keplerian disks can be highly nonisothermal. In this paper we build on our previous work and solve the full problem of the steady state nonisothermal viscous diffusion and vertical hydrostatic equilibrium. We find that the self-consistent solution departs significantly from the analytic isothermal density, with potentially large effects on the emergent spectrum. This implies that nonisothermal disk models must be used for a detailed modeling of Be star disks.

Lower semicontinuity of attractors for non-autonomous dynamical systems

This paper is concerned with the lower semicontinuity of attractors for semilinear non-autonomous differential equations in Banach spaces. We require the unperturbed attractor to be given as the union of unstable manifolds of time-dependent hyperbolic solutions, generalizing previous results valid only for gradient-like systems in which the hyperbolic solutions are equilibria. The tools employed are a study of the continuity of the local unstable manifolds of the hyperbolic solutions and results on the continuity of the exponential dichotomy of the linearization around each of these solutions.

Non-autonomous semilinear evolution equations with almost sectorial operators

Inspired by the theory of semigroups of growth a, we construct an evolution process of growth alpha. The abstract theory is applied to study semilinear singular non-autonomous parabolic problems. We prove that. under natural assumptions. a reasonable concept of solution can be given to Such semilinear singularly non-autonomous problems. Applications are considered to non-autonomous parabolic problems in space of Holder continuous functions and to a parabolic problem in a domain Omega subset of R(n) with a one dimensional handle.

The constrained compartmentalized knapsack problem: mathematical models and solution methods

The constrained compartmentalized knapsack problem can be seen as an extension of the constrained knapsack problem. However, the items are grouped into different classes so that the overall knapsack has to be divided into compartments, and each compartment is loaded with items from the same class. Moreover, building a compartment incurs a fixed cost and a fixed loss of the capacity in the original knapsack, and the compartments are lower and upper bounded. The objective is to maximize the total value of the items loaded in the overall knapsack minus the cost of the compartments. This problem has been formulated as an integer non-linear program, and in this paper, we reformulate the non-linear model as an integer linear master problem with a large number of variables. Some heuristics based on the solution of the restricted master problem are investigated. A new and more compact integer linear model is also presented, which can be solved by a branch-and-bound commercial solver that found most of the optimal solutions for the constrained compartmentalized knapsack problem. On the other hand, heuristics provide good solutions with low computational effort. (C) 2011 Elsevier BM. All rights reserved.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids. (C) 2010 Elsevier B.V. All rights reserved.

Numerical solution of the PTT constitutive equation for unsteady three-dimensional free surface flows

This work deals with the development of a numerical technique for simulating three-dimensional viscoelastic free surface flows using the PTT (Phan-Thien-Tanner) nonlinear constitutive equation. In particular, we are interested in flows possessing moving free surfaces. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are considered. The PTT equation is solved by a high order method, which requires the calculation of the extra-stress tensor on the mesh contours. To validate the numerical technique developed in this work flow predictions for fully developed pipe flow are compared with an analytic solution from the literature. Then, results of complex free surface flows using the FIT equation such as the transient extrudate swell problem and a jet flowing onto a rigid plate are presented. An investigation of the effects of the parameters epsilon and xi on the extrudate swell and jet buckling problems is reported. (C) 2010 Elsevier B.V. All rights reserved.

An unstructured CVFEM and moving interface algorithm for non-Newtonian Hele-Shaw flows in injection molding

Purpose - The purpose of this paper is to develop a novel unstructured simulation approach for injection molding processes described by the Hele-Shaw model. Design/methodology/approach - The scheme involves dual dynamic meshes with active and inactive cells determined from an initial background pointset. The quasi-static pressure solution in each timestep for this evolving unstructured mesh system is approximated using a control volume finite element method formulation coupled to a corresponding modified volume of fluid method. The flow is considered to be isothermal and non-Newtonian. Findings - Supporting numerical tests and performance studies for polystyrene described by Carreau, Cross, Ellis and Power-law fluid models are conducted. Results for the present method are shown to be comparable to those from other methods for both Newtonian fluid and polystyrene fluid injected in different mold geometries. Research limitations/implications - With respect to the methodology, the background pointset infers a mesh that is dynamically reconstructed here, and there are a number of efficiency issues and improvements that would be relevant to industrial applications. For instance, one can use the pointset to construct special bases and invoke a so-called ""meshless"" scheme using the basis. This would require some interesting strategies to deal with the dynamic point enrichment of the moving front that could benefit from the present front treatment strategy. There are also issues related to mass conservation and fill-time errors that might be addressed by introducing suitable projections. The general question of ""rate of convergence"" of these schemes requires analysis. Numerical results here suggest first-order accuracy and are consistent with the approximations made, but theoretical results are not available yet for these methods. Originality/value - This novel unstructured simulation approach involves dual meshes with active and inactive cells determined from an initial background pointset: local active dual patches are constructed ""on-the-fly"" for each ""active point"" to form a dynamic virtual mesh of active elements that evolves with the moving interface.

Numerical Solution of the Upper-Convected Maxwell Model for Three-Dimensional Free Surface Flows

This work presents a finite difference technique for simulating three-dimensional free surface flows governed by the Upper-Convected Maxwell (UCM) constitutive equation. A Marker-and-Cell approach is employed to represent the fluid free surface and formulations for calculating the non-Newtonian stress tensor on solid boundaries are developed. The complete free surface stress conditions are employed. The momentum equation is solved by an implicit technique while the UCM constitutive equation is integrated by the explicit Euler method. The resulting equations are solved by the finite difference method on a 3D-staggered grid. By using an exact solution for fully developed flow inside a pipe, validation and convergence results are provided. Numerical results include the simulation of the transient extrudate swell and the comparison between jet buckling of UCM and Newtonian fluids.

Investigation of AISI 1040 steel corrosion H(2)S solution containing chloride ions by digital image processing coupled with in electrochemical techniques

This paper presents a study of AISI 1040 steel corrosion in aqueous electrolyte of acetic acid buffer containing 3.1 and 31 x 10(-3) mol dm(-3) of Na(2)S in both the presence and absence of 3.5 wt.% NaCl. This investigation of steel corrosion was carried out using potential polarization, and open-circuit and in situ optical microscopy. The morphological analysis and classification of types of surface corrosion damage by digital image processing reveals grain boundary corrosion and shows a non-uniform sulfide film growth, which occurs preferentially over pearlitic grains through successive formation and dissolution of the film. (C) 2011 Elsevier Ltd. All rights reserved.

Structural and thermodynamic analysis of thrombin:suramin interaction in solution and crystal phases

Suramin is a hexasulfonated naphthylurea which has been recently characterized as a non-competitive inhibitor of human alpha-thrombin activity over fibrinogen, although its binding site and mode of interaction with the enzyme remain elusive. Here, we determined two X-ray structure of the thrombin: suramin complex, refined at 2.4 angstrom resolution. While a single thrombin: suramin complex was found in the asymmetric unit cell of the crystal, some of the crystallographic contacts with symmetrically related molecules are mediated by both the enzyme and the ligand. Molecular dynamics simulations with the 1:1 complex demonstrate a large rearrangement of suramin in the complex, but with the protein scaffold and the more extensive protein-ligand regions keep unchanged. Small-angle X-ray scattering measurements at high micromolar concentration demonstrate a suramin-induced dimerization of the enzyme. These data indicating a dissimilar binding mode in the monomeric and oligomeric states, with a monomeric, 1:1 complex to be more likely to exist at the thrombin physiological, nanomolar concentration range. Collectively, close understanding on the structural basis for interaction is given which might establish a basis for design of suramin analogues targeting thrombin. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.