Ab initio determination of the flexibility of 2'-aminoribonucleosides and 2'-aminoarabinonucleosides inserted in duplexes.

The sugar puckering of adenosine and uridine nucleosides with an amino group at 2' in the ribo or arabino orientations are determined using high-level quantum mechanical calculations Only the conformations that have dihedrals compatible with their insertion into a duplex are retained. The amino group has always been found to be pyramidal and its orientation governs the conformation of the sugar. The energetically most favorable conformation of the 2'-aminoribonucleosides has the south puckering but must be discarded. For another orientation of the 2'-amino group, the conformation is energetically less favorable but has the north puckering. Calculations performed in the presence of a water molecule give similar results but with a smaller energy gap. The model then explains why the insertion of a 2'-aminoribonucleotide destabilizes double-stranded RNAs and also double-stranded DNAs. In the arabino orientation, an NH(2) substituent at 2' favors north puckering. In contrast to 2'-aminoribonucleosides, deoxynucleosides inserted into a duplex remain in the most energetically favorable conformation compatible with the canonical values of the torsion angles. The whole relaxed potential map, in the amplitude/pseudorotation space, shows that for natural deoxyadenosine there is only one valley in the east running from south to north puckering.

Comparative molecular surface analysis (CoMSA) for virtual combinatorial library screening of styrylquinoline HIV-1 blocking agents.

We used comparative molecular surface analysis to design molecules for the synthesis as part of the search for new HIV-1 integrase inhibitors. We analyzed the virtual combinatorial library (VCL) constituted from various moieties of styrylquinoline and styrylquinazoline inhibitors. Since imines can be applied in a strategy of dynamic combinatorial chemistry (DCC), we also tested similar compounds in which the -C=N- or -N=C- linker connected the heteroaromatic and aromatic moieties. We then used principal component analysis (PCA) or self-organizing maps (SOM), namely, the Kohonen neural networks to obtain a clustering plot analyzing the diversity of the VCL formed. Previously synthesized compounds of known activity, used as molecular probes, were projected onto this plot, which provided a set of promising virtual drugs. Moreover, we further modified the above mentioned VCL to include the single bond linker -C-N- or -N-C-. This allowed increasing compound stability but expanded also the diversity between the available molecular probes and virtual targets. The application of the CoMSA with SOM indicated important differences between such compounds and active molecular probes. We synthesized such compounds to verify the computational predictions.

Probing of HIV-1 integrase/DNA interactions using novel analogs of viral DNA.

The specific activity of the human immunodeficiency virus, type 1 (HIV-1), integrase on the viral long terminal repeat requires the binding of the enzyme to certain sequences located in the U3 and U5 regions at the ends of viral DNA, but the determinants of this specific DNA-protein recognition are not yet completely understood. We synthesized DNA duplexes mimicking the U5 region and containing either 2'-modified nucleosides or 1,3-propanediol insertions and studied their interactions with HIV-1 integrase, using Mn2+ or Mg2+ ions as integrase cofactors. These DNA modifications had no strong effect on integrase binding to the substrate analogs but significantly affected 3'-end processing rate. The effects of nucleoside modifications at positions 5, 6, and especially 3 strongly depended on the cationic cofactor used. These effects were much more pronounced in the presence of Mg2+ than in the presence of Mn2+. Modifications of base pairs 7-9 affected 3'-end processing equally in the presence of both ions. Adenine from the 3rd bp is thought to form at least two hydrogen bonds with integrase that are crucial for specific DNA recognition. The complementary base, thymine, is not important for integrase activity. For other positions, our results suggest that integrase recognizes a fine structure of the sugar-phosphate backbone rather than heterocyclic bases. Integrase interactions with the unprocessed strand at positions 5-8 are more important than interactions with the processed strand for specific substrate recognition. Based on our results, we suggest a model for integrase interaction with the U5 substrate.

Intermolecular interactions in the crystal structures of potential HIV-1 integrase inhibitors.

2-[(2,5-dichloro-4-nitro-phenylamino)-methoxy-methyl]-8-hydroxy-quinoline 1 and 2-methyl-quinoline-5,8-dione-5-oxime 2 were obtained as potential HIV-1 integrase inhibitors and analyzed by X-ray crystallography. Semiempirical theoretical calculations of energy preferred conformations were also carried out. The crystal structures of both compounds are stabilized via hydrogen bonds and pi-pi stacking interactions. The planarity of compound 1 is caused by intramolecular hydrogen bonds.

