Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

A Global Constraint On Relative Rotation To Avoid Lumped Compliant Mechanisms In Topology Optimization

Some of the well known formulations for topology optimization of compliant mechanisms could lead to lumped compliant mechanisms. In lumped compliance, most of the elastic deformation in a mechanism occurs at few points, while rest of the mechanism remains more or less rigid. Such points are referred to as point-flexures. It has been noted in literature that high relative rotation is associated with point-flexures. In literature we also find a formulation of local constraint on relative rotations to avoid lumped compliance. However it is well known that a global constraint is easier to handle than a local constraint, by a numerical optimization algorithm. The current work presents a way of putting global constraint on relative rotations. This constraint is also simpler to implement since it uses linearized rotation at the center of finite-elements, to compute relative rotations. I show the results obtained by using this constraint oil the following benchmark problems - displacement inverter and gripper.

The nature of non-equilibrium solidification of undercooled Agsingle bondCu alloys in contact with primary copper dendrites

Liquids of silver-copper alloys with near eutectic compositions embedded in a copper matrix were undercooled. The structural and microstructural investigations of these alloys solidified from undercooled temperature indicated the absence of both the eutectic reaction and diffusionless transformation below the equal free energy curve (T0). Instead the liquid maintained local equilibrium with the copper dendrites continuously until it intersected the extended solidus of the silver rich solid solution.

Topology-Guided Tessellation Of Quadratic Elements

Topology-based methods have been successfully used for the analysis and visualization of piecewise-linear functions defined on triangle meshes. This paper describes a mechanism for extending these methods to piecewise-quadratic functions defined on triangulations of surfaces. Each triangular patch is tessellated into monotone regions, so that existing algorithms for computing topological representations of piecewise-linear functions may be applied directly to the piecewise-quadratic function. In particular, the tessellation is used for computing the Reeb graph, a topological data structure that provides a succinct representation of level sets of the function.

On The Topology Of Manifolds With Positive Isotropic Curvature

We show that a closed orientable Riemannian n-manifold, n >= 5, with positive isotropic curvature and free fundamental group is homeomorphic to the connected sum of copies of Sn-1 x S-1.

Non-equilibrium trapping in Al---In system

The comparative trapping efficiency of indium in aluminium solid solution during vapour deposition and rapid solidification at room temperature was evaluated by studying the nature of the decomposition of the metastable solid solution. A spinodal decomposition was observed in the case of rapid solidification while vapour deposited foils decompose by a nucleation and growth type of process. From this we conclude that rapid solidification effects a more efficient non-equilibrium trapping of indium in aluminium solid solution compared to vapour deposition at room temperature.

Turbulence measurements in an equilibrium axisymmetric wall jet

This paper reports measurements of turbulent quantities in an axisymmetric wall jet subjected to an adverse pressure gradient in a conical diffuser, in such a way that a suitably defined pressure-gradient parameter is everywhere small. Self-similarity is observed in the mean velocity profile, as well as the profiles of many turbulent quantities at sufficiently large distances from the injection slot. Autocorrelation measurements indicate that, in the region of turbulent production, the time scale of ν fluctuations is very much smaller than the time scale of u fluctuations. Based on the data on these time scales, a possible model is proposed for the Reynolds stress. One-dimensional energy spectra are obtained for the u, v and w components at several points in the wall jet. It is found that self-similarity is exhibited by the one-dimensional wavenumber spectrum of $\overline{q^2}(=\overline{u^2}+\overline{v^2}+\overline{w^2})$, if the half-width of the wall jet and the local mean velocity are used for forming the non-dimensional wavenumber. Both the autocorrelation curves and the spectra indicate the existence of periodicity in the flow. The rate of dissipation of turbulent energy is estimated from the $\overline{q^2}$ spectra, using a slightly modified version of a previously suggested method.

Topology-based denoising of chaos

In this article, we propose a denoising algorithm to denoise a time series y(i) = x(i) + e(i), where {x(i)} is a time series obtained from a time- T map of a uniformly hyperbolic or Anosov flow, and {e(i)} a uniformly bounded sequence of independent and identically distributed (i.i.d.) random variables. Making use of observations up to time n, we create an estimate of x(i) for i<n. We show under typical limiting behaviours of the orbit and the recurrence properties of x(i), the estimation error converges to zero as n tends to infinity with probability 1.

Determination of skin friction in strong pressure-gradient equilibrium and near-equilibrium turbulent boundary layers

The conventional Clauser-chart method for determination of local skin friction in zero or weak pressure-gradient turbulent boundary layer flows fails entirely in strong pressure-gradient situations. This failure occurs due to the large departure of the mean velocity profile from the universal logarithmic law upon which the conventional Clauser-chart method is based. It is possible to extend this method,even for strong pressure-gradient situations involving equilibrium or near-equilibrium turbulent boundary layers by making use of the so-called non-universal logarithmic laws. These non-universal log laws depend on the local strength of the pressure gradient and may be regarded as perturbations of the universal log law.The present paper shows that the modified Clauser-chart method, so developed, yields quit satisfactory results in terms of estimation of local skin friction in strongly accelerated or retarded equilibrium and near-equilibrium turbulent boundary layers that are not very close to relaminarization or separation.

CFD analysis of high frequency miniature pulse tube refrigerators for space applications with thermal non-equilibrium model

