Automated generation of MCSS-derived pharmacophoric DOCK site points for searching multiconformation databases

All docking methods employ some sort of heuristic to orient the ligand molecules into the binding site of the target structure. An automated method, MCSS2SPTS, for generating chemically labeled site points for docking is presented. MCSS2SPTS employs the program Multiple Copy Simultaneous Search (MCSS) to determine target-based theoretical pharmacophores. More specifically, chemically labeled site points are automatically extracted from selected low-energy functional-group minima and clustered together. These pharmacophoric site points can then be directly matched to the pharmacophoric features of database molecules with the use of either DOCK or PhDOCK to place the small molecules into the binding site. Several examples of the ability of MCSS2SPTS to reproduce the three-dimensional pharmacophoric features of ligands from known ligand-protein complex structures are discussed. In addition, a site-point set calculated for one human immunodeficiency virus 1 (HIV1) protease structure is used with PhDOCK to dock a set of HIV1 protease ligands; the docked poses are compared to the corresponding complex structures of the ligands. Finally, the use of an MCSS2SPTS-derived site-point set for acyl carrier protein synthase is compared to the use of atomic positions from a bound ligand as site points for a large-scale DOCK search. In general, MCSS2SPTS-generated site points focus the search on the more relevant areas and thereby allow for more effective sampling of the target site.     

MCSS functionality maps for a flexible protein

The Multiple Copy Simultaneous Search (MCSS) methodology for finding energetically favorable positions and orientations of small functional groups in a binding site is extended to include flexibility of the target. This makes possible the finding of novel minima not present in a fixed structure and so extends the diversity of inhibitors that can be constructed starting with the MCSS procedure. Quenched molecular dynamics is used to generate energetically favorable positions and orientations of the functional groups in the field of a flexible protein. The method is applied to the viral protein HIV-1 protease with methanol and methyl ammonium as a test case. If the protein is quenched with many copies of functional groups randomly distributed in the binding site, the resulting minima have ligand-protein interaction energies that are, on average, less favorable than those obtained with standard MCSS. This is a consequence of the renormalized potential function employed in the Locally Enhanced Sampling (LES) approximation. However, local optimizations of existing MCSS minima with a flexible protein results in lower energy minima in regions of the protein that are of particular interest. Their use in constructing a consensus protein model for ligand design is discussed.     

Molecular mechanisms of attenuation of the Sabin strain of poliovirus type 3

Mutations critical for the central nervous system (CNS) attenuation of the Sabin vaccine strains of poliovirus (PV) are located within the viral internal ribosome entry site (IRES). We examined the interaction of the IRESs of PV type 3 (PV3) and Sabin type 3 (Sabin3) with polypyrimidine tract-binding protein (PTB) and a neural cell-specific homologue, nPTB. PTB and nPTB were found to bind to a site directly adjacent to the attenuating mutation, and binding at this site was less efficient on the Sabin3 IRES than on the PV3 IRES. Translation mediated by the PV3 and Sabin3 IRESs in neurons of the chicken embryo spinal cord demonstrated a translation deficit for the Sabin3 IRES that could be rescued by increasing PTB expression in the CNS. These data suggest that the low levels of PTB available in the CNS, coupled to a reduced binding of PTB on the Sabin3 IRES, leads to its CNS-specific attenuation. This study also demonstrates the use of the chicken embryo to easily investigate translation of RNA within a neuron in the CNS of an intact living organism.     

Dynamic ligand design and combinatorial optimization: designing inhibitors to endothiapepsin

The dynamic ligand design (DLD) algorithm, an automated method for the creation of novel ligands, links up small functional groups that have been placed in energetically favorable positions in the binding site of a target molecule. The positions and orientations of the small functional groups can be determined using the multi-copy simultaneous search approach (MCSS) or experimental data. In this work the original DLD methodology is extended by using a modified version of the pseudo-potential energy function. A novel simulated annealing protocol is presented for optimizing the pseudo-potential energy of ligands in the binding site; the protocol is expected to be applicable to other optimization problems. The utility of the method is illustrated by designing an inhibitor for endothiapepsin. The binding affinity of the inhibitor is assessed using a thermodynamic cycle that decomposes the binding free energy into a sum of translational, rotational, configurational, hydrophobic, and electrostatic contributions. The calculations suggest that the designed molecule will bind endothiapepsin with high affinity.     

