Search PubMed⌕ Search

Biomedical subjects

Jill E Gready

Publications and source records attributed to Jill E Gready.

13 recordsLinked to original sources
PubMed2005

Computational methods for the study of enzymic reaction mechanisms III: a perturbation plus QM/MM approach for calculating relative free energies of protonation.

We describe a coupling parameter, that is, perturbation, approach to effectively create and annihilate atoms in the quantum mechanical Hamiltonian within the closed shell restricted Hartree-Fock formalism. This perturbed quantum mechanical atom (PQA) method is combined with molecular mechanics (MM) methods (PQA/MM) within a molecular dynamics simulation, to model the protein environment (MM region) effects that also make a contribution to the overall free energy change. Using the semiempirical PM3 method to model the QM region, the application of this PQA/MM method is illustrated by calculation of the relative protonation free energy of the conserved OD2 (Asp27) and the N5 (dihydrofolate) proton acceptor sites in the active site of Escherichia coli dihydrofolate reductase (DHFR) with the bound nicotinamide adenine dinucleotide phosphate (NADPH) cofactor. For a number of choices for the QM region, the relative protonation free energy was calculated as the sum of contributions from the QM region and the interaction between the QM and MM regions via the thermodynamic integration (TI) method. The results demonstrate the importance of including the whole substrate molecule in the QM region, and the overall protein (MM) environment in determining the relative stabilities of protonation sites in the enzyme active site. The PQA/MM free energies obtained by TI were also compared with those estimated by a less computationally demanding nonperturbative method based on the linear response approximation (LRA). For some choices of QM region, the total free energies calculated using the LRA method were in very close agreement with the PQA/MM values. However, the QM and QM/MM component free energies were found to differ significantly between the two methods.

Algorithms↗
PubMed2005

The prion protein gene: identifying regulatory signals using marsupial sequence.

The function of the prion protein gene (PRNP) and its normal product PrP(C) is elusive. We used comparative genomics as a strategy to understand the normal function of PRNP. As the reliability of comparisons increases with the number of species and increased evolutionary distance, we isolated and sequenced a 66.5 kb BAC containing the PRNP gene from a distantly related mammal, the model Australian marsupial Macropus eugenii (tammar wallaby). Marsupials are separated from eutherians such as human and mouse by roughly 180 million years of independent evolution. We found that tammar PRNP, like human PRNP, has two exons. Prion proteins encoded by the tammar wallaby and a distantly related marsupial, Monodelphis domestica (Brazilian opossum) PRNP contain proximal PrP repeats with a distinct, marsupial-specific composition and a variable number. Comparisons of tammar wallaby PRNP with PRNPs from human, mouse, bovine and ovine allowed us to identify non-coding gene regions conserved across the marsupial-eutherian evolutionary distance, which are candidates for regulatory regions. In the PRNP 3' UTR we found a conserved signal for nuclear-specific polyadenylation and the putative cytoplasmic polyadenylation element (CPE), indicating that post-transcriptional control of PRNP mRNA activity is important. Phylogenetic footprinting revealed conserved potential binding sites for the MZF-1 transcription factor in both upstream promoter and intron/intron 1, and for the MEF2, MyT1, Oct-1 and NFAT transcription factors in the intron(s). The presence of a conserved NFAT-binding site and CPE indicates involvement of PrP(C) in signal transduction and synaptic plasticity.

3' Untranslated Regions↗
PubMed2005

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase.

R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67*DHF*NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand-ligand and ligand-protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.

Binding Sites↗
PubMed2004

Evolution of vertebrate genes related to prion and Shadoo proteins--clues from comparative genomic analysis.

