q-Gaussian based Smoothed Functional Algorithms for Stochastic Optimization

The q-Gaussian distribution results from maximizing certain generalizations of Shannon entropy under some constraints. The importance of q-Gaussian distributions stems from the fact that they exhibit power-law behavior, and also generalize Gaussian distributions. In this paper, we propose a Smoothed Functional (SF) scheme for gradient estimation using q-Gaussian distribution, and also propose an algorithm for optimization based on the above scheme. Convergence results of the algorithm are presented. Performance of the proposed algorithm is shown by simulation results on a queuing model.

Functional characterization, homology modeling and docking studies of beta-glucosidase responsible for bioactivation of cyanogenic hydroxynitrile glucosides from Leucaena leucocephala (subabul)

Glycosyl hydrolase family 1 beta-glucosidases are important enzymes that serve many diverse functions in plants including defense, whereby hydrolyzing the defensive compounds such as hydroxynitrile glucosides. A hydroxynitrile glucoside cleaving beta-glucosidase gene (Llbglu1) was isolated from Leucaena leucocephala, cloned into pET-28a (+) and expressed in E. coli BL21 (DE3) cells. The recombinant enzyme was purified by Ni-NTA affinity chromatography. The optimal temperature and pH for this beta-glucosidase were found to be 45 A degrees C and 4.8, respectively. The purified Llbglu1 enzyme hydrolyzed the synthetic glycosides, pNPGlucoside (pNPGlc) and pNPGalactoside (pNPGal). Also, the enzyme hydrolyzed amygdalin, a hydroxynitrile glycoside and a few of the tested flavonoid and isoflavonoid glucosides. The kinetic parameters K (m) and V (max) were found to be 38.59 mu M and 0.8237 mu M/mg/min for pNPGlc, whereas for pNPGal the values were observed as 1845 mu M and 0.1037 mu M/mg/min. In the present study, a three dimensional (3D) model of the Llbglu1 was built by MODELLER software to find out the substrate binding sites and the quality of the model was examined using the program PROCHEK. Docking studies indicated that conserved active site residues are Glu 199, Glu 413, His 153, Asn 198, Val 270, Asn 340, and Trp 462. Docking of rhodiocyanoside A with the modeled Llbglu1 resulted in a binding with free energy change (Delta G) of -5.52 kcal/mol on which basis rhodiocyanoside A could be considered as a potential substrate.

Smoothed functional and Quasi-Newton algorithms for routing in multi-stage queueing network with constraints

We consider the problem of optimal routing in a multi-stage network of queues with constraints on queue lengths. We develop three algorithms for probabilistic routing for this problem using only the total end-to-end delays. These algorithms use the smoothed functional (SF) approach to optimize the routing probabilities. In our model all the queues are assumed to have constraints on the average queue length. We also propose a novel quasi-Newton based SF algorithm. Policies like Join Shortest Queue or Least Work Left work only for unconstrained routing. Besides assuming knowledge of the queue length at all the queues. If the only information available is the expected end-to-end delay as with our case such policies cannot be used. We also give simulation results showing the performance of the SF algorithms for this problem.

Vibrational spectroscopy and density functional theory analysis of 3-O-caffeoylquinic acid

Density functional theory (DFT) calculations are being performed to investigate the geometric, vibrational, and electronic properties of the chlorogenic acid isomer 3-CQA (1R,3R,4S,5R)-3-{(2E)-3-(3,4-dihydroxyphenyl)prop-2-enoyl]oxy}-1,4, 5-trihydroxycyclohexanecarboxylic acid), a major phenolic compound in coffee. DFT calculations with the 6-311G(d,p) basis set produce very good results. The electrostatic potential mapped onto an isodensity surface has been obtained. A natural bond orbital analysis (NBO) has been performed in order to study intramolecular bonding, interactions among bonds, and delocalization of unpaired electrons. HOMO-LUMO studies give insights into the interaction of the molecule with other species. The calculated HOMO and LUMO energies indicate that a charge transfer occurs within the molecule. (C) 2012 Elsevier B.V. All rights reserved.

