Variability of conformations at crystal contacts in BPTI represent true low-energy structures: correspondence among lattice packing and molecular dynamics structures.

The structures of five basic pancreatic trypsin inhibitor (BPTI) molecules are compared to establish the extent and nature of the conformational variability resulting from crystal packing effects. BPTI is an ideal system to evaluate such factors because of the availability of high resolution X-ray models of five different BPTI structures, each in a different crystal packing environment. Differences observed among the structures are found to be distributed throughout the molecule, although the regions that display most variability are associated with the loop structures (residues 14-17 and 24-29). The regions of structure that show the largest rms deviations from the mean of the five packing motifs correlate well with the presence of intermolecular contacts in the crystal lattice. For most of the molecules there is also a correspondence between a larger number of intermolecular contacts and systematically higher B-factors, although it is not apparent whether this is induced by the crystal contact or results from the fact that the contacts are made predominantly through surface loops. The conformational differences seen among the X-ray models constitute more than local shifts at the lattice contact surfaces, and in fact involve in some cases the making and breaking of intramolecular H-bonds. The magnitudes of the differences among packing models are significantly larger than those usually associated with changes induced by mutagenesis; for instance; the structural differences at the site of mutation observed on removing an internal disulfide from the molecule are significantly less than those associated with lattice contact effects. The crystal packing conformations are compared to representative structures of BPTI generated during a 96-psec molecular dynamics (MD) simulation. This comparison shows a high level of correspondence between the protein flexibility indicated by the X-ray and MD analyses, and specifically between those regions that are most variable. This suggests that the regions that show most variability among the crystal packing models are not artifacts of crystallization, but rather represent true low-energy conformers that have been preferentially selected by crystallization factors.     

Efficient conformational sampling of local side-chain flexibility

Side-chain flexibility of ligand-binding sites needs to be considered in the rational design of novel inhibitors. We have developed a method to generate conformational ensembles that efficiently sample local side-chain flexibility from a single crystal structure. The rotamer-based approach is tested here for the S1' pocket of human collagenase-1 (MMP-1), which is known to undergo conformational changes in multiple side-chains upon binding of certain inhibitors. First, a raw ensemble consisting of a large number of conformers of the S1' pocket was generated using an exhaustive search of rotamer combinations on a template crystal structure. A combination of principal component analysis and fuzzy clustering was then employed to successfully identify a core ensemble consisting of a low number of representatives from the raw ensemble. The core ensemble contained geometrically diverse conformers of stable nature, as indicated in several cases by a relative energy lower than that of the minimised template crystal structure. Through comparisons with X-ray crystallography and NMR structural data we show that the core ensemble occupied a conformational space similar to that observed under experimental conditions. The synthetic inhibitor RS-104966 is known to induce a conformational change in the side-chains of the S1' pocket of MMP-1 and could not be docked in the template crystal structure. However, the experimental binding mode was reproduced successfully using members of the core ensemble as the docking target, establishing the usefulness of the method in drug design.     

Influence of counter-ions on the crystal structures of DNA decamers: binding of [Co(NH3)6]3+ and Ba2+ to A-DNA.

A-DNA is a stable alternative right-handed double helix that is favored by certain sequences (e.g., (dG)n.(dC)n) or under low humidity conditions. Earlier A-DNA structures of several DNA oligonucleotides and RNA.DNA chimeras have revealed some conformational variation that may be the result of sequence-dependent effects or crystal packing forces. In this study, four crystal structures of three decamer oligonucleotides, d(ACCGGCCGGT), d(ACCCGCGGGT), and r(GC)d(GTATACGC) in two crystal forms (either the P6(1)22 or the P2(1)2(1)2(1) space group) have been analyzed at high resolution to provide the molecular basis of the structural difference in an experimentally consistent manner. The study reveals that molecules crystallized in the same space group have a more similar A-DNA conformation, whereas the same molecule crystallized in different space groups has different (local) conformations. This suggests that even though the local structure is influenced by the crystal packing environments, the DNA molecule adjusts to adopt an overall conformation close to canonical A-DNA. For example, the six independent CpG steps in these four structures have different base-base stacking patterns, with their helical twist angles (omega) ranging from 28 degrees to 37 degrees. Our study further reveals the structural impact of different counter-ions on the A-DNA conformers. [Co(NH3)6]3+ has three unique A-DNA binding modes. One binds at the major groove side of a GpG step at the O6/N7 sites of guanine bases via hydrogen bonds. The other two modes involve the binding of ions to phosphates, either bridging across the narrow major groove or binding between two intra-strand adjacent phosphates. Those interactions may explain the recent spectroscopic and NMR observations that [Co(NH3)6]3+ is effective in inducing the B- to A-DNA transition for DNA with (G)n sequence. Interestingly, Ba2+ binds to the same O6/N7 sites on guanine by direct coordinations.     

