Structural Mechanics of DNA Wrapping in the Nucleosome

Experimental X-ray crystal structures and a database of calculated structural parameters of DNA octamers were used in combination to analyse the mechanics of DNA bending in the nucleosome core complex. The 1kx5 X-ray crystal structure of the nucleosome core complex was used to determine the relationship between local structure at the base-step level and the global superhelical conformation observed for nucleosome-bound DNA. The superhelix is characterised by a large curvature (597°) in one plane and very little curvature (10°) in the orthogonal plane. Analysis of the curvature at the level of 10-step segments shows that there is a uniform curvature of 30° per helical turn throughout most of the structure but that there are two sharper kinks of 50° at ± 2 helical turns from the central dyad base pair. The curvature is due almost entirely to the base-step parameter roll. There are large periodic variations in roll, which are in phase with the helical twist and account for 500° of the total curvature. Although variations in the other base-step parameters perturb the local path of the DNA, they make minimal contributions to the total curvature. This implies that DNA bending in the nucleosome is achieved using the roll-slide-twist degree of freedom previously identified as the major degree of freedom in naked DNA oligomers. The energetics of bending into a nucleosome-bound conformation were therefore analysed using a database of structural parameters that we have previously developed for naked DNA oligomers. The minimum energy roll, the roll flexibility force constant and the maximum and minimum accessible roll values were obtained for each base step in the relevant octanucleotide context to account for the effects of conformational coupling that vary with sequence context. The distribution of base-step roll values and corresponding strain energy required to bend DNA into the nucleosome-bound conformation defined by the 1kx5 structure were obtained by applying a constant bending moment. When a single bending moment was applied to the entire sequence, the local details of the calculated structure did not match the experiment. However, when local 10-step bending moments were applied separately, the calculated structure showed excellent agreement with experiment. This implies that the protein applies variable bending forces along the DNA to maintain the superhelical path required for nucleosome wrapping. In particular, the 50° kinks are constraints imposed by the protein rather than a feature of the 1kx5 DNA sequence. The kinks coincide with a relatively flexible region of the sequence, and this is probably a prerequisite for high-affinity nucleosome binding, but the bending strain energy is significantly higher at these points than for the rest of the sequence. In the most rigid regions of the sequence, a higher strain energy is also required to achieve the standard 30° curvature per helical turn. We conclude that matching of the DNA sequence to the local roll periodicity required to achieve bending, together with the increased flexibility required at the kinks, determines the sequence selectivity of DNA wrapping in the nucleosome.     

Model identification for DNA sequence-structure relationships

We investigate the use of algebraic state-space models for the sequence dependent properties of DNA. By considering the DNA sequence as an input signal, rather than using an all atom physical model, computational efficiency is achieved. A challenge in deriving this type of model is obtaining its structure and estimating its parameters. Here we present two candidate model structures for the sequence dependent structural property Slide and a method of encoding the models so that a recursive least squares algorithm can be applied for parameter estimation. These models are based on the assumption that the value of Slide at a base-step is determined by the surrounding tetranucleotide sequence. The first model takes the four bases individually as inputs and has a median root mean square deviation of 0.90 A. The second model takes the four bases pairwise and has a median root mean square deviation of 0.88 A. These values indicate that the accuracy of these models is within the useful range for structure prediction. Performance is comparable to published predictions of a more physically derived model, at significantly less computational cost.     

Comparison of Crystal Field Dependent and Independent Methods to Analyse Lanthanide Induced NMR Shifts in Axially Symmetric Complexes. Part II: Systems with a C4 Symmetry Axis

Analysis of the LIS data for several series of Ln³⁺ complexes of C₄ symmetry in terms of structural changes, crystal-field effects and/or variation of hyperfine constants along the lanthanide series was undertaken using a combination of the two-nuclei and three-nuclei techniques together with the classical onenucleus technique. Isostructurality of whole series of complexes, with changes of the F_i, and B₀² parameters, was clearly defined for the complexes of L by the combination of the two first methods. Small changes, involving the three F_i, G_i and B₀² parameters, are observed for the series of complexes of L-L4, using the three data plotting methods. Some of the plots according to the two- and three-nuclei methods are accidentally linear, without necessarily implying isostructurality of the complexes, as they involve parameters, which may be insensitive to any small structural changes occurring in these systems. These parameter variations could result from a magnification, by the present graphical analysis, of the breaks expected from the gradual structural changes along the series due to the lanthanide contraction. The α and β parameters of the three-nuclei method are not diagnostic of the type of structures the complexes have in solution, due to their very indirect dependence on the geometric factors.     

