A comparison of deterministic, reliability-based and risk-based structural optimization under uncertainty

In this paper, the effects of uncertainty and expected costs of failure on optimum structural design are investigated, by comparing three distinct formulations of structural optimization problems. Deterministic Design Optimization (DDO) allows one the find the shape or configuration of a structure that is optimum in terms of mechanics, but the formulation grossly neglects parameter uncertainty and its effects on structural safety. Reliability-based Design Optimization (RBDO) has emerged as an alternative to properly model the safety-under-uncertainty part of the problem. With RBDO, one can ensure that a minimum (and measurable) level of safety is achieved by the optimum structure. However, results are dependent on the failure probabilities used as constraints in the analysis. Risk optimization (RO) increases the scope of the problem by addressing the compromising goals of economy and safety. This is accomplished by quantifying the monetary consequences of failure, as well as the costs associated with construction, operation and maintenance. RO yields the optimum topology and the optimum point of balance between economy and safety. Results are compared for some example problems. The broader RO solution is found first, and optimum results are used as constraints in DDO and RBDO. Results show that even when optimum safety coefficients are used as constraints in DDO, the formulation leads to configurations which respect these design constraints, reduce manufacturing costs but increase total expected costs (including expected costs of failure). When (optimum) system failure probability is used as a constraint in RBDO, this solution also reduces manufacturing costs but by increasing total expected costs. This happens when the costs associated with different failure modes are distinct. Hence, a general equivalence between the formulations cannot be established. Optimum structural design considering expected costs of failure cannot be controlled solely by safety factors nor by failure probability constraints, but will depend on actual structural configuration. (c) 2011 Elsevier Ltd. All rights reserved.

Optimization of anti-Trypanosoma cruzi oxadiazoles leads to identification of compounds with efficacy in infected mice

We recently showed that oxadiazoles have anti-Trypanosoma cruzi activity at micromolar concentrations. These compounds are easy to synthesize and show a number of clear and interpretable structure-activity relationships (SAR), features that make them attractive to pursue potency enhancement. We present here the structural design, synthesis, and anti-T. cruzi evaluation of new oxadiazoles denoted 5a-h and 6a-h. The design of these compounds was based on a previous model of computational docking of oxadiazoles on the T. cruzi protease cruzain. We tested the ability of these compounds to inhibit catalytic activity of cruzain, but we found no correlation between the enzyme inhibition and the antiparasitic activity of the compounds. However, we found reliable SAR data when we tested these compounds against the whole parasite. While none of these oxadiazoles showed toxicity for mammalian cells, oxadiazoles 6c (fluorine), 6d (chlorine), and 6e (bromine) reduced epimastigote proliferation and were cidal for trypomastigotes of T. cruzi Y strain. Oxadiazoles 6c and 6d have IC50 of 9.5 +/- 2.8 and 3.5 +/- 1.8 mu M for trypomastigotes, while Benznidazole, which is the currently used drug for Chagas disease treatment, showed an IC50 of 11.3 +/- 2.8 mu M. Compounds 6c and 6d impair trypomastigote development and invasion in macrophages, and also induce ultrastructural alterations in trypomastigotes. Finally, compound 6d given orally at 50 mg/kg substantially reduces the parasitemia in T. cruzi-infected BALB/c mice. Our drug design resulted in potency enhancement of oxadiazoles as anti-Chagas disease agents, and culminated with the identification of oxadiazole 6d, a trypanosomicidal compound in an animal model of infection. (C) 2012 Elsevier Ltd. All rights reserved.

Structural role of the active-site metal in the conformation of Trypanosoma brucei phosphoglycerate mutase

Phosphoglycerate mutases (PGAMs) participate in both the glycolytic and the gluconeogenic pathways in reversible isomerization of 3-phosphoglycerate and 2-phosphoglycerate. PGAMs are members of two distinct protein families: enzymes that are dependent on or independent of the 2,3-bisphosphoglycerate cofactor. We determined the X-ray structure of the monomeric Trypanosoma brucei independent PGAM (TbiPGAM) in its apoenzyme form, and confirmed this observation by small angle X-ray scattering data. Comparing the TbiPGAM structure with the Leishmania mexicana independent PGAM structure, previously reported with a phosphoglycerate molecule bound to the active site, revealed the domain movement resulting from active site occupation. The structure reported here shows the interaction between Asp319 and the metal bound to the active site, and its contribution to the domain movement. Substitution of the metal-binding residue Asp319 by Ala resulted in complete loss of independent PGAM activity, and showed for the first time its involvement in the enzymes function. As TbiPGAM is an attractive molecular target for drug development, the apoenzyme conformation described here provides opportunities for its use in structure-based drug design approaches.