Comparison of the solution structures of angiotensin I & II. Implication for structure-function relationship.

Conformational analysis of angiotensin I (AI) and II (AII) peptides has been performed through 2D 1H-NMR spectroscopy in dimethylsulfoxide and 2,2,2-trifluoroethanol/H2O. The solution structural models of AI and AII have been determined in dimethylsulfoxide using NOE distance and 3JHNHalpha coupling constants. Finally, the AI family of models resulting from restrained energy minimization (REM) refinement, exhibits pairwise rmsd values for the family ensemble 0.26 +/- 0.13 A, 1.05 +/- 0.23 A, for backbone and heavy atoms, respectively, and the distance penalty function is calculated at 0.075 +/- 0.006 A2. Comparable results have been afforded for AII ensemble (rmsd values 0.30 +/- 0.22 A, 1.38 +/- 0.48 A for backbone and heavy atoms, respectively; distance penalty function is 0.029 +/- 0.003 A2). The two peptides demonstrate similar N-terminal and different C-terminal conformation as a consequence of the presence/absence of the His9-Leu10 dipeptide, which plays an important role in the different biological function of the two peptides. Other conformational variations focused on the side-chain orientation of aromatic residues, which constitute a biologically relevant hydrophobic core and whose inter-residue contacts are strong in dimethylsulfoxide and are retained even in mixed organic-aqueous media. Detailed analysis of the peptide structural features attempts to elucidate the conformational role of the C-terminal dipeptide to the different binding affinity of AI and AII towards the AT1 receptor and sets the basis for understanding the factors that might govern free- or bound-depended AII structural differentiation.     

Flexibility and structural conservation in a c-KIT G-quadruplex

A quadruplex sequence from the promoter region of the c-KIT gene forms a stable quadruplex, as characterized by crystallographic and NMR methods. Two new crystal structures are reported here, together with molecular dynamics simulation studies on these quadruplex crystal structures and an NMR structure. The new crystal structures, each in a distinct space group and lattice packing arrangement, together with the existing structures, demonstrate that the c-KIT quadruplex fold does not change with differing environments, suggesting that quadruplex topological dynamism is not a general phenomenon. The single and dinucleotide loops in these structures show a high degree of conformational flexibility within the three crystal forms and the NMR ensemble, with no evidence of clustering to particular conformers. This is in accord with the findings of high loop flexibility from the molecular dynamics studies. It is suggested that intramolecular quadruplexes can be grouped into two broad classes (i) those with at least one single-nucleotide loop, often showing singular topologies even though loops are highly flexible, and (ii) with all loops comprising at least two nucleotides, leading to topological dynamism. The loops can have more stable and less dynamic base-stacked secondary structures.     

Assessing the native state conformational distribution of ubiquitin by peptide acidity

