From Rational Design of Drug Crystals to Understanding of Nucleic Acid Structures: Lamivudine Duplex

A DNA-like duplex of nucleosides is probable to exist even without the 5`-phosphate groups needed to assemble the chain backbone. However, double-stranded helical structures of nucleosides are unknown. Here, we report a duplex of nucleoside analogs that is spontaneously assembled due to stacking of the neutral and protonated molecules of lamivudine, a nucleoside reverse transcriptase inhibitor (NTRI) widely used in anti-HIV drug combinatory medication. The left-handed lamivudine duplex has features similar to those of i-motif DNA, as the face-to-face base stacking and the helix rise per base pair. Furthermore, the protonation pattern on alternate bases expected for it DNA-like duplex stabilized by pairing of neutral and protonated cytosine fragments was observed for the first time in the lamivudine double-stranded helix. This structure demonstrates that hydrogen bonds can substitute for covalent phosphodiester linkage in the stabilization of the duplex backbone. This interesting example of spontaneous molecular self-organization indicates that the 5`-phosphate group could not be a requirement for duplex assembly.

Crystal Engineering of an Anti-HIV Drug Based on the Recognition of Assembling Molecular Frameworks

A rational strategy was employed for design of an orthorhombic structure of lamivudine with maleic acid. On the basis of the lamivudine saccharinate structure reported in the literature, maleic acid was chosen to synthesize a salt with the anti-HIV drug because of the structural similarities between the salt formers. Maleic acid has an acid-ionization constant of the anti first proton and an arrangement of their hydrogen bonding functionalities similar to those of saccharin. Likewise, there is a saccharin-like conformational rigidity in maleic acid because of the hydrogen-bonded ring formation and the Z-configuration around the C=C double bond. As was conceivably predicted, lamivudine maleate assembles into a structure whose intermolecular architecture is related to that of saccharinate salt of the drug. Therefore, a molecular framework responsible for crystal assembly into a lamivudine saccharinate-like structure could be recognized in the salt formers. Furthermore, structural correlations and structure-solubility relationships were established for lamivudine maleate and saccharinate. Although there is a same molecular framework in maleic acid and saccharin, these salt formers are Structurally different in some aspects. When compared to saccharin, neither out-of-plane SO(2) oxygens nor a benzene group occur in maleic acid. Both features could be related to higher solubility of lamivudine maleate. Here, we also anticipate that multicomponent molecular crystals of lamivudine with other salt formers possessing the molecular framework responsible for crystal assembly can be engineered successfully.

Quantitative structure-activity relationships for a series of inhibitors of cruzain from Trypanosoma cruzi: Molecular modeling, CoMFA and CoMSIA studies

Human parasitic diseases are the foremost threat to human health and welfare around the world. Trypanosomiasis is a very serious infectious disease against which the currently available drugs are limited and not effective. Therefore, there is an urgent need for new chemotherapeutic agents. One attractive drug target is the major cysteine protease from Trypanosoma cruzi, cruzain. In the present work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted on a series of thiosemicarbazone and semicarbazone derivatives as inhibitors of cruzain. Molecular modeling studies were performed in order to identify the preferred binding mode of the inhibitors into the enzyme active site, and to generate structural alignments for the three-dimensional quantitative structure-activity relationship (3D QSAR) investigations. Statistically significant models were obtained (CoMFA. r(2) = 0.96 and q(2) = 0.78; CoMSIA, r(2) = 0.91 and q(2) = 0.73), indicating their predictive ability for untested compounds. The models were externally validated employing a test set, and the predicted values were in good agreement with the experimental results. The final QSAR models and the information gathered from the 3D CoMFA and CoMSIA contour maps provided important insights into the chemical and structural basis involved in the molecular recognition process of this family of cruzain inhibitors, and should be useful for the design of new structurally related analogs with improved potency. (C) 2009 Elsevier Inc. All rights reserved.

Conformational Polymorphism in Racemic Crystals of the Diuretic Drug Chlortalidone

