A smooth and differentiable bulk-solvent model for macromolecular diffraction.

Inclusion of low-resolution data in macromolecular crystallography requires a model for the bulk solvent. Previous methods have used a binary mask to accomplish this, which has proven to be very effective, but the mask is discontinuous at the solute-solvent boundary (i.e. the mask value jumps from zero to one) and is not differentiable with respect to atomic parameters. Here, two algorithms are introduced for computing bulk-solvent models using either a polynomial switch or a smoothly thresholded product of Gaussians, and both models are shown to be efficient and differentiable with respect to atomic coordinates. These alternative bulk-solvent models offer algorithmic improvements, while showing similar agreement of the model with the observed amplitudes relative to the binary model as monitored using R, R(free) and differences between experimental and model phases. As with the standard solvent models, the alternative models improve the agreement primarily with lower resolution (>6 A) data versus no bulk solvent. The models are easily implemented into crystallographic software packages and can be used as a general method for bulk-solvent correction in macromolecular crystallography.

Crystal structure of the cytosolic C2A-C2B domains of synaptotagmin III. Implications for Ca(+2)-independent snare complex interaction.

Synaptotagmins are synaptic vesicle-associated, phospholipid-binding proteins most commonly associated with Ca(+2)-dependent exocytotic and Ca(+2)- independent endocytotic events. Synaptotagmin III is a 63.2-kD member of the synaptotagmin homology group; one of its characteristic properties is the ability to bind divalent cations and accessory proteins promiscuously. In the cytosolic portion of this protein, a flexible seven-amino acid linker joins two homologous C2 domains. The C2A domain binds to phospholipid membranes and other accessory proteins in a divalent cation-dependent fashion. The C2B domain promotes binding to other C2B domains, as well as accessory proteins independent of divalent cations. The 3.2 A crystal structure of synaptotagmin III, residues 295-566, which includes the C2A and C2B domains, exhibits differences in the shape of the Ca(+2)-binding pocket, the electrostatic surface potential, and the stoichiometry of bound divalent cations for the two domains. These observations may explain the disparate binding properties of the two domains. The C2A and the C2B domains do not interact; synaptotagmin, therefore, covalently links two independent C2 domains, each with potentially different binding partners. A model of synaptotagmin's involvement in Ca(+2)-dependent regulation of membrane fusion through its interaction with the SNARE complex is presented.

Epsin 1 undergoes nucleocytosolic shuttling and its eps15 interactor NH(2)-terminal homology (ENTH) domain, structurally similar to Armadillo and HEAT repeats, interacts with the transcription factor promyelocytic leukemia Zn(2)+ finger protein (PLZF).

Epsin (Eps15 interactor) is a cytosolic protein involved in clathrin-mediated endocytosis via its direct interactions with clathrin, the clathrin adaptor AP-2, and Eps15. The NH(2)-terminal portion of epsin contains a phylogenetically conserved module of unknown function, known as the ENTH domain (epsin NH(2)-terminal homology domain). We have now solved the crystal structure of rat epsin 1 ENTH domain to 1.8 A resolution. This domain is structurally similar to armadillo and Heat repeats of beta-catenin and karyopherin-beta, respectively. We have also identified and characterized the interaction of epsin 1, via the ENTH domain, with the transcription factor promyelocytic leukemia Zn(2)+ finger protein (PLZF). Leptomycin B, an antifungal antibiotic, which inhibits the Crm1- dependent nuclear export pathway, induces an accumulation of epsin 1 in the nucleus. These findings suggest that epsin 1 may function in a signaling pathway connecting the endocytic machinery to the regulation of nuclear function.

Crystallographic R factor refinement by molecular dynamics.

Molecular dynamics was used to refine macromolecular structures by incorporating the difference between the observed crystallographic structure factor amplitude and that calculated from an assumed atomic model into the total energy of the system. The method has a radius of convergence that is larger than that of conventional restrained least-squares refinement. Test cases showed that the need for manual corrections during refinement of macromolecular crystal structures is reduced. In crambin, the dynamics calculation moved residues that were misplaced by more than 3 angstroms into the correct positions without human intervention.

Assessing the quality of solution nuclear magnetic resonance structures by complete cross-validation.

