Prediction methods and databases within chemoinformatics: emphasis on drugs and drug candidates

MOTIVATION: To gather information about available databases and chemoinformatics methods for prediction of properties relevant to the drug discovery and optimization process. RESULTS: We present an overview of the most important databases with 2-dimensional and 3-dimensional structural information about drugs and drug candidates, and of databases with relevant properties. Access to experimental data and numerical methods for selecting and utilizing these data is crucial for developing accurate predictive in silico models. Many interesting predictive methods for classifying the suitability of chemical compounds as potential drugs, as well as for predicting their physico-chemical and ADMET properties have been proposed in recent years. These methods are discussed, and some possible future directions in this rapidly developing field are described.     

Malaria and urbanization in sub-Saharan Africa

There are already 40 cities in Africa with over 1 million inhabitants and the United Nations Environmental Programme estimates that by 2025 over 800 million people will live in urban areas. Recognizing that malaria control can improve the health of the vulnerable and remove a major obstacle to their economic development, the Malaria Knowledge Programme of the Liverpool School of Tropical Medicine and the Systemwide Initiative on Malaria and Agriculture convened a multi-sectoral technical consultation on urban malaria in Pretoria, South Africa from 2nd to 4th December, 2004. The aim of the meeting was to identify strategies for the assessment and control of urban malaria. This commentary reflects the discussions held during the meeting and aims to inform researchers and policy makers of the potential for containing and reversing the emerging problem of urban malaria.     

Triple effect of mimetic peptides interfering with neural cell adhesion molecule homophilic cis interactions

The neural cell adhesion molecule (NCAM) is pivotal in neural development, regeneration, and learning. Here we characterize two peptides, termed P1-B and P2, derived from the homophilic binding sites in the first two N-terminal immunoglobulin (Ig) modules of NCAM, with regard to their effects on neurite extension and adhesion. To evaluate how interference of these mimetic peptides with NCAM homophilic interactions in cis influences NCAM binding in trans, we employed a coculture system in which PC12-E2 cells were grown on monolayers of fibroblasts with or without NCAM expression and the rate of neurite outgrowth subsequently was analyzed. P2, but not P1-B, induced neurite outgrowth in the absence of NCAM binding in trans. When PC12-E2 cells were grown on monolayers of NCAM-expressing fibroblasts, the effect of both P1-B and P2 on neurite outgrowth was dependent on peptide concentrations. P1-B and P2 acted as conventional antagonists, agonists, and reverse agonists of NCAM at low, intermediate, and high peptide concentrations, respectively. The demonstrated in vitro triple pharmacological effect of mimetic peptides interfering with the NCAM homophilic cis binding will be valuable for the understanding of the actions of these mimetics in vivo.     

Reversible dimerization of acid-denatured ACBP controlled by helix A4

The peptide segment corresponding to helix A4 in acyl-coenzyme-A-binding protein (ACBP) is an exceptionally stable helix in the denatured state of the protein as well as in its isolated form. Circular dichroism spectroscopy showed an alpha-helix content in the helix A4 peptide (HA4) of 45%, and under denaturing conditions at pH 2.3, helix conformations are still populated in 24% of the ensemble of molecules. The structure of HA4 at atomic resolution was assessed using nuclear magnetic resonance (NMR) spectroscopy. Long-range NOEs between remote residues at opposite peptide ends suggested the formation of an antiparallel homodimer, and the resulting structure was treated as the minimum higher-order structure. The dimerization property of helix A4 is maintained in the full-length protein under denaturing conditions. NMR diffusion studies and concentration-dependent experiments on ACBP at low pH proved the formation of dimers and revealed a cooperative stabilization of helix A4 in this process. This emphasizes its special role in the structure formation in the denatured state of ACBP. No dimers are formed in the presence of guanidine hydrochloride, which underlines the fundamental difference between the nature of these two denatured states.     

Structure of the plasminogen kringle 4 binding calcium-free form of the C-type lectin-like domain of tetranectin