Chlortalidone (HIGROTON) is a diuretic drug widely used in antihypertensive therapy. Thus far, only two solid-state polymorphs of chlortalidone have been reported. We elucidated the structure of chlortalidone form I and a new polymorph. This new phase, namely, chlortalidone form III, was also entirely characterized. It was possible to conclude that it is a conformer with a different orientation of the chlorobenzenesulfonamide moiety. Compared to form I, it has a rotation of about 90 degrees on the axis of the C-C bond bridging the substituted phenyl and isoindolinyl rings. This conformational feature is related to the crystal packing patterns of the chlortalidone forms. Furthermore, certain intermolecular hydrogen bonds are present in both polymorphs, giving rise to ribbons with chlortalidone enantiomers alternately placed into them. The chlortalidone form I and form III crystallize in the triclinic space group P (1) over bar as racemic mixtures. Additional conformational details also differentiate the chlortalidone conformers. Slight twists on the isoindolinyl and sulfamyl groups exist. Considering all structural relationships, the fingerprint plots derived from the Hirshfeld surfaces exhibited the characteristics of the chlortalidone form I and form III crystal structures.

Structure-Based Approach for the Study of Estrogen Receptor Binding Affinity and Subtype Selectivity

Estrogens exert important physiological effects through the modulation of two human estrogen receptor (hER) subtypes, alpa (hER alpha) and beta (hER beta). Because the levels and relative proportion of hER alpha and hER beta differ significantly in different target cells, selective hER ligands could target specific tissues or pathways regulated by one receptor subtype without affecting the other. To understand the structural and chemical basis by which small molecule modulators are able to discriminate between the two subtypes, we have applied three-dimensional target-based approaches employing a series of potent hER-ligands. Comparative molecular field analysis (CoMFA) studies were applied to a data set of 81 hER modulators, for which binding affinity values were collected for both hER alpha and hER beta. Significant statistical coefficients were obtained (hER alpha, q(2) = 0.76; hER beta, q(2) = 0.70), indicating the internal consistency of the models. The generated models were validated using external test sets, and the predicted values were in good agreement with the experimental results. Five hER crystal structures were used in GRID/PCA investigations to generate molecular interaction fields (MIF) maps. hER alpha and hER beta were separated using one factor. The resulting 3D information was integrated with the aim of revealing the most relevant structural features involved in hER subtype selectivity. The final QSAR and GRID/PCA models and the information gathered from 3D contour maps should be useful for the design or novel hER modulators with improved selectivity.

Three-dimensional quantitative structure-activity relationships for a large series of potent antitubercular agents

Comparative molecular field analysis (CoMFA) studies were conducted on a series of 100 isoniazid derivatives as anti-tuberculosis agents using two receptor-independent structural data set alignment strategies: (1) rigid-body fit, and (2) pharmacophore-based. Significant cross-validated correlation coefficients were obtained (CoMFA(1), q(2) = 0,75 and CoMFA(2), q(2) = 0.74), indicating the potential of the models for untested compounds. The models were then used to predict the inhibitory potency of 20 test set compounds that were not included in the training set, and the predicted values were in good agreement with the experimental results.

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.

Triple Structural Transition below Room Temperature in the Antifilarial Drug Diethylcarbamazine Citrate

A very unusual triple structural transition pattern below room temperature was observed for the antifilarial drug diethylcarbamazine citrate. Besides the first thermal, crystallographic, and vibrational investigations of this first-line drug used in clinical treatment for lymphatic filariasis, a noteworthy behavior with three structural transformations as a function of temperature was demonstrated by differential scanning calorimetry, Raman spectroscopy, and single-crystal X-ray diffractometry. Our X-ray data on single crystals allow for a complete featuring and understanding of all transitions, since the four structures associated with the three solid-solid phase transformations were accurately determined. Two of three structural transitions show an order-disorder mechanism and temperature hysteresis with exothermic peaks at 224 K (T(1)`) and 213 K (T(2)`) upon cooling and endothermic ones at 248 K (T(1)) and 226 K (T(2)) upon heating. The other transition occurs at 108 K (T(3)) and it is temperature-rate sensitive. Molecular displacements onto the (010) plane and conformational changes of the diethylcarbamazine backbone as a consequence of the C-H center dot center dot center dot N hydrogen bonding formation/cleavage between drug molecules explain the mechanism of the transitions at T(1)`/T(2). However, such changes are observed only on alternate columns of the drug intercalated by citrate chains, which leads to a doubling of the lattice period along the a axis of the 235 K structure with respect to the 150 and 293 K structures. At T(2)`/T(1), these structural alterations occur in all columns of the drug. At T(3), there is a rotation on the axis of the N-C bond between the carbamoyl moiety and an ethyl group of one crystallographically independent diethylcarbamazine molecule besides molecular shifts and other conformational alterations. The impact of this study is based on the fascinating finding in which the versatile capability of structural adaptation dependent on the thermal history was observed for a relatively simple organic salt, diethylcarbamazine citrate.