Ligand binding and activation of UTP-activated G protein-coupled P2Y2 and P2Y4 receptors elucidated by mutagenesis,pharmacological and computational studies

The nucleotide receptors P2Y₂ and P2Y₄ are the most closely related G protein-coupled receptors (GPCRs) of the P2Y receptor (P2YR) family. Both subtypes couple to G_q proteins and are activated by the pyrimidine nucleotide UTP, but only P2Y₂R is also activated by the purine nucleotide ATP. Agonists and antagonists of both receptor subtypes have potential as drugs e.g. for neurodegenerative and inflammatory diseases. So far, potent and selective, “drug-like” ligands for both receptors are scarce, but would be required for target validation and as lead structures for drug development. Structural information on the receptors is lacking since no X-ray structures or cryo-electron microscopy images are available. Thus, we performed receptor homology modeling and docking studies combined with mutagenesis experiments on both receptors to address the question how ligand binding selectivity for these closely related P2YR subtypes can be achieved. The orthosteric binding site of P2Y₂R appeared to be more spacious than that of P2Y₄R. Mutation of Y197 to alanine in P2Y₄R resulted in a gain of ATP sensitivity. Anthraquinone-derived antagonists are likely to bind to the orthosteric or an allosteric site depending on their substitution pattern and the nature of the orthosteric binding site of the respective P2YR subtype. These insights into the architecture of P2Y₂- and P2Y₄Rs and their interactions with structurally diverse agonists and antagonist provide a solid basis for the future design of potent and selective ligands.     

Utilization of a mammalian cell-based RNA binding assay to characterize the RNA binding properties of picornavirus 3C proteinases.

Using an assay capable of detecting sequence-specific RNA/protein interactions in mammalian cells, we demonstrate that the poliovirus and rhinovirus 3C proteinases are able to bind structured target RNA sequences derived from their respective 5' noncoding regions in vivo. Specific RNA binding by poliovirus 3C was found to be dependent on the integrity of stem-loop d of the RNA cloverleaf structure located at the 5' end of poliovirus genomic RNA. In contrast, mutation of stem-loop b did not prevent this in vivo interaction. However, mutation of stem-loop b, which serves as the RNA binding site for a cellular co-factor important for efficient poliovirus replication, did significantly attenuate the efficiency of 3C RNA binding in vivo and 3CD RNA binding in vitro. This in vivo protein:RNA binding assay was also used to identify several residues in 3C that are critical for RNA binding, but dispensable for 3C proteinase activity. The mammalian cell-based RNA binding assay described in this study may have considerable potential utility in the future detection or analysis of in vivo RNA/protein interactions unrelated to the 3C/RNA interaction described here.     

Genetic determinants of cell type-specific poliovirus propagation in HEK 293 cells

The ability of poliovirus to propagate in neuronal cells can be reduced by introducing appropriate nucleotide substitutions into the viral genome. Specific mutations scattered throughout the poliovirus genome yielded the live attenuated vaccine strains of poliovirus. Neuron-specific propagation deficits of the Sabin strains are partially encrypted within a confined region of the internal ribosomal entry site (IRES), which carries attenuating point mutations in all three serotypes. Recently, high levels of neurovirulence attenuation were achieved with genetically engineered polioviruses containing heterologous IRES elements. This is exemplified with poliovirus recombinants replicating under control of a human rhinovirus type 2 (HRV2) IRES element. We have carried out experiments delineating the genetic basis for neuronal IRES function. Neuronal dysfunction of the HRV2 IRES is determined mainly by IRES stem-loop domain V, the locus for attenuating point mutations within the Sabin strains. Neuronal incompetence associated with HRV2 IRES domain V is substantially more pronounced than that observed with the attenuating IRES point mutation of the Sabin serotype 1 vaccine strain. Mix-and-match recombination of polio and HRV2 IRES domain V suggests that the attenuation phenotype correlates with overall structural features rather than primary sequence. Our experiments have identified HEK 293 cells as a novel system for the study of neuron-specific replication phenotypes of poliovirus. This cell line, originally derived from embryonic human kidney, has recently been described to display neuronal characteristics. We report propagation properties in HEK 293 cells for poliovirus recombinants with attenuated neurovirulence in experimental animals that corroborate this observation.     