A study of the complex-formation of phenylalanyl-tRNA-synthetase from Escherichia coli using the tRNA-Phe method of small-angle x-ray scattering

The method of small-angle X-ray scattering was employed to analyse the equilibrium enzyme-substrate complexes in solution. A new approach of analysis of the experimental data was developed. This type of analysis provides the determination of dissociation constants and structural parameters of enzyme-substrate complexes. The radius of gyration (Rg) and dimensions of half-axis of the equivalent prolonged ellipsoid (a, b) of E. coliphenylalanyl-tRNA synthetase and its complexes with one or two tRNA(Phe) molecules have been determined. The values of these parameters speak in favour of structural rearrangements due to the interaction of the enzyme with tRNA(Phe). The thermodynamic characteristics of phenylalanyl-tRNA synthetase complexes with tRNA(Phe) testify to the negative cooperativity in binding of two tRNA molecules with the enzyme.     

Structural analysis of multiprotein complexes by cross-linking, mass spectrometry, and database searching

Most protein complexes are inaccessible to high resolution structural analysis. We report the results of a combined approach of cross-linking, mass spectrometry, and bioinformatics to two human complexes containing large coiled-coil segments, the NDEL1 homodimer and the NDC80 heterotetramer. An important limitation of the cross-linking approach, so far, was the identification of cross-linked peptides from fragmentation spectra. Our novel approach overcomes the data analysis bottleneck of cross-linking and mass spectrometry. We constructed a purpose-built database to match spectra with cross-linked peptides, define a score that expresses the quality of our identification, and estimate false positive rates. We show that our analysis sheds light on critical structural parameters such as the directionality of the homodimeric coiled coil of NDEL1, the register of the heterodimeric coiled coils of the NDC80 complex, and the organization of a tetramerization region in the NDC80 complex. Our approach is especially useful to address complexes that are difficult in addressing by standard structural methods.     

Exploring overlapping functional units with various structure in protein interaction networks

Revealing functional units in protein-protein interaction (PPI) networks are important for understanding cellular functional organization. Current algorithms for identifying functional units mainly focus on cohesive protein complexes which have more internal interactions than external interactions. Most of these approaches do not handle overlaps among complexes since they usually allow a protein to belong to only one complex. Moreover, recent studies have shown that other non-cohesive structural functional units beyond complexes also exist in PPI networks. Thus previous algorithms that just focus on non-overlapping cohesive complexes are not able to present the biological reality fully. Here, we develop a new regularized sparse random graph model (RSRGM) to explore overlapping and various structural functional units in PPI networks. RSRGM is principally dominated by two model parameters. One is used to define the functional units as groups of proteins that have similar patterns of connections to others, which allows RSRGM to detect non-cohesive structural functional units. The other one is used to represent the degree of proteins belonging to the units, which supports a protein belonging to more than one revealed unit. We also propose a regularizer to control the smoothness between the estimators of these two parameters. Experimental results on four S. cerevisiae PPI networks show that the performance of RSRGM on detecting cohesive complexes and overlapping complexes is superior to that of previous competing algorithms. Moreover, RSRGM has the ability to discover biological significant functional units besides complexes.     

A generalized conformational energy function of DNA derived from molecular dynamics simulations

Proteins recognize DNA sequences by two different mechanisms. The first is direct readout, in which recognition is mediated by direct interactions between the protein and the DNA bases. The second is indirect readout, which is caused by the dependence of conformation and the deformability of the DNA structure on the sequence. Various energy functions have been proposed to evaluate the contribution of indirect readout to the free-energy changes in complex formations. We developed a new generalized energy function to estimate the dependence of the deformability of DNA on the sequence. This function was derived from molecular dynamics simulations previously conducted on B-DNA dodecamers, each of which had one possible tetramer sequence embedded at its center. By taking the logarithm of the probability distribution function (PDF) for the base-step parameters of the central base-pair step of the tetramer, its ability to distinguish the native sequence from random ones was superior to that with the previous method that approximated the energy function in harmonic form. From a comparison of the energy profiles calculated with these two methods, we found that the harmonic approximation caused significant errors in the conformational energies of the tetramers that adopted multiple stable conformations.     