Extensive structural change of the envelope protein of dengue virus induced by a tuned ionic strength: conformational and energetic analyses

The Dengue has become a global public health threat, with over 100 million infections annually; to date there is no specific vaccine or any antiviral drug. The structures of the envelope (E) proteins of the four known serotype of the dengue virus (DENV) are already known, but there are insufficient molecular details of their structural behavior in solution in the distinct environmental conditions in which the DENVs are submitted, from the digestive tract of the mosquito up to its replication inside the host cell. Such detailed knowledge becomes important because of the multifunctional character of the E protein: it mediates the early events in cell entry, via receptor endocytosis and, as a class II protein, participates determinately in the process of membrane fusion. The proposed infection mechanism asserts that once in the endosome, at low pH, the E homodimers dissociate and insert into the endosomal lipid membrane, after an extensive conformational change, mainly on the relative arrangement of its three domains. In this work we employ all-atom explicit solvent Molecular Dynamics simulations to specify the thermodynamic conditions in that the E proteins are induced to experience extensive structural changes, such as during the process of reducing pH. We study the structural behavior of the E protein monomer at acid pH solution of distinct ionic strength. Extensive simulations are carried out with all the histidine residues in its full protonated form at four distinct ionic strengths. The results are analyzed in detail from structural and energetic perspectives, and the virtual protein movements are described by means of the principal component analyses. As the main result, we found that at acid pH and physiological ionic strength, the E protein suffers a major structural change; for lower or higher ionic strengths, the crystal structure is essentially maintained along of all extensive simulations. On the other hand, at basic pH, when all histidine residues are in the unprotonated form, the protein structure is very stable for ionic strengths ranging from 0 to 225 mM. Therefore, our findings support the hypothesis that the histidines constitute the hot points that induce configurational changes of E protein in acid pH, and give extra motivation to the development of new ideas for antivirus compound design.

Structure- and ligand-based drug design approaches for neglected tropical diseases

Drug discovery has moved toward more rational strategies based on our increasing understanding of the fundamental principles of protein-ligand interactions. Structure( SBDD) and ligand-based drug design (LBDD) approaches bring together the most powerful concepts in modern chemistry and biology, linking medicinal chemistry with structural biology. The definition and assessment of both chemical and biological space have revitalized the importance of exploring the intrinsic complementary nature of experimental and computational methods in drug design. Major challenges in this field include the identification of promising hits and the development of high-quality leads for further development into clinical candidates. It becomes particularly important in the case of neglected tropical diseases (NTDs) that affect disproportionately poor people living in rural and remote regions worldwide, and for which there is an insufficient number of new chemical entities being evaluated owing to the lack of innovation and R&D investment by the pharmaceutical industry. This perspective paper outlines the utility and applications of SBDD and LBDD approaches for the identification and design of new small-molecule agents for NTDs.

Structural Insights into Human Peroxisome Proliferator Activated Receptor Delta (PPAR-Delta) Selective Ligand Binding

Peroxisome proliferator activated receptors (PPARs delta, alpha and gamma) are closely related transcription factors that exert distinct effects on fatty acid and glucose metabolism, cardiac disease, inflammatory response and other processes. Several groups developed PPAR subtype specific modulators to trigger desirable effects of particular PPARs without harmful side effects associated with activation of other subtypes. Presently, however, many compounds that bind to one of the PPARs cross-react with others and rational strategies to obtain highly selective PPAR modulators are far from clear. GW0742 is a synthetic ligand that binds PPAR delta more than 300-fold more tightly than PPAR alpha or PPAR gamma but the structural basis of PPAR delta: GW0742 interactions and reasons for strong selectivity are not clear. Here we report the crystal structure of the PPAR delta:GW0742 complex. Comparisons of the PPAR delta:GW0742 complex with published structures of PPARs in complex with alpha and gamma selective agonists and pan agonists suggests that two residues (Val312 and Ile328) in the buried hormone binding pocket play special roles in PPAR delta selective binding and experimental and computational analysis of effects of mutations in these residues confirms this and suggests that bulky substituents that line the PPAR alpha and gamma ligand binding pockets as structural barriers for GW0742 binding. This analysis suggests general strategies for selective PPAR delta ligand design.

