Molecular Modeling Suggests Cruzain Specificity for Peptide Primaquine Prodrugs

Chagas` disease, infection caused by the protozoan Trypanosoma cruzi, is an important, social and medical ailment in the Latin America. This disease is endemic in 21 countries, mostly Latin America countries, with more than 300,000 new cases every year and about 16-18 million infected people. Current therapy is not effective in the chronic phase of the disease. Thus, new and better drugs are urgently needed. In this sense, the in vitro activity of primaquine (PQ) was reported. Based on this, peptide prodrugs of primaquine containing dipeptides - lysine-arginine (LysArg), phenylalanine-alanine (PheAla) and phenylalanine-arginine (PheArg) -- as carriers, were designed to be selectively cleaved by cruzain, a specific cysteine protease of T. cruzi. The prodrugs have shown to be active against tripomastigote forms according to this order: LysArg-PQ> PheAla-PQ> PheArg-PQ. The molecular mechanism of action considered a probable nucleophilic attack of the catalytic residue of cruzain (Cys25) on the respective prodrug amide carbonyl carbon, releasing PQ. In order to test this hypothesis, molecular modeling studies were performed, physicochemical parameters and stereoelectronic features calculated by using the AM1 semi-empirical method suggest that the amide carbonyl carbon is favorable for cleavage, where the LysArg showed the most electronic reactive and sterically disposable, leading to the prodrug release and action. In addition, the docking study indicates the occurrence of specific interactions between prodrugs and the pockets S1 and S2 of cruzain through the dipeptides carriers, being the distance between cruzain Cys25 and the amide carbonyl group related to the biological activity of the prodrugs.

QSAR Modeling of a Set of Pyrazinoate Esters as Antituberculosis Prodrugs

Tuberculosis is an infection caused mainly by Mycobacterium tuberculosis. A first-line antimycobacterial drug is pyrazinamide (PZA), which acts partially as a prodrug activated by a pyrazinamidase releasing the active agent, pyrazinoic acid (POA). As pyrazinoic acid presents some difficulty to cross the mycobacterial cell wall, and also the pyrazinamide-resistant strains do not express the pyrazinamidase, a set of pyrazinoic acid esters have been evaluated as antimycobacterial agents. In this work, a QSAR approach was applied to a set of forty-three pyrazinoates against M. tuberculosis ATCC 27294, using genetic algorithm function and partial least squares regression (WOLF 5.5 program). The independent variables selected were the Balaban index (I), calculated n-octanol/water partition coefficient (ClogP), van-der-Waals surface area, dipole moment, and stretching-energy contribution. The final QSAR model (N = 32, r(2) = 0.68, q(2) = 0.59, LOF = 0.25, and LSE = 0.19) was fully validated employing leave-N-out cross-validation and y-scrambling techniques. The test set (N = 11) presented an external prediction power of 73%. In conclusion, the QSAR model generated can be used as a valuable tool to optimize the activity of future pyrazinoic acid esters in the designing of new antituberculosis agents.

Molecular modeling as a promising tool to study dendrimer prodrugs delivery

Molecular modeling methodologies were applied to perform preliminary studies concerning the release of active agents from potentially antichagasic and antileishmanial dendrimer prodrugs. The dendrimer was designed having myo-inositol as a core, L-malic acid as a spacer group, and hydroxymethylnitrofurazone (NFOH), 3-hydroxyflavone or quercetin, as active compounds. Each dendrimer presented a particular behavior concerning to the following investigated properties: spatial hindrance, map of electrostatic potential (MEP), and the lowest unoccupied molecular orbital energy (E(LUMO)). Additionally, the findings suggested that the carbonyl group next to the active agent seems to be the most promising ester breaking point. (C) 2009 Elsevier B.V. All rights reserved.

Molecular modeling and QSAR studies of a set of indole and benzimidazole derivatives as H(4) receptor antagonists

Histamine is an important biogenic amine, which acts with a group of four G-protein coupled receptors (GPCRs), namely H(1) to H(4) (H(1)R - H(4)R) receptors. The actions of histamine at H(4)R are related to immunological and inflammatory processes, particularly in pathophysiology of asthma, and H(4)R ligands having antagonistic properties could be helpful as antiinflammatory agents. In this work, molecular modeling and QSAR studies of a set of 30 compounds, indole and benzimidazole derivatives, as H(4)R antagonists were performed. The QSAR models were built and optimized using a genetic algorithm function and partial least squares regression (WOLF 5.5 program). The best QSAR model constructed with training set (N = 25) presented the following statistical measures: r (2) = 0.76, q (2) = 0.62, LOF = 0.15, and LSE = 0.07, and was validated using the LNO and y-randomization techniques. Four of five compounds of test set were well predicted by the selected QSAR model, which presented an external prediction power of 80%. These findings can be quite useful to aid the designing of new anti-H(4) compounds with improved biological response.

Host-guest complexation of a nitroheterocyclic compound with cyclodextrins: a spectrofluorimetric and molecular modeling study

Nitroheterocyclic compounds (NC) were candidate drugs proposed for Chagas disease chemotherapy. In this study, we investigated the complexation of hydroxymethylnitrofurazone (NFOH), a potential antichagasic compound, with alpha-cyclodextrin (alpha-CD), beta-cyclodextrin (beta-CD), Hydroxypropyl-beta-cyclodextrin (HP-beta-CD), Dimethyl-beta-cyclodextrin (DM-beta-CD) and gamma-cyclodextrin (gamma-CD) by fluorescence spectroscopy and molecular modeling studies. Hildebrand-Benesi equation was used to calculate the formation constants of NFOH with cyclodextrins based on the fluorescence differences in the CDs solution. The complexing capacity of NFOH with different CDs was compared through the results of association constant according to the following order: DM-beta-CD > beta-CD > alpha-CD > HP-beta-CD > gamma-CD. Molecular modeling studies give support for the experimental assignments, in favor of the formation of an inclusion complex between cyclodextrins with NFOH. This is an important study to investigate the effects of different kinds of cyclodextrins on the inclusion complex formation with NFOH and to better characterize a potential formulations to be used as therapeutic options for the oral treatment of Chagas disease.

Design, synthesis, antimicrobial activity and molecular modeling studies of novel benzofuroxan derivatives against Staphylococcus aureus

Molecular modi. cation is a quite promising strategy in the design and development of drug analogs with better bioavailability, higher intrinsic activity and less toxicity. In the search of new leads with potential antimicrobial activity, a new series of 14 4-substituted [N`-(benzofuroxan-5-yl) methylene] benzohydrazides, nifuroxazide derivatives, were synthesized and tested against standard and multidrug-resistant Staphylococcus aureus strains. The selection of the substituent groups was based on physicochemical properties, such as hydrophobicity and electronic effect. These properties were also evaluated through the lipophilic and electrostatic potential maps, respectively, considering the compounds with better biological pro. le. Twelve compounds exhibited similar bacteriostatic activity against standard and multidrug-resistant strains. The most active compound was the 4-CF(3) substituted derivative, which presented a minimum inhibitory concentration (MIC) value of 14.6-13.1 mu g/mL, and a ClogP value of 1.87. The results highlight the benzofuroxan derivatives as potential leads for designing new future antimicrobial drug candidates. (C) 2009 Elsevier Ltd. All rights reserved.