Differential interactions of antitumor antibiotics chromomycin A(3) and mithramycin with d(TATGCATA)(2) in presence of Mg(2+)

The antitumor antibiotics chromomycin A(3) (CHR) and mithramycin (MTR) are known to inhibit macromolecular biosynthesis by reversibly binding to double stranded DNA with a GC base specificity via the minor groove in the presence of a divalent cation such as Mg(2+). Earlier reports from our laboratory showed that the antibiotics form two types of complexes with Mg(2+): complex I with 1:1 stoichiometry and complex II with 2:1 stoichiometry in terms of the antibiotic and Mg(2+). The binding potential of an octanucleotide, d(TATGCATA)(2), which contains one potential site of association with the above complexes of the two antibiotics, was examined using spectroscopic techniques such as absorption, fluorescence, and circular dichroism. We also evaluated thermodynamic parameters for the interaction. In spite of the presence of two structural moieties of the antibiotic in complex II, a major characteristic feature was the association of a single ligand molecule per molecule of octameric duplex in all cases. This indicated that the modes of association for the two types of complexes with the oligomeric DNA were different. The association was dependent on the nature of the antibiotics. Spectroscopic characterization along with analysis of binding and thermodynamic parameters showed that differences in the mode of recognition by complexes I and II of the antibiotics with polymeric DNA existed at the oligomeric level. Analysis of the thermodynamic parameters led us to propose a partial accommodation of the ligand in the groove without the displacement of bound water molecules and supported earlier results on the DNA structural transition from B --> A type geometry as an obligatory requirement for the accommodation of the bulkier complex II of the two drugs. The role of the carbohydrate moieties of the antibiotics in the DNA recognition process was indicated when we compared the DNA binding properties with the same type of Mg(2+) complex for the two antibiotics.     

The intrinsic curvature of DNA in solution

We propose a detailed quantitative scheme for explaining the anomalous electrophoretic mobility in polyacrylamide gels of repeating sequence DNA. We assume that such DNA adopts a superhelical configuration in these circumstances, and migrates less quickly than straight DNA of the same length because it can only pass through larger holes. The retardation is maximal when the length of the DNA reaches one superhelical turn, but is less for shorter pieces. We attribute the curvature of the superhelix to different angles of roll at each kind of dinucleotide step, i.e. an opening up of an angle by an increased separation on the minor-groove side. The main effect is due to a difference of about 3 degrees in roll values between AA/TT and other steps, together with a difference of about 1 degree in the angle of helical twist: we deduce these values explicitly from some of the available data on gel-running. The scheme involves a simple calculation of the superhelical parameters for any given repeating sequence, and it gives a good correlation with all of the available data. We argue that these same base-step angular parameters are also consistent with observations from X-ray diffraction of crystallized oligomers, and particularly with the recent data on CGCA6GCG from Nelson et al. We are concerned here with the intrinsic curvature of unconstrained DNA, as distinct from the curvature of DNA in association with protein molecules; and this paper represents a first attempt at an absolute determination.     

Rings, bracelet or snaps: fashionable alternatives for Smc complexes

The mechanism of higher order chromosome organization has eluded researchers for over 100 years. A breakthrough occurred with the discovery of multi-subunit protein complexes that contain a core of two molecules from the structural maintenance of chromosome (Smc) family. Smc complexes are important structural components of chromosome organization in diverse aspects of DNA metabolism, including sister chromatid cohesion, condensation, global gene repression, DNA repair and homologous recombination. In these different processes, Smc complexes may facilitate chromosome organization by tethering together two parts of the same or different chromatin strands. The mechanism of tethering by Smc complexes remains to be elucidated, but a number of intriguing topological alternatives are suggested by the unusual structural features of Smc complexes, including their large coiled-coil domains and ATPase activities. Distinguishing between these possibilities will require innovative new approaches.     

Simulated dipeptide recognition by vancomycin