Solvothermal Preparation of Drug Crystals: Didanosine

For the first time, crystals of suitable size for X-ray diffractometry structure determination (Dian important anti-HI V drug were prepared under solvothermal conditions. In this study, the crystal structure of didanosine (2`,3`-dideoxyinosine, ddI) in the form of a hydrate was determined using single-crystal X-ray diffractometry. Powder X-ray diffraction analysis revealed that the solid-state phase of the drug incorporated into pharmaceutical solid dosage forms is isostructural to the solvothermally prepared ddI material, even though they do not exhibit an identical chemical composition due to different water fractions occupying hydrophobic channels formed within the crystal lattice. Two ddI conformers are present in the structure, in agreement with a previous structure elucidation attempt. Concerning the keto enol equilibrium of ddI, our crystal data and vibrational characterizations by Fourier transform infrared (FTIR) and FT-Raman spectroscopy techniques were conclusive to state that both conformers exist in the keto form, contrary to solid-state NMR spectroscopic assignments that suggested ddI molecules occur as enol tautomers. In addition, characterizations by thermal (differential scanning calorimetry) and spectroscopic techniques allowed us to understand the structural similarities and the differences related to the hydration pattern of the nonstoichiometric hydrates.

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.

CDK9 a potential target for drug development

The family of Cyclin-Dependent Kinases (CDKs) can be subdivided into two major functional groups based on their roles in cell cycle and/or transcriptional control. CDK9 is the catalytic subunit of positive transcription elongation factor b (P-TEFb). CDK9 is the kinase of the TAK complex (Tat-associated kinase complex), and binds to Tat protein of HIV, suggesting a possible role for CDK9 in AIDS progression. CDK9 complexed with its regulatory partner cyclin T1, serves as a cellular mediator of the transactivation function of the HIV Tat protein. P-TEFb is responsible for the phosphorylation of the carboxyl-terminal domain of RNA Pol II, resulting in stimulation of transcription. Furthermore, the complexes containing CDK9 induce the differentiation in distinct tissue. The CDK9/cyclin T1 complex is expressed at higher level in more differentiated primary neuroectodermal and neuroblastoma tumors, showing a correlation between the kinase expression and tumor differentiation grade. This may have clinical and therapeutical implications for these tumor types. Among the CDK inhibitors two have shown to be effective against CDK9: Roscovitine and Flavopiridol. These two inhibitors prevented the replication of human immunodeficiency virus (HIV) type 1 by blocking Tat transactivation of the HIV type 1 promoter. These compounds inhibit CDKs by binding to the catalytic domain in place of ATP, preventing transfer of a phosphate group to the substrate. More sensitive therapeutic agents of CDK9 can be designed, and structural studies can add information in the understanding of this kinase. The major features related to CDK9 inhibition will be reviewed in this article.

Crystal Structure Determination for Eugenyl Acetate

Essential oils are good candidates for the substitution of conventional medicinal treatments. Many articles and patents for their use have been published in recent years. The most attractive aspects of using essential oils as medicaments are their natural source and rapid permeability. Besides permeability, the solubility behavior of a drug is a key determinant of its oral bioavailability. Based on these characteristics, the aim of this study was to synthesize an essential oil derivative compound, using the raw oil extracted from Syzygium aromaticum L., without previous purification. The Eugenol molecular modification may diminish the problems of water solubility and bioavailability. The Eugenyl acetate molecule was characterized and its molecular modification investigated, including its structural properties and stereochemistry. This study was performed applying techniques, such as carbon-13 nuclear magnetic resonance spectroscopy (C-13 NMR), X-ray crystallographic analysis (XRD), powder X-ray diffraction (PXRD) and microscopic recording.