Search PubMed⌕ Search

Biomedical subjects

J D Hirst

Publications and source records attributed to J D Hirst.

13 recordsLinked to original sources

The evolutionary landscape of functional model proteins.

To study the distinct influences of structure and function on evolution, we propose a minimalist model for proteins with binding pockets, called functional model proteins, based on a shifted-HP model on a two-dimensional square lattice. These model proteins are not maximally compact and contain an empty lattice site surrounded by at least three nearest neighbors, thus providing a binding pocket. Functional model proteins possess a unique native state, cooperative folding and tolerance to mutation. Due to the explicit functionality in these models (by design), we have been able to explore their fitness or evolutionary landscapes, as characterized by the size and distribution of homologous families and by the complexity of the inter-relatedness of the functional model proteins. Mindful that these minimalist models are highly idealized and two-dimensional, functional model proteins should nevertheless provide a useful means for exploring the constraints of maintaining structure and function on the evolution of proteins.

Amino Acid Sequence↗

Do active site conformations of small ligands correspond to low free-energy solution structures?

We compare the low free energy structures of ten small, polar ligands in solution to their conformations in their respective receptor active sites. The solution conformations are generated by a systematic search and the free energies of representative structures are computed with a continuum solvation model. Based on the values of torsion angles, we find little similarity between low energy solution structures of small ligands and their active site conformations. However, in nine out of ten cases, the positions of 'anchor points' (key atoms responsible for tight binding) in the lowest energy solution structures are very similar to the positions of these atoms in the active site conformations. A metric that more closely captures the essentials of binding supports the basic premise underlying pharmacophore mapping, namely that active site conformations of small flexible ligands correspond to their low energy structures in solution. This work supports the efforts of building pharmacophore models based on the information present in solution structures of small isolated ligands.

Binding Sites↗

Nonlinear quantitative structure-activity relationship for the inhibition of dihydrofolate reductase by pyrimidines.

A novel method for quantitative structure-activity relationship (QSAR) analysis is presented. The method, which does not assume any particular functional form for the QSAR, develops nonlinear relationships between parameters describing a set of molecules and the activity of the molecules. For the QSAR of the inhibition of Escherichia coli dihydrofolate reductase by 2,4-diamino-5-(substituted benzyl)pyrimidines, the method compares favorably to other nonlinear methods. Cross-validation trials demonstrate that the predictive ability is as accurate as other methods, and the method is simpler and faster than neural network and machine-learning methods. Consequently, its implementation is much easier, and interpretation of the generated QSAR is more straightforward.

Escherichia coli↗

Predicting leucine zipper structures from sequence.

The leucine zipper structure is adopted by one family of the coiled coil proteins. Leucine zippers have a characteristic leucine repeat: Leu-X6-Leu-X6-Leu-X6-Liu (where X may be any residue). However, many sequences have the leucine repeat, but do not adopt the leucine zipper structure (we shall refer to these as non-zippers). We have found and analyzed residue pair patterns that allow one to identify correctly 90% of leucine zippers and 97% of non-zippers. Simpler analyses, based on the frequency of occurrence of residues at certain positions, specify, at most, 65% of zippers and 80-90% of non-zippers. Both short and long patterns contribute to the successful discrimination of leucine zippers from non-zippers. A number of these patterns involve hydrophobic residues that would be placed on the solvent-exposed surface of the helix, were the sequence to adopt a leucine zipper structure. Thus, an analysis of protein sequences has allowed us to improve discrimination between leucine zippers and non-zippers, and has provided some further insight into the physical factors influencing the leucine zipper structure.

Amino Acid Sequence↗

Molecular dynamics simulations of isolated helices of myoglobin.

The apo form of myoglobin has two non-native stable states that have been experimentally characterized. Investigation of these states has suggested possible folding pathways for myoglobin. We have performed molecular dynamics simulations on solvated isolated helices of myoglobin to investigate the relationship between the intrinsic stabilities of the isolated helices and the structure and folding pathway of apomyoglobin. Analyses of hydrogen bonding and fluctuations from simulations at 298 and 368 K are used to explore the relative stabilities of the helices of myoglobin. The ordering observed is A approximately G approximately H > B > E > F, which mirrors both the experimental equilibrium and kinetic data available for apomyoglobin. The experimental observation that a subdomain comprising helices A, G, and H is an important early intermediate and our result that these helices are the most stable suggest that the intrinsically more stable helices form early in the folding process and that this significantly influences the folding pathway.

