11β-Hydroxysteroid dehydrogenase inhibitors and selectivity modeling of 11β-HSDs

11β-hydroksisteroididehydrogenaasientsyymit (11β-HSD) 1 ja 2 säätelevät kortisonin ja kortisolin määrää kudoksissa. 11β-HSD1 -entsyymin ylimäärä erityisesti viskeraalisessa rasvakudoksessa aiheuttaa metaboliseen oireyhtymän klassisia oireita, mikä tarjoaa mahdollisuuden metabolisen oireyhtymän hoitoon 11β-HSD1 -entsyymin selektiivisellä estämisellä. 11β-HSD2 -entsyymin inhibitio aiheuttaa kortisonivälitteisen mineralokortikoidireseptorien aktivoitumisen, mikä puolestaan johtaa hypertensiivisiin haittavaikutuksiin. Haittavaikutuksista huolimatta 11β-HSD2 -entsyymin estäminen saattaa olla hyödyllistä tilanteissa, joissa halutaan nostaa kortisolin määrä elimistössä. Lukuisia selektiivisiä 11β-HSD1 inhibiittoreita on kehitetty, mutta 11β-HSD2-inhibiittoreita on raportoitu vähemmän. Ero näiden kahden isotsyymin aktiivisen kohdan välillä on myös tuntematon, mikä vaikeuttaa selektiivisten inhibiittoreiden kehittämistä kummallekin entsyymille. Tällä työllä oli kaksi tarkoitusta: (1) löytää ero 11β-HSD entsyymien välillä ja (2) kehittää farmakoforimalli, jota voitaisiin käyttää selektiivisten 11β-HSD2 -inhibiittoreiden virtuaaliseulontaan. Ongelmaa lähestyttiin tietokoneavusteisesti: homologimallinnuksella, pienmolekyylien telakoinnilla proteiiniin, ligandipohjaisella farmakoforimallinnuksella ja virtuaaliseulonnalla. Homologimallinnukseen käytettiin SwissModeler -ohjelmaa, ja luotu malli oli hyvin päällekäinaseteltavissa niin templaattinsa (17β-HSD1) kuin 11β-HSD1 -entsyymin kanssa. Eroa entsyymien välillä ei löytynyt tarkastelemalla päällekäinaseteltuja entsyymejä. Seitsemän yhdistettä, joista kuusi on 11β-HSD2 -selektiivisiä, telakoitiin molempiin entsyymeihin käyttäen ohjelmaa GOLD. 11β-HSD1 -entsyymiin yhdisteet kiinnittyivät kuten suurin osa 11β-HSD1 -selektiivisistä tai epäselektiivisistä inhibiittoreista, kun taas 11β-HSD2 -entsyymiin kaikki yhdisteet olivat telakoituneet käänteisesti. Tällainen sitoutumistapa mahdollistaa vetysidokset Ser310:een ja Asn171:een, aminohappoihin, jotka olivat nähtävissä vain 11β-HSD2 -entsyymissä. Farmakoforimallinnukseen käytettiin ohjelmaa LigandScout3.0, jolla ajettiin myös virtuaaliseulonnat. Luodut kaksi farmakoforimallia, jotka perustuivat aiemmin telakointiinkin käytettyihin kuuteen 11β-HSD2 -selektiiviseen yhdisteeseen, koostuivat kuudesta ominaisuudesta (vetysidosakseptori, vetysidosdonori ja hydrofobinen), ja kieltoalueista. 11β-HSD2 -selektiivisyyden kannalta tärkeimmät ominaisuudet ovat vetysidosakseptori, joka voi muodostaa sidoksen Ser310 kanssa ja vetysidosdonori sen vieressä. Tälle vetysidosdonorille ei löytynyt vuorovaikutusparia 11β-HSD2-mallista. Sopivasti proteiiniin orientoitunut vesimolekyyli voisi kuitenkin olla sopiva ratkaisu puuttuvalle vuorovaikutusparille. Koska molemmat farmakoforimallit löysivät 11β-HSD2 -selektiivisiä yhdisteitä ja jättivät epäselektiivisiä pois testiseulonnassa, käytettiin molempia malleja Innsbruckin yliopistossa säilytettävistä yhdisteistä (2700 kappaletta) koostetun tietokannan seulontaan. Molemmista seulonnoista löytyneistä hiteistä valittiin yhteensä kymmenen kappaletta, jotka lähetettiin biologisiin testeihin. Biologisien testien tulokset vahvistavat lopullisesti sen kuinka hyvin luodut mallit edustavat todellisuudessa 11β-HSD2 -selektiivisyyttä.