Multi-component, Self-assembled, Functional Soft Materials

Supramolecular chemistry is an emerging tool for devising materials that can perform specified functions. The self-assembly of facially amphiphilic bile acid molecules has been extensively utilized for the development of functional soft materials. Supramolecular hydrogels derived from the bile acid backbone act as useful templates for the intercalation of multiple components. Based on this, synthesis of gel-nanoparticle hybrid materials, photoluminescent coating materials, development of a new enzyme assay technique, etc. were achieved in the author's laboratory. The present account highlights some of these achievements.

A practical approach for cross-functional vehicle body weight optimization

The goal of optimization in vehicle design is often blurred by the myriads of requirements belonging to attributes that may not be quite related. If solutions are sought by optimizing attribute performance-related objectives separately starting with a common baseline design configuration as in a traditional design environment, it becomes an arduous task to integrate the potentially conflicting solutions into one satisfactory design. It may be thus more desirable to carry out a combined multi-disciplinary design optimization (MDO) with vehicle weight as an objective function and cross-functional attribute performance targets as constraints. For the particular case of vehicle body structure design, the initial design is likely to be arrived at taking into account styling, packaging and market-driven requirements. The problem with performing a combined cross-functional optimization is the time associated with running such CAE algorithms that can provide a single optimal solution for heterogeneous areas such as NVH and crash safety. In the present paper, a practical MDO methodology is suggested that can be applied to weight optimization of automotive body structures by specifying constraints on frequency and crash performance. Because of the reduced number of cases to be analyzed for crash safety in comparison with other MDO approaches, the present methodology can generate a single size-optimized solution without having to take recourse to empirical techniques such as response surface-based prediction of crash performance and associated successive response surface updating for convergence. An example of weight optimization of spaceframe-based BIW of an aluminum-intensive vehicle is given to illustrate the steps involved in the current optimization process.

Phase synchronization in brain networks derived from correlation between probabilities of recurrences in functional MRI data

It is increasingly being recognized that resting state brain connectivity derived from functional magnetic resonance imaging (fMRI) data is an important marker of brain function both in healthy and clinical populations. Though linear correlation has been extensively used to characterize brain connectivity, it is limited to detecting first order dependencies. In this study, we propose a framework where in phase synchronization (PS) between brain regions is characterized using a new metric ``correlation between probabilities of recurrence'' (CPR) and subsequent graph-theoretic analysis of the ensuing networks. We applied this method to resting state fMRI data obtained from human subjects with and without administration of propofol anesthetic. Our results showed decreased PS during anesthesia and a biologically more plausible community structure using CPR rather than linear correlation. We conclude that CPR provides an attractive nonparametric method for modeling interactions in brain networks as compared to standard correlation for obtaining physiologically meaningful insights about brain function.

Cytochrome P450-mediated metabolism in brain: functional roles and their implications

Introduction: Cytochromes P450 (P450) and associated monooxygenases are a family of heme proteins involved in metabolism of endogenous compounds (arachidonic acid, eicosanoids and prostaglandins) as also xenobiotics including drugs and environmental chemicals. Liver is the major organ involved in P450-mediated metabolism and hepatic enzymes have been characterized. Extrahepatic organs, such as lung, kidney and brain have the capability for biotransformation through P450 enzymes. Brain, including human brain, expresses P450 enzymes that metabolize xenobiotics and endogenous compounds. Areas covered: An overview of P450-mediated metabolism in brain is presented focusing on distinct differences seen in expression of P450 enzymes, generation of unique P450 enzymes in brain through alternate splicing and their consequences in terms of metabolism of psychoactive drugs and inflammatory prompts, such as leukotrienes, thus modulating inflammatory response. Expert opinion: The brain possesses unique P450s that metabolize drugs and endogenous compounds through pathways that are markedly different from that seen in liver indicating that extrapolation directly from liver to brain is not appropriate. It is therefore necessary to characterize the unique brain P450s and their ability to metabolize xenobiotics and endogenous compounds to better understand the functions of this important class of enzymes in brain, especially human brain.

