Search PubMed⌕ Search

Biomedical subjects

Michael S Chapman

Publications and source records attributed to Michael S Chapman.

At least 19 recordsLinked to original sources

Bias in cross-validated free R factors: mitigation of the effects of non-crystallographic symmetry.

Current methods of free R factor cross-validation assume that the structure factors of the test and working sets are independent of one another. This assumption is only an approximation when the modeled structure occupies anything less than the full asymmetric unit. Through progressive elimination of reflections from the working set, starting with those expected to be most correlated to the test set, small biases in free R can be measured, presumably because of over-sampling of the Fourier transform owing to bulk solvent in the crystal. This level of bias may be of little practical importance, but it rises to significant levels with increasing non-crystallographic symmetry owing to wider correlations between structure factors than hitherto appreciated. In the presence of 15-fold non-crystallographic symmetry, with resolutions commonly attainable in macromolecular crystallography, it may not be possible to calculate an unbiased free R factor. Methods are developed for the calculation of reduced-bias free R factors through elimination of the strongest correlations between test and working sets. With 180-fold non-crystallographic symmetry they may not be an accurate indicator of absolute quality, but they do yield the correct optimal weighting for stereochemical restraints.

Algorithms↗

Generalized anomeric interpretation of the "high-energy" N-P bond in N-methyl-N'-phosphorylguanidine: importance of reinforcing stereoelectronic effects in "high-energy" phosphoester bonds.

Electronic structure calculations have been performed on a model N-phosphorylguanidine, or phosphagen, to understand the stereoelectronic factors contributing to the lability of the "high-energy" N-P bond. The lability of the N-P bond is central to the physiological role of phosphagens involving phosphoryl transfer reactions important in cellular energy buffering and metabolism. Eight protonated forms of N-methyl-N'-phosphorylguanidine have been energy minimized at levels of theory ranging up to B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) to investigate the correlation between protonation state and N-P bond length. Selected forms have also been minimized using the CCSD/6-311++G(d,p) and QCISD/6-311++G(d,p) levels of theory. Bulk solvation energies using the polarized continuum model (PCM) with B3LYP/6-311++G(d,p) test the influence of the surroundings on computed structures and energies. The N-P bond length depends on the overall protonation state where increased protonation at the phosphoryl group or deprotonation at the unsubstituted N'' nitrogen results in shorter, stronger N-P bonds. Natural bond orbital analysis shows that the protonation state affects the N-P bond length by altering the magnitude of stabilizing n(O) --> sigma*(N-P) stereoelectronic interactions and to a lesser extent the sigma(N-P) --> sigma*(C-N'') and sigma(N-P) --> sigma*(C-N) interactions. The computations do not provide evidence of a competition between the phosphoryl and guanidinium groups for the same lone pair on the bridging nitrogen, as previously suggested by opposing resonance theory. The computed n(O) --> sigma*(N-P) anomeric effect provides a novel explanation of "high-energy" N-P bond lability. This offers new mechanistic insight into phosphoryl transfer reactions involving both phosphagens and other biochemically important "high-energy" phosphoester bonds.

Guanidines↗

Fitting of high-resolution structures into electron microscopy reconstruction images.

Dynamic macromolecular assemblies, such as ribosomes, viruses, and muscle protein complexes, are often more amenable to visualization by electron microscopy than by high-resolution X-ray crystallography or NMR. When high-resolution structures of component structures are available, it is possible to build an atomic model that gives information about the molecular interactions at greater detail than the experimental resolution, due to constraints of modeling placed upon the interpretation. There are now several competing computational methods to search systematically for orientations and positions of components that match the experimental image density, and continuing developments will be reviewed. Attention is now also moving toward the related task of optimization, with flexible and/or multifragment models and sometimes with stereochemically restrained refinement methods. This paper will review the various approaches and describe advances in the authors' methods and applications of real-space refinement.

Animals↗

Large-scale production, purification and crystallization of wild-type adeno-associated virus-2.

Adeno-associated virus-2 (AAV-2) has long been recognized as a potential vector for human gene therapy. Although much progress has been made in the molecular virology of AAV-2, structural studies of AAV-2 have been hampered by the low efficiency of virus production in culture, the low purity of preparations, and the low solubility of pure virus particles in solution. Methods of larger scale AAV-2 production have been developed through adaptation to suspension culture and re-optimization of the times of infection and transfection with respect to particle production. The methods allow the purification of 10mg ( approximately 10(15) particles) of AAV-2 per preparation at approximately 99% purity as judged by SDS-PAGE. This was sufficient for the screening of conditions for the formation of diffraction-grade crystals, ultimately leading to an atomic structure for AAV-2.