New computational approaches for wrinkled and slack membranes

The static response of thin, wrinkled membranes is studied using both a tension field approximation based on plane stress conditions and a 3D nonlinear elasticityformulation, discretized through 8-noded Cosserat point elements. While the tension field approach only obtains the wrinkled/slack regions and at best a measure of the extent of wrinkliness, the 3D elasticity solution provides, in principle, the deformed shape of a wrinkled/slack membrane. However, since membranes barely resist compression, the discretized and linearized system equations via both the approaches are ill-conditioned and solutions could thus be sensitive to discretizations errors as well as other sources of noises/imperfections. We propose a regularized, pseudo-dynamical recursion scheme that provides a sequence of updates, which are almost insensitive to theregularizing term as well as the time step size used for integrating the pseudo-dynamical form. This is borne out through several numerical examples wherein the relative performance of the proposed recursion scheme vis-a-vis a regularized Newton strategy is compared. The pseudo-time marching strategy, when implemented using 3D Cosserat point elements, also provides a computationally cheaper, numerically accurate and simpler alternative to that using geometrically exact shell theories for computing large deformations of membranes in the presence of wrinkles. (C) 2010 Elsevier Ltd. All rights reserved.

Experimental and Computational Characterization of Alloy--Spinel--Corundum Equilibrium in the System Fe--Co--Al--O at 1873K

The three-phase equilibrium between alloy, spinel solid solution and alpha -Al sub 2 O sub 3 in the Fe--Co--Al--O system at 1873k was fully characterized as a function of alloy composition using both experimental and computational methods. The equilibrium oxygen content of the liquid alloy was measured by suction sampling and inert gas fusion analysis. The O potential corresponding to the three-phase equilibrium was determined by emf measurements on a solid state galvanic cell incorporating (Y sub 2 O sub 3 )ThO sub 2 as the solid electrolyte and Cr + Cr sub 2 O sub 3 as the reference electrode. The equilibrium composition of the spinel phase formed at the interface between the alloy and alumina crucible was measured by electron probe microanalysis (EPMA). The experimental results were compared with the values computed using a thermodynamic model. The model used values for standard Gibbs energies of formation of pure end-member spinels and Gibbs energies of solution of gaseous O in liquid Fe and cobalt available in the literature. The activity--composition relationship in the spinel solid solution was computed using a cation distribution model. The variation of the activity coefficient of O with alloy composition in the Fe--Co--O system was estimated using both the quasichemical model of Jacob and Alcock and Wagner's model along with the correlations of Chiang and Chang and Kuo and Chang. The computed results of spinel composition and O potential are in excellent agreement with the experimental data. Graphs. 29 ref.--AA

sRNAscanner: A Computational Tool for Intergenic Small RNA Detection in Bacterial Genomes

