A structural model for the retroviral proteases

The CARD8 inflammasome in HIV infection

The biggest challenge to immune control of HIV infection is the rapid within-host viral evolution, which allows selection of viral variants that escape from T cell and antibody recognition. Thus, it is impossible to clear HIV infection without targeting "immutable" components of the virus. Unlike the adaptive immune system that recognizes cognate epitopes, the CARD8 inflammasome senses the essential enzymatic activity of the HIV-1 protease, which is immutable for the virus. Hence, all subtypes of HIV clinical isolates can be recognized by CARD8. In HIV-infected cells, the viral protease is expressed as a subunit of the viral Gag-Pol polyprotein and remains functionally inactive prior to viral budding. A class of anti-HIV drugs, the non-nucleoside reverse transcriptase inhibitors (NNRTIs), can promote Gag-pol dimerization and subsequent premature intracellular activation of the viral protease. NNRTI treatment triggers CARD8 inflammasome activation, which leads to pyroptosis of HIV-infected CD4⁺ T cells and macrophages. Targeting the CARD8 inflammasome can be a potent and broadly effective strategy for HIV eradication.

Beyond Inhibition: A Novel Strategy of Targeting HIV-1 Protease to Eliminate Viral Reservoirs

HIV-1 protease (PR) is a viral enzyme that cleaves the Gag and Gag-Pol polyprotein precursors to convert them into their functional forms, a process which is essential to generate infectious viral particles. Due to its broad substrate specificity, HIV-1 PR can also cleave certain host cell proteins. Several studies have identified host cell substrates of HIV-1 PR and described the potential impact of their cleavage on HIV-1-infected cells. Of particular interest is the interaction between PR and the caspase recruitment domain-containing protein 8 (CARD8) inflammasome. A recent study demonstrated that CARD8 can sense HIV-1 PR activity and induce cell death. While PR typically has low levels of intracellular activity prior to viral budding, premature PR activation can be achieved using certain non-nucleoside reverse transcriptase inhibitors (NNRTIs), resulting in CARD8 cleavage and downstream pyroptosis. Used together with latency reversal agents, the induction of premature PR activation to trigger CARD8-mediated cell killing may help eliminate latent reservoirs in people living with HIV. This represents a novel strategy of utilizing PR as an antiviral target through premature activation rather than inhibition. In this review, we discuss the viral and host substrates of HIV-1 protease and highlight potential applications and advantages of targeting CARD8 sensing of HIV-1 PR.

Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations

HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.

Structures of Substrate Complexes of Foamy Viral Protease-Reverse Transcriptase

Reverse transcriptases (RTs) use their DNA polymerase and RNase H activities to catalyze the conversion of single-stranded RNA to double-stranded DNA (dsDNA), a crucial process for the replication of retroviruses. Foamy viruses (FVs) possess a unique RT, which is a fusion with the protease (PR) domain. The mechanism of substrate binding by this enzyme has been unknown. Here, we report a crystal structure of monomeric full-length marmoset FV (MFV) PR-RT in complex with an RNA/DNA hybrid substrate. We also describe a structure of MFV PR-RT with an RNase H deletion in complex with a dsDNA substrate in which the enzyme forms an asymmetric homodimer. Cryo-electron microscopy reconstruction of the full-length MFV PR-RT–dsDNA complex confirmed the dimeric architecture. These findings represent the first structural description of nucleic acid binding by a foamy viral RT and demonstrate its ability to change its oligomeric state depending on the type of bound nucleic acid.

IMPORTANCE Reverse transcriptases (RTs) are intriguing enzymes converting single-stranded RNA to dsDNA. Their activity is essential for retroviruses, which are divided into two subfamilies differing significantly in their life cycles: Orthoretrovirinae and Spumaretrovirinae. The latter family is much more ancient and comprises five genera. A unique feature of foamy viral RTs is that they contain N-terminal protease (PR) domains, which are not present in orthoretroviral enzymes. So far, no structural information for full-length foamy viral PR-RT interacting with nucleic substrates has been reported. Here, we present crystal and cryo-electron microscopy structures of marmoset foamy virus (MFV) PR-RT. These structures revealed the mode of binding of RNA/DNA and dsDNA substrates. Moreover, unexpectedly, the structures and biochemical data showed that foamy viral PR-RT can adopt both a monomeric configuration, which is observed in our structures in the presence of an RNA/DNA hybrid, and an asymmetric dimer arrangement, which we observed in the presence of dsDNA.

A novel class III endogenous retrovirus with a class I envelope gene in African frogs with an intact genome and developmentally regulated transcripts in Xenopus tropicalis

BACKGROUND

Retroviruses exist as exogenous infectious agents and as endogenous retroviruses (ERVs) integrated into host chromosomes. Such endogenous retroviruses (ERVs) are grouped into three classes roughly corresponding to the seven genera of infectious retroviruses: class I (gamma-, epsilonretroviruses), class II (alpha-, beta-, delta-, lentiretroviruses) and class III (spumaretroviruses). Some ERVs have counterparts among the known infectious retroviruses, while others represent paleovirological relics of extinct or undiscovered retroviruses.

RESULTS

Here we identify an intact ERV in the Anuran amphibian, Xenopus tropicalis. XtERV-S has open reading frames (ORFs) for gag, pol (polymerase) and env (envelope) genes, with a small additional ORF in pol and a serine tRNA primer binding site. It has unusual features and domain relationships to known retroviruses. Analyses based on phylogeny and functional motifs establish that XtERV-S gag and pol genes are related to the ancient env-less class III ERV-L family but the surface subunit of env is unrelated to known retroviruses while its transmembrane subunit is class I-like. LTR constructs show transcriptional activity, and XtERV-S transcripts are detected in embryos after the maternal to zygotic mid-blastula transition and before the late tailbud stage. Tagged Gag protein shows typical subcellular localization. The presence of ORFs in all three protein-coding regions along with identical 5' and 3' LTRs (long terminal repeats) indicate this is a very recent germline acquisition. There are older, full-length, nonorthologous, defective copies in Xenopus laevis and the distantly related African bullfrog, Pyxicephalus adspersus. Additional older, internally deleted copies in X. tropicalis carry a 300 bp LTR substitution.

CONCLUSIONS