Structural Investigation of Anti-Trypanosoma cruzi 2-Iminothiazolidin-4-ones Allows the Identification of Agents with Efficacy in Infected Mice

We modified the thiazolidinic ring at positions N3, C4, and C5, yielding compounds 6-24. Compounds with a phenyl at position N3, 15-19, 22-24, exhibited better inhibitory properties for cruzain and against the parasite than 2-iminothiazolidin-4-one S. We were able to identify one high-efficacy trypanocidal compound, 2-minothiazolidin-4-one 18, which inhibited the activity of cruzain and the proliferation of epirnastigotes and was cidal for trypomastigotes but was not toxic for splenocytes. Having located some of the structural determinants of the trypanocidal properties, we subsequently wished to determine if the exchange of the thiazolidine for a thiazole ring leaves the functional properties unaffected. We therefore tested thiazoles 26-45 and observed that they did not inhibit cruzain, but they exhibited trypanocidal effects. Parasite development was severely impaired when treated with 18, thus reinforcing the notion that this class of heterocycles can lead to useful cidal agents for Chagas disease.

Purification, characterization and structural determination of chitinases produced by Moniliophthora perniciosa

The enzyme chitinase from Moniliophthora perniciosa the causative agent of the witches' broom disease in Theobroma cacao, was partially purified with ammonium sulfate and filtration by Sephacryl S-200 using sodium phosphate as an extraction buffer. Response surface methodology (RSM) was used to determine the optimum pH and temperature conditions. Four different isoenzymes were obtained: ChitMp I, ChitMp II, ChitMp III and ChitMp IV. ChitMp I had an optimum temperature at 44-73ºC and an optimum pH at 7.0-8.4. ChitMp II had an optimum temperature at 45-73ºC and an optimum pH at 7.0-8.4. ChitMp III had an optimum temperature at 54-67ºC and an optimum pH at 7.3-8.8. ChitMp IV had an optimum temperature at 60ºC and an optimum pH at 7.0. For the computational biology, the primary sequence was determined in silico from the database of the Genome/Proteome Project of M. perniciosa, yielding a sequence with 564 bp and 188 amino acids that was used for the three-dimensional design in a comparative modeling methodology. The generated models were submitted to validation using Procheck 3.0 and ANOLEA. The model proposed for the chitinase was subjected to a dynamic analysis over a 1 ns interval, resulting in a model with 91.7% of the residues occupying favorable places on the Ramachandran plot and an RMS of 2.68.

Morphogenetic and structural characteristics of guinea grass tillers at different ages under intermittent stocking

The objective of this research was to assess morphogenetic and structural characteristics of tillers of guinea grass cv. Tanzania at different ages. The pastures of guinea grass were managed in six pasture conditions related to the combination of three frequencies (90, 95, and 99% light interception) and two post-grazing heights (25 and 50 cm). In these six pastures conditions, three tiller ages were evaluated (young, mature, and old). The design was of completely randomized block with three replications. Young tillers exhibited higher leaf appearance rate and leaf elongation rate and, consequently, higher final leaf length and number of live leaves than mature and old tillers, regardless of the pasture condition. On pastures managed with 90 or 95% light interception associated with a post-grazing height of 25 cm, old tillers presented longer leaf lifespan than young and mature ones. There is a progressive reduction in the vigor of growth of pastures of guinea grass cv. Tanzania with advancing tiller age.

Reliability-based design optimization strategies based on FORM: a review

