Quantitative structure-activity studies on a series of migrastatin analogs as inhibitors of cancer cell metastasis

Migrastatin, a macrolide natural product, and its structurally related analogs are potent inhibitors of cancer cell metastasis, invasion and migration. In the present work, a specialized fragment-based method was employed to develop QSAR models for a series of migrastatin and isomigrastatin analogs. Significant correlation coefficients were obtained (best model, q(2) = 0.76 and r(2) = 0.91) indicating that the QSAR models possess high internal consistency. The best model was then used to predict the potency of an external test set, and the predicted values were in good agreement with the experimental results (R(2) (pred) = 0.85). The final model and the corresponding contribution maps, combined with molecular modeling studies, provided important insights into the key structural features for the anticancer activity of this family of synthetic compounds based on natural products.

Fragment-based and classical quantitative structure-activity relationships for a series of hydrazides as antituberculosis agents

Worldwide, tuberculosis (TB) is the leading cause of death among curable infectious diseases. Multidrug-resistant Mycobacterium tuberculosis is an emerging problem of great importance to public health, and there is an urgent need for new anti-TB drugs. In the present work, classical 2D quantitative structure-activity relationships (QSAR) and hologram QSAR (HQSAR) studies were performed on a training set of 91 isoniazid derivatives. Significant statistical models (classical QSAR, q(2) = 0.68 and r(2) = 0.72; HQSAR, q(2) = 0.63 and r(2) = 0.86) were obtained, indicating their consistency for untested compounds. The models were then used to evaluate an external test set containing 24 compounds which were not included in the training set, and the predicted values were in good agreement with the experimental results (HQSAR, r(pred)(2) = 0.87; classical QSAR, r(pred)(2) = 0.75).

Classical and fragment-based hologram structure-activity relationships for a series of analgesic cyclic imides

Cyclic imides have been widely employed in drug design research due to their multiple pharmacological and biological properties. In the present study, two-dimensional quantitative structure-activity relationship (2D QSAR) studies were conducted on a series of potent analgesic cyclic imides using both classical and hologram QSAR (HQSAR) methods, yielding significant statistical models (classical QSAR, q(2) = 0.80; HQSAR, q(2) = 0.84). The models were then used to evaluate an external data test, and the predicted values were in good agreement with the experimental results, indicating their consistency for untested compounds.

Structure-activity relationships for a class of selective inhibitors of the major cysteine protease from Trypanosoma cruzi