Recent findings of new genes in fish related to the prion protein (PrP) gene PRNP, including our recent report of SPRN coding for Shadoo (Sho) protein found also in mammals, raise issues of their function and evolution. Here we report additional novel fish genes found in public databases, including a duplicated SPRN gene, SPRNB, in Fugu, Tetraodon, carp, and zebrafish encoding the Sho2 protein, and we use comparative genomic analysis to analyze the evolutionary relationships and to infer evolutionary trajectories of the complete data set. Phylogenetic footprinting performed on aligned human, mouse, and Fugu SPRN genes to define candidate regulatory promoter regions, detected 16 conserved motifs, three of which are known transcription factor-binding sites for a receptor and transcription factors specific to or associated with expression in brain. This result and other homology-based (VISTA global genomic alignment; protein sequence alignment and phylogenetics) and context-dependent (genomic context; relative gene order and orientation) criteria indicate fish and mammalian SPRN genes are orthologous and suggest a strongly conserved basic function in brain. Whereas tetrapod PRNPs share context with the analogous stPrP-2-coding gene in fish, their sequences are diverged, suggesting that the tetrapod and fish genes are likely to have significantly different functions. Phylogenetic analysis predicts the SPRN/SPRNB duplication occurred before divergence of fish from tetrapods, whereas that of stPrP-1 and stPrP-2 occurred in fish. Whereas Sho appears to have a conserved function in vertebrate brain, PrP seems to have an adaptive role fine-tuned in a lineage-specific fashion. An evolutionary model consistent with our findings and literature knowledge is proposed that has an ancestral prevertebrate SPRN-like gene leading to all vertebrate PrP-related and Sho-related genes. This provides a new framework for exploring the evolution of this unusual family of proteins and for searching for members in other fish branches and intermediate vertebrate groups.

Animals↗
PubMed2004

C-type lectin-like domains in Fugu rubripes.

BACKGROUND: Members of the C-type lectin domain (CTLD) superfamily are metazoan proteins functionally important in glycoprotein metabolism, mechanisms of multicellular integration and immunity. Three genome-level studies on human, C. elegans and D. melanogaster reported previously demonstrated almost complete divergence among invertebrate and mammalian families of CTLD-containing proteins (CTLDcps). RESULTS: We have performed an analysis of CTLD family composition in Fugu rubripes using the draft genome sequence. The results show that all but two groups of CTLDcps identified in mammals are also found in fish, and that most of the groups have the same members as in mammals. We failed to detect representatives for CTLD groups V (NK cell receptors) and VII (lithostathine), while the DC-SIGN subgroup of group II is overrepresented in Fugu. Several new CTLD-containing genes, highly conserved between Fugu and human, were discovered using the Fugu genome sequence as a reference, including a CSPG family member and an SCP-domain-containing soluble protein. A distinct group of soluble dual-CTLD proteins has been identified, which may be the first reported CTLDcp group shared by invertebrates and vertebrates. We show that CTLDcp-encoding genes are selectively duplicated in Fugu, in a manner that suggests an ancient large-scale duplication event. We have verified 32 gene structures and predicted 63 new ones, and make our annotations available through a distributed annotation system (DAS) server http://anz.anu.edu.au:8080/Fugu_rubripes/ and their sequences as additional files with this paper. CONCLUSIONS: The vertebrate CTLDcp family was essentially formed early in vertebrate evolution and is completely different from the invertebrate families. Comparison of fish and mammalian genomes revealed three groups of CTLDcps and several new members of the known groups, which are highly conserved between fish and mammals, but were not identified in the study using only mammalian genomes. Despite limitations of the draft sequence, the Fugu rubripes genome is a powerful instrument for gene discovery and vertebrate evolutionary analysis. The composition of the CTLDcp superfamily in fish and mammals suggests that large-scale duplication events played an important role in the evolution of vertebrates.

Amino Acid Sequence↗
PubMed2004

Conformation of prion protein repeat peptides probed by FRET measurements and molecular dynamics simulations.