At equilibrium, every energetically feasible conformation of a protein occurs with a non-zero probability. Quantitative analysis of protein flexibility is thus synonymous with determining the proper Boltzmann-weighting of this conformational distribution. The exchange reactivity of solvent-exposed amide hydrogens greatly varies with conformation, while the short-lived peptide anion intermediate implies an insensitivity to the dynamics of conformational motion. Amides that are well-exposed in model conformational ensembles of ubiquitin vary a million-fold in exchange rates which continuum dielectric methods can predict with an rmsd of 3. However, the exchange rates for many of the more rarely exposed amides are markedly overestimated in the PDB-deposited 2K39 and 2KN5 ubiquitin ensembles, while the 2NR2 ensemble predictions are largely consistent with those of the Boltzmann-weighted conformational distribution sampled at the level of 1%. The correlation between the fraction of solvent-accessible conformations for a given amide hydrogen and the exchange rate constant for that residue provides a useful monitor of the degree of completeness with which a given ensemble has sampled the energetically accessible conformational space. These exchange predictions correlate with the degree to which each ensemble deviates from a set of 46 ubiquitin X-ray structures. Kolmogorov-Smirnov analysis for the distribution of intra- and inter-ensemble pairwise structural rmsd values assisted the identification of a subensemble of 2K39 that eliminates the overestimations of hydrogen exchange rates observed for the full ensemble. The relative merits of incorporating experimental restraints into the conformational sampling process are compared to using these restraints as filters to select subpopulations consistent with the experimental data.     

Local RNA structural changes induced by crystallization are revealed by SHAPE

We present a simple approach to locate sites that undergo conformational changes upon crystallization by comparative structural mapping of the same RNA in three different environments. As a proof of principle, we probed the readily crystallized P4-P6DeltaC209 domain from the Tetrahymena thermophila group I intron in a native solution, in a solution mimicking the crystallization drop, and in crystals. We chose the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry, which monitors the flexibility and the conformation of each nucleotide. First, SHAPE successfully revealed the structural changes that occur during the crystallization process. Specifically, 64% of the nucleotides implicated in packing contacts and present in the portion of the molecule analyzed were identified. Second, reactivity differences for some of these nucleotides were already observed in the crystallization solution, suggesting that the crystallization buffer locked down a particular structure that was favorable to crystal formation. Third, the probing of a known structure extends our understanding of the structural basis for the SHAPE reaction by suggesting that reactivity is enhanced by a C2'-endo sugar pucker. Furthermore, by identifying local conformational changes of the RNA that take place during crystallization, SHAPE could be combined with the in vitro selection of stable mutants to rationalize the design of RNA candidates for crystallization.     

Crystallization of substrate and product analog complexes of tryptophanyl-tRNA synthetase

We have prepared crystals of tryptophanyl-tRNA synthetase from Bacillus stearothermophilus complexed to tryptophan (type II*), and to tryptophanyl-3'(2')-ATP (type IV). The latter compound is a product analog, enzymatically synthesized by acyl transfer of tryptophan from the tryptophanyl-5'-AMP intermediate to a second molecule of ATP. It resembles the 3'-terminal fragment, tryptophanyl-3'(2')-adenosine, of Trp-tRNATrp. Both crystal forms diffract to high resolution. Although both forms are grown from 2 M K2HPO4, they are dramatically different in the shape of the unit cell and in space group symmetry. Type II* crystals are monoclinic (space group P21). However, low-resolution reflections obey the symmetry of space group P321, which indicates both the existence and the location of noncrystallographic symmetry in the monoclinic unit cell. Type IV crystals belong to space group P41212 (or its enantiomorph) and the unit cell is elongated along the fourfold screw axis. Analysis of molecular packing suggests that intermolecular contacts in the two crystal types are very different. Thus, the two structures may exhibit conformational differences related to catalysis by this enzyme. Solution of type II* and type IV crystal structures may provide representations resembling a Michaelis complex and an acyl transfer product complex.     

Plasminogen activator inhibitor 1. Structure of the native serpin, comparison to its other conformers and implications for serpin inactivation