Structure determination of macromolecules in solution by nuclear magnetic resonance (NMR) spectroscopy involves the fitting of atomic models to the observed nuclear Overhauser effect (NOE) data. Complete cross-validation has been used to define reliable and unbiased criteria for the quality of solution NMR structures. The method is based on the partitioning of NOE data into test sets and the cross-validation of statistical quantities for each of the test sets. A high correlation between cross-validated measures of fit, such as distance bound violations and NMR R values, and the quality of solution NMR structures was observed. Less complete data resulted in poorer satisfaction of the cross-validated measures of fit. Optimization of cross-validated measures of fit will likely produce solution NMR structures with maximal information content.

Direct observation of protein solvation and discrete disorder with experimental crystallographic phases.

A complete and accurate set of experimental crystallographic phases to a resolution of 1.8 angstroms was obtained for a 230-residue dimeric fragment of rat mannose-binding protein A with the use of multiwavelength anomalous dispersion (MAD) phasing. An accurate image of the crystal structure could thus be obtained without resort to phases calculated from a model. Partially reduced disulfide bonds, local disorder, and differences in the mobility of chemically equivalent molecules are apparent in the experimental electron density map. A solvation layer is visible that includes well-ordered sites of hydration around polar and charged protein atoms, as well as diffuse, partially disordered solvent shells around exposed hydrophobic groups. Because the experimental phases and the resulting electron density map are free from the influence of a model, they provide a stringent test of theoretical models of macromolecular solvation, motion, and conformational heterogeneity.

Solution of a protein crystal structure with a model obtained from NMR interproton distance restraints.

Model calculations were performed to test the possibility of solving crystal structures of proteins by Patterson search techniques with three-dimensional structures obtained from nuclear magnetic resonance (NMR) interproton distance restraints. Structures for crambin obtained from simulated NMR data were used as the test system; the root-mean-square deviations of the NMR structures from the x-ray structure were 1.5 to 2.2 A for backbone atoms and 2.0 to 2.8 A for side-chain atoms. Patterson searches were made to determine the orientation and position of the NMR structures in the unit cell. The correct solution was obtained by comparing the rotation function results of several of the NMR structures and the average structure derived from them. Conventional refinement techniques reduced the R factor from 0.43 at 4 A resolution to 0.27 at 2 A resolution without inclusion of water molecules. The partially refined structure has root-mean-square backbone and side-chain atom deviations from the x-ray structure of 0.5 and 1.3 A, respectively.

Three-dimensional structure of an angiotensin II-Fab complex at 3 A: hormone recognition by an anti-idiotypic antibody.

The elucidation of bioactive conformations of small peptide hormones remains an elusive goal to structural chemists because of the inherent flexibility of these molecules. Angiotensin II (AII), the major effector of the renin-angiotensin system, is an octapeptide hormone for which no clear structural models exist. Peptide hormones such as AII share the property that they bind to their receptors with high affinities, in spite of the fact that they must overcome an extremely large conformational entropy barrier to bind in one conformation. A "surrogate system" that consists of a high-affinity monoclonal antibody (MAb) and AII has been used to study a bound conformation of AII. The crystallographic structure of the complex reveals a structure of AII that is compatible with predicted bioactive conformations of AII derived from structure-activity studies and theoretical calculations. In the complex, the deeply bound hormone is folded into a compact structure in which two turns bring the amino and carboxyl termini close together. The antibody of this complex (MAb 131) has the unusual property that it was not generated against AII, but rather against an anti-idiotypic antibody reactive with a MAb to AII, which renders this antibody an anti-anti-idiotypic antibody. The high affinity for AII of the original MAb to AII was passed on to MAb 131 through a structural determinant on the anti-idiotypic antibody. Strikingly, the conformation of AII in this complex is highly similar to complementarity determining region loops of antibodies, possibly indicating that a true molecular mimic of bound AII was present on the anti-idiotypic antibody against which MAb 131 was elicited.

Structural insights into the molecular mechanism of calcium-dependent vesicle-membrane fusion.

The fusion of vesicles with target membranes is controlled by a complex network of protein-protein and protein-lipid interactions. Recently determined structures of the SNARE complex, synaptotagmin III, nSec1, domains of the NSF chaperone and its adaptor (SNAP), and Rab3 and some of its effectors provide the framework for developing molecular models of vesicle fusion and for designing experiments to test these models. Ultimately, knowledge of the structures of higher-order complexes and their dynamic behavior will be required to obtain a full understanding of the vesicle fusion protein machinery.