Gap renormalization of molecular crystals from density-functional theory

Fundamental gap renormalization due to electronic polarization is a basic phenomenon in molecular crystals. Despite its ubiquity and importance, all conventional approaches within density-functional theory completely fail to capture it, even qualitatively. Here, we present a new screened range-separated hybrid functional, which, through judicious introduction of the scalar dielectric constant, quantitatively captures polarization-induced gap renormalization, as demonstrated on the prototypical organic molecular crystals of benzene, pentacene, and C-60. This functional is predictive, as it contains system-specific adjustable parameters that are determined from first principles, rather than from empirical considerations.

Single-step synthesis of internally functionalizable hyperbranched polyethers

Radical catalyzed thiol-ene reaction has become a useful alternative to the Huisgen-type azide-yne click reaction as it helps expand the variability in reaction conditions as well as the range of clickable entities. In this study, the direct generation of a hyperbranched polyether (HBPE) having decyl units at the periphery and a pendant allyl group on every repeat unit of the polymer backbone is described; the allyl groups serve as a reactive handle for postpolymerization modifications and permits the generation of a variety of internally functionalized HBPEs. In this design, the AB(2) monomer carries two decylbenzyl ether units (B-functionality), an aliphatic OH (A-functionality) and a pendant allyl group within the spacer segment; polymerization of the monomer readily occurs at 150 degrees C via melt transetherification process by continuous removal of 1-decanol under reduced pressure. The resulting HBPE has a hydrophobic periphery due to the presence of numerous decyl chains, while the allyl groups that remain unaffected during the melt polymerization provides an opportunity to install a variety of functional groups within the interior; thiol-ene click reaction with two different thiols, namely 3-mercaptopropionic acid and mercaptosuccinic acid, generated interesting amphiphilic structures. Preliminary field emission scanning electron microscope (FESEM) and Atomic Force Microscopy (AFM) imaging studies reveal the formation of fairly uniform spherical aggregates in water with sizes ranging from 200 to 400 nm; this suggests that these amphiphilic HBPs is able to reconfigure to generate jellyfish-like conformations that subsequently aggregate in an alkaline medium. The internal allyl functional groups were also used to generate intramolecularly core-crosslinked HBPEs, by the use of dithiol crosslinkers; gel permeation chromatography traces provided clear evidence for reduction in the size after crosslinking. In summary, we have developed a simple route to prepare core-clickable HBPEs and have demonstrated the quantitative reaction of the allyl groups present within the interior of the polymers; such HB polymeric systems that carry numerous functional groups within the core could have interesting applications in analyte sequestration and possibly sensing, especially from organic media. (c) 2013 Wiley Periodicals, Inc. J. Polym. Sci., Part A: Polym. Chem. 2013, 51, 4125-4135

Energy Hyperspace for Stacking Interaction in AU/AU Dinucleotide Step: Dispersion-Corrected Density Functional Theory Study

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence-directed features, such as small roll values for the purine-pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near-zero slide to be most favorable for the purine-pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar-phosphate backbone with C3-endo sugars and this demands C1-C1 distance of about 5.4 angstrom along the chains. Consideration of an energy penalty term for deviation of C1-C1 distance from the mean value, to the recent DFT-D functionals, specifically B97X-D appears to predict reliable energy contour for AU/AU step. Such distance-based penalty improves energy contours for the other purine-pyrimidine sequences also. (c) 2013 Wiley Periodicals, Inc. Biopolymers 101: 107-120, 2014.

Functional specificity among yeast mitochondrial Hsp70s paralogs and orthologs