The crystal structure of a constitutively active multiple site mutant of plasminogen activator inhibitor 1 (PAI-1) was determined and refined at a resolution of 2.7 A.The present structure comprises a dimer of two crystallographically independent PAI-1 molecules that pack by association of the residues P6 to P3 of the reactive centre loop of one molecule (A) with the edge of the main beta-sheet A of the other molecule (B).Thus, the reactive centre loop is ordered for molecule A by crystal packing forces, while for molecule B it is unconstrained by crystal packing contacts and is disordered.The overall structure of active PAI-1 is similar to the structures of other active inhibitory serpins exhibiting as the major secondary structural feature a five-stranded beta-sheet A and an intact proteinase-binding loop protruding from the one end of the elongated molecule. No preinsertion of the reactive centre loop is observed in this structure.A comparison of the present structure with the previously determined crystal structures of PAI-1 in its alternative conformations reveals that, upon cleavage of an intact form of PAI-1 or formation of latent PAI-1, the well-characterised rearrangements of the serpin secondary structural elements are accompanied by dramatic and partly unexpected conformational changes of helical and loop structures proximal to beta-sheet A.The present structure explains the stabilising effects of the mutated residues, reveals the structural cause for the observed spectroscopic differences between active and latent PAI-1, and provides new insights into possible mechanisms of stabilisation by its natural binding partner, vitronectin.     

Packing analysis of carbohydrates and polysaccharides. IV. A new method for detailed crystal structure refinement of polysaccharides and its application to V-amylose

A method is described for predicting and solving crystal structures of linear homopolysaccharides. The method is based on the refinement of the structure with respect to either stereochemical constraints or x-ray diffraction intensities. In the refinement process, all conformational and packing features of the molecule, such as bond lengths, bond angles, conformational angles, nonbonded contacts, hydrogen bonds, etc., can be allowed to vary until the structure reaches both a conformation and crystalline packing that are in minimum disagreement with the stereochemical restraints and the diffraction data. In this fashion, both packing and conformational features of the structure can be simultaneously refined, and not separately as has been the custom in the past. The refinement procedure is based on a method of constrained optimization which possesses improved characteristics of reaching a solution and avoiding false minima, in comparison with least squares methods. The procedure is, in addition, capable of easily finding molecules of solvent of crystallization. The method was applied to further refining the previously solved crystal structure of V-amylose. The results indicated that contrary to the previously found six-fold molecular symmetry in the P2₁2₁2₁ space group, the V-amylose molecule exhibits only two-fold symmetry with the asymmetric unit consisting of three glucose residues in one-half turn of the helix. The three residues are nonequivalent principally due to unequal rotational positions of the hydroxymethyl groups. The crystal structure of V-amylose predicted from stereochemical refinement was identical in all details with that obtained from refining against X-ray data. The excellent agreement with the diffraction data was indicated by the crystallographic disagreement index R = 0.25.     

Non-sequential structure-based alignments reveal topology-independent core packing arrangements in proteins

MOTIVATION: Proteins of the same class often share a secondary structure packing arrangement but differ in how the secondary structure units are ordered in the sequence. We find that proteins that share a common core also share local sequence-structure similarities, and these can be exploited to align structures with different topologies. In this study, segments from a library of local sequence-structure alignments were assembled hierarchically, enforcing the compactness and conserved inter-residue contacts but not sequential ordering. Previous structure-based alignment methods often ignore sequence similarity, local structural equivalence and compactness. RESULTS: The new program, SCALI (Structural Core ALIgnment), can efficiently find conserved packing arrangements, even if they are non-sequentially ordered in space. SCALI alignments conserve remote sequence similarity and contain fewer alignment errors. Clustering of our pairwise non-sequential alignments shows that recurrent packing arrangements exist in topologically different structures. For example, the three-layer sandwich domain architecture may be divided into four structural subclasses based on internal packing arrangements. These subclasses represent an intermediate level of structure classification, more general than topology, but more specific than architecture as defined in CATH. A strategy is presented for developing a set of predictive hidden Markov models based on multiple SCALI alignments.     

Force field based conformational analysis of RNA structural motifs: GNRA tetraloops and their pyrimidine relatives