Determinants of attenuation and temperature sensitivity in the type 1 poliovirus Sabin vaccine.

To identify determinants of attenuation in the poliovirus type 1 Sabin vaccine strain, a series of recombinant viruses were constructed by using infectious cDNA clones of the virulent type 1 poliovirus P1/Mahoney and the attenuated type 1 vaccine strain P1/Sabin. Intracerebral inoculation of these viruses into transgenic mice which express the human receptor for poliovirus identified regions of the genome that conferred reduced neurovirulence. Exchange of smaller restriction fragments and site-directed mutagenesis were used to identify the nucleotide changes responsible for attenuation. P1/Sabin mutations at nucleotides 935 of VP4, 2438 of VP3, and 2795 and 2879 of VP1 were all shown to be determinants of attenuation. The recombinant viruses and site-directed mutants were also used to identify the nucleotide changes which are involved in the temperature sensitivity of P1/Sabin. Determinants of this phenotype in HeLa cells were mapped to changes at nucleotides 935 of VP4, 2438 of VP3, and 2741 of VP1. The 3Dpol gene of P1/Sabin, which contains three amino acid differences from its parent P1/Mahoney, also contributes to the temperature sensitivity of P1/Sabin; however, mutants containing individual amino acid changes grew as well as P1/Mahoney at elevated temperatures, suggesting that either some combination or all three changes are required for temperature sensitivity. In addition, the 3'-noncoding region of P1/Sabin augments the temperature-sensitive phenotype conferred by 3Dpol. Although nucleotide 2741, 3Dpol, and the 3'-noncoding region of P1/Sabin contribute to the temperature sensitivity of P1/Sabin, they do not contribute to attenuation in transgenic mice expressing the poliovirus receptor, demonstrating that determinants of attenuation and temperature sensitivity can be genetically separated.     

Identifying Ligand Binding Conformations of the β2-Adrenergic Receptor by Using Its Agonists as Computational Probes

Recently available G-protein coupled receptor (GPCR) structures and biophysical studies suggest that the difference between the effects of various agonists and antagonists cannot be explained by single structures alone, but rather that the conformational ensembles of the proteins need to be considered. Here we use an elastic network model-guided molecular dynamics simulation protocol to generate an ensemble of conformers of a prototypical GPCR, β₂-adrenergic receptor (β₂AR). The resulting conformers are clustered into groups based on the conformations of the ligand binding site, and distinct conformers from each group are assessed for their binding to known agonists of β₂AR. We show that the select ligands bind preferentially to different predicted conformers of β₂AR, and identify a role of β₂AR extracellular region as an allosteric binding site for larger drugs such as salmeterol. Thus, drugs and ligands can be used as “computational probes” to systematically identify protein conformers with likely biological significance.     

Improved chromatographic performance of a modified human albumin based stationary phase

Derivatization of the free cys₃₄ in human serum albumin (HSA) anchored to a silica matrix has been performed by in situ reaction with ethacrynic acid. This modification, which is reported to occur under physiological conditions, gives rise in practice to a new column with different binding properties with respect to the column based on the native protein. Significant differences were observed in the binding of drugs known to bind to site I, (R)-(S)-warfarin and phenylbutazone, and to site II, 1,4-benzodiazepin-2-ones and nonsteroidal anti-inflammatory agents. In particular, the chromatographic retentions markedly decreased for most of the drugs, and, in the case of chiral compounds, significant differences were often observed in the behavior of the two enantiomers, with higher values of enantioselectivity obtained for some of the examined compounds. Furthermore, the noncovalent binding of ethacrynic acid to the protein modifies the binding properties of the albumin. Chirality 9:335–340, 1997. © 1997 Wiley-Liss, Inc.     