We report the combined use of steady-state fluorescence resonance energy transfer (FRET) experiments and molecular dynamics (MD) simulations to investigate conformational distributions of the prion protein (PrP) repeat system. FRET was used for the first time to probe the distance, as a function of temperature and pH, between a donor Trp residue and an acceptor dansyl group attached to the N-terminus in seven model peptides containing one to three repeats of the second decarepeat of PrP from marsupial possum (PHPGGSNWGQ)nG, and one and two human PrP consensus octarepeats (PHGGGWGQ)nG. In multirepeat peptides, single-Trp mutants were made by replacing other Trp(s) with Phe. As previous work has shown PrP repeats do not adopt a single preferred stable conformation, the FRET values are averages reflecting heterogeneity in the donor-acceptor distances. The T-dependence of the conformational distributions, and derived average dansyl-Trp distances, were obtained directly from MD simulation of the marsupial dansyl-PHPGGSNWGQG peptide. The results show excellent agreement between the FRET and MD T-dependent distances, and demonstrate the remarkable sensitivity and reproducibility of the FRET method in this first-time use for a set of disordered peptides. Based on the results, we propose a model involving cation-pi or pi-pi His-Trp interactions to explain the T- (5-85 degrees C) and pH- (6.0, 7.2) dependencies on distance, with HW i, i + 4 or WH i, i + 4 separations in sequence being more stable than HW i, i + 6 or WH i, i + 6 separations. The model has peptides adopting loosely folded conformations, with dansyl-Trp distances very much less than estimates for fully extended conformations, for example, approximately 16 vs. 33, approximately 21 vs. 69, and approximately 22 vs. 106 A for 1-3 decarepeats, and approximately 14 vs. 25 and approximately 19 vs. 54 A for 1-2 octarepeats, respectively. The study demonstrates the usefulness of combining FRET with MD, a combination reported only once previously. Initial "mapping" of the conformational distribution of flexible peptides by simulation can assist in designing and interpreting experiments using steady-state intensity methods, and indicating how time-resolved or anisotropy methods might be used.

Amino Acid Sequence↗
PubMed2004

An efficient Z-score algorithm for assessing sequence alignments.

We describe an alternative method for scoring of the pairwise alignment of two biological sequences. Designed to overcome the bias due to the composition of the alignment, it measures the distance (in standard deviations) between the given alignment and the mean value of all other alignments that can be obtained by a permutation of either sequence. We demonstrate that the standard deviation can be calculated efficiently. By concentrating upon the ungapped case, the mean and standard deviation can be calculated exactly and in two steps, the first being O(N) time, where N is the length of the sequence, the second in a fixed number of calculations, i.e., in O(1) time. We argue that this statistic is a more consistent measure than a similarity score based upon a standard scoring matrix. Even in the ungapped case, the statistic proves in many cases to be more accurate than the commonly used (FASTA) (Pearson and Lipman, 1988) gapped Z-score in which the sequence is matched against a random sample of the database. We demonstrate the use of the POZ-score as a secondary filter which screens out several well-known types of false positive, reducing the amount of manual screening to be done by the biologist.

Algorithms↗
PubMed2003

Shadoo, a new protein highly conserved from fish to mammals and with similarity to prion protein.

We report evidence from cDNA isolation and expression analysis as well as analyses of genome, expressed sequence tag (EST), cDNA and expression databases for a new gene named SPRN (shadow of prion protein). SPRN comprises two exons, with the open reading frame (ORF) contained within exon 2, and codes for a protein of 130-150 amino acids named Shadoo (Japanese shadow), predicted to be extracellular and GPI-anchored. The SPRN gene was found in fish (zebrafish, Fugu) and mammals (mouse, rat, human). Conservation of order and transcription orientation of two proximal genes between fishes and mammals strongly indicates gene orthology. Sequence comparison shows: a highly conserved N-terminal signal sequence; Arg-rich basic region containing up to six tetrarepeats of consensus XXRG (where X is G, A or S); a hydrophobic region of 20 residues with strong homology to PrP; a less conserved C-terminal domain containing a conserved glycosylation motif; and a C-terminal peptide predicted to be a signal sequence for glycophosphotidylinositol (GPI)-anchor attachment. Fish Shadoos (Sho) show well conserved sequences (identity 54%) over 106 amino acids (zebrafish length), and conservation among the mammalian sequences is very high (identity 81-96%). The fish and mammalian sequences are also well conserved, particularly for zebrafish, to beyond the end of the hydrophobic sequence (identity 41-53%, 78 amino acids, zebrafish length). The overall structure appears closely related to prion proteins (PrPs), although the C-terminal domains of Shos are quite different from those of PrPs, for which conformational changes in mammals are implicated in disease. The structural similarity is particularly interesting given recent reports of three new genes with similarities to PrPs found in Fugu (PrP-like, PrP-461/stPrP-1 and stPrP-2) and other fish, but for which direct evolution to higher vertebrate PrPs is unlikely and for which no other mammalian homologues have been found. Database information indicates expression of SPRN in embryo, brain and retina of mouse and rat, hippocampus of human, and in embryo and retina of zebrafish, and we directly confirmed a strikingly specific expression of the mammalian (human, mouse, rat) transcripts in whole brain. This result together with some common structural features led to the suggestive hypothesis of a possible functional link between mammalian PrP and Sho proteins.