Recent developments for the efficient crystallographic refinement of macromolecular structures.

Macromolecular crystallographic refinement has recently been made more efficient by the use of cross-validated maximum likelihood targets and torsion-angle molecular dynamics simulated annealing. In combination with automated model building methods, the amount of manual intervention required to complete and refine a structure is greatly reduced.

Structural insights into the molecular mechanism of Ca(2+)-dependent exocytosis.

The fusion of vesicles with target membranes is controlled by a complex network of protein-protein and protein-lipid interactions. Recent structures of the SNARE complex, synaptotagmin III, nSec1, domains of NSF and its adaptor SNAP, along with Rab3 and some of its effectors, provide the framework for developing molecular models of vesicle fusion and for designing experiments to test these models. Ultimately, this knowledge of the structures of higher-order complexes and their dynamic behavior will allow us to obtain a full understanding of the vesicle fusion protein machinery.

Human ornithine aminotransferase complexed with L-canaline and gabaculine: structural basis for substrate recognition.

BACKGROUND: Ornithine aminotransferase (OAT) is a 45 kDa pyridoxal-5'-phosphate (PLP)-dependent enzyme that catalyzes the conversion of L-ornithine and 2-oxoglutarate to glutamate-delta-semialdehyde and glutamic acid, respectively. In humans, loss of OAT function causes an accumulation of ornithine that results in gyrate atrophy of the choroid and retina, a disease that progressively leads to blindness. In an effort to learn more about the structural basis of this enzyme's function, we have determined the X-ray structures of OAT in complex with two enzyme-activated suicide substrates: L-canaline, an ornithine analog, and gabaculine, an irreversible inhibitor of several related aminotransferases. RESULTS: The structures of human OAT bound to the inhibitors gabaculine and L-canaline were solved to 2.3 A at 110K by difference Fourier techniques. Both inhibitors coordinate similarly in the active site, binding covalently to the PLP cofactor and causing a 20 degrees rotation in the cofactor tilt relative to the ligand-free form. Aromatic-aromatic interactions occur between the bound gabaculine molecule and active-site residues Tyr85 and Phe177, whereas Tyr55 and Arg180 provide specific contacts to the alpha-amino and carboxyl groups of L-canaline. CONCLUSIONS: The OAT-L-canaline complex structure implicates Tyr55 and Arg180 as the residues involved in coordinating with the natural substrate ornithine during normal enzyme turnover. This correlates well with two enzyme-inactivating point mutations associated with gyrate atrophy, Tyr55-->His and Arg180-->Thr. The OAT-gabaculine complex provides the first structural evidence that the potency of the inhibitor is due to energetically favourable aromatic interactions with residues in the active site. This aromatic-binding mode may be relevant to structure-based drug design efforts against other omega-aminotransferase targets, such as GABA aminotransferase.

Dramatic structural and thermodynamic consequences of repacking a protein's hydrophobic core.

BACKGROUND: Rop is an RNA binding, dimeric, four-helix bundle protein with a well-defined, regular hydrophobic core ideally suited for redesign studies. A family of Rop variants in which the hydrophobic core was systematically redesigned has previously been created and characterized. RESULTS: We present a structural and thermodynamic analysis of Ala2Ile2-6, a variant of Rop with an extensively redesigned hydrophobic core. The structure of Ala2Ile2-6 reveals a completely new fold formed by a conformational "flip" of the two protomers around the dimeric interface. The free-energy profile of Ala2Ile2-6 is also very different from that of wild-type Rop. Ala2Ile2-6 has a higher melting temperature than Rop, but undergoes a slightly smaller free-energy change on unfolding. CONCLUSIONS: The structure of Ala2Ile2-6, along with molecular modeling results, demonstrate the importance of tight packing of core residues and the adoption of favorable core side chain rotamer values in determining helix-helix interactions in the four-helix bundle fold. Structural disorder at the N and C termini of Ala2Ile2-6 provides a basis for the large differences in the enthalpy and entropy of Ala2Ile2-6 folding compared with wildtype Rop.

Statistical analysis of predicted transmembrane alpha-helices.