The antimicrobial activity of vancomycin and related glycopeptide antibiotics is due to stereospecific recognition of polypeptide components in bacterial cell walls. To better understand how these antibiotics recognize polypeptide determinants, we have developed dynamic models of the complexes formed by the vancomycin aglycon and two different dipeptide ligands, Ac-D-ala-D-ala and Ac-D-ala-gly. Molecular dynamics simulations of the two complexes, initially conditioned with distance constraints derived from two-dimensional nuclear magnetic resonance (NMR) studies, are conformationally stable and propagate in a manner consistent with the NMR-derived constraints after the constraints are removed. Free energy calculations accurately predict the relative binding affinity of these two complexes and help validate the simulation models for detailed structural analysis. Although the two ligands adopt similar conformations when bound to the antibiotic, there are clear differences in the configuration of intermolecular hydrogen bonds, the overall shape of the antibiotic, and other structural features of the two complexes. This analysis illustrates how complex structural and dynamic factors interrelate and contribute to differences in binding affinity. © 1997 John Wiley & Sons, Ltd.     

Application of the small-angle x-ray scattering technique for the study of two-step equilibrium enzyme-substrate interactions

The small-angle x-ray scattering (SAXS) technique is used for the investigation of two-stage equilibrium macromolecular interactions of the enzyme-substrate type in solution. Experimental procedures and methods of analyzing the data obtained from SAXS have been elaborated. The algorithm for the data analysis allows one to determine the stoichiometric, equilibrium, and structural parameters of the enzyme-substrate complexes obtained. The thermodynamic characteristics for the formation of complexes of double-stranded oligonucleotide with Eco dam methyltransferase (MTase) have been determined and demonstrate a high cooperativity of MTase binding when the ternary complex containing the dimeric enzyme is formed. The structural parameters (R_g, R_c, semiaxes) have been determined for free enzyme and polynucleotides and of enzyme-substrate complexes, indicating structural rearrangements of the enzyme in the interaction with substrates. © 1996 John Wiley & Sons, Inc.     

Indirect readout in drug-DNA recognition: role of sequence-dependent DNA conformation

DNA-binding drugs have numerous applications in the engineered gene regulation. However, the drug-DNA recognition mechanism is poorly understood. Drugs can recognize specific DNA sequences not only through direct contacts but also indirectly through sequence-dependent conformation, in a similar manner to the indirect readout mechanism in protein-DNA recognition. We used a knowledge-based technique that takes advantage of known DNA structures to evaluate the conformational energies. We built a dataset of non-redundant free B-DNA crystal structures to calculate the distributions of adjacent base-step and base-pair conformations, and estimated the effective harmonic potentials of mean force (PMF). These PMFs were used to calculate the conformational energy of drug-DNA complexes, and the Z-score as a measure of the binding specificity. Comparing the Z-scores for drug-DNA complexes with those for free DNA structures with the same sequence, we observed that in several cases the Z-scores became more negative upon drug binding. Furthermore, the specificity is position-dependent within the drug-bound region of DNA. These results suggest that DNA conformation plays an important role in the drug-DNA recognition. The presented method provides a tool for the analysis of drug-DNA recognition and can facilitate the development of drugs for targeting a specific DNA sequence.     

Mass spectrometry on hydrogen/deuterium exchange of dihydrofolate reductase: effects of ligand binding

To address the effects of ligand binding on the structural fluctuations of Escherichia coli dihydrofolate reductase (DHFR), the hydrogen/deuterium (H/D) exchange kinetics of its binary and ternary complexes formed with various ligands (folate, dihydrofolate, tetrahydrofolate, NADPH, NADP(+), and methotrexate) were examined using electrospray ionization mass spectrometry. The kinetic parameters of H/D exchange reactions, which consisted of two phases with fast and slow rates, were sensitively influenced by ligand binding, indicating that changes in the structural fluctuation of the DHFR molecule are associated with the alternating binding and release of the cofactor and substrate. No additivity was observed in the kinetic parameters between a ternary complex and its constitutive binary complexes, indicating that ligand binding cooperatively affects the structural fluctuation of the DHFR molecule via long-range interactions. The local H/D exchange profile of pepsin digestion fragments was determined by matrix-assisted laser desorption/ionization mass spectrometry, and the helix and loop regions that appear to participate in substrate binding, largely fluctuating in the apo-form, are dominantly influenced by ligand binding. These results demonstrate that the structural fluctuation of kinetic intermediates plays an important role in enzyme function, and that mass spectrometry on H/D exchange coupled with ligand binding and protease digestion provide new insight into the structure-fluctuation-function relationship of enzymes.     

Effects of amino group modification in discoidal apolipoprotein A-I-egg phosphatidylcholine-cholesterol complexes on their reactions with lecithin:cholesterol acyltransferase