XtERV-S represents a genera-spanning member of the largely env-less class III ERV that has ancient and modern copies in Anurans. This provirus has an env ORF with a surface subunit unrelated to known retroviruses and a transmembrane subunit related to class I gammaretroviruses in sequence and organization, and is expressed in early embryogenesis. Additional XtERV-S-related but defective copies are present in X. tropicalis and other African frog taxa. XtERV-S is an unusual class III ERV variant, and it may represent an important transitional retroviral form that has been spreading in African frogs for tens of millions of years.

The impact of structural bioinformatics tools and resources on SARS-CoV-2 research and therapeutic strategies

SARS-CoV-2 is the causative agent of COVID-19, the ongoing global pandemic. It has posed a worldwide challenge to human health as no effective treatment is currently available to combat the disease. Its severity has led to unprecedented collaborative initiatives for therapeutic solutions against COVID-19. Studies resorting to structure-based drug design for COVID-19 are plethoric and show good promise. Structural biology provides key insights into 3D structures, critical residues/mutations in SARS-CoV-2 proteins, implicated in infectivity, molecular recognition and susceptibility to a broad range of host species. The detailed understanding of viral proteins and their complexes with host receptors and candidate epitope/lead compounds is the key to developing a structure-guided therapeutic design. Since the discovery of SARS-CoV-2, several structures of its proteins have been determined experimentally at an unprecedented speed and deposited in the Protein Data Bank. Further, specialized structural bioinformatics tools and resources have been developed for theoretical models, data on protein dynamics from computer simulations, impact of variants/mutations and molecular therapeutics. Here, we provide an overview of ongoing efforts on developing structural bioinformatics tools and resources for COVID-19 research. We also discuss the impact of these resources and structure-based studies, to understand various aspects of SARS-CoV-2 infection and therapeutic development. These include (i) understanding differences between SARS-CoV-2 and SARS-CoV, leading to increased infectivity of SARS-CoV-2, (ii) deciphering key residues in the SARS-CoV-2 involved in receptor-antibody recognition, (iii) analysis of variants in host proteins that affect host susceptibility to infection and (iv) analyses facilitating structure-based drug and vaccine design against SARS-CoV-2.

A synergy of activity, stability, and inhibitor-interaction of HIV-1 protease mutants evolved under drug-pressure

A clinically-relevant, drug-resistant mutant of HIV-1 protease (PR), termed Flap+_(I54V) and containing L10I, G48V, I54V and V82A mutations, is known to produce significant changes in the entropy and enthalpy balance of drug-PR interactions, compared to wild-type PR. A similar mutant, Flap+_(I54A) , which evolves from Flap+_(I54V) and contains the single change at residue 54 relative to Flap+_(I54V) , does not. Yet, how Flap+_(I54A) behaves in solution is not known. To understand the molecular basis of V54A evolution, we compared nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy, isothermal titration calorimetry, and enzymatic assay data from four PR proteins: PR (pWT), Flap+_(I54V) , Flap+_(I54A) , and Flap+_(I54) , a control mutant that contains only L10I, G48V and V82A mutations. Our data consistently show that selection to the smaller side chain at residue 54, not only decreases inhibitor affinity, but also restores the catalytic activity.

A Personal History of Using Crystals and Crystallography to Understand Biology and Advanced Drug Discovery

Over the past 60 years, the use of crystals to define structures of complexes using X-ray analysis has contributed to the discovery of new medicines in a very significant way. This has been in understanding not only small-molecule inhibitors of proteins, such as enzymes, but also protein or peptide hormones or growth factors that bind to cell surface receptors. Experimental structures from crystallography have also been exploited in software to allow prediction of structures of important targets based on knowledge of homologues. Crystals and crystallography continue to contribute to drug design and provide a successful example of academia–industry collaboration.

Mexican HIV-1 Protease Sequence Diversity

Protease is one of three enzymes encoded within HIV's pol gene, responsible for the cleavage of viral Gag-Pol polypeptide into mature viral proteins and a target of current anti-retroviral therapy. Protease diversity analysis in Latin America has been lacking in spite of extensive studies of protease-inhibitor resistance mutations. We studied the diversity of 777 Mexican protease sequences and found that all were subtype B except one (CRF02_AG). Phylogenetic analysis suggested the existence of six different clades with geospecific contributions. Thirty-three percent of sites were conserved, 25% had conservative substitutions, and 41% exhibited physicochemical changes. The most conserved regions surrounded the active site, most of the flap domain, and a region between the 60's loop and C-terminal triad. A single sequence exhibited an active site mutation (T26S). Variable sites were mapped to a crystallographic structure, providing further insight into the distribution and functional relevance of variable sites among Mexican isolates.

Chemical Synthesis of an Enzyme Containing an Artificial Catalytic Apparatus

With the goal of investigating electronic aspects of the catalysis of peptide bond hydrolysis, an analogue of HIV-1 protease was designed in which a non-peptide hydroxy-isoquinolinone artificial catalytic apparatus replaced the conserved Asp25–Thr26–Gly27 sequence in each 99-residue polypeptide chain of the homodimeric enzyme molecule. The enzyme analogue was prepared by total chemical synthesis and had detectable catalytic activity on known HIV-1 protease peptide substrates. Compared with uncatalyzed hydrolysis, the analogue enzyme increased the rate of peptide bond hydrolysis by ∼108-fold. Extensions of this unique approach to the study of enzyme catalysis in HIV-1 protease are discussed.

Structural and Functional Aspects of Foamy Virus Protease-Reverse Transcriptase

Reverse transcription describes the process of the transformation of single-stranded RNA into double-stranded DNA via an RNA/DNA duplex intermediate, and is catalyzed by the viral enzyme reverse transcriptase (RT). This event is a pivotal step in the life cycle of all retroviruses. In contrast to orthoretroviruses, the domain structure of the mature RT of foamy viruses is different, i.e., it harbors the protease (PR) domain at its N-terminus, thus being a PR-RT. This structural feature has consequences on PR activation, since the enzyme is monomeric in solution and retroviral PRs are only active as dimers. This review focuses on the structural and functional aspects of simian and prototype foamy virus reverse transcription and reverse transcriptase, as well as special features of reverse transcription that deviate from orthoretroviral processes, e.g., PR activation.

NMR and MD studies combined to elucidate inhibitor and water interactions of HIV-1 protease and their modulations with resistance mutations