Molecular modeling of wild-type and D816V c-Kit inhibition based on ATP-competitive binding of ellipticine derivatives to tyrosine kinases.

The D816V activating mutation of the c-Kit kinase domain often causes human mastocytosis. Although inhibitors of wild-type c-Kit are known (e.g. STI-571), they are at least 10 times less active against the c-Kit mutant. Several derivatives of ellipticine (5,11-dimethyl-6H-pyrido[3,4-b]carbazole), substituted at positions 1, 2, 9, and 11, were found to inhibit purified D816V and wild-type c-Kit kinase domains with comparable potencies by competing with ATP binding. We investigated the difference between these inhibitors by modeling the D816V mutation in crystal structures of inactive and active c-Kit. Molecular dynamics simulations strongly suggested that the D816V point mutation shifts the conformational equilibrium of c-Kit kinase domain toward the active conformation. All ellipticine compounds were subsequently docked to the D816V mutant c-Kit model. The model provides possible explanations for the structure-activity relationships observed among ellipticine compounds, resulting in new insights into D816V c-Kit mutant inhibition.

Functional analysis of novel sonic hedgehog gene mutations identified in basal cell carcinomas from xeroderma pigmentosum patients.

Altered sonic hedgehog (SHH) signaling is crucial in the development of basal cell carcinomas (BCC), the most common human cancer. Mutations in SHH signal transducers, PATCHED and SMOOTHENED, have already been identified, but SHH mutations are extremely rare; only 1 was detected in 74 sporadic BCCs. We present data showing unique SHH mutations in BCCs from repair-deficient, skin cancer-prone xeroderma pigmentosum (XP) patients, which are characterized by high levels of UV-specific mutations in key genes involved in skin carcinogenesis, including PATCHED and SMOOTHENED. Thus, 6 UV-specific SHH mutations were detected in 5 of 33 XP BCCs. These missense SHH alterations are not activating mutations for its postulated proto-oncogene function, as the mutant SHH proteins do not show transforming activity and induce differentiation or stimulate proliferation to the same level as the wild-type protein. Structural modeling studies of the 4 proteins altered at the surface residues, G57S, G64K, D147N, and R155C, show that they do not effect the protein conformation. Interestingly, they are all located on one face of the compact SHH protein suggesting that they may have altered affinity for different partners, which may be important in altering other functions. Additional functional analysis of the SHH mutations found in vivo in XP BCCs will help shed light on the role of SHH in skin carcinogenesis. In conclusion, we report for the first time, significant levels of SHH mutations found only in XP BCCs and none in squamous cell carcinomas, indicating their importance in the specific development of BCCs.

Linker-modified quinoline derivatives targeting HIV-1 integrase: synthesis and biological activity.

A novel series of HIV-1 integrase inhibitors was synthesized and tested in both in vitro and ex vivo assays. These inhibitors are featured by the presence of a quinoline subunit and an ancillary aromatic ring linked by functionalized spacers such as amide, hydrazide, urea and 1-hydroxyprop-1-en-3-one moiety. Amide derivatives are the most promising ones and could serve as leads for further developments.

Mechanism of HIV-1 integrase inhibition by styrylquinoline derivatives in vitro.

