We present the first receptor-based pharmacophore model for HIV-1 integrase. The development of "dynamic" pharmacophore models is a new method that accounts for the inherent flexibility of the active site and aims to reduce the entropic penalties associated with binding a ligand. Furthermore, this new drug discovery method overcomes the limitation of an incomplete crystal structure of the target protein. A molecular dynamics (MD) simulation describes the flexibility of the uncomplexed protein. Many conformational models of the protein are saved from the MD simulations and used in a series of multi-unit search for interacting conformers (MUSIC) simulations. MUSIC is a multiple-copy minimization method, available in the BOSS program; it is used to determine binding regions for probe molecules containing functional groups that complement the active site. All protein conformations from the MD are overlaid, and conserved binding regions for the probe molecules are identified. Those conserved binding regions define the dynamic pharmacophore model. Here, the dynamic model is compared to known inhibitors of the integrase as well as a three-point, ligand-based pharmacophore model from the literature. Also, a "static" pharmacophore model was determined in the standard fashion, using a single crystal structure. Inhibitors thought to bind in the active site of HIV-1 integrase fit the dynamic model but not the static model. Finally, we have identified a set of compounds from the Available Chemicals Directory that fit the dynamic pharmacophore model, and experimental testing of the compounds has confirmed several new inhibitors.

No abstract

We report the crystal structure of alanine racemase from Mycobacterium tuberculosis (Alr(Mtb)) at 1.9 A resolution. In our structure, Alr(Mtb) is found to be a dimer formed by two crystallographically different monomers, each comprising 384 residues. The domain makeup of each monomer is similar to that of Bacillus and Pseudomonas alanine racemases and includes both an alpha/beta-barrel at the N-terminus and a C-terminus primarily made of beta-strands. The hinge angle between these two domains is unique for Alr(Mtb), but the active site geometry is conserved. In Alr(Mtb), the PLP cofactor is covalently bound to the protein via an internal aldimine bond with Lys42. No guest substrate is noted in its active site, although some residual electron density is observed in the enzyme's active site pocket. Analysis of the active site pocket, in the context of other known alanine racemases, allows us to propose the inclusion of conserved residues found at the entrance to the binding pocket as additional targets in ongoing structure-aided drug design efforts. Also, as observed in other alanine racemase structures, PLP adopts a conformation that significantly distorts the planarity of the extended conjugated system between the PLP ring and the internal aldimine bond.

No abstract

Structure-based computational methods continue to enhance progress in the discovery and refinement of therapeutic agents. Several such methods and their applications are described. These include molecular visualization and molecular modeling, docking, fragment methods, 3-D database techniques, and free-energy perturbation. Related issues that are discussed include the use of simplified potential energy functions and the determination of the positions of tightly bound waters. Strengths and weaknesses of the various methods are described.

Breast cancer brain metastasis is resistant to therapy and a particularly poor prognostic feature in patient survival. Altered metabolism is a common feature of cancer cells but little is known as to what metabolic changes benefit breast cancer brain metastases. We found that brain-metastatic breast cancer cells evolved the ability to survive and proliferate independent of glucose due to enhanced gluconeogenesis and oxidations of glutamine and branched chain amino acids, which together sustain the non-oxidative pentose pathway for purine synthesis. Silencing expression of fructose-1,6-bisphosphatases (FBPs) in brain metastatic cells reduced their viability and improved the survival of metastasis-bearing immunocompetent hosts. Clinically, we showed that brain metastases from human breast cancer patients expressed higher levels of FBP and glycogen than the corresponding primary tumors. Together, our findings identify a critical metabolic condition required to sustain brain metastasis, and suggest that targeting gluconeogenesis may help eradicate this deadly feature in advanced breast cancer patients.

The HIV-1 integrase, which is essential for viral replication, catalyzes the insertion of viral DNA into the host chromosome, thereby recruiting host cell machinery into making viral proteins. It represents the third main HIV enzyme target for inhibitor design, the first two being the reverse transcriptase and the protease. Two 1-ns molecular dynamics simulations have been carried out on completely hydrated models of the HIV-1 integrase catalytic domain, one with no metal ions and another with one magnesium ion in the catalytic site. The simulations predict that the region of the active site that is missing in the published crystal structures has (at the time of this work) more secondary structure than previously thought. The flexibility of this region has been discussed with respect to the mechanistic function of the enzyme. The results of these simulations will be used as part of inhibitor design projects directed against the catalytic domain of the enzyme.

The human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an essential enzyme in the life cycle of the virus and is an attractive target for the development of new drugs useful in acquired immunodeficiency syndrome multidrug therapy. Starting from the crystal structure of the 5CITEP inhibitor bound to the active site in the catalytic domain of the HIV-1 IN, two different molecular dynamics simulations in water have been carried out. In the first simulation the wild-type IN was used, whereas in the second one the double mutation T66I/M154I, described to lead to drug resistance, was introduced in the protein. Compelling differences have been observed in these two structures during analyses of the molecular dynamics trajectories, particularly in the inhibitor binding modes and in the conformational flexibility of the loop (residues 138-149) located near the three catalytic residues in the active site (Asp(64), Asp(116), Glu(152)). Because the conformational flexibility of this region is important for efficient biological activity and its behavior is quite different in the two models, we suggest a hypothetical mechanism for the inhibition and drug resistance of HIV-1 IN. These results can be useful for the rational design of more potent and selective integrase inhibitors and may allow for the design of inhibitors that will be more robust against known resistance mutations.