Chagas` disease is a parasitic infection widely distributed throughout Latin America, with devastating consequences in terms of human morbidity and mortality. Cruzain, the major cysteine protease from Trypanosoma cruzi, is an attractive target for antitrypanosomal chemotherapy. In the present work, classical two-dimensional quantitative structure-activity relationships (2D QSAR) and hologram QSAR (HQSAR) studies were performed on a training set of 45 thiosemicarbazone and semicarbazone derivatives as inhibitors of T. cruzi cruzain. Significant statistical models (HQSAR, q2=0.75 and r2=0.96; classical QSAR, q2=0.72 and r2=0.83) were obtained, indicating their consistency for untested compounds. The models were then used to evaluate an external test set containing 10 compounds which were not included in the training set, and the predicted values were in good agreement with the experimental results (HQSAR, [image omitted]=0.95; classical QSAR, [image omitted]=0.91), indicating the existence of complementary between the two ligand-based drug design techniques.

Crystal structure and statistical coupling analysis of highly glycosylated peroxidase from royal palm tree (Roystonea regia)

Royal palm tree peroxidase (RPTP) is a very stable enzyme in regards to acidity, temperature, H(2)O(2), and organic solvents. Thus, RPTP is a promising candidate for developing H(2)O(2)-sensitive biosensors for diverse applications in industry and analytical chemistry. RPTP belongs to the family of class III secretory plant peroxidases, which include horseradish peroxidase isozyme C, soybean and peanut peroxidases. Here we report the X-ray structure of native RPTP isolated from royal palm tree (Roystonea regia) refined to a resolution of 1.85 angstrom. RPTP has the same overall folding pattern of the plant peroxidase superfamily, and it contains one heme group and two calcium-binding sites in similar locations. The three-dimensional structure of RPTP was solved for a hydroperoxide complex state, and it revealed a bound 2-(N-morpholino) ethanesulfonic acid molecule (MES) positioned at a putative substrate-binding secondary site. Nine N-glycosylation sites are clearly defined in the RPTP electron-density maps, revealing for the first time conformations of the glycan chains of this highly glycosylated enzyme. Furthermore, statistical coupling analysis (SCA) of the plant peroxidase superfamily was performed. This sequence-based method identified a set of evolutionarily conserved sites that mapped to regions surrounding the heme prosthetic group. The SCA matrix also predicted a set of energetically coupled residues that are involved in the maintenance of the structural folding of plant peroxidases. The combination of crystallographic data and SCA analysis provides information about the key structural elements that could contribute to explaining the unique stability of RPTP. (C) 2009 Elsevier Inc. All rights reserved.

Structure and thermal reactivity of Zn(II) salts of isocinchomeronic acid (2,5-pyridinedicarboxylic acid)

The synthesis, an improved refined crystal and molecular structure re-determination, and the thermal decomposition behavior of two Zn(II) derivatives of isocinchomeronic acid (2,5-pyridinedicarboxylic acid or H(2)2,5-pydc) are presented. [Zn(2,5-pydc)(H(2)O)(3)Zn(2,5-pydc)(H(2)O)(2)](2) (1) crystallizes in the triclinic P-1 space group with a = 7.106(2), b = 11.450(2), c = 11.869(1) angstrom, alpha = 107.29(1), beta = 104.08(1), gamma = 90.32(2)degrees, and Z = 2. [Zn(2,5-pydc)(H(2)O)(2)] center dot H(2)O (2) is orthorhombic (P2(1)2(1)2(1) space group), with a = 7.342(1), b = 9.430(1), c = 13.834(2) angstrom, and Z = 4. The structures were refined to agreement R(1)-factors of 0.0315 (1) and 0.0336 (2). Complex (1) is arranged as molecular Zn(4)(2,5-pydc)(4)(H(2)O)(10) tetramers, the cages of which define channels that remain unblocked by anions. Compound (2) is polymeric with Zn(2,5-pydc)(H(2)O)(2) and Zn(2,5-pydc)(H(2)O)(3) units linked through bridging ligands. Both compounds were synthesized under mild conditions in aqueous media, without need to resort to hydrothermal media. Changing the pH from 4.51 to 5.75 suffices to direct the chemical processes toward the orthorhombic compound rather than to the triclinic one.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

Structure and Calcium-Binding Activity of LipL32, the Major Surface Antigen of Pathogenic Leptospira sp.

Leptospixosis, a spirochaetal zoonotic disease caused by Leptospira, has been recognized as an important emerging infectious disease. LipL32 is the major exposed outer membrane protein found exclusively in pathogenic leptospires, where it accounts for up to 75% of the total outer membrane proteins. It is highly immunogenic, and recent studies have implicated LipL32 as an extracellular matrix binding protein, interacting with collagens, fibronectin, and laminin. In order to better understand the biological role and the structural requirements for the function of this important lipoprotein, we have determined the 2.25-angstrom-resolution structure of recombinant LipL32 protein corresponding to residues 21-272 of the wild-type protein (LipL32(21-272)). The LipL32(21-272) monomer is made of a jelly-roll fold core from which several peripheral secondary structures protrude. LipL32(21-272) is structurally similar to several other jelly-roll proteins, some of which bind calcium ions and extracellular matrix proteins. Indeed, spectroscopic data (circular dichroism, intrinsic tryptophan fluorescence, and extrinsic 1-amino-2-naphthol-4-sulfonic acid fluorescence) confirmed the calcium-binding properties of LipL32(21-272). Ca(2+) binding resulted in a significant increase in the thermal stability of the protein, and binding was specific for Ca(2+) as no structural or stability perturbations were observed for Mg(2+), Zn(2+), or Cu(2+). Careful examination of the crystal lographic structure suggests the locations of putative regions that could mediate Ca(2+) binding as well as binding to other interacting host proteins, such as collagens, fibronectin, and lamixidn. (C) 2009 Elsevier Ltd. All rights reserved.

The Aedes aegypti larval transcriptome: a comparative perspective with emphasis on trypsins and the domain structure of peritrophins

The genome sequence of Aedes aegypti was recently reported. A significant amount of Expressed Sequence Tags (ESTs) were sequenced to aid in the gene prediction process. In the present work we describe an integrated analysis of the genomic and EST data, focusing on genes with preferential expression in larvae (LG), adults (AG) and in both stages (SG). A total of 913 genes (5.4% of the transcript complement) are LG, including ion transporters and cuticle proteins that are important for ion homeostasis and defense. From a starting set of 245 genes encoding the trypsin domain, we identified 66 putative LG, AG, and SG trypsins by manual curation. Phylogenetic analyses showed that AG trypsins are divergent from their larval counterparts (LG), grouping with blood-induced trypsins from Anopheles gambiae and Simulium vittatum. These results support the hypothesis that blood-feeding arose only once, in the ancestral Culicomorpha. Peritrophins are proteins that interlock chitin fibrils to form the peritrophic membrane (PM) that compartmentalizes the food in the midgut. These proteins are recognized by having chitin-binding domains with 6 conserved Cys and may also present mucin-like domains (regions expected to be highly O-glycosylated). PM may be formed by a ring of cells (type 2, seen in Ae. aegypti larvae and Drosophila melanogaster) or by most midgut cells (type 1, found in Ae. aegypti adult and Tribolium castaneum). LG and D. melanogaster peritrophins have more complex domain structures than AG and T. castaneum peritrophins. Furthermore, mucin-like domains of peritrophins from T. castaneum (feeding on rough food) are lengthier than those of adult Ae. aegypti (blood-feeding). This suggests, for the first time, that type 1 and type 2 PM may have variable molecular architectures determined by different peritrophins and/or ancillary proteins, which may be partly modulated by diet.

Designing Novel Antitrypanosomal Agents from a Mixed Graph-Theoretical Substructural Approach

Chagas disease is nowadays the most serious parasitic health problem. This disease is caused by Trypanosoma cruzi. The great number of deaths and the insufficient effectiveness of drugs against this parasite have alarmed the scientific community worldwide. In an attempt to overcome this problem, a model for the design and prediction of new antitrypanosomal agents was obtained. This used a mixed approach, containing simple descriptors based on fragments and topological substructural molecular design descriptors. A data set was made up of 188 compounds, 99 of them characterized an antitrypanosomal activity and 88 compounds that belong to other pharmaceutical categories. The model showed sensitivity, specificity and accuracy values above 85%. Quantitative fragmental contributions were also calculated. Then, and to confirm the quality of the model, 15 structures of molecules tested as antitrypanosomal compounds (that we did not include in this study) were predicted, taking into account the information on the abovementioned calculated fragmental contributions. The model showed an accuracy of 100% which means that the ""in silico"" methodology developed by our team is promising for the rational design of new antitrypanosomal drugs. (C) 2009 Wiley Periodicals, Inc. J Comput Chem 31: 882-894. 2010

Three-Dimensional Quantitative Structure-Activity Relationship Study of Antitumor 2-Formylpyridine Thiosemicarbazones Derivatives as Inhibitors of Ribonucleotide Reductase

In order to extend previous SAR and QSAR studies, 3D-QSAR analysis has been performed using CoMFA and CoMSIA approaches applied to a set of 39 alpha-(N)-heterocyclic carboxaldehydes thiosemicarbazones with their inhibitory activity values (IC(50)) evaluated against ribonucleotide reductase (RNR) of H.Ep.-2 cells (human epidermoid carcinoma), taken from selected literature. Both rigid and field alignment methods, taking the unsubstituted 2-formylpyridine thiosemicarbazone in its syn conformation as template, have been used to generate multiple predictive CoMFA and CoMSIA models derived from training sets and validated with the corresponding test sets. Acceptable predictive correlation coefficients (Q(cv)(2) from 0.360 to 0.609 for CoMFA and Q(cv)(2) from 0.394 to 0.580 for CoMSIA models) with high fitted correlation coefficients (r` from 0.881 to 0.981 for CoMFA and r(2) from 0.938 to 0.993 for CoMSIA models) and low standard errors (s from 0.135 to 0.383 for CoMFA and s from 0.098 to 0.240 for CoMSIA models) were obtained. More precise CoMFA and CoMSIA models have been derived considering the subset of thiosemicarbazones (TSC) substituted only at 5-position of the pyridine ring (n=22). Reasonable predictive correlation coefficients (Q(cv)(2) from 0.486 to 0.683 for CoMFA and Q(cv)(2) from 0.565 to 0.791 for CoMSIA models) with high fitted correlation coefficients (r(2) from 0.896 to 0.997 for CoMFA and r(2) from 0.991 to 0.998 for CoMSIA models) and very low standard errors (s from 0.040 to 0.179 for CoMFA and s from 0.029 to 0.068 for CoMSIA models) were obtained. The stability of each CoMFA and CoMSIA models was further assessed by performing bootstrapping analysis. For the two sets the generated CoMSIA models showed, in general, better statistics than the corresponding CoMFA models. The analysis of CoMFA and CoMSIA contour maps suggest that a hydrogen bond acceptor near the nitrogen of the pyridine ring can enhance inhibitory activity values. This observation agrees with literature data, which suggests that the nitrogen pyridine lone pairs can complex with the iron ion leading to species that inhibits RNR. The derived CoMFA and CoMSIA models contribute to understand the structural features of this class of TSC as antitumor agents in terms of steric, electrostatic, hydrophobic and hydrogen bond donor and hydrogen bond acceptor fields as well as to the rational design of this key enzyme inhibitors.