High frequency, miniature, pulse tube cryocoolers are extensively used in space applications because of their simplicity. Parametric studies of inertance type pulse tube cooler are performed with different length-to-diameter ratios of the pulse tube with the help of the FLUENT (R) package. The local thermal non-equilibrium of the gas and the matrix is taken into account for the modeling of porous zones, in addition to the wall thickness of the components. Dynamic characteristics and the actual mechanism of energy transfer in pulse are examined with the help of the pulse tube wall time constant. The heat interaction between pulse tube wall and the oscillating gas, leading to surface heat pumping, is quantified. The axial heat conduction is found to reduce the performance of the pulse tube refrigerator. The thermal non-equilibrium predicts a higher cold heat exchanger temperature compared to thermal equilibrium. The pressure drop through the porous medium has a strong non-linear effect due to the dominating influence of Forchheimer term over that of the linear Darcy term at high operating frequencies. The phase angle relationships among the pressure, temperature and the mass flow rate in the porous zones are also important in determining the performance of pulse tuberefrigerator.

On the coalescence of twisted flux tubes

We study the problem of the coalescence of twisted flux tubes by assuming that the azimuthal field lines reconnect at a current sheet during the coalescence process and everywhere else the magnetic field is frozen in the fluid. We derive relations connecting the topology of the coalesced flux tube with the topologies of the initial flux tubes, and then obtain a structure equation for calculating the field configuration of the coalesced flux tube from the given topology. Some solutions for the two extreme cases of low-β plasma and high-β plasma are discussed. The coalesced flux tube has less twist than the initial flux tube. Magnetic helicity is found to be exactly conserved during the coalescence, but the assumptions in the model put a constraint on the energy dissipation so that we do not get a relaxation to the minimum-energy Taylor state in the low-β case. It is pointed out that the structure equation connecting the topology and the equilibrium configuration is quite general and can be of use in many two-dimensional flux tube problems.

Electrical switching and thermal studies on Ge22Te78-xIx chalcohalide glasses: The effect of iodine on network-topology

Investigations on the electrical switching behavior and thermal studies using Alternating Differential Scanning Calorimetry have been undertaken on bulk, melt-quenched Ge22Te78-,Is (3 <= x <= 10) chalcohalide glasses. All the glasses studied have been found to exhibit memory-type electrical switching. The threshold voltages of Ge22Te78-I-x(x) glasses have been found to increase with the addition of iodine and the composition dependence of threshold voltages of Ge22Te78-xIx glasses exhibits a cusp at 5 at.% of iodine. Also, the variation with composition of the glass transition temperature (Tg) of Ge22Te78-I-x(x) glasses, exhibits a broad hump around this composition. Based on the present results, the composition x = 5 has been identified as the inverse rigidity percolation threshold at which Ge22Te78-I-x(x) glassy system exhibits a change from a stressed rigid amorphous solid to a flexible polymeric glass. Further, a sharp minimum is seen in the composition dependence of non-reversing enthalpy (Delta H-nr) of Ge22Te78-I-x(x) glasses at x = 5, which is suggestive of a thermally reversing window at this composition. (C) 2007 Elsevier Ltd. All rights reserved.

Solution structures of conformationally equilibrium forms of holo-acyl carrier protein (PfACP) from Plasmodium falciparum provides insight into the mechanism of activation of ACPs

Acyl Carrier Protein (ACP) from the malaria parasite, Plasmodium falciparum (PfACP) in its holo form is found to exist in two conformational states in solution. Unique 3D solution structures of holo-PfACP have been determined for both equilibrium conformations, using high-resolution NMR methods. Twenty high-resolution solution structures for each of the two forms of holo-PfACP have been determined on the basis of 1226 and 1218 unambiguously assigned NOEs (including NOEs between 4 '-phosphopantetheine prosthetic group (4 '-PP) and protein), 55 backbone dihedral angles and 26 hydrogen bonds. The atomic rmsd values of the determined structures of two equilibrium forms, about the mean coordinates of the backbone and heavy atoms, are 0.48 +/- 0.09 and 0.92 +/- 0.10 and 0.49 +/- 0.08 and 0.97 +/- 0.11 angstrom, respectively. The interaction of 4 '-PP with the polypeptide backbone is reported here for the first time for any of the ACPs. The structures of holo-PfACP consist of three well-defined helices that are tightly packed. The structured regions of the molecule are stabilized by extensive hydrophobic interactions. The difference between the two forms arises from a reorientation of the 4 '-PP group. The enthalpy difference between the two forms, although small, implies that a conformational switch is essential for the activation of holo-ACP. Sequence and structures of holo-PfACP have been compared with those of the ACPs from type I and type II fatty acid biosynthesis pathways (FAS), in particular with the ACP from rat and the butyryl-ACP from E. coli. The PfACP structure, thus determined has several novel features hitherto not seen in other ACPs.

Composition of aqueous solutions in equilibrium with sulfides and oxides of iron at 3500°C

Solutions of potassium chloride (pH-buffered and 1-molat) equilibrated at 350°C with pyrrhotite, pyrite, and magnetite contained approximately 1 millimole of reduced sulfur and less than 0.1 millimole of oxidized sulfur per kilogram. Similar solutions equilibrated with pyrite, magnetite, and hematite contained approximately 1 millimole of reduced sulfur, but 3 to 6 millimoles of oxidized sulfur per kilogram. Both types of solutions contained less than 0.1 millimole of iron per kilogram at pH ≥ 6 and approximately 100 millimoles per kilogram at pH 2.