Over the last two decades, both the sensitivity of NMR and the time scale of molecular dynamics (MD) simulation have increased tremendously and have advanced the field of protein dynamics. HIV-1 protease has been extensively studied using these two methods, and has presented a framework for cross-evaluation of structural ensembles and internal dynamics by integrating the two methods. Here, we review studies from our laboratories over the last several years, to understand the mechanistic basis of protease drug-resistance mutations and inhibitor responses, using NMR and crystal structure-based parallel MD simulations. Our studies demonstrate that NMR relaxation experiments, together with crystal structures and MD simulations, significantly contributed to the current understanding of structural/dynamic changes due to HIV-1 protease drug resistance mutations.

An NMR strategy to detect conformational differences in a protein complexed with highly analogous inhibitors in solution

This manuscript presents an NMR strategy to investigate conformational differences in protein-inhibitor complexes, when the inhibitors tightly bind to a protein at sub-nanomolar dissociation constants and are highly analogous to each other. Using HIV-1 protease (PR), we previously evaluated amide chemical shift differences, ΔCSPs, of PR bound to darunavir (DRV) compared to PR bound to several DRV analogue inhibitors, to investigate subtle but significant long-distance conformation changes caused by the inhibitor's chemical moiety variation [Khan, S. N., Persons, J. D. Paulsen, J. L., Guerrero, M., Schiffer, C. A., Kurt-Yilmaz, N., and Ishima, R., Biochemistry, (2018), 57, 1652-1662]. However, ΔCSPs are not ideal for investigating subtle PR-inhibitor interface differences because intrinsic differences in the electron shielding of the inhibitors affect protein ΔCSPs. NMR relaxation is also not suitable as it is not sensitive enough to detect small conformational differences in rigid regions among similar PR-inhibitor complexes. Thus, to gain insight into conformational differences at the inhibitor-protein interface, we recorded ¹⁵N-half filtered NOESY spectra of PR bound to two highly analogous inhibitors and assessed NOEs between PR amide protons and inhibitor protons, between PR amide protons and hydroxyl side chains, and between PR amide protons and water protons. We also verified the PR amide-water NOEs using 2D water-NOE/ROE experiments. Differences in water-amide proton NOE peaks, possibly due to amide-protein hydrogen bonds, were observed between subunit A and subunit B, and between the DRV-bound form and an analogous inhibitor-bound form, which may contribute to remote conformational changes.

An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant

The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host-guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.

Structural Biology and the Design of New Therapeutics: From HIV and Cancer to Mycobacterial Infections: A Paper Dedicated to John Kendrew

Interest in applications of protein crystallography to medicine was evident, as the first high-resolution structures emerged in the 50s and 60s. In Cambridge, Max Perutz and John Kendrew sought to understand mutations in sickle cell and other genetic diseases related to hemoglobin, while in Oxford, the group of Dorothy Hodgkin became interested in long-lasting zinc-insulin crystals for treatment of diabetes and later considered insulin redesign, as synthetic insulins became possible. The use of protein crystallography in structure-guided drug discovery emerged as enzyme structures allowed the identification of potential inhibitor-binding sites and optimization of interactions of hits using the structure of the target protein. Early examples of this approach were the use of the structure of renin to design antihypertensives and the structure of HIV protease in design of AIDS antivirals. More recently, use of structure-guided design with fragment-based drug discovery, which reduces the size of screening libraries by decreasing complexity, has improved ligand efficiency in drug design and has been used to progress three oncology drugs through clinical trials to FDA approval. We exemplify current developments in structure-guided target identification and fragment-based lead discovery with efforts to develop new antimicrobials for mycobacterial infections.

Protein crystallography and drug discovery: recollections of knowledge exchange between academia and industry

The development of structure-guided drug discovery is a story of knowledge exchange where new ideas originate from all parts of the research ecosystem. Dorothy Crowfoot Hodgkin obtained insulin from Boots Pure Drug Company in the 1930s and insulin crystallization was optimized in the company Novo in the 1950s, allowing the structure to be determined at Oxford University. The structure of renin was developed in academia, on this occasion in London, in response to a need to develop antihypertensives in pharma. The idea of a dimeric aspartic protease came from an international academic team and was discovered in HIV; it eventually led to new HIV antivirals being developed in industry. Structure-guided fragment-based discovery was developed in large pharma and biotechs, but has been exploited in academia for the development of new inhibitors targeting protein-protein interactions and also antimicrobials to combat mycobacterial infections such as tuberculosis. These observations provide a strong argument against the so-called 'linear model', where ideas flow only in one direction from academic institutions to industry. Structure-guided drug discovery is a story of applications of protein crystallography and knowledge exhange between academia and industry that has led to new drug approvals for cancer and other common medical conditions by the Food and Drug Administration in the USA, as well as hope for the treatment of rare genetic diseases and infectious diseases that are a particular challenge in the developing world.

Altered Plasmodium falciparum Sensitivity to the Antiretroviral Protease Inhibitor Lopinavir Associated with Polymorphisms in pfmdr1

The HIV protease inhibitor lopinavir inhibits Plasmodium falciparum aspartic proteases (plasmepsins) and parasite development, and children receiving lopinavir-ritonavir experienced fewer episodes of malaria than those receiving other antiretroviral regimens. Resistance to lopinavir was selected in vitro over ∼9 months, with ∼4-fold decreased sensitivity. Whole-genome sequencing of resistant parasites showed a mutation and increased copy number in pfmdr1 and a mutation in a protein of unknown function, but no polymorphisms in plasmepsin genes.

Strategies for Antiviral Drug Discovery

The need for antiviral drugs is growing rapidly as more viral diseases are recognized. The methods used to discover these drugs have evolved considerably over the past 40 years and the overall process of discovery can be broken down into sub-processes which include lead generation, lead optimization and lead development. Various methods are now employed to ensure these processes are carried out efficiently. For lead generation, screening methodologies have developed to the extent where hundreds of thousands of compounds can be screened against a particular target. An alternative approach is to use the structures of enzyme substrates as a starting point for drug discovery. Much use is now made of X-ray crystallographic data of target–inhibitor complexes for the optimization of lead structures, and methods for preparing libraries of compounds to assist both generation and optimization of leads are welldeveloped. The methods used to predict and improve the pharmacokinetic properties of compounds are also changing rapidly. Finally, novel approaches to antiviral therapy using oligonucleotide-based compounds or modulating the host immune response are also being explored. This review discusses these approaches, provides examples of where their application has been successful and sets them against a historical background.