In deterministic optimization, the uncertainties of the structural system (i.e. dimension, model, material, loads, etc) are not explicitly taken into account. Hence, resulting optimal solutions may lead to reduced reliability levels. The objective of reliability based design optimization (RBDO) is to optimize structures guaranteeing that a minimum level of reliability, chosen a priori by the designer, is maintained. Since reliability analysis using the First Order Reliability Method (FORM) is an optimization procedure itself, RBDO (in its classical version) is a double-loop strategy: the reliability analysis (inner loop) and the structural optimization (outer loop). The coupling of these two loops leads to very high computational costs. To reduce the computational burden of RBDO based on FORM, several authors propose decoupling the structural optimization and the reliability analysis. These procedures may be divided in two groups: (i) serial single loop methods and (ii) unilevel methods. The basic idea of serial single loop methods is to decouple the two loops and solve them sequentially, until some convergence criterion is achieved. On the other hand, uni-level methods employ different strategies to obtain a single loop of optimization to solve the RBDO problem. This paper presents a review of such RBDO strategies. A comparison of the performance (computational cost) of the main strategies is presented for several variants of two benchmark problems from the literature and for a structure modeled using the finite element method.

Design methodology of piezoelectric energy-harvesting skin using topology optimization

Abstract This paper describes a design methodology for piezoelectric energy harvester s that thinly encapsulate the mechanical devices and expl oit resonances from higher- order vibrational modes. The direction of polarization determines the sign of the pi ezoelectric tensor to avoid cancellations of electric fields from opposite polarizations in the same circuit. The resultant modified equations of state are solved by finite element method (FEM). Com- bining this method with the solid isotropic material with penalization (SIMP) method for piezoelectric material, we have developed an optimization methodology that optimizes the piezoelectric material layout and polarization direc- tion. Updating the density function of the SIMP method is performed based on sensitivity analysis, the sequen- tial linear programming on the early stage of the opti- mization, and the phase field method on the latter stage

Reactive sputter magnetron reactor for preparation of thin films and simultaneous in-situ structural study by X-ray diffraction.

Reactive Sputter Magnetron (RSM) is a widely used technique to thin films growing of compounds both, in research laboratories and in industrial processes. The nature of the deposited compound will depend then on the nature of the magnetron target and the nature of the ions generated in the plasma. One important aspect of the problem is the knowledge of the evolution of the film during the process of growing itself. In this work, we present the design, construction of a chamber to be installed in the Huber goniometer in the XRD2 line of LNLS in Campinas, which allows in situ growing kinetic studies of thin films.

Structure-based drug design studies on a series of aldolase inhibitors

Human African trypanosomiasis, also known as sleeping sickness, is a major cause of death in Africa, and for which there are no safe and effective treatments available. The enzyme aldolase from Trypanosoma brucei is an attractive, validated target for drug development. A series of alkyl‑glycolamido and alkyl-monoglycolate derivatives was studied employing a combination of drug design approaches. Three-dimensional quantitative structure-activity relationships (3D QSAR) models were generated using the comparative molecular field analysis (CoMFA). Significant results were obtained for the best QSAR model (r2 = 0.95, non-cross-validated correlation coefficient, and q2 = 0.80, cross-validated correlation coefficient), indicating its predictive ability for untested compounds. The model was then used to predict values of the dependent variables (pKi) of an external test set,the predicted values were in good agreement with the experimental results. The integration of 3D QSAR, molecular docking and molecular dynamics simulations provided further insight into the structural basis for selective inhibition of the target enzyme.

Inhibitors of Trypanosoma brucei trypanothione reductase: comparative molecular field analysis modeling and structural basis for selective inhibition

Background: Sleeping sickness is a major cause of death in Africa. Since no secure treatment is available, the development of novel therapeutic agents is urgent. In this context, the enzyme trypanothione reductase (TR) is a prominent molecular target that has been investigated in drug design for sleeping sickness. Results: In this study, comparative molecular field analysis models were generated for a series of Trypanosoma brucei TR inhibitors. Statistically significant results were obtained and the models were applied to predict the activity of external test sets, with good correlation between predicted and experimental results. We have also investigated the structural requirements for the selective inhibition of the parasite's enzyme over the human glutathione reductase. Conclusion: The quantitative structure-activity relationship models provided valuable information regarding the essential molecular requirements for the inhibitory activity upon the target protein, providing important insights into the design of more potent and selective TR inhibitors.