Statistical analyses were undertaken for putative transmembrane alpha-helices obtained from a database representing the subset of membrane proteins available in Swiss-Prot. The average length of a transmembrane alpha-helix was found to be 22-21 amino acids with a large variation around the mean. The transfer free energy from water to oil of a transmembrane alpha-helix in bitopic proteins, -48 kcal/mol, is higher than that in polytopic proteins, -39 kcal/mol, and is nearly identical to that obtained by assuming a random distribution of solely hydrophobic amino acids in the alpha-helix. The amino acid composition of hydrophobic residues is similar in bitopic and polytopic proteins. In contrast, the more polar the amino acids are, the less likely they are to be found in bitopic proteins compared to polytopic ones. This most likely reflects the ability of alpha-helical bundles to shield the polarity of residues from the hydrophobic bilayer. One half of all amino acids were distributed nonrandomly in both bitopic and polytopic proteins. A preference was found for tyrosine and tryptophan residues to be at the ends of transmembrane alpha-helices. Correlated distribution analysis of amino acid pairs indicated that most amino acids are independently distributed in each helix. Exceptions are cysteine, tyrosine, and tryptophan which appear to cluster closely to one another and glycines which are preferentially found on the same side of alpha-helices.

Helicity, membrane incorporation, orientation and thermal stability of the large conductance mechanosensitive ion channel from E. coli.

In this report, we present structural studies on the large conductance mechanosensitive ion channel (MscL) from E. coli in detergent micelles and lipid vesicles. Both transmission Fourier transform infrared spectroscopy and circular dichroism (CD) spectra indicate that the protein is highly helical in detergents as well as liposomes. The secondary structure of the proteins was shown to be highly resistant towards denaturation (25-95 degrees C) based on an ellipticity thermal profile. Amide H+/D+ exchange was shown to be extensive (ca. 66%), implying that two thirds of the protein are water accessible. MscL, reconstituted in oriented lipid bilayers, was shown to possess a net bilayer orientation using dichroic ratios measured by attenuated total-reflection Fourier transform infrared spectroscopy. Here, we present and discuss this initial set of structural data on this new family of ion-channel proteins.

Crystallographic refinement by simulated annealing. Application to a 2.8 A resolution structure of aspartate aminotransferase.

Crystallographic refinement by simulated annealing with molecular dynamics has been applied to a 2.8 A (1 A = 0.1 nm) resolution X-ray structure of aspartate aminotransferase. Comparison of the refined structure and a structure obtained by combined restrained least-squares refinement and manual re-fitting shows a similar R factor, stereochemistry, and mean difference from the isomorphous replacement phase centroids. Crystallographic refinement by simulated annealing accomplished structural changes and improvements of the electron density maps that were not possible by using restrained least-squares refinement without manual re-fitting. Crystallographic refinement by simulated annealing can generate an ensemble of structures, each of which agrees with the diffraction information. Regions of large variations of the ensemble indicate either erroneously fitted or disordered segments of the macromolecule.

The 1.8 A crystal structure of a statically disordered 17 base-pair RNA duplex: principles of RNA crystal packing and its effect on nucleic acid structure.

The crystal structure of a 17 base RNA oligomer, r(CACCGGAUG GUUCGGUG), has been solved to a resolution of 1.8 A through a combination of molecular replacement, multiple isomorphous replacement phasing, and analysis of observed intensity distributions. The oligomer, which forms a stem-loop in solution, crystallized as a pseudo-infinite duplex in spacegroup P321. The asymmetric unit of the crystal contains four superimposed orientations of the duplex that are out of register, such that backbones superimpose, but base identity differs. This static disorder was initially discovered by brominating a single residue per strand in the sequence, and observing four peaks per strand in difference maps phased with a native molecular replacement solution. The presence of four superimposed duplex "motifs" related by non-crystallographic hypersymmetry was detected by computing /2 and Wilson ratios for the observed intensities. The observed ratios matched those produced from calculated intensities of a 4-fold statically disordered model. Multi-conformer simulated annealing refinement against a maximum-likelihood target incorporating experimental phase information was used to refine the 4-fold disordered model to an Rfree and R of 29.35% and 25.5%, respectively. The resulting structure reveals four distinct conformations of the duplex, with an average pairwise backbone rmsd of 2.35 A. The structural differences between the four conformations, which can be attributed to differences in packing environment, highlight the possible influence of crystal packing forces on nucleic acid X-ray structures. Analysis of inter-helical packing between symmetry-related molecules reveals an RNA "zipper" that mediates direct phosphate oxygen-2' hydroxyl interactions between close-packed phosphate-sugar backbones. This may be a general mode for RNA tertiary interaction that does not depend on metal ions or primary sequence.