Crystallization↗

Impact of a Poisson-Boltzmann electrostatic restraint on protein structures refined at medium resolution.

The Poisson-Boltzmann formulism has been developed as a restraint for electrostatic interactions during the crystallographic refinement of macromolecules. It accounts implicitly for the effects of solvent and mobile ions, which are usually not included as restraints in the refinement of experimental structures. The electrostatic restraint has been implemented by combining software for numerically solving the three-dimensional Poisson-Boltzmann equation with a package for stereochemically restrained refinement. Its application to medium-resolution protein structures leads to a reduced free R factor, overfitting and to improved interactions in salt bridges and between polar or charged amino acids and the solvent. In contrast, Coulombic and screened Coulombic treatments did not lead to significant gains. The work leads to a modest improvement in refinement methods, confirmation that the Poisson-Boltzmann formulism is more consistent with experimental structure than the Coulombic approach, and to a reduction in the discrepancy between experimental and electrostatically optimized atomic models.

Crystallography↗

The active site cysteine of arginine kinase: structural and functional analysis of partially active mutants.

Arginine kinase buffers cellular ATP levels by catalyzing reversible phosphoryl transfer between ATP and arginine. A conserved cysteine has long been thought important in catalysis. Here, cysteine 271 of horseshoe crab arginine kinase has been mutated to serine, alanine, asparagine, or aspartate. Catalytic turnover rates were 0.02-1.0% of wild type, but the activity of uncharged mutations could be partially rescued with chloride. Steady-state binding constants were slightly increased, more so for phospho-L-arginine than ADP. Substrate binding synergy observed in many phosphagen kinases was reduced or eliminated in mutant enzymes. The crystallographic structure of the alanine mutant at 2.3 A resolution, determined as a transition state analogue complex with arginine, nitrate, and MgADP, was nearly identical to wild type. Enzyme-substrate interactions are maintained as in wild type, and substrates remain at least roughly aligned for in-line phosphoryl transfer. Homology models with serine, asparagine, or aspartate replacing the active site cysteine similarly show only minor structural changes. Most striking, however, is the presence in the C271A mutant crystallographic structure of a chloride ion within 3.5 A of the nonreactive N(eta) substrate nitrogen, approximating the position of the sulfur in the wild-type's cysteine. Together, the results contradict prevailing speculation that the cysteine mediates a substrate-induced conformational change, confirm that it is the thiolate form that is relevant to catalysis, and suggest that one of its roles is to help to enhance the catalytic rate through electrostatic stabilization of the transition state.

Arginine Kinase↗

Mathematical aspects of protein structure determination with NMR orientational restraints.

The field of structural biology is becoming increasingly important as new technological developments facilitate the collection of data on the atomic structures of proteins and nucleic acids. The solid-state NMR method is a relatively new biophysical technique that holds particular promise for determining the structures of peptides and proteins that are located within the cell membrane. This method provides information on the orientation of the peptide planes relative to an external magnetic field. In this article, we discuss some of the mathematical methods and tools that are useful in deriving the atomic structure from these orientational data. We first discuss how the data are viewed as tensors, and how these tensors can be used to construct an initial atomic model, assuming ideal stereochemistry. We then discuss methods for refining the models using global optimization, with stereochemistry constraints treated as penalty functions. These two processes, initial model building followed by refinement, are the two crucial steps between data collection and the final atomic model.

Magnetic Resonance Spectroscopy↗

The role of phosphagen specificity loops in arginine kinase.

Phosphagen kinases catalyze the reversible transfer of a phosphate between ATP and guanidino substrates, a reaction that is central to cellular energy homeostasis. Members of this conserved family include creatine and arginine kinases and have similar reaction mechanisms, but they have distinct specificities for different guanidino substrates. There has not been a full structural rationalization of specificity, but two loops have been implicated repeatedly. A small domain loop is of length that complements the size of the guanidino substrate, and is located where it could mediate a lock-and-key mechanism. The second loop contacts the substrate with a valine in the methyl-substituted guanidinium of creatine, and with a glutamate in the unsubstituted arginine substrate, leading to the proposal of a discriminating hydrophobic/hydrophilic minipocket. In the present work, chimeric mutants were constructed with creatine kinase loop elements inserted into arginine kinase. Contrary to the prior rationalizations of specificity, most had measurable arginine kinase activity but no creatine kinase activity or enhanced phosphocreatine binding. Guided by structure, additional mutations were introduced in each loop, recovering arginine kinase activities as high as 15% and 64% of wild type, respectively, even though little activity would be expected in the constructs if the implicated sites had dominant roles in specificity. An atomic structure of the mismatched complex of arginine kinase with creatine and ADP indicates that specificity can also be mediated by an active site that allows substrate prealignment that is optimal for reactivity only with cognate substrates and not with close homologs that bind but do not react.