Natural type 3/type 2 intertypic vaccine-related poliovirus recombinants with the first crossover sites within the VP1 capsid coding region

Background

Ten uncommon natural type 3/type 2 intertypic poliovirus recombinants were isolated from stool specimens from nine acute flaccid paralysis case patients and one healthy vaccinee in China from 2001 to 2008.

Principal Findings

Complete genomic sequences revealed their vaccine-related genomic features and showed that their first crossover sites were randomly distributed in the 3′ end of the VP1 coding region. The length of donor Sabin 2 sequences ranged from 55 to 136 nucleotides, which is the longest donor sequence reported in the literature for this type of poliovirus recombination. The recombination resulted in the introduction of Sabin 2 neutralizing antigenic site 3a (NAg3a) into a Sabin 3 genomic background in the VP1 coding region, which may have been altered by some of the type 3-specific antigenic properties, but had not acquired any type 2-specific characterizations. NAg3a of the Sabin 3 strain seems atypical; other wild-type poliovirus isolates that have circulated in recent years have sequences of NAg3a more like the Sabin 2 strain.

Conclusions

10 natural type 3/type 2 intertypic VP1 capsid-recombinant polioviruses, in which the first crossover sites were found to be in the VP1 coding region, were isolated and characterized. In spite of the complete replacement of NAg3a by type 2-specific amino acids, the serotypes of the recombinants were not altered, and they were totally neutralized by polyclonal type 3 antisera but not at all by type 2 antisera. It is possible that recent type 3 wild poliovirus isolates may be a recombinant having NAg3a sequences derived from another strain during between 1967 and 1980, and the type 3/type 2 recombination events in the 3′ end of the VP1 coding region may result in a higher fitness.     

Synergistic effects in the designs of neuraminidase ligands: Analysis from docking and molecular dynamics studies

Docking and molecular dynamics were used to study the nine ligands (see Scheme 1) at the neuraminidase (NA) active sites. Their binding modes are structurally and energetically different, with details given in the text. Compared with 1A (oseltamivir carboxylate), the changes of core template or/and functional groups in the other ligands cause the reductions of interaction energies and numbers of H-bonds with the NA proteins. Nonetheless, all these ligands occupy the proximity space at the NA active sites and share some commonness in their binding modes. The fragment approach was then used to analyze and understand the binding specificities of the nine ligands. The contributions of each core template and functional group were evaluated. It was found that the core templates rather than functional groups play a larger role during the binding processes; in addition, the binding qualities are determined by the synergistic effects of the core templates and functional groups. Among the nine ligands, 1A (oseltamivir carboxylate) has the largest synergistic energy and its functional groups fit perfectly with the NA active site, consistent with the largest interaction energy, numerous H-bonds with the NA active-site residues as well as experimentally lowest IC₅₀ value. Owing to the poorer metabolizability than oseltamivir, large contribution of the benzene core template and fine synergistic effects of the functional groups, the 4-(N-acetylamino)-5-guanidino-3-(3-pentyloxy)benzoic acid should be an ideal lead compound for optimizing NA drugs.     

Impaired binding of standard initiation factors mediates poliovirus translation attenuation

In the oral poliovirus vaccine, three attenuated virus strains generated by Albert Sabin are used. However, insufficient genetic stability of these strains causes major problems in poliovirus eradication. In infected cells, translation of the plus-strand poliovirus RNA genome is directed by the internal ribosome entry site (IRES), a cis-acting RNA element that facilitates the cap-independent binding of ribosomes to an internal site of the viral RNA. In each Sabin vaccine strain, a single point mutation in the IRES secondary-structure domain V is a major determinant of neurovirulence attenuation. Here we report how these decisive mutations in the IRES confer a reduction in poliovirus translation efficiency. These single-nucleotide exchanges impair the interaction of the standard translation initiation factor eIF4G with the IRES domain V. Moreover, binding of eIF4B and the polypyrimidine tract-binding protein and the association of ribosomes with the viral RNA are affected by these mutations. However, the negative effects of the IRES mutations are completely relieved by addition of purified eIF4F. This indicates that eIF4G is the crucial factor that initially binds to the poliovirus IRES and recruits the IRES to the other components of the translational apparatus, while impaired binding of eIF4G plays a key role in attenuation of poliovirus neurovirulence.     