Tetranectin is a homotrimeric protein containing a C-type lectin-like domain. This domain (TN3) can bind calcium, but in the absence of calcium, the domain binds a number of kringle-type protein ligands. Two of the calcium-coordinating residues are also critical for binding plasminogen kringle 4 (K4). The structure of the calcium free-form of TN3 (apoTN3) has been determined by NMR. Compared to the structure of the calcium-bound form of TN3 (holoTN3), the core region of secondary structural elements is conserved, while large displacements occur in the loops involved in calcium or K4 binding. A conserved proline, which was found to be in the cis conformation in holoTN3, is in apoTN3 predominantly in the trans conformation. Backbone dynamics indicate that, in apoTN3 especially, two of the three calcium-binding loops and two of the three K4-binding residues exhibit increased flexibility, whereas no such flexibility is observed in holoTN3. In the 20 best nuclear magnetic resonance structures of apoTN3, the residues critical for K4 binding span a large conformational space. Together with the relaxation data, this indicates that the K4-ligand-binding site in apoTN3 is not preformed.     

AICAR stimulates adiponectin and inhibits cytokines in adipose tissue

5-Aminoimidazole-4-carboxamide ribonucleoside (AICAR) can be used as an experimental tool to activate 5'-AMP-activated protein kinase (AMPK) and has been shown to improve insulin sensitivity. In parallel adiponectin also seems to activate AMPK and to improve insulin sensitivity. We have investigated the effects of AICAR on the gene expression of adiponectin and on gene expression and release of cytokines in human adipose tissue in vitro. AICAR stimulated AMPK alpha1 activity 3-4-fold (p<0.001), and dose-dependently increased adiponectin mRNA levels with significant stimulation (2-4-fold) at AICAR concentrations of 0.5-2mM (p<0.05). The adipose tissue protein release of tumor necrosis factor-alpha (TNF- alpha) and interleukin-6 (IL-6) was decreased by AICAR (p<0.05). In conclusion, AICAR stimulated adipose tissue AMPK alpha1 activity and adiponectin gene expression, while attenuating the release of TNF-alpha and IL-6. Reduced concentrations of these cytokines and increased levels of adiponectin might play a role for the insulin sensitizing effects of AICAR.     

Early Selection in Gag by Protective HLA Alleles Contributes to Reduced HIV-1 Replication Capacity That May Be Largely Compensated for in Chronic Infection