Adenosine Diphosphate↗

Study of the structural dynamics of the E coli 70S ribosome using real-space refinement.

Cryo-EM density maps showing the 70S ribosome of E. coli in two different functional states related by a ratchet-like motion were analyzed using real-space refinement. Comparison of the two resulting atomic models shows that the ribosome changes from a compact structure to a looser one, coupled with the rearrangement of many of the proteins. Furthermore, in contrast to the unchanged inter-subunit bridges formed wholly by RNA, the bridges involving proteins undergo large conformational changes following the ratchet-like motion, suggesting an important role of ribosomal proteins in facilitating the dynamics of translation.

Bacterial Proteins↗

Structure determination of adeno-associated virus 2: three complete virus particles per asymmetric unit.

The atomic structure of adeno-associated virus 2 (AAV-2) has been determined to 3.0 A resolution. AAV-2 crystallized in space group P1, with unit-cell parameters a = 249.7, b = 249.7, c = 644.8 A, alpha = 90.0, beta = 101.2, gamma = 120.0 degrees. The crystals contained three full virus particles in the asymmetric unit, allowing 180-fold non-crystallographic symmetry averaging. The particle orientations were determined using the self-rotation function and found to have similar but resolvably different orientations. Approximate alignment of icosahedral and interparticle threefold screw symmetry led to a native Patterson that was interpretable in terms of approximate particle positions. Accurate positions required a Patterson correlation search that was constrained to be consistent with non-crystallographic threefold projection symmetry evident in the diffraction intensities. Initial phases to 15.0 A resolution were calculated by molecular replacement using the known structure of a distantly related homolog (23% sequence identity). Real-space averaging was performed and phases were extended from 15.0 to 3.0 A. An atomic model was fitted and refined using a simulated-annealing real-space procedure.

Crystallization↗

The putative catalytic bases have, at most, an accessory role in the mechanism of arginine kinase.

Arginine kinase is a member of the phosphagen kinase family that includes creatine kinase and likely shares a common reaction mechanism in catalyzing the buffering of cellular ATP energy levels. Abstraction of a proton from the substrate guanidinium by a catalytic base has long been thought to be an early mechanistic step. The structure of arginine kinase as a transition state analog complex (Zhou, G., Somasundaram, T., Blanc, E., Parthasarathy, G., Ellington, W. R., and Chapman, M. S. (1998) Proc. Natl. Acad. Sci. U. S. A. 95, 8449-8454) showed that Glu-225 and Glu-314 were the only potential catalytic residues contacting the phosphorylated nitrogen. In the present study, these residues were changed to Asp, Gln, and Val or Ala in several single and multisite mutant enzymes. These mutations had little impact on the substrate binding constants. The effect upon activity varied with reductions in kcat between 3000-fold and less than 2-fold. The retention of significant activity in some mutants contrasts with published studies of homologues and suggests that acid-base catalysis by these residues may enhance the rate but is not absolutely essential. Crystal structures of mutant enzymes E314D at 1.9 A and E225Q at 2.8 A resolution showed that the precise alignment of substrates is subtly distorted. Thus, pre-ordering of substrates might be just as important as acid-base chemistry, electrostatics, or other potential effects in the modest impact of these residues upon catalysis.

Amino Acid Substitution↗

Atomic refinement with correlated solid-state NMR restraints.

The orientation data provided by solid-state NMR can provide a great deal of structural information about membrane proteins. The quality of the information provided is, however, somewhat degraded by sign degeneracies in measurements of the dipolar coupling tensor. This is reflected in the dipolar coupling penalty function used in atomic refinement, which is less capable of properly restraining atoms when dipolar sign degeneracies are present. In this report we generate simulated solid-state NMR data using a variety of procedures, including back-calculation from crystal structures of alpha-helical and beta-sheet membrane proteins. We demonstrate that a large fraction of the dipolar sign degeneracies are resolved if anisotropic dipolar coupling measurements are correlated with anisotropic chemical shift measurements, and that all sign degeneracies can be resolved if three data types are correlated. The advantages of correlating data are demonstrated with atomic refinement of two test membrane proteins. When refinement is performed using correlated dipolar couplings and chemical shifts, perturbed structures converge to conformations with a larger fraction of correct dipolar signs than when data are uncorrelated. In addition, the final structures are closer to the original unperturbed structures when correlated data are used in the refinement. Thus, refinement with correlated data leads to improved atomic structures. The software used to correlate dipolar coupling and chemical shift data and to set up energy functions and their derivatives for refinement, CNS-SS02, is available at our web site.