Discoidal complexes of human apolipoprotein A-I-egg phosphatidylcholine-cholesterol were prepared by the sodium cholate dialysis procedure and were reacted to varying extents with the amino group reagents citraconic anhydride, diketene, and formaldehyde in the presence of sodium borohydride. Modification of positive lysine residues with negative or neutral groups (citraconic anhydride and diketene, respectively) resulted, for extensively reacted complexes (90%), in structural alterations and in a marked decrease in reactivity with purified human lecithin:cholesterol acyltransferase. The structural and kinetic effects were partially reversible by removal of the modifying groups or by increased ionic strength. Similar extents of modification (84%) with retention of positive charge and introduction of two methyl groups (reductive methylation) had no effect on the structure or the reactivity of the complexes. These results, together with kinetic data at variable complex concentrations or at variable temperatures, indicate that specific lysine residues of apolipoprotein A-I are not involved in the lecithin:cholesterol acyltransferase activation process; instead, charge interactions and structural changes are responsible for the observed decrease in activating capacity. In terms of kinetic parameters, intrinsic K*m values and probably enzyme-substrate particle dissociation constants are affected, but the activation energies remain the same upon chemical modification.     

Structure of recombinant rabies virus nucleoprotein-RNA complex and identification of the phosphoprotein binding site

Rabies virus nucleoprotein (N) was produced in insect cells, in which it forms nucleoprotein-RNA (N-RNA) complexes that are biochemically and biophysically indistinguishable from rabies virus N-RNA. We selected recombinant N-RNA complexes that were bound to short insect cellular RNAs which formed small rings containing 9 to 11 N monomers. We also produced recombinant N-RNA rings and viral N-RNA that were treated with trypsin and that had lost the C-terminal quarter of the nucleoprotein. Trypsin-treated N-RNA no longer bound to recombinant rabies virus phosphoprotein (the viral polymerase cofactor), so the presence of the C-terminal part of N is needed for binding of the phosphoprotein. Both intact and trypsin-treated recombinant N-RNA rings were analyzed with cryoelectron microscopy, and three-dimensional models were calculated from single-particle image analysis combined with back projection. Nucleoprotein has a bilobed shape, and each monomer has two sites of interaction with each neighbor. Trypsin treatment cuts off part of one of the lobes without shortening the protein or changing other structural parameters. Using negative-stain electron microscopy, we visualized phosphoprotein bound to the tips of the N-RNA rings, most likely at the site that can be removed by trypsin. Based on the shape of N determined here and on structural parameters derived from electron microscopy on free rabies virus N-RNA and from nucleocapsid in virus, we propose a low-resolution model for rabies virus N-RNA in the virus.     

On the interface of light-harvesting antenna complexes and reaction centers in oxygenic photosynthesis

Photosynthetic pigment-protein complexes (PPCs) accomplish light-energy capture and photochemistry in natural photosynthesis. In this review, we examine three pigment protein complexes in oxygenic photosynthesis: light-harvesting antenna complexes and two reaction centers: Photosystem II (PSII), and Photosystem I (PSI). Recent technological developments promise unprecedented insights into how these multi–component protein complexes are assembled into higher order structures and thereby execute their function. Furthermore, the interfacial domain between light-harvesting antenna complexes and PSII, especially the potential roles of the structural loops from CP29 and the PB–loop of ApcE in higher plant and cyanobacteria, respectively, are discussed. It is emphasized that the structural nuances are required for the structural dynamics and consequently for functional regulation in response to an ever–changing and challenging environment.     

Biological insights from hydrogen exchange mass spectrometry

Over the past two decades, hydrogen exchange mass spectrometry (HXMS) has achieved the status of a widespread and routine approach in the structural biology toolbox. The ability of hydrogen exchange to detect a range of protein dynamics coupled with the accessibility of mass spectrometry to mixtures and large complexes at low concentrations result in an unmatched tool for investigating proteins challenging to many other structural techniques. Recent advances in methodology and data analysis are helping HXMS deliver on its potential to uncover the connection between conformation, dynamics and the biological function of proteins and complexes. This review provides a brief overview of the HXMS method and focuses on four recent reports to highlight applications that monitor structure and dynamics of proteins and complexes, track protein folding, and map the thermodynamics and kinetics of protein unfolding at equilibrium. These case studies illustrate typical data, analysis and results for each application and demonstrate a range of biological systems for which the interpretation of HXMS in terms of structure and conformational parameters provides unique insights into function. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.