Volsurf Descriptors to Analyse Anti-HCV and Cytotoxic Activities of Sesquiterpene Lactones from Asteraceae Family

Various significant anti-HCV and cytotoxic sesquiterpene lactones (SLs) have been characterized. In this work, the chemometric tool Principal Component Analysis (PCA) was applied to two sets of SLs and the variance of the biological activity was explored. The first principal component accounts for as much of the variability in the data as possible, and each succeeding component accounts for as much of the remaining variability as possible. The calculations were performed using VolSurf program. For anti-HCV activity, PC1 (First Principal Component) explained 30.3% and PC2 (Second Principal Component) explained 26.5% of matrix total variance, while for cytotoxic activity, PC1 explained 30.9% and PC2 explained 15.6% of the total variance. The formalism employed generated good exploratory and predictive results and we identified some structural features, for both sets, important to the suitable biological activity and pharmacokinetic profile.

Sisal Fibers Treated with NaOH and Benzophenonetetracarboxylic Dianhydride as Reinforcement of Phenolic Matrix

In this work, composites based on a phenolic matrix and untreated- and treated sisal fibers were prepared. The treated sisal fibers used were those reacted with NaOH 2% solution and esterified using benzophenonetetracarboxylic dianhydride (BTDA). These treated fibers were modified with the objective of improving the adhesion of the fiber-matrix interface, which in turn influences the properties of the composites. BTDA was chosen as the esterifying agent to take advantage of the possibility of introducing; the polar and aromatic groups that are also present in the matrix structure into the surface of the fiber, which could then intensify the interactions occurring in the fiber-matrix interface. The fibers were then analyzed by SEM and FTIR to ascertain their chemical composition. The results showed that the fibers had been successfully modified. The composites (reinforced with 15%, w/w of 3.0 cm length sisal fiber randomly distributed) were characterized by SEM, impact strength, and water absorption capacity. In the tests conducted, the response of the composites was affected both by properties of the matrix and the fibers, besides the interfacial properties of the fiber-matrix. Overall, the results showed that the fiber treatment resulted in a composite that was less hygroscopic although with somewhat lower impact strength, when compared with the composite reinforced with untreated sisal fibers. (C) 2009 Wiley Periodicals, Inc. J Appl Polym Sci 115: 269-276, 2010