Algorithms↗

Induced fit in guanidino kinases--comparison of substrate-free and transition state analog structures of arginine kinase.

Arginine kinase (AK) is a member of the guanidino kinase family that plays an important role in buffering ATP concentration in cells with high and fluctuating energy demands. The AK specifically catalyzes the reversible phosphoryl transfer between ATP and arginine. We have determined the crystal structure of AK from the horseshoe crab (Limulus polyphemus) in its open (substrate-free) form. The final model has been refined at 2.35 A with a final R of 22.3% (R(free) = 23.7%). The structure of the open form is compared to the previously determined structure of the transition state analog complex in the closed form. Classically, the protein would be considered two domain, but dynamic domain (DynDom) analysis shows that most of the differences between the two structures can be considered as the motion between four rigid groups of nonsequential residues. ATP binds near a cluster of positively charged residues of a fixed dynamic domain. The other three dynamic domains close the active site with separate hinge rotations relative to the fixed domain. Several residues of key importance for the induced motion are conserved within the phosphagen kinase family, including creatine kinase. Substantial conformational changes are induced in different parts of the enzyme as intimate interactions are formed with both substrates. Thus, although induced fit occurs in a number of phosphoryl transfer enzymes, the conformational changes in phosphagen kinases appear to be more complicated than in prior examples.

Adenosine Diphosphate↗

Refinement of the arginine kinase transition-state analogue complex at 1.2 A resolution: mechanistic insights.

The three-dimensional crystal structure of an arginine kinase transition-state analogue complex has been refined at 1.2 A resolution, with an overall R factor of 12.3%. The current model provides a unique opportunity to analyze the structure of a bimolecular (phosphagen kinase) enzyme in its transition state. This atomic resolution structure confirms in-line transfer of the phosphoryl group and the catalytic importance of the precise alignment of the substrates. The structure is consistent with a concerted proton transfer that has been proposed for an unrelated kinase. Refinement of anisotropic temperature factors and translation-libration-screw (TLS) analyses led to the identification of four rigid groups and their prevalent modes of motion in the transition state. The relative magnitudes of the mobility of rigid groups are consistent with their proposed roles in catalysis.

Animals↗

The atomic structure of adeno-associated virus (AAV-2), a vector for human gene therapy.

The structure of the adeno-associated virus (AAV-2) has been determined to 3-A resolution by x-ray crystallography. AAV is being developed as a vector for gene therapy to treat diseases including hemophilia, cancer, and cystic fibrosis. As in the distantly related autonomous parvoviruses, the capsid protein has a beta-barrel fold, but long loops between the beta-strands share little structural homology with other parvoviruses, leading to unique surface features. Most prominent are groups of threefold-related peaks, each an intimate association of loops from two neighboring subunits. Mutations affecting cell entry and receptor binding are clustered near the positively charged side of each peak, implicating the region in attachment to the cellular receptor, heparan sulfate proteoglycan. Amino acids involved in antibody binding are in the same general vicinity. The structure will guide rational engineering of vector capsids to tailor cellular targeting and to avoid immediate neutralization by an immune system sensitized by prior exposure to AAV.

Amino Acid Sequence↗

Simulated-annealing real-space refinement as a tool in model building.

Methods have been developed that further automate the building of macromolecular structures into electron-density maps. The software supports molecular-dynamics real-space refinement of an atomic model to local regions of a map within the context of a popular molecular-modeling program, O [Jones et al. (1991), Acta Cryst. A47, 110-119]. It is implemented as a module to the CNS refinement package [Brünger et al. (1998), Acta Cryst. D54, 905-921], controlled by a graphical user interface and commands executed directly through the molecular-graphics package. The method is illustrated with examples of the building and rebuilding of protein and nucleic acid models in which laborious manual adjustments are avoided. The resulting models show improved convergence during subsequent reciprocal-space refinement. The novel feature of the RSRef2000 software is the combination of simulated-annealing optimization with local real-space refinement, allowing several local minima to be explored quickly and automatically within the context of interactive model building.

Arginine Kinase↗