Molecular modeling studies and in vitro bioactivity evaluation of a set of novel 5-nitro-heterocyclic derivatives as anti-T. cruzi agents

In this study, in vitro anti-T. cruzi activity assays of nifuroxazide (NX) analogues, such as 5-nitro-2-furfuryliden and 5-nitro-2-theniliden derivatives, were performed. A molecular modeling approach was also carried out to relate the lipophilicity potential ( LP) property and biological activity data. The majority of the NX derivatives showed increased anti-T. cruzi activity in comparison to the reference drug, benznidazole (BZN). Additionally, the 5-nitro-2-furfuryliden derivatives presented better pharmacological profile than the 5-nitro-2-theniliden analogues. The LP maps and corresponding ClogP values indicate that there is an optimum lipophilicity value, which must be observed in the design of new potential anti-T. cruzi agents. (c) 2009 Elsevier Ltd. All rights reserved.

Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

Synthesis, antichagasic in vitro evaluation, cytotoxicity assays, molecular modeling and SAR/QSAR studies of a 2-phenyl-3-(1-phenyl-1H-pyrazol-4-yl)-acrylic acid benzylidene-carbohydrazide series

Chagas disease (American trypanosomiasis) is one of the most important parasitic diseases with serious social and economic impacts mainly on Latin America. This work reports the synthesis, in vitro trypanocidal evaluation, cytotoxicity assays, and molecular modeling and SAR/QSAR studies of a new series of N-phenylpyrazole benzylidene-carbohydrazides. The results pointed 6k (X = H, Y = p-NO(2), pIC(50) = 4.55 M) and 6l (X = F, Y = p-CN, pIC(50) = 4.27 M) as the most potent derivatives compared to crystal violet (pIC(50) = 3.77 M). The halogen-benzylidene-carbohydrazide presented the lowest potency whereas 6l showed the most promising pro. le with low toxicity (0% of cell death). The best equation from the 4D-QSAR analysis (Model 1) was able to explain 85% of the activity variability. The QSAR graphical representation revealed that bulky X-substituents decreased the potency whereas hydrophobic and hydrogen bond acceptor Y-substituents increased it. (C) 2008 Elsevier Ltd. All rights reserved.