Background:Bacterial non-coding small RNAs (sRNAs) have attracted considerable attention due to their ubiquitous nature and contribution to numerous cellular processes including survival, adaptation and pathogenesis. Existing computational approaches for identifying bacterial sRNAs demonstrate varying levels of success and there remains considerable room for improvement. Methodology/Principal Findings: Here we have proposed a transcriptional signal-based computational method to identify intergenic sRNA transcriptional units (TUs) in completely sequenced bacterial genomes. Our sRNAscanner tool uses position weight matrices derived from experimentally defined E. coli K-12 MG1655 sRNA promoter and rho-independent terminator signals to identify intergenic sRNA TUs through sliding window based genome scans. Analysis of genomes representative of twelve species suggested that sRNAscanner demonstrated equivalent sensitivity to sRNAPredict2, the best performing bioinformatics tool available presently. However, each algorithm yielded substantial numbers of known and uncharacterized hits that were unique to one or the other tool only. sRNAscanner identified 118 novel putative intergenic sRNA genes in Salmonella enterica Typhimurium LT2, none of which were flagged by sRNAPredict2. Candidate sRNA locations were compared with available deep sequencing libraries derived from Hfq-co-immunoprecipitated RNA purified from a second Typhimurium strain (Sittka et al. (2008) PLoS Genetics 4: e1000163). Sixteen potential novel sRNAs computationally predicted and detected in deep sequencing libraries were selected for experimental validation by Northern analysis using total RNA isolated from bacteria grown under eleven different growth conditions. RNA bands of expected sizes were detected in Northern blots for six of the examined candidates. Furthermore, the 5'-ends of these six Northern-supported sRNA candidates were successfully mapped using 5'-RACE analysis. Conclusions/Significance: We have developed, computationally examined and experimentally validated the sRNAscanner algorithm. Data derived from this study has successfully identified six novel S. Typhimurium sRNA genes. In addition, the computational specificity analysis we have undertaken suggests that similar to 40% of sRNAscanner hits with high cumulative sum of scores represent genuine, undiscovered sRNA genes. Collectively, these data strongly support the utility of sRNAscanner and offer a glimpse of its potential to reveal large numbers of sRNA genes that have to date defied identification. sRNAscanner is available from: http://bicmku.in:8081/sRNAscanner or http://cluster.physics.iisc.ernet.in/sRNAscanner/.

New computational approaches for inelastic x-ray scattering

Inelastic x-ray scattering spectroscopy is a versatile experimental technique for probing the electronic structure of materials. It provides a wealth of information on the sample's atomic-scale structure, but extracting this information from the experimental data can be challenging because there is no direct relation between the structure and the measured spectrum. Theoretical calculations can bridge this gap by explaining the structural origins of the spectral features. Reliable methods for modeling inelastic x-ray scattering require accurate electronic structure calculations. This work presents the development and implementation of new schemes for modeling the inelastic scattering of x-rays from non-periodic systems. The methods are based on density functional theory and are applicable for a wide variety of molecular materials. Applications are presented in this work for amorphous silicon monoxide and several gas phase systems. Valuable new information on their structure and properties could be extracted with the combination of experimental and computational methods.

Computational Methods for Studies of Quantum Dots and Boson Condensates

This thesis presents ab initio studies of two kinds of physical systems, quantum dots and bosons, using two program packages of which the bosonic one has mainly been developed by the author. The implemented models, \emph{i.e.}, configuration interaction (CI) and coupled cluster (CC) take the correlated motion of the particles into account, and provide a hierarchy of computational schemes, on top of which the exact solution, within the limit of the single-particle basis set, is obtained. The theory underlying the models is presented in some detail, in order to provide insight into the approximations made and the circumstances under which they hold. Some of the computational methods are also highlighted. In the final sections the results are summarized. The CI and CC calculations on multiexciton complexes in self-assembled semiconductor quantum dots are presented and compared, along with radiative and non-radiative transition rates. Full CI calculations on quantum rings and double quantum rings are also presented. In the latter case, experimental and theoretical results from the literature are re-examined and an alternative explanation for the reported photoluminescence spectra is found. The boson program is first applied on a fictitious model system consisting of bosonic electrons in a central Coulomb field for which CI at the singles and doubles level is found to account for almost all of the correlation energy. Finally, the boson program is employed to study Bose-Einstein condensates confined in different anisotropic trap potentials. The effects of the anisotropy on the relative correlation energy is examined, as well as the effect of varying the interaction potential.}

Computational Fluid Dynamics in Combustion

Computational fluid dynamics has reached a stage where flow field in practical situation can be predicted to aid the design and to probe into the fundamental flow physics to understand and resolve the issues in fundamental fluid mechanics The study examines the computation of reacting flows After exploring the conservation equations for species and energy, the methods of closing the reaction rate terms in turbulent flow have been examined briefly Two cases of computation where combustion-flow interaction plays important role, have been discussed to illustrate the computational aspects and the physical insight that can be gained by the reacting flow computation