Local Structure and Photochromic Response in Ormosils Containing Dodecatungstophosphoric Acid

Two hybrid materials based on dodecatungstophosphoric acid (HPW) dispersed in ormosils modified with 3-aminopropiltrietoxysilane (APTS) or with N-(3-(trimethoxysilyl)-propyl)-ethylene-diamine (TSPEN) show reversible photochromic response induced by irradiation in the 200-390 nm UV range. A set of solid-state nuclear magnetic resonance (NMR) techniques was used to analyze the structural properties of the main components of these hybrids (the HPW polyanion, the inorganic matrix, and the organic functionalities). For the ormosils, the use of (29)Si NMR, {(1)H}-(29)Si cross-polarization, and {(1)H}-(29)Si HETCOR revealed a homogeneous distribution of silicon species Q ``, T(2), and T(3) for the APTS hybrid, contrasting with the separation of T(3) species in the TSPEN hybrid. The combination of (31)P NMR, {(1)H}-(31)P cross-polarization and (31)P-{(1)H} spin-echo double resonance (SEDOR) revealed the dispersion of the HPW ions in the ormosil, occupying sites with a high number of close protons (>50). Differences in the molecular dynamics at room temperature, inferred from SEDOR experiments, indicate a state of restricted mobility of the HPW ion and the surrounding molecular groups in the TSPEN hybrid. This behavior is consistent with the presence of more amino groups in the TSPEN, acting as chelating groups to the HPW ion. This hybrid, with the strong chelate interaction of the diamine group, shows the most intense photochromic response, in agreement with the charge transfer models proposed to explain the photochromic effect. Electronic reflectance spectroscopy in irradiated samples revealed the presence of one-electron and two-electron reduced polyanions. The one-electron reduced species could be detected also by (31)P NMR spectroscopy immediately after UV irradiation.