Determination of coalescence kernels for high-shear granulation using DEM simulations

Discrete element method (DEM) modeling is used in parallel with a model for coalescence of deformable surface wet granules. This produces a method capable of predicting both collision rates and coalescence efficiencies for use in derivation of an overall coalescence kernel. These coalescence kernels can then be used in computationally efficient meso-scale models such as population balance equation (PBE) models. A soft-sphere DEM model using periodic boundary conditions and a unique boxing scheme was utilized to simulate particle flow inside a high-shear mixer. Analysis of the simulation results provided collision frequency, aggregation frequency, kinetic energy, coalescence efficiency and compaction rates for the granulation process. This information can be used to bridge the gap in multi-scale modeling of granulation processes between the micro-scale DEM/coalescence modeling approach and a meso-scale PBE modeling approach.

XSophe-Sophe-XeprView (R). A computer simulation software suite (v. 1.1.3) for the analysis of continuous wave EPR spectra

The XSophe-Sophe-XeprView((R)) computer simulation software suite enables scientists to easily determine spin Hamiltonian parameters from isotropic, randomly oriented and single crystal continuous wave electron paramagnetic resonance (CW EPR) spectra from radicals and isolated paramagnetic metal ion centers or clusters found in metalloproteins, chemical systems and materials science. XSophe provides an X-windows graphical user interface to the Sophe programme and allows: creation of multiple input files, local and remote execution of Sophe, the display of sophelog (output from Sophe) and input parameters/files. Sophe is a sophisticated computer simulation software programme employing a number of innovative technologies including; the Sydney OPera HousE (SOPHE) partition and interpolation schemes, a field segmentation algorithm, the mosaic misorientation linewidth model, parallelization and spectral optimisation. In conjunction with the SOPHE partition scheme and the field segmentation algorithm, the SOPHE interpolation scheme and the mosaic misorientation linewidth model greatly increase the speed of simulations for most spin systems. Employing brute force matrix diagonalization in the simulation of an EPR spectrum from a high spin Cr(III) complex with the spin Hamiltonian parameters g(e) = 2.00, D = 0.10 cm(-1), E/D = 0.25, A(x) = 120.0, A(y) = 120.0, A(z) = 240.0 x 10(-4) cm(-1) requires a SOPHE grid size of N = 400 (to produce a good signal to noise ratio) and takes 229.47 s. In contrast the use of either the SOPHE interpolation scheme or the mosaic misorientation linewidth model requires a SOPHE grid size of only N = 18 and takes 44.08 and 0.79 s, respectively. Results from Sophe are transferred via the Common Object Request Broker Architecture (CORBA) to XSophe and subsequently to XeprView((R)) where the simulated CW EPR spectra (1D and 2D) can be compared to the experimental spectra. Energy level diagrams, transition roadmaps and transition surfaces aid the interpretation of complicated randomly oriented CW EPR spectra and can be viewed with a web browser and an OpenInventor scene graph viewer.

Modeling control objectives for business process compliance

Business process design is primarily driven by process improvement objectives. However, the role of control objectives stemming from regulations and standards is becoming increasingly important for businesses in light of recent events that led to some of the largest scandals in corporate history. As organizations strive to meet compliance agendas, there is an evident need to provide systematic approaches that assist in the understanding of the interplay between (often conflicting) business and control objectives during business process design. In this paper, our objective is twofold. We will firstly present a research agenda in the space of business process compliance, identifying major technical and organizational challenges. We then tackle a part of the overall problem space, which deals with the effective modeling of control objectives and subsequently their propagation onto business process models. Control objective modeling is proposed through a specialized modal logic based on normative systems theory, and the visualization of control objectives on business process models is achieved procedurally. The proposed approach is demonstrated in the context of a purchase-to-pay scenario.

User guide for shock and blast simulation with the OctVCE code (version 3.5+), Department of Mechanical Engineering Report 2007/13

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel as it is simple to code and sufficient for practical engineering design problems. This also makes the code much more ‘user-friendly’ than structured grid approaches as the gridding process is done automatically. The CFD methodology relies on a finite-volume formulation of the unsteady Euler equations and is solved using a standard explicit Godonov (MUSCL) scheme. Both octree-based adaptive mesh refinement and shared-memory parallel processing capability have also been incorporated. For further details on the theory behind the code, see the companion report 2007/12.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.