A study of symbolic processing and computational aspects in helicopter dynamics

Even research models of helicopter dynamics often lead to a large number of equations of motion with periodic coefficients; and Floquet theory is a widely used mathematical tool for dynamic analysis. Presently, three approaches are used in generating the equations of motion. These are (1) general-purpose symbolic processors such as REDUCE and MACSYMA, (2) a special-purpose symbolic processor, DEHIM (Dynamic Equations for Helicopter Interpretive Models), and (3) completely numerical approaches. In this paper, comparative aspects of the first two purely algebraic approaches are studied by applying REDUCE and DEHIM to the same set of problems. These problems range from a linear model with one degree of freedom to a mildly non-linear multi-bladed rotor model with several degrees of freedom. Further, computational issues in applying Floquet theory are also studied, which refer to (1) the equilibrium solution for periodic forced response together with the transition matrix for perturbations about that response and (2) a small number of eigenvalues and eigenvectors of the unsymmetric transition matrix. The study showed the following: (1) compared to REDUCE, DEHIM is far more portable and economical, but it is also less user-friendly, particularly during learning phases; (2) the problems of finding the periodic response and eigenvalues are well conditioned.

A simple computational model for predicting .pi.-facial selectivity in reductions of sterically unbiased ketones. Relative importance of electrostatic and orbital interactions

Various factore controlling the preferred facial selectivity in the reductions of a number of sterically unbiased ketones have been evaluated using a semiempirical MO procedure. MNDO optimized geometries do not reveal any significant ground-state distortions which can be correlated with the observed face selectivities. Electrostatic effecta due to an approaching reagent were modeled by placing a test negative charge at a fixed distance from the carbonyl carbon on each of the two faces. A second series of calculations was carried out using the hydride ion as a test nucleophile. The latter calculations effectively include orbital interactions involving the u and u* orbitals of the newly formed bond in the reaction. The computed energy differences with the charge model are generally much larger compared to those with the hydride ion. However, both models lead to predictions which are qualitatively consistent with the experimentally determined facial preferences for most of the systems. Thus, electrostatic interactions between the nucleophile and the substrate seem to effectively determine the face selectivities in these molecules. However, there are a few exceptions in which orbital interactions are found to contribute significantly and occasionally reverse the preference dictated by electrostatic effecta. The remarkable succew of the hydride model calculations, in spite of retaining the unperturbed geometries of the substrates, points to the unimportance of torsional effeds and orbital distortions associated with the pyramidalized carbonyl unit in the transition state in most of the substrates considered. Additional experimental results are reported which provide useful calibration for the present computational approach.

Selectivity in Garratt-Braverman Cyclization: An Experimental and Computational Study

Bispropargyl sulfones equipped with aromatic rings of dissimilar nature were synthesized. Under basic conditions, these sulfones isomerized to the bisallenic sulfones, creating a competitive scenario between two alternate Garratt-Braverman (GB) cyclization pathways. The observed product distribution ruled out the involvement of any ionic intermediate and supported the diradical mechanism with greater involvement of the electron-rich aromatic ring via the more nucleophilic radical. DFT-based calculations supported the diradical mechanism along with the observed selectivity.

Theoretical and Computational Analysis of Static and Dynamic Anomalies in Water-DMSO Binary Mixture at Low DMSO Concentrations