Structural biology in antiviral drug discovery

Structural biology has emerged during the last thirty years as a powerful tool for rational drug discovery. Crystal structures of biological targets alone and in complex with ligands and inhibitors provide essential insights into the mechanisms of actions of enzymes, their conformational changes upon ligand binding, the architectures and interactions of binding pockets. Structure-based methods such as crystallographic fragment screening represent nowadays invaluable instruments for the identification of new biologically active compounds. In this context, three-dimensional protein structures have played essential roles for the understanding of the activity and for the design of novel antiviral agents against several different viruses. In this review, the evolution in the resolution of viral structures is analysed, along with the role of crystal structures in the discovery and optimisation of new antivirals.

Effects of Hinge-region Natural Polymorphisms on Human Immunodeficiency Virus-Type 1 Protease Structure, Dynamics, and Drug Pressure Evolution

Multidrug resistance to current Food and Drug Administration-approved HIV-1 protease (PR) inhibitors drives the need to understand the fundamental mechanisms of how drug pressure-selected mutations, which are oftentimes natural polymorphisms, elicit their effect on enzyme function and resistance. Here, the impacts of the hinge-region natural polymorphism at residue 35, glutamate to aspartate (E35D), alone and in conjunction with residue 57, arginine to lysine (R57K), are characterized with the goal of understanding how altered salt bridge interactions between the hinge and flap regions are associated with changes in structure, motional dynamics, conformational sampling, kinetic parameters, and inhibitor affinity. The combined results reveal that the single E35D substitution leads to diminished salt bridge interactions between residues 35 and 57 and gives rise to the stabilization of open-like conformational states with overall increased backbone dynamics. In HIV-1 PR constructs where sites 35 and 57 are both mutated (e.g. E35D and R57K), x-ray structures reveal an altered network of interactions that replace the salt bridge thus stabilizing the structural integrity between the flap and hinge regions. Despite the altered conformational sampling and dynamics when the salt bridge is disrupted, enzyme kinetic parameters and inhibition constants are similar to those obtained for subtype B PR. Results demonstrate that these hinge-region natural polymorphisms, which may arise as drug pressure secondary mutations, alter protein dynamics and the conformational landscape, which are important thermodynamic parameters to consider for development of inhibitors that target for non-subtype B PR.

Effects of a Specific Inhibitor of HIV Proteinase (Ro 31-8959) on Virus Maturation in a Chronically Infected Promonocytic Cell Line (U1)

The human immunodeficiency virus (HIV) proteinase inhibitor Ro 31-8959 prevents the maturation of virus in phorbol 12-myristate 13-acetate (PMA)-stimulated U1 cells, a chronically infected promonocytic cell line. Inhibition of both the morphological maturation of virions and the enzymic processing of gag polyprotein (p56) to produce capsid protein p24 was demonstrated at nanomolar concentrations of the compound. Furthermore, prolonged inhibition of the processing of p56 antigen was confirmed in pulse-chase experiments. The conclusion is that Ro 31-8959 can inhibit production of mature virions in a promonocyte cell line which is infected chronically/latently with HIV.

In vitro Sequential Selection and Characterization of Human Immunodeficiency Virus Type 1 Variants with Reduced Sensitivity to Hydroxyethylurea Protease Inhibitors

In vitro resistance to the human immunodeficiency virus (HIV) protease inhibitors SC-52151 and SC-55389A was evaluated in an in vitro sequential selection scheme. HIVRF variants were selected for reduced sensitivity to SC-52151 and subsequently passaged in both SC-52151 and a structurally different hydroxyethylurea protease inhibitor, SC-55389A, to select for dual-resistant virus. SC-52151 selection alone resulted in a 23-fold reduction in virus sensitivity whereas selection in both inhibitors resulted in 34- and eightfold reductions in virus sensitivity to SC-52151 and SC-55389A, respectively. Sequence analysis of the protease gene revealed that SC-52151 -resistant virus had a Gly to Val substitution at residue 48 (G48V) and, in 58% of subclones, an accompanying Val to Ala substitution at residue 82 (V82A). Dual-resistant virus had both G48V and V82A substitutions present and, in the majority of subclones, an lle to Thr and/or Leu to Pro substitution at residues 54 and 63, respectively. Drug susceptibility assays with limiting dilution-cloned HIVRFR (G48V/V82A) and HIVRFRR (G48V/154T/L63P/V82A) viruses demonstrated moderate to high-level cross-resistance to additional structurally non-related protease inhibitors. Recombinant HIVHXB2 proviral clones with G48V, L63P and V82A substitutions showed that one active site mutation was permissible, but the presence of both G48V and V82A substitutions together significantly reduced infectious virus production. Insight into the contributions of the observed substitutions to drug resistance is presented in molecular modelling studies.

Identification and characterization of a LTR retrotransposon from the genome of Cyprinus carpio var. Jian

A Ty3/gypsy-retrotransposon-type transposon was found in the genome of the Jian carp (Cyprinus carpio var. Jian) in a previous study (unpublished), and was designated a JRE retrotransposon (Jian retrotransposon). The full-length JRE retrotransposon is 5126 bp, which includes two long terminal repeats of 470 bp at the 5' end and 453 bp at the 3' end, and two open reading frames between them: 4203 bp encoding the group-specific antigen (GAG) and polyprotein (POL). The pol gene has a typical Ty3/gypsy retrotransposon structure, and the gene order is protease, reverse transcriptase, RNase H, and integrase (PR-RT-RH-IN). A phylogenetic analysis of the pol gene showed that it has similarities of 40.7, 40, and 32.8 %, to retrotransposons of Azumapecten farreri, Mizuhopecten yessoensis, and Xiphophorus maculatus, respectively. Therefore, JRE might belong to the JULE retrotransposon family. The copy number of the JRE transposon in the genome of the Jian carp is 124, determined with real-time quantitative PCR. The mRNA of the JRE retrotransposon is expressed in five Jian carp tissues, the liver, kidney, blood, muscle, and gonad, and slightly higher in the kidney and liver than in the other tissues.

Elucidation of the structure of retroviral proteases: a reminiscence