Amino Acid Sequence↗
PubMed2003

Comparative analysis of structural properties of the C-type-lectin-like domain (CTLD).

The superfamily of proteins containing the C-type-lectin-like domain (CTLD) is a group of abundant extracellular metazoan proteins characterized by evolutionary flexibility and functional versatility. Several CTLDs are also found in parasitic prokaryotes and viruses. The 37 distinct currently available CTLD structures demonstrate significant structural conservation despite low or undetectable sequence similarity. Our aim in this study was to perform an extensive comparative analysis of all available CTLD structures to establish the most conserved structural features of the fold, and to test and extend the early analysis of Drickamer. By implication, these features should be those critical for maintenance of integrity of the fold. By analyzing CTLD structures superimposed by several methods, we have established groups of conserved structural positions involved in fold maintenance but not in ligand binding; these are consistent with the fold's known functional flexibility. In addition to the well-recognized disulfide bridges, groups of conserved residues are involved in hydrophobic interactions stabilizing the core of the fold and the long loop region, and in an alpha2-beta1-beta5 polar interaction. Evaluation of the conclusions of the structure comparison study compared with alignments of all available human, mouse and Caenorhabditis elegans CTLD sequences showed that conservation patterns are preserved throughout the whole CTLD sequence space. Our observations provide an improved understanding of CTLD structure, and will help in identification of new CTLDs and the mechanisms that drive and constrain the coevolution of the structure and function of the fold.

Algorithms↗
PubMed2002

Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis.

QM/MM methods have been developed as a computationally feasible solution to QM simulation of chemical processes, such as enzyme-catalyzed reactions, within a more approximate MM representation of the condensed-phase environment. However, there has been no independent method for checking the quality of this representation, especially for highly nonisotropic protein environments such as those surrounding enzyme active sites. Hence, the validity of QM/MM methods is largely untested. Here we use the possibility of performing all-QM calculations at the semiempirical PM3 level with a linear-scaling method (MOZYME) to assess the performance of a QM/MM method (PM3/AMBER94 force field). Using two model pathways for the hydride-ion transfer reaction of the enzyme dihydrofolate reductase studied previously (Titmuss et al., Chem Phys Lett 2000, 320, 169-176), we have analyzed the reaction energy contributions (QM, QM/MM, and MM) from the QM/MM results and compared them with analogous-region components calculated via an energy partitioning scheme implemented into MOZYME. This analysis further divided the MOZYME components into Coulomb, resonance and exchange energy terms. For the model in which the MM coordinates are kept fixed during the reaction, we find that the MOZYME and QM/MM total energy profiles agree very well, but that there are significant differences in the energy components. Most significantly there is a large change (approximately 16 kcal/mol) in the MOZYME MM component due to polarization of the MM region surrounding the active site, and which arises mostly from MM atoms close to (<10 A) the active-site QM region, which is not modelled explicitly by our QM/MM method. However, for the model where the MM coordinates are allowed to vary during the reaction, we find large differences in the MOZYME and QM/MM total energy profiles, with a discrepancy of 52 kcal/mol between the relative reaction (product-reactant) energies. This is largely due to a difference in the MM energies of 58 kcal/mol, of which we can attribute approximately 40 kcal/mol to geometry effects in the MM region and the remainder, as before, to MM region polarization. Contrary to the fixed-geometry model, there is no correlation of the MM energy changes with distance from the QM region, nor are they contributed by only a few residues. Overall, the results suggest that merely extending the size of the QM region in the QM/MM calculation is not a universal solution to the MOZYME- and QM/MM-method differences. They also suggest that attaching physical significance to MOZYME Coulomb, resonance and exchange components is problematic. Although we conclude that it would be possible to reparameterize the QM/MM force field to reproduce MOZYME energies, a better way to account for both the effects of the protein environment and known deficiencies in semiempirical methods would be to parameterize the force field based on data from DFT or ab initio QM linear-scaling calculations. Such a force field could be used efficiently in MD simulations to calculate free energies.