Conformational variability of solution nuclear magnetic resonance structures.

In structure determination by X-ray crystallography and solution NMR spectroscopy, experimental data are collected as time and ensemble-averages. Thus, in principle, appropriate time and ensemble-averaged models should be used. Refinement of an ensemble of conformers rather than one single structure against the experimental NMR data could, however, result in overfitting the data because of the significantly increased number of parameters. To avoid overfitting, complete cross-validation, which provides an unbiased measure of the fit, has been applied to nuclear Overhauser effect derived distance refinement. Using two synthetic test cases, a correlation was demonstrated between the cross-validated measure to the fit (defined in terms of root-mean-square deviations from the distance restraints and number of violations) and the number of models that best reproduce the conformational variability in solution. A new method, based on a probability map, has been used to generate good representations of the resulting ensembles of structures. The method has also been applied to observed NMR data for two proteins, interleukin 4 and interleukin 8. For interleukin 4, cross-validation indicates that a single-conformer model gives the most accurate representation of the structure, whereas conventional measures of fit between the experimental data and those calculated from the model decrease when increasing the number of conformers, indicating overfitting. For interleukin 8, complete cross-validation predicts a twin-conformer model to be the most faithful representation of the experimental data. Two distinct conformations for the loop formed by residues 16 to 22 emerge from the family of twin-conformer structures. The putative alternate conformation of the loop is not observed in the crystal structure of interleukin 8. However, because of crystal packing contacts in this region this does not necessarily exclude the presence of the alternate conformation in solution. The twin-conformer model is supported by observed chemical exchange line broadening for the amide of His18 obtained by 15N relaxation studies. This region has also been implied to be involved in receptor binding.

Crystal structures of a Rab protein in its inactive and active conformations.

We have determined crystal structures of Sec4, a member of the Rab family in the G protein superfamily, in two states: bound to GDP, and to a non-hydrolyzable GTP analog, guanosine-5'-(beta, gamma)-imidotriphosphate (GppNHp). This represents the first structure of a Rab protein bound to GDP. Sec4 in both states grossly resembles other G proteins bound to GDP and GppNHp. In Sec4-GppNHp, structural features common to active Rab proteins are observed. In Sec4-GDP, the switch I region is highly disordered and displaced relative to the switch I region of Ras-GDP. In two of the four molecules of Sec4-GDP in the asymmetric unit of the Sec4-GDP crystals, the switch II region adopts a conformation similar to that seen in the structure of the small G protein Ran bound to GDP. This allows residues threonine 76, glutamate 80, and arginine 81 of Sec4 to make contacts with other conserved residues and water molecules important for nucleotide binding. In the other two molecules in the asymmetric unit, these interactions do not take place. This structural variability in both the switch I and switch II regions of GDP-bound Sec4 provides a possible explanation for the high off-rate of GDP bound to Sec4, and suggests a mechanism for regulation of the GTPase cycle of Rab proteins by GDI proteins.

The 1.0 A crystal structure of Ca(2+)-bound calmodulin: an analysis of disorder and implications for functionally relevant plasticity.

Calmodulin (CaM) is a highly conserved 17 kDa eukaryotic protein that can bind specifically to over 100 protein targets in response to a Ca(2+) signal. Ca(2+)-CaM requires a considerable degree of structural plasticity to accomplish this physiological role; however, the nature and extent of this plasticity remain poorly characterized. Here, we present the 1.0 A crystal structure of Paramecium tetraurelia Ca(2+)-CaM, including 36 discretely disordered residues and a fifth Ca(2+) that mediates a crystal contact. The 36 discretely disordered residues are located primarily in the central helix and the two hydrophobic binding pockets, and reveal correlated side-chain disorder that may assist target-specific deformation of the binding pockets. Evidence of domain displacements and discrete backbone disorder is provided by translation-libration-screw (TLS) analysis and multiconformer models of protein disorder, respectively. In total, the evidence for disorder at every accessible length-scale in Ca(2+)-CaM suggests that the protein occupies a large number of hierarchically arranged conformational substates in the crystalline environment and may sample a quasi-continuous spectrum of conformations in solution. Therefore, we propose that the functionally distinct forms of CaM are less structurally distinct than previously believed, and that the different activities of CaM in response to Ca(2+) may result primarily from Ca(2+)-mediated alterations in the dynamics of the protein.