Determinations of only a very few protein structures had consequences comparable to the impact exerted by the structure of the protease encoded by HIV-1, published just over 25 years ago. The structure of this relatively small protein and its cousins from other retroviruses provided a clear target for a spectacularly successful structure-assisted drug design effort that offered new hope for controlling the then-escalating AIDS epidemic. This reminiscence is limited primarily to work conducted at the National Cancer Institute, and is not meant to be a comprehensive history of the field, but is rather an attempt to provide a very personal account of how the structures of this most thoroughly studied crystallographic target were determined.

Molecular characterization of HIV-1 genome in fission yeast Schizosaccharomyces pombe

Background

The human immunodeficiency virus type 1 (HIV-1) genome (~9 kb RNA) is flanked by two long terminal repeats (LTR) promoter regions with nine open reading frames, which encode Gag, Pol and Env polyproteins, four accessory proteins (Vpu, Vif, Vpr, Nef) and two regulatory proteins (Rev, Tat). In this study, we carried out a genome-wide and functional analysis of the HIV-1 genome in fission yeast (Schizosaccharomyces pombe).

Results

Each one of the HIV-1 genes was cloned and expressed individually in fission yeast. Subcellular localization of each viral protein was first examined. The effect of protein expression on cellular proliferation and colony formations, an indication of cytotoxicity, were observed. Overall, there is a general correlation of subcellular localization of each viral protein between fission yeast and mammalian cells. Three viral proteins, viral protein R (Vpr), protease (PR) and regulator of expression of viral protein (Rev), were found to inhibit cellular proliferation. Rev was chosen for further analysis in fission yeast and mammalian cells. Consistent with the observation in fission yeast, expression of HIV-1 rev gene also caused growth retardation in mammalian cells. However, the observed growth delay was neither due to the cytotoxic effect nor due to alterations in cell cycling. Mechanistic testing of the Rev effect suggests it triggers transient induction of cellular oxidative stress.

Conclusions

Some of the behavioral and functional similarities of Rev between fission yeast and mammalian cells suggest fission yeast might be a useful model system for further studies of molecular functions of Rev and other HIV-1 viral proteins.

Mutations Proximal to Sites of Autoproteolysis and the α-Helix That Co-evolve under Drug Pressure Modulate the Autoprocessing and Vitality of HIV-1 Protease

N-Terminal self-cleavage (autoprocessing) of the HIV-1 protease precursor is crucial for liberating the active dimer. Under drug pressure, evolving mutations are predicted to modulate autoprocessing, and the reduced catalytic activity of the mature protease (PR) is likely compensated by enhanced conformational/dimer stability and reduced susceptibility to self-degradation (autoproteolysis). One such highly evolved, multidrug resistant protease, PR20, bears 19 mutations contiguous to sites of autoproteolysis in retroviral proteases, namely clusters 1-3 comprising residues 30-37, 60-67, and 88-95, respectively, accounting for 11 of the 19 mutations. By systematically replacing corresponding clusters in PR with those of PR20, and vice versa, we assess their influence on the properties mentioned above and observe no strict correlation. A 10-35-fold decrease in the cleavage efficiency of peptide substrates by PR20, relative to PR, is reflected by an only ∼4-fold decrease in the rate of Gag processing with no change in cleavage order. Importantly, optimal N-terminal autoprocessing requires all 19 PR20 mutations as evaluated in vitro using the model precursor TFR-PR20 in which PR is flanked by the transframe region. Substituting PR20 cluster 3 into TFR-PR (TFR-PR(PR20-3)) requires the presence of PR20 cluster 1 and/or 2 for autoprocessing. In accordance, substituting PR clusters 1 and 2 into TFR-PR20 affects the rate of autoprocessing more drastically (>300-fold) compared to that of TFR-PR(PR20-3) because of the cumulative effect of eight noncluster mutations present in TFR-PR20(PR-12). Overall, these studies imply that drug resistance involves a complex synchronized selection of mutations modulating all of the properties mentioned above governing PR regulation and function.

A Direct Interaction with RNA Dramatically Enhances the Catalytic Activity of the HIV-1 Protease In Vitro

Though the steps of human immunodeficiency virus type 1 (HIV-1) virion maturation are well documented, the mechanisms regulating the proteolysis of the Gag and Gag-Pro-Pol polyproteins by the HIV-1 protease (PR) remain obscure. One proposed mechanism argues that the maturation intermediate p15NC must interact with RNA for efficient cleavage by the PR. We investigated this phenomenon and found that processing of multiple substrates by the HIV-1 PR was enhanced in the presence of RNA. The acceleration of proteolysis occurred independently from the substrate's ability to interact with nucleic acid, indicating that a direct interaction between substrate and RNA is not necessary for enhancement. Gel-shift assays demonstrated the HIV-1 PR is capable of interacting with nucleic acids, suggesting that RNA accelerates processing reactions by interacting with the PR rather than the substrate. All HIV-1 PRs examined have this ability; however, the HIV-2 PR does not interact with RNA and does not exhibit enhanced catalytic activity in the presence of RNA. No specific sequence or structure was required in the RNA for a productive interaction with the HIV-1 PR, which appears to be principally, though not exclusively, driven by electrostatic forces. For a peptide substrate, RNA increased the kinetic efficiency of the HIV-1 PR by an order of magnitude, affecting both turnover rate (k(cat)) and substrate affinity (K(m)). These results suggest that an allosteric binding site exists on the HIV-1 PR and that HIV-1 PR activity during maturation could be regulated in part by the juxtaposition of the enzyme with virion-packaged RNA.

The triple threat of HIV-1 protease inhibitors

Newly released human immunodeficiency virus type 1 (HIV-1) particles obligatorily undergo a maturation process to become infectious. The HIV-1 protease (PR) initiates this step, catalyzing the cleavage of the Gag and Gag-Pro-Pol structural polyproteins. Proper organization of the mature virus core requires that cleavage of these polyprotein substrates proceeds in a highly regulated, specific series of events. The vital role the HIV-1 PR plays in the viral life cycle has made it an extremely attractive target for inhibition and has accordingly fostered the development of a number of highly potent substrate-analog inhibitors. Though the PR inhibitors (PIs) inhibit only the HIV-1 PR, their effects manifest at multiple different stages in the life cycle due to the critical importance of the PR in preparing the virus for these subsequent events. Effectively, PIs masquerade as entry inhibitors, reverse transcription inhibitors, and potentially even inhibitors of post-reverse transcription steps. In this chapter, we review the triple threat of PIs: the intermolecular cooperativity in the form of a cooperative dose-response for inhibition in which the apparent potency increases with increasing inhibition; the pleiotropic effects of HIV-1 PR inhibition on entry, reverse transcription, and post-reverse transcription steps; and their potency as transition state analogs that have the potential for further improvement that could lead to an inability of the virus to evolve resistance in the context of single drug therapy.