Styrylquinoline derivatives (SQ) efficiently inhibit the 3'-processing activity of integrase (IN) with IC50 values of between 0.5 and 5 microM. We studied the mechanism of action of these compounds in vitro. First, we used steady-state fluorescence anisotropy to assay the effects of the SQ derivatives on the formation of IN-viral DNA complexes independently of the catalytic process. The IC50 values obtained in activity and DNA-binding tests were similar, suggesting that the inhibition of 3'-processing can be fully explained by the prevention of IN-DNA recognition. SQ compounds act in a competitive manner, with Ki values of between 400 and 900 nM. In contrast, SQs did not inhibit 3'-processing when IN-DNA complexes were preassembled. Computational docking followed or not by molecular dynamics using the catalytic core of HIV-1 IN suggested a competitive inhibition mechanism, which is consistent with our previous data obtained with the corresponding Rous sarcoma virus domain. Second, we used preassembled IN-preprocessed DNA complexes to assay the potency of SQs against the strand transfer reaction, independently of 3'-processing. Inhibition occurred even if the efficiency was decreased by about 5- to 10-fold. Our results suggest that two inhibitor-binding modes exist: the first one prevents the binding of the viral DNA and then the two subsequent reactions (i.e., 3'-processing and strand transfer), whereas the second one prevents the binding of target DNA, thus inhibiting strand transfer. SQ derivatives have a higher affinity for the first site, in contrast to that observed for the diketo acids, which preferentially bind to the second one.

Differential effects of two growth inhibitory K vitamin analogs on cell cycle regulating proteins in human hepatoma cells.

A comparison was made between two K vitamin analogs. Growth in vitro of Hep G2 hepatoma cells was inhibited both by Compound 5 (Cpd 5), a recently synthesized thioalkyl analog of vitamin K or 2-(2-mercaptoethanol)-3-methyl-1, 4-naphthoquinone, as well as by synthetic vitamin K3 (menadione). Using synchronized Hep G2 hepatoma cells, the actions of both Cpd 5 and vitamin K3 on cell cycle regulating proteins were examined. Cpd 5 decreased the levels of cyclin D1, Cdk4, p16, p21 and cyclin B1. By contrast, VK3 only decreased the level of cyclin D1, but had no effect on the levels of Cdk4, p16 or p21. Interestingly, both VK3 and VK2 increased the levels of p21. The naturally occurring K vitamins had little effect on cell growth and none on the cyclins or Cdks. Amounts and activity of the G1/S phase controlling Cdc25A were measured. We found that Cpd 5 directly inhibited both Cdc25A activity and its protein expression, whereas VK3 did not. Thus, the main effects of Cpd 5 were on G1 and S phase proteins, especially Cdk4 and Cdc25A amounts in contrast to VK3. Computer docking studies of Cpd 5 and VK3 to Cdc25A phosphatase showed three binding sites. In the best conformation, Cpd 5 was found to be closer to the enzyme active site than VK3. These findings show that Cpd 5 represents a new class of anticancer agent, being a protein tyrosine phosphatase (PTP) antagonist, that binds to Cdc25A with suppression of its activity. Tumors expressing high levels of oncogenic Cdc25A phosphatase may thus be susceptible to the growth inhibitory activities of this class of compound.

Use of the Kohonen neural network for rapid screening of ex vivo anti-HIV activity of styrylquinolines.

Using the Kohonen neural network, the electrostatic potentials on the molecular surfaces of 14 styrylquinoline derivatives were drawn as comparative two-dimensional maps and compared with their known human immunodeficiency virus (HIV)-1 replication blocking potency in cells. A feature of the potential map was discovered to be related with the HIV-1 blocking activity and was used to unmask the activity of further five analogues, previously described but whose cytotoxicity precluded an estimation of their activity, and to predict the activity of 10 new compounds while the experimental data were unknown. The measurements performed later turned out to agree with the predictions.

Hamming distance geometry of a protein conformational space: application to the clustering of a 4-ns molecular dynamics trajectory of the HIV-1 integrase catalytic core.

Protein structures can be encoded into binary sequences (Gabarro-Arpa et al., Comput Chem 2000;24:693-698) these are used to define a Hamming distance in conformational space: the distance between two different molecular conformations is the number of different bits in their sequences. Each bit in the sequence arises from a partition of conformational space in two halves. Thus, the information encoded in the binary sequences is also used to characterize the regions of conformational space visited by the system. We apply this distance and their associated geometric structures to the clustering and analysis of conformations sampled during a 4-ns molecular dynamics simulation of the HIV-1 integrase catalytic core. The cluster analysis of the simulation shows a division of the trajectory into two segments of 2.6 and 1.4 ns length, which are qualitatively different: the data points to the fact that equilibration is only reached at the end of the first segment. The Hamming distance is compared also to the r.m.s. deviation measure. The analysis of the cases studied so far shows that under the same conditions the two measures behave quite differently, and that the Hamming distance appears to be more robust than the r.m.s. deviation.