A novel approach for large-scale polypeptide folding based on elastic networks using continuous optimization

We present a new computationally efficient method for large-scale polypeptide folding using coarse-grained elastic networks and gradient-based continuous optimization techniques. The folding is governed by minimization of energy based on Miyazawa–Jernigan contact potentials. Using this method we are able to substantially reduce the computation time on ordinary desktop computers for simulation of polypeptide folding starting from a fully unfolded state. We compare our results with available native state structures from Protein Data Bank (PDB) for a few de-novo proteins and two natural proteins, Ubiquitin and Lysozyme. Based on our simulations we are able to draw the energy landscape for a small de-novo protein, Chignolin. We also use two well known protein structure prediction software, MODELLER and GROMACS to compare our results. In the end, we show how a modification of normal elastic network model can lead to higher accuracy and lower time required for simulation.

Low Power continuous time common mode sensing for common mode feedback circuits

Continuous common mode feedback (CMFB) circuits having high input impedance and low distortion are proposed. The proposed circuits are characterized for 0.18 mu m CMOS process with 1.8 V supply. Simulation results indicate that the proposed common mode detector consumes no standby power and CMFB circuit consumes 27-34% less power than previous high swing CMFB circuits.

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.

Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework

Partition of unity methods, such as the extended finite element method, allows discontinuities to be simulated independently of the mesh (Int. J. Numer. Meth. Engng. 1999; 45:601-620). This eliminates the need for the mesh to be aligned with the discontinuity or cumbersome re-meshing, as the discontinuity evolves. However, to compute the stiffness matrix of the elements intersected by the discontinuity, a subdivision of the elements into quadrature subcells aligned with the discontinuity is commonly adopted. In this paper, we use a simple integration technique, proposed for polygonal domains (Int. J. Nuttier Meth. Engng 2009; 80(1):103-134. DOI: 10.1002/nme.2589) to suppress the need for element subdivision. Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics and a multi-material problem show that the proposed method yields accurate results. Owing to its simplicity, the proposed integration technique can be easily integrated in any existing code. Copyright (C) 2010 John Wiley & Sons, Ltd.

Three-Dimensional Unsteady Flow With Heat and Mass Transfer Over a Continuous Stretching Surface

A numerical solution of the unsteady boundary layer equations under similarity assumptions is obtained. The solution represents the three-dimensional unsteady fluid motion caused by the time-dependent stretching of a flat boundary. It has been shown that a self-similar solution exists when either the rate of stretching is decreasing with time or it is constant. Three different numerical techniques are applied and a comparison is made among them as well as with earlier results. Analysis is made for various situations like deceleration in stretching of the boundary, mass transfer at the surface, saddle and nodal point flows, and the effect of a magnetic field. Both the constant temperature and constant heat flux conditions at the wall have been studied.

Continuous-cooling method of specific heat measurement in the temperature range 100-300 K

The authors have developed a simple continuous-cooling method to determine specific heat of liquids and solids in the temperature range 100-300 K. The technique employs very simple instrumentation and continuously records the sample temperature as it cools to the bath temperature through a calibrated heat link. They have obtained specific heat values which agree with the reported data to within 3% for the samples investigated. This method also facilitates easy detection of abrupt changes in specific heat, as demonstrated in the observation of glass transition in some organic glass-forming systems. The method is sensitive to the study of relaxing heat capacity in supercooled liquids.

Integration-free algorithms for fixed end-point regulator problems

A new formulation is suggested for the fixed end-point regulator problem, which, in conjunction with the recently developed integration-free algorithms, provides an efficient means of obtaining numerical solutions to such problems.

An accurate numerical integration scheme for finite rotations using rotation vector parametrization