The evolutionary diversity of the HSP70 gene family at the genetic level has generated complex structural variations leading to altered functional specificity and mode of regulation in different cellular compartments. By utilizing Saccharomyces cerevisiae as a model system for better understanding the global functional cooperativity between Hsp70 paralogs, we have dissected the differences in functional properties at the biochemical level between mitochondrial heat shock protein 70 (mtHsp70) Ssc1 and an uncharacterized Ssc3 paralog. Based on the evolutionary origin of Ssc3 and a high degree of sequence homology with Ssc1, it has been proposed that both have a close functional overlap in the mitochondrial matrix. Surprisingly, our results demonstrate that there is no functional cross-talk between Ssc1 and Ssc3 paralogs. The lack of in vivo functional overlap is due to altered conformation and significant lower stability associated with Ssc3. The substrate-binding domain of Ssc3 showed poor affinity toward mitochondrial client proteins and Tim44 due to the open conformation in ADP-bound state. In addition to that, the nucleotide-binding domain of Ssc3 showed an altered regulation by the Mge1 co-chaperone due to a high degree of conformational plasticity, which strongly promotes aggregation. Besides, Ssc3 possesses a dysfunctional inter-domain interface thus rendering it unable to perform functions similar to generic Hsp70s. Moreover, we have identified the critical amino acid sequence of Ssc1 and Ssc3 that can “make or break” mtHsp70 chaperone function. Together, our analysis provides the first evidence to show that the nucleotide-binding domain of mtHsp70s plays a critical role in determining the functional specificity among paralogs and orthologs across kingdoms.

Molecular flexibility of Mycobacterium tuberculosis ribosome recycling factor and its functional consequences: An exploration involving mutants

Internal mobility of the two domain molecule of ribosome recycling factor (RRF) is known to be important for its action. Mycobacterium tuberculosis RRF does not complement E. coli for its deficiency of RRF (in the presence of E. coli EF-G alone). Crystal structure had revealed higher rigidity of the M. tuberculosis RRF due to the presence of additional salt bridges between domains. Two inter-domain salt bridges and one between the linker region and the domain containing C-terminal residues were disrupted by appropriate mutations. Except for a C-terminal deletion mutant, all mutants showed RRF activity in E. coli when M. tuberculosis EF-G was also co-expressed. The crystal structures of the point mutants, that of the C-terminal deletion mutant and that of the protein grown in the presence of a detergent, were determined. The increased mobility resulting from the disruption of the salt bridge involving the hinge region allows the appropriate mutant to weakly complement E. coli for its deficiency of RRF even in the absence of simultaneous expression of the mycobacterial EF-G. The loss of activity of the C-terminal deletion mutant appears to be partly due to the rigidification of the molecule consequent to changes in the hinge region.

Functional analysis of ultra high information rates conveyed by rat vibrissal primary afferents

Sensory receptors determine the type and the quantity of information available for perception. Here, we quantified and characterized the information transferred by primary afferents in the rat whisker system using neural system identification. Quantification of ``how much'' information is conveyed by primary afferents, using the direct method (DM), a classical information theoretic tool, revealed that primary afferents transfer huge amounts of information (up to 529 bits/s). Information theoretic analysis of instantaneous spike-triggered kinematic stimulus features was used to gain functional insight on ``what'' is coded by primary afferents. Amongst the kinematic variables tested-position, velocity, and acceleration-primary afferent spikes encoded velocity best. The other two variables contributed to information transfer, but only if combined with velocity. We further revealed three additional characteristics that play a role in information transfer by primary afferents. Firstly, primary afferent spikes show preference for well separated multiple stimuli (i.e., well separated sets of combinations of the three instantaneous kinematic variables). Secondly, neurons are sensitive to short strips of the stimulus trajectory (up to 10 ms pre-spike time), and thirdly, they show spike patterns (precise doublet and triplet spiking). In order to deal with these complexities, we used a flexible probabilistic neuron model fitting mixtures of Gaussians to the spike triggered stimulus distributions, which quantitatively captured the contribution of the mentioned features and allowed us to achieve a full functional analysis of the total information rate indicated by the DM. We found that instantaneous position, velocity, and acceleration explained about 50% of the total information rate. Adding a 10 ms pre-spike interval of stimulus trajectory achieved 80-90%. The final 10-20% were found to be due to non-linear coding by spike bursts.