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.

Transfer of saline tolerance from one strain of rice to another by injection of DNA

Direct injection of genomic DNA from salt tolerant cv. Pokkali into developing floral tillers on IR20 produced transgenic seeds similar to Pokkali in husk colour and which germinated well in 0.2 M NaCl and had a 4-6-fold higher proline content.

Synthesis, X-ray structure and catalytic properties of a copper(II) Schiff base complex modeling the activity of the Cu-B site of dopamine beta-hydroxylase

The copper(II) complex [Cu(salgly) (bpy)] . 4H(2)O (1), where salgly is a tridentate glycinatosalicylaldimine Schiffbase Ligand, is prepared and structurally characterized. The complex is found to be catalytically active in the oxidation of ascorbic acid by dioxygen and the process is also effective in the presence of benzylamine giving benzaldehyde as a product, thus modeling the activity of the Cu-B site of dopamine beta-hydroxylase. (C) 2000 Elsevier Science S.A. All rights reserved.

Modeling server-unreliability in closed queuing-networks

A method is presented to model server unreliability in closed queuing networks. Breakdowns and repairs of servers, assumed to be time-dependent, are modeled using virtual customers and virtual servers in the system. The problem is thus converted into a closed queue with all reliable servers and preemptive resume priority centers. Several recent preemptive priority approximations and an approximation of the one proposed are used in the analysis. This method has approximately the same computational requirements as that of mean-value analysis for a network of identical dimensions and is therefore very efficient

Natural Convection on a Horizontal Cone in a Porous Medium with Non-Uniform Wall Temperature/Concentration or Heat/Mass Flux and Suction/Injection

The steady natural convection flow on a horizontal cone embedded in a saturated porous medium with non-uniform wall temperature/concentration or heat/mass flux and suction/injection has been investigated. Non-similar solutions have been obtained. The nonlinear couple differential equations under boundary layer approximations governing the flow have been numerically solved. The Nusselt and Sherwood numbers are found to depend on the buoyancy forces, suction/injection rates, variation of wall temperature/concentration or heat/mass flux, Lewis number and the non-Darcy parameter.

Modeling of the sea surface for daylight imagery studies

Plastic-coated paper is shown to possess reflectivity characteristics quite similar to those of the surface of water. This correspondence has been used with a conversion factor to model a sea surface by means of plastic-coated paper. Such a paper model is then suitably illuminated and photographed, yielding physically simulated daylight imagery of the sea surface under controlled conditions. A simple example of sinusoidal surface simulation is described.

Studies on follicular growth in the immature rat and hamster: Effect of a single injection of gonadotropin or estrogen on the rate of3H-thymidine incorporation into ovarian DNAin vitro

Initiation of follicular growth by specific hormonal stimuli in ovaries of immature rats and hamsters was studied by determining the rate of incorporation of3H-thymidine into ovarian DNAin vitro. Incorporation was considered as an index of DNA synthesis and cell multiplication. A single injection of pregnant mare serum gonadotropin could thus maximally stimulate by 18 hr3H-thymidine incorporation into DNA of the ovary of immature hamsters. Neutralization of pregnant mare serum gonadotropin by an antiserum to ovine follicle stimulating hormone only during the initial 8–10 hr and not later could inhibit the increase in3H-thymidine incorporationin vitro observed at 18 hr, suggesting that the continued presence of gonadotropin stimulus was not necessary for this response. The other indices of follicular growth monitored such as ovarian weight, serum estradiol and uterine weight showed discernible increase at periods only after the above initial event. A single injection of estrogen (diethyl stilbesterol or estradiol-l7β) could similarly cause 18 hr later, a stimulation in the rate of incorporation of3H-thymidine into DNAin vitro in ovaries of immature rats. The presence of endogenous gonadotropins, however, was obligatory for observing this response to estrogen. Evidence in support of the above was two-fold: (i) administration of antiserum to follicle stimulating hormone or luteinizing hormone along with estrogen completely inhibited the increase in3H-thymidine incorporation into ovarian DNAin vitro; (ii) a radioimmunological measurement revealed following estrogen treatment, the presence of a higher concentration of endogenous follicle stimulating hormone in the ovary. Finally, administration of varying doses of ovine follicle stimulating hormone along with a constant dose of estrogen to immature rats produced a dose-dependent increment in the incorporation of3H-thymidine into ovarian DNAin vitro. These observations suggested the potentiality of this system for developing a sensitive bioassay for follicle stimulating hormone.

Experimental studies and phase field modeling of microstructure evolution during solidification with electromagnetic stirring