Amino Acid Sequence↗

Helicity, circular dichroism and molecular dynamics of proteins.

In protein unfolding studies, reduction in circular dichroism (CD) at 222 nm has been interpreted as loss of helicity. Estimates of the helicity of a protein from its CD spectrum are calibrated by reference to X-ray crystal structures, based on the assumption that the mean residue ellipticity at 222 nm is directly proportional to the number of residues in a helix. We have examined various influences on the CD at 222 nm, using molecular dynamics simulations to provide the structural detail required. We have found that the fragmentation of long helices, without a reduction in the number of helical residues, significantly reduces the mean residue ellipticity at 222 nm. The dynamical motion of the protein and the precise conformation of helical residues also play an important role. We discuss the implications of these factors to the interpretation of CD for the partial unfolding of apomyoglobin.

Apoproteins↗

Quantitative structure-activity relationships by neural networks and inductive logic programming. I. The inhibition of dihydrofolate reductase by pyrimidines.

Neural networks and inductive logic programming (ILP) have been compared to linear regression for modelling the QSAR of the inhibition of E. coli dihydrofolate reductase (DHFR) by 2,4-diamino-5-(substituted benzyl)pyrimidines, and, in the subsequent paper [Hirst, J.D., King, R.D. and Sternberg, M.J.E. J. Comput.-Aided Mol. Design, 8 (1994) 421], the inhibition of rodent DHFR by 2,4-diamino-6,6-dimethyl-5-phenyl-dihydrotriazines. Cross-validation trials provide a statistically rigorous assessment of the predictive capabilities of the methods, with training and testing data selected randomly and all the methods developed using identical training data. For the ILP analysis, molecules are represented by attributes other than Hansch parameters. Neural networks and ILP perform better than linear regression using the attribute representation, but the difference is not statistically significant. The major benefit from the ILP analysis is the formulation of understandable rules relating the activity of the inhibitors to their chemical structure.

Animals↗

Quantitative structure-activity relationships by neural networks and inductive logic programming. II. The inhibition of dihydrofolate reductase by triazines.

One of the largest available data sets for developing a quantitative structure-activity relationship (QSAR)--the inhibition of dihydrofolate reductase (DHFR) by 2,4-diamino-6,6-dimethyl-5-phenyl-dihydrotriazine derivatives--has been used for a sixfold cross-validation trial of neural networks, inductive logic programming (ILP) and linear regression. No statistically significant difference was found between the predictive capabilities of the methods. However, the representation of molecules by attributes, which is integral to the ILP approach, provides understandable rules about drug-receptor interactions.

Animals↗

Prediction of structural and functional features of protein and nucleic acid sequences by artificial neural networks.

The applications of artificial neural networks to the prediction of structural and functional features of protein and nucleic acid sequences are reviewed. A brief introduction to neural networks is given, including a discussion of learning algorithms and sequence encoding. The protein applications mostly involve the prediction of secondary and tertiary structure from sequence. The problems in nucleic acid analysis tackled by neural networks are the prediction of translation initiation sites in Escherichia coli, the recognition of splice junctions in human mRNA, and the prediction of promoter sites in E. coli. The performance of the approach is compared with other current statistical methods.

Algorithms↗

Prediction of ATP/GTP-binding motif: a comparison of a perceptron type neural network and a consensus sequence method [corrected].

Neural networks have been applied to a number of protein structure problems. In some applications their success has not been substantiated by a comparison with the performance of a suitable alternative statistical method on the same data. In this paper, a two-layer feed-forward neural network has been trained to recognize ATP/GTP-binding [corrected] local sequence motifs. The neural network correctly classified 78% of the 349 sequences used. This was much better than a simple motif-searching program. A more sophisticated statistical method was developed, however, which performed marginally better (80% correct classification) than the neural network. The neural network and the statistical method performed similarly on sequences of varying degrees of homology. These results do not imply that neural networks, especially those with hidden layers, are not useful tools, but they do suggest that two-layer networks in particular should be carefully tested against other statistical methods.

Adenosine Triphosphate↗