The effect of a tightly bound water molecule on scaffold diversity in the computer-aided de novo ligand design of CDK2 inhibitors

We have determined the effects that tightly bound water molecules have on the de novo design of cyclin-dependent kinase-2 (CDK2) ligands. In particular, we have analyzed the impact of a specific structural water molecule on the chemical diversity and binding mode of ligands generated through a de novo structure-based ligand generation method in the binding site of CDK2. The tightly bound water molecule modifies the size and shape of the binding site and we have found that it also imposed constraints on the observed binding modes of the generated ligands. This in turn had the indirect effect of reducing the chemical diversity of the underlying molecular scaffolds that were able to bind to the enzyme satisfactorily.      

Regulatory interactions in the dimeric cytochrome bc1 complex: The advantages of being a twin

The dimeric cytochrome bc₁ complex catalyzes the oxidation-reduction of quinol and quinone at sites located in opposite sides of the membrane in which it resides. We review the kinetics of electron transfer and inhibitor binding that reveal functional interactions between the quinol oxidation site at center P and quinone reduction site at center N in opposite monomers in conjunction with electron equilibration between the cytochrome b subunits of the dimer. A model for the mechanism of the bc₁ complex has emerged from these studies in which binding of ligands that mimic semiquinone at center N regulates half-of-the-sites reactivity at center P and binding of ligands that mimic catalytically competent binding of ubiquinol at center P regulates half-of-the-sites reactivity at center N. An additional feature of this model is that inhibition of quinol oxidation at the quinone reduction site is avoided by allowing catalysis in only one monomer at a time, which maximizes the number of redox acceptor centers available in cytochrome b for electrons coming from quinol oxidation reactions at center P and minimizes the leakage of electrons that would result in the generation of damaging oxygen radicals.     

Interaction between the 5'-terminal cloverleaf and 3AB/3CDpro of poliovirus is essential for RNA replication.

On the basis of sequence alignments and secondary structure comparisons of the first 100 nucleotides of enterovirus and rhinovirus RNAs, chimeric constructs in which this region of poliovirus type 1 Mahoney [PV1(M)] is replaced with that of human rhinovirus type 2 (HRV2) or HRV14 have been engineered. These chimeric constructs contain the internal ribosomal entry site of either poliovirus or encephalomyocarditis virus. Independent of the internal ribosomal entry site elements, only the constructs containing either the PV1(M) or HRV2 cloverleaf sequences yielded viable viruses. The secondary structures of all three cloverleaves are quite similar. However, highly purified polioviral proteins 3CDpro and 3AB together bound to the PV1(M) and HRV2 cloverleaves, albeit with different affinities, whereas the HRV14 homolog did not interact with these proteins to any appreciable extent. These results support a mechanism of poliovirus genomic replication in which the formation of a complex between the cloverleaf structure and the 3CDpro/3AB proteins of poliovirus plays an essential role.     

Ligand binding study of human PEBP1/RKIP: interaction with nucleotides and Raf-1 peptides evidenced by NMR and mass spectrometry

Background

Human Phosphatidylethanolamine binding protein 1 (hPEBP1) also known as Raf kinase inhibitory protein (RKIP), affects various cellular processes, and is implicated in metastasis formation and Alzheimer''s disease. Human PEBP1 has also been shown to inhibit the Raf/MEK/ERK pathway. Numerous reports concern various mammalian PEBP1 binding ligands. However, since PEBP1 proteins from many different species were investigated, drawing general conclusions regarding human PEBP1 binding properties is rather difficult. Moreover, the binding site of Raf-1 on hPEBP1 is still unknown.

Methods/Findings