The protocol of conformational analysis applied here to ribonucleotide oligomers combines conformational search in the space of torsion angles and energy minimization using the AMBER4.1 force field with a continuum treatment of electrostatic solute-solvent interactions. RNA fragments with 5′-GGGCGNNAGCCU-3′ sequences commonly fold into hairpins with four-membered loops. The combinatorial search for acceptable conformations using the MC-SYM program was restricted to loop nucleotides and yielded roughly 1500 structures being compatible with a double-stranded stem. After energy minimization by the JUMNA program (without applying any experimental constraints), these structures converged into an ensemble of 74 different conformers including 26 structures which contained the sheared G-A base pair observed in experimental studies of GNRA tetraloops. Energetic analysis shows that inclusion of solvent electrostatic effects is critically important for the selection of conformers that agree with experimentally determined structures. The continuum model accounts for solvent polarization by means of the electrostatic reaction field. In the case of GNRA loop sequences, the contributions of the reaction field shift relative stabilities towards conformations showing most of the structural features derived from NMR studies. The agreement of computed conformations with the experimental structures of GAAA, GCAA, and GAGA tetraloops suggests that the continuum treatment of the solvent represents a definitive improvement over methods using simple damping models in electrostatic energy calculations. Application of the procedure described here to the evaluation of the relative stabilities of conformers resulting from searching the conformational space of RNA structural motifs provides some progress in (non-homology based) RNA 3D-structure prediction. Received: 20 January 1999 / Revised version: 4 June 1999 / Accepted: 10 June 1999     

Conformational ensemble of an intrinsically flexible loop in mitochondrial import protein Tim21 studied by modeling and molecular dynamics simulations

BackgroundTim21, a subunit of a highly dynamic translocase of the inner mitochondrial membrane (TIM23) complex, translocates proteins by interacting with subunits in the translocase of the outer membrane (TOM) complex and Tim23 channel in the TIM23 complex. A loop segment in Tim21, which is in close proximity of the binding site of Tim23, has different conformations in X-ray, NMR and new crystal contact-free space (CCFS) structures. MD simulations can provide information on the structure and dynamics of the loop in solution.MethodsThe conformational ensemble of the loop was characterized using loop modeling and molecular dynamics (MD) simulations.ResultsMD simulations confirmed mobility of the loop. Multidimensional scaling and clustering were used to characterize the dynamic conformational ensemble of the loop. Free energy landscape showed that the CCFS crystal structure occupied a low energy region as compared to the conventional X-ray crystal structure. Analysis of crystal packing indicates that the CCFS provides larger conformational space for the motions of the loop.ConclusionsOur work reported the conformational ensemble of the loop in solution, which is in agreement with the structure obtained from CCFS approach. The combination of the experimental techniques and computational methods is beneficial for studying highly flexible regions of proteins.General significanceComputational methods, such as loop modeling and MD simulations, have proved to be useful for studying conformational flexibility of proteins. These methods in integration with experimental techniques such as CCFS has the potential to transform the studies on flexible regions of proteins.     

Thermal movements of proteins: small-scale fluctuation and conformation substates

A comparison is made of the published data on different physical methods concerning small-scale (small-amplitude) thermal motions, conformers and conformational substates of myoglobin and some other proteins in crystal and solution. Because of the underestimation of the statical and dynamical rotational disorder of the crystal lattice, the X-ray diffraction method may somewhat overestimate the amplitude values of the intramolecular thermal motions, which leads to the underestimation of the extent of conformational rigidity of the protein molecule. In solution the rigid-body rotation of the macromolecules makes a significant contribution to the apparent amplitude of atomic displacement. The damping of these motions upon the increase of viscosity or freezing of solution, crystallization or intramolecular crosslinking can partially simulate the depression of intramolecular thermal motions. According to infrared spectroscopy data the crystallization of carboxymyoglobin as well as temperature decrease down to 5K are not accompanied by preferential stabilization of one of the several discrete conformers of the active site or by diminished distribution breadth of conformational substates within individual conformers. The structural differences between these discrete conformers and the more so between conformational substates are not sizable and are confined within a limit of the amplitude of thermal motion in the polymeric monocrystal including their paracrystalline regions. The dehydration of the hemoglobin and myoglobin films leads to the loss of conformational rigidity of the active site and to the distortion of its native conformation, i.e. to the inactivation of protein.     

High-resolution structure of the OmpA membrane domain