Thixocasting requires manufacturing of billets with non-dendritic microstructure. Aluminum alloy A356 billets were produced by rheocasting in a mould placed inside a linear electromagnetic stirrer. Subsequent heat treatment was used to produce a transition from rosette to globular microstructure. The current and the duration of stirring were explored as control parameters. Simultaneous induction heating of the billet during stirring was quantified using experimentally determined thermal profiles. The effect of processing parameters on the dendrite fragmentation was discussed. Corresponding computational modeling of the process was performed using phase-field modeling of alloy solidification in order to gain insight into the process of morphological changes of a solid during this process. A non-isothermal alloy solidification model was used for simulations. The morphological evolution under such imposed thermal cycles was simulated and compared with experimentally determined one. Suitable scaling using the thermosolutal diffusion distances was used to overcome computational difficulties in quantitative comparison at system scale. The results were interpreted in the light of existing theories of microstructure refinement and globularisation.

Complexation of siRNA with Dendrimer: A Molecular Modeling Approach

This paper reports the structural behavior and thermodynamics of the complexation of siRNA with poly(amidoamine) (PAMAM) dendrimers of generation 3 (G3) and 4 (G4) through fully atomistic molecular dynamics (MD) simulations accompanied by free energy calculations and inherent structure determination. We have also done simulation with one siRNA and two dendrimers (2 x G3 or 2xG4) to get the microscopic picture of various binding modes. Our simulation results reveal the formation of stable siRNA-dendrimer complex over nanosecond time scale. With the increase in dendrimcr generation, the charge ratio increases and hence the binding energy between siRNA and dendrimer also increases in accordance with available experimental measurements. Calculated radial distribution functions of amines groups of various subgenerations in a given generation of dendrimer and phosphate in backbone of siRNA reveals that one dendrimer of generation 4 shows better binding with siRNA almost wrapping the dendrimer when compared to the binding with lower generation dendrimer like G3. In contrast, two dendrimers of generation 4 show binding without siRNA wrapping the den-rimer because of repulsion between two dendrimers. The counterion distribution around the complex and the water molecules in the hydration shell of siRNA give microscopic picture of the binding dynamics. We see a clear correlation between water. counterions motions and the complexation i.e. the water molecules and counterions which condensed around siRNA are moved away from the siRNA backbone when dendrimer start binding to the siRNA back hone. As siRNA wraps/bind to the dendrimer counterions originally condensed onto siRNA (Na-1) and dendrimer (Cl-) get released. We give a quantitative estimate of the entropy of counterions and show that there is gain in entropy due to counterions release during the complexation. Furthermore, the free energy of complexation of IG3 and IG4 at two different salt concentrations shows that increase in salt concentration leads to the weakening of the binding affinity of siRNA and dendrimer.

Analysis of performance of a hot gas injection thrust vector control system

The complex three-dimensional flowfield produced by secondary injection of hot gases in a rocket nozzle for thrust vector control is analyzed by solving unsteady three-dimensional Euler equations with appropriate boundary conditions. Various system performance parameters like secondary jet amplification factor and axial thrust augmentation are deduced by integrating the nozzle wall pressure distributions obtained as part of the flowfield solution and compared with measurements taken in actual static tests. The agreement is good within the practical range of secondary injectant flow rates for thrust vector control applications.

Performance study of high operating temperature HgCdTe mid wave infrared detector through numerical modeling

The design of present generation uncooled Hg1-xCdxTe infrared photon detectors relies on complex heterostructures with a basic unit cell of type (n) under bar (+)/pi/(p) under bar (+). We present an analysis of double barrier (n) under bar (+)/pi/(p) under bar (+) mid wave infrared (x = 0.3) HgCdTe detector for near room temperature operation using numerical computations. The present work proposes an accurate and generalized methodology in terms of the device design, material properties, and operation temperature to study the effects of position dependence of carrier concentration, electrostatic potential, and generation-recombination (g-r) rates on detector performance. Position dependent profiles of electrostatic potential, carrier concentration, and g-r rates were simulated numerically. Performance of detector was studied as function of doping concentration of absorber and contact layers, width of both layers and minority carrier lifetime. Responsivity similar to 0.38 A W-1, noise current similar to 6 x 10(-14) A/Hz(1/2) and D* similar to 3.1 x 10(10)cm Hz(1/2) W-1 at 0.1 V reverse bias have been calculated using optimized values of doping concentration, absorber width and carrier lifetime. The suitability of the method has been illustrated by demonstrating the feasibility of achieving the optimum device performance by carefully selecting the device design and other parameters. (C) 2010 American Institute of Physics. doi:10.1063/1.3463379]