In the present study, we investigated human PEBP1 by NMR to determine the binding site of four different ligands: GTP, FMN, and one Raf-1 peptide in tri-phosphorylated and non-phosphorylated forms. The study was carried out by NMR in near physiological conditions, allowing for the identification of the binding site and the determination of the affinity constants K_D for different ligands. Native mass spectrometry was used as an alternative method for measuring K_D values.

Conclusions/Significance

Our study demonstrates and/or confirms the binding of hPEBP1 to the four studied ligands. All of them bind to the same region centered on the conserved ligand-binding pocket of hPEBP1. Although the affinities for GTP and FMN decrease as pH, salt concentration and temperature increase from pH 6.5/NaCl 0 mM/20°C to pH 7.5/NaCl 100 mM/30°C, both ligands clearly do bind under conditions similar to what is found in cells regarding pH, salt concentration and temperature. In addition, our work confirms that residues in the vicinity of the pocket rather than those within the pocket seem to be required for interaction with Raf-1.     

Structural factors that control conformational transitions and serotype specificity in type 3 poliovirus.

The three-dimensional structure of the Sabin strain of type 3 poliovirus has been determined at 2.4 A resolution. Significant structural differences with the Mahoney strain of type 1 poliovirus are confined to loops and terminal extensions of the capsid proteins, occur in all of the major antigenic sites of the virion and typically involve insertions, deletions or the replacement of prolines. Several newly identified components of the structure participate in assembly-dependent interactions which are relevant to the biologically important processes of viral assembly and uncoating. These include two sites of lipid substitution, two putative nucleotides and a beta sheet formed by the N-termini of capsid proteins VP4 and VP1. The structure provides an explanation for the temperature sensitive phenotype of the P3/Sabin strain. Amino acids that regulate temperature sensitivity in type 3 poliovirus are located in the interfaces between promoters, in the binding site for a lipid substituent and in an assembly-dependent extended beta sheet that stabilizes the association of pentamers. Several lines of evidence indicate that these structural components also control conformational transitions at various stages of the viral life cycle.     

Trypsin sensitivity of the Sabin strain of type 1 poliovirus: cleavage sites in virions and related particles.

Treatment of the Sabin strain of type 1 poliovirus with trypsin produced two stable fragments of capsid protein VP1 which remained associated with the virions. Trypsinized virus was fully infectious and was neutralized by type-specific antisera. The susceptible site in the Sabin 1 strain was between the lysine at position 99 and the asparagine at position 100. A similar tryptic cleavage occurred in the Leon and Sabin strains of type 3 poliovirus, probably at the arginine at position 100, but not in the type 1 Mahoney strain, which lacks a basic residue at either position 99 or position 100. Tryptic treatment of heat-treated virus and 14S assembly intermediates produced unique stable fragments which were different from those produced in virions. The implications of our results for future characterization of the surface structures of these particles and structural rearrangements in the poliovirus capsid are discussed.     

Substrate requirements of human rhinovirus 3C protease for peptide cleavage in vitro

A series of synthetic peptides representing authentic proteolytic cleavage sites of human rhinovirus type 14 were assayed as substrates for purified 3C protease. Competition cleavage assays were employed to determine the relative specificity constants (Kcat/Km) for substrates with sequences related to the viral 2C-3A cleavage site. Variable length peptides representing the 2C-3A cleavage site were cleaved with comparable efficiency. These studies defined a minimum substrate of 6 amino acids (TLFQ/GP), although retention of the residue at position P5 (ETLFQ/GP) resulted in a better substrate by an order of magnitude. Amino acid substitutions at position P5, P4, P1', or P2' indicated that the identity of the residue at position P5 was not critical, whereas substitutions at position P4, P1' or P2' resulted in substrates with Kcat/Km values varying over 2 orders of magnitude. In contrast to the 2C-3A cleavage site, small peptide derivatives representative of the 3A-3B cleavage site were relatively poor substrates, which suggested that residues flanking the minimum core sequence may influence susceptibility to cleavage. The 3C protease of rhinovirus type 14 was also capable of cleaving peptides representing comparable cleavage sites predicted for coxsackie B virus and poliovirus.