A numerical integration procedure for rotational motion using a rotation vector parametrization is explored from an engineering perspective by using rudimentary vector analysis. The incremental rotation vector, angular velocity and acceleration correspond to different tangent spaces of the rotation manifold at different times and have a non-vectorial character. We rewrite the equation of motion in terms of vectors lying in the same tangent space, facilitating vector space operations consistent with the underlying geometric structure. While any integration algorithm (that works within a vector space setting) may be used, we presently employ a family of explicit Runge-Kutta algorithms to solve this equation. While this work is primarily motivated out of a need for highly accurate numerical solutions of dissipative rotational systems of engineering interest, we also compare the numerical performance of the present scheme with some of the invariant preserving schemes, namely ALGO-C1, STW, LIEMIDEA] and SUBCYC-M. Numerical results show better local accuracy via the present approach vis-a-vis the preserving algorithms. It is also noted that the preserving algorithms do not simultaneously preserve all constants of motion. We incorporate adaptive time-stepping within the present scheme and this in turn enables still higher accuracy and a `near preservation' of constants of motion over significantly longer intervals. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

Solution concepts in continuous-kernel multicriteria games

This paper considers nonzero-sum multicriteria games with continuous kernels. Solution concepts based on the notions of Pareto optimality, equilibrium, and security are extended to these games. Separate necessary and sufficient conditions and existence results are presented for equilibrium, Pareto-optimal response, and Pareto-optimal security strategies of the players.

Study Of Fluid-Flow And Residence-Time Distribution In A Continuous Slab Casting Mold

The fluid-flow pattern and residence-time distribution (r.t.d.) of the fluid in a continuous casting mould have been studied using a water model. The two recirculating zones below the discharge ports have been found to be asymmetric. The effect of casting speed, discharge port diameter, shroud well depth and the immersion depth on r.t.d. have been investigated. The r.t.d. curve has been well represented by a model of two backmix cells of equal volume in series. The exist of the fluid has been found to be non-uniform across the cross-section of the mould. The fluid-flow pattern has been observed to change with time in a random fashion. Dead volume of upto 31.8% has been found with smaller discharge ports.

Performance comparison of batch and continuous flow surface aeration systems

The oxygen transfer rate and the corresponding power requirement to operate the rotor are vital for design and scale-up of surface aerators. The aeration process can be analyzed in two ways such as batch and continuous systems. The process behaviors of batch and continuous flow systems are different from each other. The experimental and numerical results obtained through the batch systems cannot be relied on and applied for the designing of the continuous aeration tank. Based on the experimentation on batch and continuous type systems, the present work compares the performance of both the batch and continuous surface aeration systems in terms of their oxygen transfer capacity and power consumption. A simulation equation developed through experimentation has shown that continuous flow surface aeration systems are taking more energy than the batch systems. It has been found that batch systems are economical and better for the field application but not feasible where large quantity of wastewater is produced.

A Deterministic Optimization Approach to Protein Sequence Design Using Continuous Models

Determining the sequence of amino acid residues in a heteropolymer chain of a protein with a given conformation is a discrete combinatorial problem that is not generally amenable for gradient-based continuous optimization algorithms. In this paper we present a new approach to this problem using continuous models. In this modeling, continuous "state functions" are proposed to designate the type of each residue in the chain. Such a continuous model helps define a continuous sequence space in which a chosen criterion is optimized to find the most appropriate sequence. Searching a continuous sequence space using a deterministic optimization algorithm makes it possible to find the optimal sequences with much less computation than many other approaches. The computational efficiency of this method is further improved by combining it with a graph spectral method, which explicitly takes into account the topology of the desired conformation and also helps make the combined method more robust. The continuous modeling used here appears to have additional advantages in mimicking the folding pathways and in creating the energy landscapes that help find sequences with high stability and kinetic accessibility. To illustrate the new approach, a widely used simplifying assumption is made by considering only two types of residues: hydrophobic (H) and polar (P). Self-avoiding compact lattice models are used to validate the method with known results in the literature and data that can be practically obtained by exhaustive enumeration on a desktop computer. We also present examples of sequence design for the HP models of some real proteins, which are solved in less than five minutes on a single-processor desktop computer Some open issues and future extensions are noted.

Continuous ordering below the tricritical point: A critical test for the kinetic pathway in Fe-12at.% Si alloy

A critical test has been presented to establish the nature of the kinetic pathways for the decomposition of Fe-12 at.% Si alloy below the metastable tricritical point. The results, based on the measurements of saturation magnetization, establish that a congruent ordering from B2 --> D0(3) precedes the development of a B2 + D0(3) two-phase field, consistent with the predictions in 1976 of Allen and Cahn.

Supermacroporous Polymer-Based Cryogel Bioreactor for Monoclonal Antibody Production in Continuous Culture Using Hybridoma Cells

Cryogel matrices composed of different polymeric blends were synthesized, yielding a unique combination of hydrophilicity and hydrophobicity with the presence or absence of charged surface. Four such cryogel matrices composed of polyacrylamide-chitosan (PAAC), poly(N-isopropylacrylamide)-chitosan, polyacrylonitrile (PAN), and poly(N-isopropylacrylamide) were tested for growth of different hybridoma cell lines and production of antibody in static culture. All the matrices were capable for the adherence of hybridoma cell lines 6A4D7, B7B10, and H9E10 to the polymeric surfaces as well as for the efficient monoclonal antibody (mAb) production. PAAC proved to be relatively better in terms of both mAb production and cell growth. Further, PAAC cryogel was designed into three different formats, monolith, disks, and beads, and used as packing material for packed-bed bioreactor. Longterm cultivation of 6A4D7 cell line on PAAC cryogel scaffold in all the three formats could be successfully done for a period of 6 weeks under static conditions. Continuous packed-bed bioreactor was setup using 6A4D7 hybridoma cell line in the three reactor formats. The reactors ran continuously for a period of 60 days during which mAb production and metabolism of cells in the bioreactors were monitored periodically. The monolith bioreactor performed most efficiently over a period of 60 days and produced a total of 57.5 mg of antibody in the first 30 days (in 500 mL) with a highest concentration of 115 mu g mL(-1), which is fourfold higher than t-flask culture. The results demonstrate that appropriate chemistry and geometry of the bioreactor matrix for cell growth and immobilization can enhance the reactor productivity. (C) 2010 American Institute of Chemical Engineers Biotechnol. Prog., 27: 170-180, 2011