Mutations that allow escape from CD8 T-cell responses are common in HIV-1 and may attenuate pathogenesis by reducing viral fitness. While this has been demonstrated for individual cases, a systematic investigation of the consequence of HLA class I-mediated selection on HIV-1 in vitro replication capacity (RC) has not been undertaken. We examined this question by generating recombinant viruses expressing plasma HIV-1 RNA-derived Gag-Protease sequences from 66 acute/early and 803 chronic untreated subtype B-infected individuals in an NL4-3 background and measuring their RCs using a green fluorescent protein (GFP) reporter CD4 T-cell assay. In acute/early infection, viruses derived from individuals expressing the protective alleles HLA-B*57, -B*5801, and/or -B*13 displayed significantly lower RCs than did viruses from individuals lacking these alleles (P < 0.05). Furthermore, acute/early RC inversely correlated with the presence of HLA-B-associated Gag polymorphisms (R = −0.27; P = 0.03), suggesting a cumulative effect of primary escape mutations on fitness during the first months of infection. At the chronic stage of infection, no strong correlations were observed between RC and protective HLA-B alleles or with the presence of HLA-B-associated polymorphisms restricted by protective alleles despite increased statistical power to detect these associations. However, RC correlated positively with the presence of known compensatory mutations in chronic viruses from B*57-expressing individuals harboring the Gag T242N mutation (n = 50; R = 0.36; P = 0.01), suggesting that the rescue of fitness defects occurred through mutations at secondary sites. Additional mutations in Gag that may modulate the impact of the T242N mutation on RC were identified. A modest inverse correlation was observed between RC and CD4 cell count in chronic infection (R = −0.17; P < 0.0001), suggesting that Gag-Protease RC could increase over the disease course. Notably, this association was stronger for individuals who expressed B*57, B*58, or B*13 (R = −0.27; P = 0.004). Taken together, these data indicate that certain protective HLA alleles contribute to early defects in HIV-1 fitness through the selection of detrimental mutations in Gag; however, these effects wane as compensatory mutations accumulate in chronic infection. The long-term control of HIV-1 in some persons who express protective alleles suggests that early fitness hits may provide lasting benefits.The host immune response elicited by CD8⁺ cytotoxic T lymphocytes (CTLs) is a major contributor to viral control following human immunodeficiency virus type 1 (HIV-1) infection (6, 39), but antiviral pressure exerted by CTLs is diminished by the selection of escape mutations in targeted regions throughout the viral proteome (7, 18, 29, 35, 41, 45, 57). A comprehensive identification of HLA-associated viral polymorphisms has recently been achieved through population-based analyses of HIV-1 sequences and HLA class I types from different cohorts worldwide (3, 8, 13-15, 34, 43, 50, 56, 63). However, despite improved characterization of the sites and pathways of immune escape, effective ways to incorporate these findings into immunogen design remain an area of debate. A better understanding of the impact of escape mutations on viral fitness may provide novel directions for HIV-1 vaccines that are designed to attenuate pathogenesis.The development of innovative vaccine strategies that can overcome the extreme diversity of HIV is a key priority (4). One proposed approach is to target the most conserved T-cell epitopes, which presumably cannot escape from CTL pressure easily due to structural or functional constraints on the viral protein (55). Complementary approaches include the design of polyvalent and/or mosaic immunogens that incorporate commonly observed viral diversity (4, 38) or the specific targeting of vulnerable regions of the viral proteome that do escape but only at a substantial cost to viral replication capacity (RC) (1, 40). A chief target of such vaccine approaches is the major HIV-1 structural protein Gag, which is known to be highly immunogenic and to elicit CTL responses that correlate with the natural control of infection (22, 36, 66). Indeed, several lines of evidence support a relationship between the selection of CTL escape mutations and reduced HIV-1 fitness. These include the reversion of escape mutations following transmission to an HLA-mismatched recipient who cannot target the epitope (19, 24, 41) as well as reduced plasma viral load (pVL) set point following the transmission of certain escape variants from donors who expressed protective HLA alleles (17, 27). Notably, these in vivo observations have been made most often for variations within Gag that are attributed to CTL responses restricted by the protective alleles HLA-B*57 and -B*5801 (17, 19, 27, 41). Most recently, reduced in vitro RCs of clinical isolates and/or engineered strains encoding single or multiple escape mutations in Gag selected in the context of certain protective HLA alleles, including B*57, B*5801, B*27, and B*13, have been demonstrated (9, 10, 42, 53, 59, 62). Despite these efforts, the goal of a T-cell vaccine that targets highly conserved and attenuation-inducing sites is hampered by a lack of knowledge concerning the contribution of most escape mutations to HIV-1 fitness as well as a poor understanding of the relative influence of HLA on the viral RC at different stages of infection.The mutability of HIV-1 permits the generation of progeny viruses encoding compensatory mutations that restore normal protein function and/or viral fitness. Detailed studies have demonstrated that the in vitro RC of escape variants in human and primate immunodeficiency viruses can be enhanced by the addition of secondary mutations outside the targeted epitope (10, 20, 52, 59, 65). Thus, vaccine strategies aimed at attenuating HIV-1 must also consider, among other factors, the frequency, time course, and extent to which compensation might overcome attenuation mediated by CTL-induced escape. Despite its anticipated utility for HIV-1 vaccine design, systematic studies to examine the consequences of naturally occurring CTL escape and compensatory mutations on viral RC have not been undertaken.We have described previously an in vitro recombinant viral assay to examine the impact of Gag-Protease mutations on HIV-1 RC (47, 49). Gag and protease have been included in each virus to minimize the impact of sequence polymorphisms at Gag cleavage sites, which coevolve with changes in protease (5, 37). Using this approach, we have demonstrated that viruses derived from HIV-1 controllers replicated significantly less well than those derived from noncontrollers and that these differences were detectable at both the acute/early (49) and chronic (47) stages. Escape mutations in Gag associated with the protective HLA-B*57 allele, as well as putative compensatory mutations outside known CTL epitopes, contributed to this difference in RC (47). However, substantial variability was observed for viruses from controllers and noncontrollers, indicating that additional factors were likely to be involved. Benefits of this assay include its relatively high-throughput capacity as well as the fact that clinically derived HIV-1 sequences are used in their entirety. Thus, it is possible to examine a large number of “real-world” Gag-Protease sequences, to define an RC value for each one, and to identify sequences within the population of recombinant strains that are responsible for RC differences.Here, we use this recombinant virus approach to examine the contribution of HLA-associated immune pressure on Gag-Protease RC during acute/early (n = 66) and chronic (n = 803) infections in the context of naturally occurring HIV-1 subtype B isolates from untreated individuals. In a recent report (64), we employed this system to examine the Gag-Protease RC in a similar cohort of chronic HIV-1 subtype C-infected individuals. The results of these studies provide important insights into the roles of immune pressure and fitness constraints on HIV-1 evolution that may contribute to the rational design of an effective vaccine.