Viruses and viral proteins

For more than 30 years X-ray crystallography has been by far the most powerful approach for determining the structures of viruses and viral proteins at atomic resolution. The information provided by these structures, which covers many important aspects of the viral life cycle such as cell-receptor recognition, viral entry, nucleic acid transfer and genome replication, has extensively enriched our vision of the virus world. Many of the structures available correspond to potential targets for antiviral drugs against important human pathogens. This article provides an overview of the current knowledge of different structural aspects of the above-mentioned processes.

Recent patents and emerging therapeutics for HIV infections: a focus on protease inhibitors

The inclusion of protease inhibitors (PIs) in highly active antiretroviral therapy has significantly improved clinical outcomes in HIV-1-infected patients. To date, PIs are considered to be the most important therapeutic agents for the treatment of HIV infections. Despite high anti-HIV-1 potency, poor oral bioavailability of PIs has been a major concern. For achieving therapeutic concentrations, large doses of PIs are administered, which results in unacceptable systemic toxicities. Such severe and long-term toxicities necessitate the development of safer and potentially promising PIs. Recently, considerable attention has been paid to the development of newer compounds capable of inhibiting wild-type and resistant HIV-1 protease. Some of these PIs have displayed potent HIV-1 protease inhibitory activity. In this review, we have made an attempt to provide an overview on clinically approved and newly developing PIs, and related recent patents in the development of novel PIs.

Investigation on the mechanism for the binding and drug resistance of wild type and mutations of G86 residue in HIV-1 protease complexed with Darunavir by molecular dynamic simulation and free energy calculation

Residue Gly86 is considered as the highly conversed residue in the HIV-1 protease. In our work, the detailed binding free energies for the wild-type (WT) and mutated proteases binding to the TMC-114 are estimated to investigate the protein-inhibitor binding and drug resistance mechanism by molecule dynamic simulations and molecular mechanics Poisson Boltzmann surface area (MM-PBSA) method. The binding affinities between the mutants and inhibitor are different than that in the wild-type complex and the major resistance to Darunavir (DRV) of G86A and G86S originate from the electrostatic energy and entropy, respectively. Furthermore, free energy decomposition analysis for the WT and mutated complexes on the basis of per-residue indicates that the mutagenesis influences the energy contribution of the residue located at three regions: active site region (residue 24-32), the flap region, and the region around the mutated residue G86 (residue 79-88), especially the flap region. Finally, further hydrogen bonds and structure analysis are carried out to detect the relationship between the energy and conformation. In all, the G86 mutations change the flap region's conformation. The experimental results are in good agreement with available results.

Understanding HIV-1 protease autoprocessing for novel therapeutic development

In the infected cell, HIV-1 protease (PR) is initially synthesized as part of the GagPol polyprotein. PR autoprocessing is a virus-specific process by which the PR domain embedded in the precursor catalyzes proteolytic reactions responsible for liberation of free mature PRs, which then recognize and cleave at least ten different peptide sequences in the Gag and GagPol polyproteins. Despite extensive structure and function studies of the mature PRs as well as the successful development of ten US FDA-approved catalytic-site inhibitors, the precursor autoprocessing mechanism remains an intriguing yet-to-be-solved puzzle. This article discusses current understanding of the autoprocessing mechanism, in an effort to prompt the development of novel anti-HIV drugs that selectively target precursor autoprocessing.

Identification of the intracellular gate for a member of the equilibrative nucleoside transporter (ENT) family

Equilibrative nucleoside transporters of the SLC29 family play important roles in many physiological and pharmacological processes, including import of drugs for treatment of cancer, AIDS, cardiovascular, and parasitic diseases. However, no crystal structure is available for any member of this family. In previous studies we generated a computational model of the Leishmania donovani nucleoside transporter 1.1 (LdNT1.1) that captured this permease in the outward-closed conformation, and we identified the extracellular gate. In the present study we have modeled the inward-closed conformation of LdNT1.1 using the crystal structure of the Escherichia coli fucose transporter FucP and have identified four transmembrane helices whose ends close to form a predicted intracellular gate. We have tested this prediction by site-directed mutagenesis of relevant helix residues and by cross-linking of introduced cysteine pairs. The results are consistent with the predictions of the computational model and suggest that a similarly constituted gate operates in other members of the equilibrative nucleoside transporter family.

Nuclear localization signals in prototype foamy viral integrase for successive infection and replication in dividing cells

We identified four basic amino acid residues as nuclear localization signals (NLS) in the C-terminal domain of the prototype foamy viral (PFV) integrase (IN) protein that were essential for viral replication. We constructed seven point mutants in the C-terminal domain by changing the lysine and arginine at residues 305, 308, 313, 315, 318, 324, and 329 to threonine or proline, respectively, to identify residues conferring NLS activity. Our results showed that mutation of these residues had no effect on expression assembly, release of viral particles, or in vitro recombinant IN enzymatic activity. However, mutations at residues 305 (R → T), 313(R → T), 315(R → P), and 329(R → T) lead to the production of defective viral particles with loss of infectivity, whereas non-defective mutations at residues 308(R → T), 318(K → T), and 324(K → T) did not show any adverse effects on subsequent production or release of viral particles. Sub-cellular fractionation and immunostaining for viral protein PFV-IN and PFV-Gag localization revealed predominant cytoplasmic localization of PFV-IN in defective mutants, whereas cytoplasmic and nuclear localization of PFV-IN was observed in wild type and non-defective mutants. However sub-cellular localization of PFV-Gag resulted in predominant nuclear localization and less presence in the cytoplasm of the wild type and non-defective mutants. But defective mutants showed only nuclear localization of Gag. Therefore, we postulate that four basic arginine residues at 305, 313, 315 and 329 confer the karyoplilic properties of PFV-IN and are essential for successful viral integration and replication.

Enhanced stability of monomer fold correlates with extreme drug resistance of HIV-1 protease