Experiments have repeatedly observed both thermodynamic and dynamic anomalies in aqueous binary mixtures, surprisingly at low solute concentration. Examples of such binary mixtures include water-DMSO, water-ethanol, water-tertiary butyl alcohol (TBA), and water-dioxane, to name a few. The anomalies have often been attributed to the onset of a structural transition, whose nature, however, has been left rather unclear. Here we study the origin of such anomalies using large scale computer simulations and theoretical analysis in water-DMSO binary mixture. At very low DMSO concentration (below 10%), small aggregates of DMSO are solvated by water through the formation of DMSO-(H2O)(2) moieties. As the concentration is increased beyond 10-12% of DMSO, spanning clusters comprising the same moieties appear in the system. Those clusters are formed and stabilized not only through H-bonding but also through the association of CH3 groups of DMSO. We attribute the experimentally observed anomalies to a continuum percolation-like transition at DMSO concentration X-DMSO approximate to 12-15%. The largest cluster size of CH3-CH3 aggregation clearly indicates the formation of such percolating clusters. As a result, a significant slowing down is observed in the decay of associated rotational auto time correlation functions (of the S = O bond vector of DMSO and O-H bond vector of water). Markedly unusual behavior in the mean square fluctuation of total dipole moment again suggests a structural transition around the same concentration range. Furthermore, we map our findings to an interacting lattice model which substantiates the continuum percolation model as the reason for low concentration anomalies in binary mixtures where the solutes involved have both hydrophilic and hydrophobic moieties.

Recognition of Polycyclic Aromatic Hydrocarbons and Their Derivatives by the 1,3-Dinaphthalimide Conjugate of Calix4]arene: Emission, Absorption, Crystal Structures, and Computational Studies

Polycyclic aromatic hydrocarbons (PAHs) are environmental pollutants as well as well-known carcinogens. Therefore, it is important to develop an effective receptor for the detection and quantification of such molecules in solution. In view of this, a 1,3-dinaphthalimide derivative of calix4]arene (L) has been synthesized and characterized, and the structure has been established by single crystal XRD. In the crystal lattice, intermolecular arm-to-arm pi center dot center dot center dot pi overlap dominates and thus L becomes a promising receptor for providing interactions with the aromatic species in solution, which can be monitored by following the changes that occur in its fluorescence and absorption spectra. On the basis of the solution studies carried out with about 17 derivatives of the aromatic guest molecular systems, it may be concluded that the changes that occur in the fluorescence intensity seem to be proportional to the number of aromatic rings present and thus proportional to the extent of pi center dot center dot center dot pi interaction present between the naphthalimide moieties and the aromatic portion of the guest molecule. Though the nonaromatic portion of the guest species affects the fluorescence quenching, the trend is still based on the number of rings present in these. Four guest aldehydes are bound to L with K-ass of 2000-6000 M-1 and their minimum detection limit is in the range of 8-35 mu M. The crystal structure of a naphthaldehyde complex, L.2b, exhibits intermolecular arm-to-arm as well as arm-to-naphthaldehyde pi center dot center dot center dot pi interactions. Molecular dynamics studies of L carried out in the presence of aromatic aldehydes under vacuum as well as in acetonitrile resulted in exhibiting interactions observed in the solid state and hence the changes observed in the fluorescence and absorption spectra are attributable for such interactions. Complex formation has also been delineated through ESI MS studies. Thus L is a promising receptor that can recognize PAHs by providing spectral changes proportional to the aromatic conjugation of the guest and the extent of aromatic pi center dot center dot center dot pi interactions present between L and the guest.

Diastereospecific Coupling of Imines by Low-Valent Titanium: An Experimental and Computational Study

The reaction of phenylsilane (PhSiH3) and titanium(IV) isopropoxide Ti(OiPr)(4)] generates low-valent titanium alkoxides that reduce and reductively couple imines. The C-C coupling reaction is diastereospecific, with exclusive formation of the (+/-)-isomer. The yield is dependent on the concentration of titanium(IV) isopropoxide used. Using imines with varying electronic demand revealed that the coupling is most efficient for unsubstituted imines. The involvement of a trimeric titanium biradical species is invoked to explain the observed concentration dependence and diastereospecificity in the reaction. Multilayer ONIOM (''our own n-layered integrated molecular orbital and molecular mechanics'') calculations were carried out on the plausible intermediates involved by using the Guassian 09 suite of programs. Formation of the trimeric biradical species leading to the formation of the (+/-)-isomer is more favored than formation of the trimeric biradical species, which leads to the meso-isomer. Similar calculations with dimeric intermediates do not predict the (+/-)-selectivity. Mass spectral analysis of the reaction mixture indicates the presence of a trimeric titanium alkoxide species.