Biomedical subjects

S C Lovell

Publications and source records attributed to S C Lovell.

8 recordsLinked to original sources

PubMed2001

Sequence-structure homology recognition by iterative alignment refinement and comparative modeling.

Our approach to fold recognition for the fourth critical assessment of techniques for protein structure prediction (CASP4) experiment involved the use of the FUGUE sequence-structure homology recognition program (http://www-cryst.bioc.cam.ac.uk/fugue), followed by model building. We treat models as hypotheses and examine these to determine whether they explain the available data. Our method depends heavily on environment-specific substitution tables derived from our database of structural alignments of homologous proteins (HOMSTRAD, http://www-cryst.bioc.cam.ac.uk/homstrad/). FUGUE uses these tables to incorporate structural information into profiles created from HOMSTRAD alignments that are matched against a profile created for the target from multiple sequence alignment. In addition, environment-specific substitution tables are used throughout the modeling procedure and as part of the model evaluation. Annotation of sequence alignments with JOY, to reflect local structural features, proved valuable, both for modifying hypotheses, and for rejecting predictions when the expected pattern of conservation is not observed. Our stringency in rejecting incorrect predictions led us to submit a relatively small number of models, including only a low number of false positives, resulting in a high average score.

Amino Acid Sequence↗

PubMed2000

The penultimate rotamer library.

All published rotamer libraries contain some rotamers that exhibit impossible internal atomic overlaps if built in ideal geometry with all hydrogen atoms. Removal of uncertain residues (mainly those with B-factors >/=40 or van der Waals overlaps >/=0.4 A) greatly improves the clustering of rotamer populations. Asn, Gln, or His side chains additionally benefit from flipping of their planar terminal groups when required by atomic overlaps or H-bonding. Sensitivity to skew and to the boundaries of chi angle bins is avoided by using modes rather than traditional mean values. Rotamer definitions are listed both as the modal values and in a preferred version that maximizes common atoms between related rotamers. The resulting library shows significant differences from previous ones, differences validated by considering the likelihood of systematic misfitting of models to electron density maps and by plotting changes in rotamer frequency with B-factor. Few rotamers now show atomic overlaps in ideal geometry; those overlaps are relatively small and can be understood in terms of bond angle distortions compensated by favorable interactions. The new library covers 94.5% of examples in the highest quality protein data with 153 rotamers and can make a significant contribution to improving the accuracy of new structures. Proteins 2000;40:389-408.

Amino Acids↗

PubMed2000

Exploring steric constraints on protein mutations using MAGE/PROBE.

When planning a mutation to test some hypothesis, one crucial question is whether the new side chain is compatible with the existing structure; only if it is compatible can the interpretation of mutational results be straightforward. This paper presents a simple way of using the sensitive geometry of all-atom contacts (including hydrogens) to answer that question. The interactive MAGE/PROBE system lets the biologist explore conformational space for the mutant side chain, with an interactively updated kinemage display of its all-atom contacts to the original structure. The Autobondrot function in PROBE systematically explores that same conformational space, outputting contact scores at each point, which are then contoured and displayed. These procedures are applied here in two types of test cases, with known mutant structures. In ricin A chain, the ability of a neighboring glutamate to rescue activity of an active-site mutant is modeled successfully. In T4 lysozyme, six mutations to Leu are analyzed within the wild-type background structure, and their Autobondrot score maps correctly predict whether or not their surroundings must shift significantly in the actual mutant structures; interactive examination of contacts for the conformations involved explains which clashes are relieved by the motions. These programs are easy to use, are available free for UNIX or Microsoft Windows operating systems, and should be of significant help in choosing good mutation experiments or in understanding puzzling results.

Bacteriophage T4↗

PubMed1999

Visualizing and quantifying molecular goodness-of-fit: small-probe contact dots with explicit hydrogen atoms.

The technique of small-probe contact dot surfaces is described as a method for calculating and displaying the detailed atomic contacts inside or between molecules. It allows one both to measure and to visualize directly the goodness-of-fit of packing interactions. It requires both highly accurate structures and also the explicit inclusion of all hydrogen atoms and their van der Waals interactions. A reference dataset of 100 protein structures was chosen on the basis of resolution (1.7 A or better), crystallographic R-value, non-homology, and the absence of any unusual problems. Hydrogen atoms were added in standard geometry and, where needed, with rotational optimization of OH, SH, and NH+3 positions. Side-chain amide orientations were corrected where required by NH van der Waals clashes, as described in the accompanying paper. It was determined that, in general, methyl groups pack well in the default staggered conformation, except for the terminal methyl groups of methionine residues, which required rotational optimization. The distribution of serious clashes (i.e. non-H-bond overlap of >/=0.4 A) was studied as a function of resolution, alternate conformations, and temperature factor (B), leading to the decision that packing and other structural features would not be analyzed for residues in 'b' alternate conformations or with B-factors of 40 or above. At the level of the fine details analyzed here, structural accuracy improves quite significantly over the range from 1.7 to 1.0 A resolution. These high-resolution structures show impressively well-fitted packing interactions, with some regions thoroughly interdigitated and other regions somewhat sparser. Lower-resolution structures or model structures could undoubtedly be improved in accuracy by the incorporation of this additional information: for example, nucleic acid structures in non-canonical conformations are often very accurate for the bases and much less reliable for the backbone, whose conformation could be specified better by including explicit H atom geometry and contacts. The contact dots are an extremely sensitive method of finding problem areas, and often they can suggest how to make improvements. They can also provide explanations for structural features that have been described only as empirical regularities, which is illustrated by showing that the commonest rotamer of methionine (a left-handed spiral, with all chi values near -60 degrees) is preferred because it provides up to five good H atom van der Waals contacts. This methodology is thus applicable in two different ways: (1) for finding and correcting errors in structure models (either experimental or theoretical); and (2) for analyzing interaction patterns in the molecules themselves.

Crystallography, X-Ray↗

PubMed1999

Asparagine and glutamine: using hydrogen atom contacts in the choice of side-chain amide orientation.

Small-probe contact dot surface analysis, with all explicit hydrogen atoms added and their van der Waals contacts included, was used to choose between the two possible orientations for each of 1554 asparagine (Asn) and glutamine (Gln) side-chain amide groups in a dataset of 100 unrelated, high-quality protein crystal structures at 0.9 to 1.7 A resolution. For the movable-H groups, each connected, closed set of local H-bonds was optimized for both H-bonds and van der Waals overlaps. In addition to the Asn/Gln "flips", this process included rotation of OH, SH, NH3+, and methionine methyl H atoms, flip and protonation state of histidine rings, interaction with bound ligands, and a simple model of water interactions. However, except for switching N and O identity for amide flips (or N and C identity for His flips), no non-H atoms were shifted. Even in these very high-quality structures, about 20 % of the Asn/Gln side-chains required a 180 degrees flip to optimize H-bonding and/or to avoid NH2 clashes with neighboring atoms (incorporating a conservative score penalty which, for marginal cases, favors the assignment in the original coordinate file). The programs Reduce, Probe, and Mage provide not only a suggested amide orientation, but also a numerical score comparison, a categorization of the marginal cases, and a direct visualization of all relevant interactions in both orientations. Visual examination allowed confirmation of the raw score assignment for about 40 % of those Asn/Gln flips placed within the "marginal" penalty range by the automated algorithm, while uncovering only a small number of cases whose automated assignment was incorrect because of special circumstances not yet handled by the algorithm. It seems that the H-bond and the atomic-clash criteria independently look at the same structural realities: when both criteria gave a clear answer they agreed every time. But consideration of van der Waals clashes settled many additional cases for which H-bonding was either absent or approximately equivalent for the two main alternatives. With this extra information, 86 % of all side-chain amide groups could be oriented quite unambiguously. In the absence of further experimental data, it would probably be inappropriate to assign many more than this. Some of the remaining 14 % are ambiguous because of coordinate error or inadequacy of the theoretical model, but the great majority of ambiguous cases probably occur as a dynamic mix of both flip states in the actual protein molecule. The software and the 100 coordinate files with all H atoms added and optimized and with amide flips corrected are publicly available.

Algorithms↗

PubMed1999

Asparagine and glutamine rotamers: B-factor cutoff and correction of amide flips yield distinct clustering.

Previous rotamer libraries showed little significant clustering for asparagine chi2 or glutamine chi3 values, but none of those studies corrected amide orientations or omitted disordered side chains. The current survey used 240 proteins at </=1.7 A resolution with <50% homology and <30 clashes per thousand atoms (atomic overlap >/=0.4 A). All H atoms were added and optimized, and amide orientation was flipped by 180 degrees if required by H bonding or atomic clashes. A side chain was included only if its amide orientation was clearly determined and if no atom had a B factor >/=40, alternate conformation, or severe clash; that selection process yielded 1,490 Asn and 863 Gln side chains. Clear clustering was observed for Asn chi2 and Gln chi3 (except when Gln chi2 is trans). For Gln, five major and four minor rotamers cover 87% of examples. For Asn, there are seven backbone-independent rotamers covering 94% of examples plus rotamers specified for strictly alpha-helical, beta, and left-handed (+phi) Asn. Although the strongest influence on chi angles is avoidance of atomic clashes (especially with the NH2 hydrogens), some Asn or Gln rotamers are influenced by favorable van der Waals contacts and others by specific local H-bond patterns.

Asparagine↗

PubMed1998

Cooperativity in Bacillus stearothermophilus pyruvate kinase.

The enzyme pyruvate kinase (PK) from the moderate thermophile Bacillus stearothermophilus has been used as a model system with which to investigate the homotropic and heterotropic cooperative interactions of the enzyme. Cooperative ligand binding by the wild-type enzyme was measured using pre-steady-state and steady-state fluorescence spectroscopy, and steady-state kinetics. The results suggest that the cooperative structural changes induced by the substrate phosphoenolpyruvate (PEP) are distinct from those induced by the allosteric activator ribose- 5-phosphate (R5P). Furthermore the structural transition induced by the binding of saturating amounts of both PEP and R5P is itself distinct. This conclusion was further substantiated by the production of five mutant proteins in which the R5P- and PEP-induced homotropic cooperative transitions were separated. These results suggest that the cooperativity exhibited by pyruvate kinase from B. stearothermophilus does not conform to a simple two-state model. A putative four-state model is proposed.

Amino Acid Sequence↗

PubMed1992

Cyclic AMP potentiates receptor-stimulated phosphoinositide hydrolysis in human neuroepithelioma cells.

A stimulatory role for cAMP in the regulation of receptor-activated phosphoinositide hydrolysis has been examined in human SK-N-MCIXC and SK-N-MCIIE neuroepithelioma cells. The addition of optimal concentrations of oxotremorine-M, norepinephrine, endothelin-1, and ATP enhanced the release of inositol phosphates by 2-9-fold after activation of muscarinic, alpha 1-adrenergic, endothelin, and P2 nucleotide receptors, respectively. All combinations of these agonists elicited a release of inositol phosphates that was at least additive. However, the combined presence of oxotremorine-M and norepinephrine resulted in a phosphoinositide hydrolysis that was 30% greater than additive. This potentiation of inositol lipid hydrolysis resulted from an increased activity of the muscarinic receptor after the addition or norepinephrine and persisted after alpha 1-adrenergic receptor blockade. The enhancement of muscarinic receptor-stimulated inositol phosphate release could be quantitatively mimicked by inclusion of the beta-adrenergic agonist isoproterenol (EC50 approximately 0.1 microM), but not by alpha 1- or alpha 2-adrenergic agonists. Potentiation of oxotremorine-M-stimulated inositol lipid hydrolysis observed in the presence of either norepinephrine or isoproterenol was reduced in the absence of added Ca2+. Addition of either norepinephrine or isoproterenol to SK-N-MCIXC cells also resulted in a 16-fold increase in cAMP concentration. Although the cell-permeant 8-chloro-4-phenylthio-cAMP had a small inhibitory effect on basal inositol phosphate release, its inclusion resulted in a 19-31% enhancement of muscarinic, endothelin, ATP, and alpha 1-adrenergic receptor-stimulated phosphoinositide hydrolysis. We conclude 1) that, in SK-N-MCIXC cells, the addition of beta-adrenergic agonists selectively enhances muscarinic receptor-stimulated phosphoinositide hydrolysis through a cAMP-dependent process and 2) that the ability of exogenously added cAMP to enhance the activation of all four inositol lipid-linked receptors indicates that the effects of cAMP on inositol lipid hydrolysis are compartmentalized in these cells.

Adenosine Triphosphate↗