During treatment, mutations in HIV-1 protease (PR) are selected rapidly that confer resistance by decreasing affinity to clinical protease inhibitors (PIs). As these unique drug resistance mutations can compromise the fitness of the virus to replicate, mutations that restore conformational stability and activity while retaining drug resistance are selected on further evolution. Here we identify several compensating mechanisms by which an extreme drug-resistant mutant bearing 20 mutations (PR20) with >5-fold increased Kd and >4000-fold decreased affinity to the PI darunavir functions. (1) PR20 cleaves, albeit poorly, Gag polyprotein substrates essential for viral maturation. (2) PR20 dimer, which exhibits distinctly enhanced thermal stability, has highly attenuated autoproteolysis, thus likely prolonging its lifetime in vivo. (3) The enhanced stability of PR20 results from stabilization of the monomer fold. Both monomeric PR20(T26A) and dimeric PR20 exhibit Tm values 6-7.5 °C higher than those for their PR counterparts. Two specific mutations in PR20, L33F and L63P at sites of autoproteolysis, increase the Tm of monomeric PR(T26A) by ~8 °C, similar to PR20(T26A). However, without other compensatory mutations as seen in PR20, L33F and L63P substitutions, together, neither restrict autoproteolysis nor significantly reduce binding affinity to darunavir. To determine whether dimer stability contributes to binding affinity for inhibitors, we examined single-chain dimers of PR and PR(D25N) in which the corresponding identical monomer units were covalently linked by GGSSG sequence. Linking of the subunits did not appreciably change the ΔTm on inhibitor binding; thus stabilization by tethering appears to have little direct effect on enhancing inhibitor affinity.

Whiskers-less HIV-protease: a possible way for HIV-1 deactivation

BACKGROUND

Among viral enzymes, the human HIV-1 protease comprises the most interesting target for drug discovery. There are increasing efforts focused on designing more effective inhibitors for HIV-1 protease in order to prevent viral replication in AIDS patients. The frequent and continuous mutation of HIV-1 protease gene creates a formidable obstacle for enzyme inhibition which could not be overcome by the traditional single drug therapy. Nowadays, in vitro and in silico studies of protease inhibition constitute an advanced field in biological researches. In this article, we tried to simulate protease-substrate complexes in different states; a native state and states with whiskers deleted from one and two subunits. Molecular dynamic simulations were carried out in a cubic box filled with explicit water at 37°C and in 1atomsphere of pressure.

RESULTS

Our results showed that whisker truncation of protease subunits causes the dimer structure to decrease in compactness, disrupts substrate-binding site interactions and changes in flap status simultaneously.

CONCLUSIONS

Based on our findings we claim that whisker truncation even when applied to a single subunit, threats dimer association which probably leads to enzyme inactivation. We may postulate that inserting a gene to express truncated protease inside infected cells can interfere with protease dimerization. The resulted proteases would presumably have a combination of native and truncated subunits in their structures which exert no enzyme activities as evidenced by the present work. Our finding may create a new field of research in HIV gene therapy for protease inhibition, circumventing problems of drug resistance.

Crystal structure of a putative aspartic proteinase domain of the Mycobacterium tuberculosis cell surface antigen PE_PGRS16

We report the crystal structure of the first prokaryotic aspartic proteinase-like domain identified in the genome of Mycobacterium tuberculosis. A search in the genomes of Mycobacterium species showed that the C-terminal domains of some of the PE family proteins contain two classic DT/SG motifs of aspartic proteinases with a low overall sequence similarity to HIV proteinase. The three-dimensional structure of one of them, Rv0977 (PE_PGRS16) of M. tuberculosis revealed the characteristic pepsin-fold and catalytic site architecture. However, the active site was completely blocked by the N-terminal His-tag. Surprisingly, the enzyme was found to be inactive even after the removal of the N-terminal His-tag. A comparison of the structure with pepsins showed significant differences in the critical substrate binding residues and in the flap tyrosine conformation that could contribute to the lack of proteolytic activity of Rv0977.

Foamy virus assembly with emphasis on pol encapsidation

Foamy viruses (FVs) differ from all other genera of retroviruses (orthoretroviruses) in many aspects of viral replication. In this review, we discuss FV assembly, with special emphasis on Pol incorporation. FV assembly takes place intracellularly, near the pericentriolar region, at a site similar to that used by betaretroviruses. The regions of Gag, Pol and genomic RNA required for viral assembly are described. In contrast to orthoretroviral Pol, which is synthesized as a Gag-Pol fusion protein and packaged through Gag-Gag interactions, FV Pol is synthesized from a spliced mRNA lacking all Gag sequences. Thus, encapsidation of FV Pol requires a different mechanism. We detail how WT Pol lacking Gag sequences is incorporated into virus particles. In addition, a mutant in which Pol is expressed as an orthoretroviral-like Gag-Pol fusion protein is discussed. We also discuss temporal regulation of the protease, reverse transcriptase and integrase activities of WT FV Pol.

Effects of HIV-1 protease on cellular functions and their potential applications in antiretroviral therapy

Human Immunodeficiency Virus Type 1 (HIV-1) protease inhibitors (PIs) are the most potent class of drugs in antiretroviral therapies. However, viral drug resistance to PIs could emerge rapidly thus reducing the effectiveness of those drugs. Of note, all current FDA-approved PIs are competitive inhibitors, i.e., inhibitors that compete with substrates for the active enzymatic site. This common inhibitory approach increases the likelihood of developing drug resistant HIV-1 strains that are resistant to many or all current PIs. Hence, new PIs that move away from the current target of the active enzymatic site are needed. Specifically, allosteric inhibitors, inhibitors that prohibit PR enzymatic activities through non-competitive binding to PR, should be sought. Another common feature of current PIs is they were all developed based on the structure-based design. Drugs derived from a structure-based strategy may generate target specific and potent inhibitors. However, this type of drug design can only target one site at a time and drugs discovered by this method are often associated with strong side effects such as cellular toxicity, limiting its number of target choices, efficacy, and applicability. In contrast, a cell-based system may provide a useful alternative strategy that can overcome many of the inherited shortcomings associated with structure-based drug designs. For example, allosteric PIs can be sought using a cell-based system without considering the site or mechanism of inhibition. In addition, a cell-based system can eliminate those PIs that have strong cytotoxic effect. Most importantly, a simple, economical, and easy-to-maintained eukaryotic cellular system such as yeast will allow us to search for potential PIs in a large-scaled high throughput screening (HTS) system, thus increasing the chances of success. Based on our many years of experience in using fission yeast as a model system to study HIV-1 Vpr, we propose the use of fission yeast as a possible surrogate system to study the effects of HIV-1 protease on cellular functions and to explore its utility as a HTS system to search for new PIs to battle HIV-1 resistant strains.