The membrane domain of OmpA consists of an eight-stranded all-next-neighbor antiparallel beta-barrel with short turns at the periplasmic barrel end and long flexible loops at the external end. The structure analysis has been extended from medium resolution to 1. 65 A (1 A=0.1 nm), and the molecular model has been refined anisotropically to show oriented mobilities of the structural elements. The improved data allowed us to locate five further detergent molecules and 11 more water molecules. Moreover, the two large non-polar packing contacts have now been defined in detail. The analysis indicates that the beta-barrel constitutes a solid scaffold such that the long external loops need not contribute to stability. These loops are highly mobile and thus cause a major problem during the crystallization process. The beta-barrel was related to those of lipocalins. Two further crystal forms with exceptionally dense packing arrangements were established at medium resolution.     

Conformer selection and induced fit in flexible backbone protein-protein docking using computational and NMR ensembles

Accommodating backbone flexibility continues to be the most difficult challenge in computational docking of protein-protein complexes. Towards that end, we simulate four distinct biophysical models of protein binding in RosettaDock, a multiscale Monte-Carlo-based algorithm that uses a quasi-kinetic search process to emulate the diffusional encounter of two proteins and to identify low-energy complexes. The four binding models are as follows: (1) key-lock (KL) model, using rigid-backbone docking; (2) conformer selection (CS) model, using a novel ensemble docking algorithm; (3) induced fit (IF) model, using energy-gradient-based backbone minimization; and (4) combined conformer selection/induced fit (CS/IF) model. Backbone flexibility was limited to the smaller partner of the complex, structural ensembles were generated using Rosetta refinement methods, and docking consisted of local perturbations around the complexed conformation using unbound component crystal structures for a set of 21 target complexes. The lowest-energy structure contained > 30% of the native residue-residue contacts for 9, 13, 13, and 14 targets for KL, CS, IF, and CS/IF docking, respectively. When applied to 15 targets using nuclear magnetic resonance ensembles of the smaller protein, the lowest-energy structure recovered at least 30% native residue contacts in 3, 8, 4, and 8 targets for KL, CS, IF, and CS/IF docking, respectively. CS/IF docking of the nuclear magnetic resonance ensemble performed equally well or better than KL docking with the unbound crystal structure in 10 of 15 cases. The marked success of CS and CS/IF docking shows that ensemble docking can be a versatile and effective method for accommodating conformational plasticity in docking and serves as a demonstration for the CS theory—that binding-competent conformers exist in the unbound ensemble and can be selected based on their favorable binding energies.     

Validation of NMR side-chain conformations by packing calculations

The precision and accuracy of protein structures determined by nuclear magnetic resonance (NMR) spectroscopy depend on the completeness of input experimental data set. Typically, rather than a single structure, an ensemble of up to 20 equally representative conformers is generated and routinely deposited in the Protein Database. There are substantially more experimentally derived restraints available to define the main-chain coordinates than those of the side chains. Consequently, the side-chain conformations among the conformers are more variable and less well defined than those of the backbone. Even when a side chain is determined with high precision and is found to adopt very similar orientations among all the conformers in the ensemble, it is possible that its orientation might still be incorrect. Thus, it would be helpful if there were a method to assess independently the side-chain orientations determined by NMR. Recently, homology modeling by side-chain packing algorithms has been shown to be successful in predicting the side-chain conformations of the buried residues for a protein when the main-chain coordinates and sequence information are given. Since the main-chain coordinates determined by NMR are consistently more reliable than those of the side-chains, we have applied the side-chain packing algorithms to predict side-chain conformations that are compatible with the NMR-derived backbone. Using four test cases where the NMR solution structures and the X-ray crystal structure of the same protein are available, we demonstrate that the side-chain packing method can provide independent validation for the side-chain conformations of NMR structures. Comparison of the side-chain conformations derived by side-chain packing prediction and by NMR spectroscopy demonstrates that when there is agreement between the NMR model and the predicted model, on average 78% of the time the X-ray structure also concurs. While the side-chain packing method can confirm the reliable residue conformations in NMR models, more importantly, it can also identify the questionable residue conformations with an accuracy of 60%. This validation method can serve to increase the confidence level for potential users of structural models determined by NMR.