Binding Sites↗
PubMed2002

Copper-induced conformational change in a marsupial prion protein repeat peptide probed using FTIR spectroscopy.

We report the first Fourier transform infrared analysis of prion protein (PrP) repeats and the first study of PrP repeats of marsupial origin. Large changes in the secondary structure and an increase in hydrogen bonding within the peptide groups were evident from a red shift of the amide I band by >7 cm(-1) and an approximately five-fold reduction in amide hydrogen-deuterium exchange for peptide interacting with Cu(2+) ions. Changes in the tertiary structure upon copper binding were also evident from the appearance of a new band at 1564 cm(-1), which arises from the ring vibration of histidine. The copper-induced conformational change is pH dependent, and occurs at pH >7.

Animals↗
PubMedDate not supplied

Synthesis of quaternised 2-aminopyrimido[4,5-d]pyrimidin-4(3H)-ones and their biological activity with dihydrofolate reductase.

In a program to design and develop mechanism-based compounds active as substrates and inhibitors of dihydrofolate reductase (DHFR), we report the synthesis and physical properties of the 6-methyl- (7), 8-methyl- (8a), and 8-ethyl- (8b) derivatives of the parent 2-aminopyrimido[4,5-d]pyrimidin-4-(3H)-one (6). These compounds are the first members of a class of heterocycles related to 8-alkylpterins (N8-alkyl-2-aminopteridin-4(8H)-ones) (2a-2c), which have been shown to be novel substrates for DHFR. Three methods were developed for the synthesis of target compounds 7, 8a and 8b; however, the optimum yields (1-8%) could not be improved because the products decomposed by ring opening (e.g. to 2,4-diamino-5-methyliminomethylpyrimidin-6(1H)-one (9)) under the reaction conditions. The marked pi-electron deficiency of compounds 7, 8a and 8b is the likely cause for the susceptibility of the quaternised pyrimidine ring in the related cations 10, 15a and 15b, respectively, to add nucleophiles, thus promoting the opening of the pyrimidopyrimidine ring system. 1H-NMR spectroscopic studies of compounds 7, 8a and 8b revealed a fast and reversible covalent hydration of the associated cations across the C7z.sbnd;N8 bond for the N6-methyl derivative 7 and across the N6z.sbnd;C7 bond for the N8-methyl derivative 8a. UV spectroscopic studies of methyl derivatives 7 and 8a as well as the parent heterocycle 6 showed that protonation of the latter occurred at N1, while methylation with iodomethane proceeded at N6 and N8. The basicities of the N-methyl derivatives 7 and 8a (pK(a) ca. 5.5) are similar to those of 8-alkylpterins 2; this is an essential element of the design to promote binding to DHFR in their protonated form. Enzyme kinetics of 7, 8a and 8b with chicken DHFR confirmed our predictions that they are substrates, with apparent K(m) values of 3.8, 0.08, and 0.65 mM, and apparent V(max) values of 0.47, 2.27, and 0.30 nmol L(-1) min(-1) (for enzyme concentration 0.122 micro M), respectively. The parent compound 6 was not a substrate.

Animals↗