Single point mutation induced alterations in the equilibrium structural transitions on the folding landscape of HIV-1 protease

Equilibrium folding-unfolding transitions are hard to study in HIV-1 protease (PR) because of its autolytic properties. Further, the protease exhibits many tolerant point mutations some of which also impart drug resistance to the protein. It is conceivable that the mutations affect protein's function by altering its folding characteristics; these would clearly depend on the nature of the mutations themselves. In this background, we report here NMR studies on the effects of D25 N mutation, which removes one negative charge from the protein at the active site, on the equilibrium folding behaviour of PR starting from its acetic acid denatured state. It is observed that in PRD25N two slowly exchanging conformations are present at the N-terminal. One of them is similar to that of PR. Though the conformational and dynamics preferences of PR and PRD25N are fairly similar in 9 M acetic acid, they seem to undergo different folding transitions when acetic acid concentration is reduced. The differences are seen in the active site, in the flap, and in the hinge of the flap regions. The present study suggests that such differences, though different in detail, would occur for other mutations as well, and also for different initial denatured states. These would have significant regulatory implications for the efficacy of protease function.

Context surrounding processing sites is crucial in determining cleavage rate of a subset of processing sites in HIV-1 Gag and Gag-Pro-Pol polyprotein precursors by viral protease

Processing of the human immunodeficiency virus type 1 (HIV-1) Gag and Gag-Pro-Pol polyproteins by the HIV-1 protease (PR) is essential for the production of infectious particles. However, the determinants governing the rates of processing of these substrates are not clearly understood. We studied the effect of substrate context on processing by utilizing a novel protease assay in which a substrate containing HIV-1 matrix (MA) and the N-terminal domain of capsid (CA) is labeled with a FlAsH (fluorescein arsenical hairpin) reagent. When the seven cleavage sites within the Gag and Gag-Pro-Pol polyproteins were placed at the MA/CA site, the rates of cleavage changed dramatically compared with that of the cognate sites in the natural context reported previously. The rate of processing was affected the most for three sites: CA/spacer peptide 1 (SP1) (≈10-fold increase), SP1/nucleocapsid (NC) (≈10-30-fold decrease), and SP2/p6 (≈30-fold decrease). One of two multidrug-resistant (MDR) PR variants altered the pattern of processing rates significantly. Cleavage sites within the Pro-Pol region were cleaved in a context-independent manner, suggesting for these sites that the sequence itself was the determinant of rate. In addition, a chimera consisting of SP1/NC P4-P1 and MA/CA P1'-P4' residues (ATIM↓PIVQ) abolished processing by wild type and MDR proteases, and the reciprocal chimera consisting of MA/CA P4-P1 and SP1/NC P1'-4' (SQNY↓IQKG) was cleaved only by one of the MDR proteases. These results suggest that complex substrate interactions both beyond the active site of the enzyme and across the scissile bond contribute to defining the rate of processing by the HIV-1 PR.

Biomolecularmodeling and simulation: a field coming of age

We assess the progress in biomolecular modeling and simulation, focusing on structure prediction and dynamics, by presenting the field’s history, metrics for its rise in popularity, early expressed expectations, and current significant applications. The increases in computational power combined with improvements in algorithms and force fields have led to considerable success, especially in protein folding, specificity of ligand/biomolecule interactions, and interpretation of complex experimental phenomena (e.g. NMR relaxation, protein-folding kinetics and multiple conformational states) through the generation of structural hypotheses and pathway mechanisms. Although far from a general automated tool, structure prediction is notable for proteins and RNA that preceded the experiment, especially by knowledge-based approaches. Thus, despite early unrealistic expectations and the realization that computer technology alone will not quickly bridge the gap between experimental and theoretical time frames, ongoing improvements to enhance the accuracy and scope of modeling and simulation are propelling the field onto a productive trajectory to become full partner with experiment and a field on its own right.

Molecular analysis of the complete genomic sequences of four isolates of Gooseberry vein banding associated virus

The presence of Gooseberry vein banding associated virus (GVBaV), a badnavirus in the family Caulimoviridae, is strongly correlated with gooseberry vein banding disease in Ribes spp. In this study, full-length genomic sequences of four GVBaV isolates from different hosts and geographic regions were determined to be 7649-7663 nucleotides. These isolates share identities of 96.4-97.3% for the complete genomic sequence, indicating low genetic diversity among them. The GVBaV genome contains three open reading frames (ORFs) on the plus strand that potentially encode proteins of 26, 16, and 216 kDa. The size and organization of GVBaV ORFs 1-3 are similar to those of most other badnaviruses. The putative amino acid sequence of GVBaV ORF 3 contained motifs that are conserved among badnavirus proteins including aspartic protease, reverse transcriptase, and ribonuclease H. The highly conserved putative plant tRNA(met)-binding site is also present in the 935-bp intergenic region of GVBaV. The identities of the genomic sequences of GVBaV and other badnaviruses range from 49.1% (Sugarcane bacilliform Mor virus) to 51.7% (Pelargonium vein banding virus, PVBV). Phylogenetic analysis using the amino acid sequence of the ORF 3 putative protein shows that GVBaV groups most closely to Dioscorea bacilliform virus, PVBV, and Taro bacilliform virus. These results confirm that GVBaV is a pararetrovirus of the genus Badnavirus in the family Caulimoviridae.

Protein functional sites prediction using modified bio-basis function and quantitative indices

The prediction of functional sites in proteins is an important issue in protein function studies and drug design. To apply the kernel based pattern recognition algorithms such as support vector machines for protein functional sites prediction, a new string kernel function, termed as the modified bio-basis function, is proposed recently. The bio-basis strings for the new kernel function are selected by an efficient method that integrates the Fisher ratio and the concept of degree of resemblance. In this regard, this paper introduces some quantitative indices for evaluating the quality of selected bio-basis strings. Moreover, the effectiveness of the new string kernel function and bio-basis string selection method, along with a comparison with existing bio-basis function and related bio-basis string selection methods, is demonstrated on different protein data sets using the proposed quantitative indices and support vector machines.