Systematic study of nuclear Overhauser effects vis-à-vis local helical parameters, sugar puckers, and glycosidic torsions in B DNA: insensitivity of NOE to local transitions in B DNA oligonucleotides due to internal structural compensations.

A method has been developed to solve structures of DNA oligomers in solution from the experimental NOE data. The method is a combination of two approaches: (1) full matrix NOESY simulations and (2) conformational calculations of DNA double helix based on generalized helical parameters. The process of the refinement of a solution structure does not involve NMR-derived interproton distance constraints; rather it consists of a direct fitting of a structure to the experimental NOE data, a weighted sum of energy, and R factor being under minimization. A helical parameters-based generation of DNA forms makes it possible to organize the search for the optimal structure more effectively, systematically varying starting conformations. The method has been used to calculate a structure for the self-complementary DNA hexamer GGATCC, which is consistent with the available experimental data. The structure belongs to the B family of forms, although the local structural heterogeneity is very strong. Sugar puckers vary from O4'-exo to C3'-exo; helical steps are open with different magnitudes toward the minor groove. Next, we have addressed the question of how uniquely the structure is defined by the existing NMR data. Different structural parameters have been systematically varied, and their effect on individual NOE's and the R factor has been studied. Two energetically conjugated parameters, sugar puckers and glycosidic angles, can be determined very reliably, because of the strong dependences of the intraresidue H6/H8 to H2'/H2'/H3' NOE's. In contrast, the local helical conformation of DNA and the geometry of base pairs proved to be underdetermined by the existing NOE information, because the effect of any helical parameter on interproton distances can be compensated by the concerted changes in other parameters.     

Studies on the structure of HCN oligomers

The structure of the water-insoluble fraction of HCN oligomers (azulmic acid) was studied by IR, NMR, gel permeation chromatography, and chemical methods. The results show that nearly half of the nitrogen atoms contained in the oligomers are of the primary amino type and the other half are involved in -C = N- type bonding. It was found that the oligomers are easily and almost quantitatively acetylated, and the acetylated oligomers show characteristic acetyl amide IR absorption and NMR spectra. Owing to the greatly improved solubility due to acetylation, the molecular weight distribution was determined for the acetylated oligomers by gel permeation chromatography. Composite peaks were obtained ranging from 300 to 900 in molecular weight. Our results are essentially consistent with the structure proposed by Völker, but we point out there may be other possible structures also consistent with our experimental results.     

Characterization of the complex between mannose-binding lectin trimer and mannose-binding lectin-associated serine proteases

Mannose-binding lectin (MBL) is an oligomeric serum lectin involved in innate immunity. Human MBL is complexed with three types of serine proteases (MASP-1, MASP-2 and MASP-3) and two types of their truncated forms (sMAP and MAp44). When an MBL complex binds to carbohydrates of pathogens, the complement system is activated via the lectin pathway. Human MBL is a mixture of different sized oligomers that range mainly from trimers to hexamers. It has been suggested that different MBL oligomers may have distinct MASP compositions. In the present study, an MBL trimer (MBL-I) exclusive of other oligomers was isolated from human serum by chromatography. Immunoblot analysis of MBL-I revealed that it had been co-purified with MASP-1 and sMAP. This suggests that MASP-1 and sMAP are bound to each other in MBL-I. The MBL-I complex was found to activate C2, but to lack the ability to activate C4 due to the absence of MASP-2.     

Heterogeneous triangular structures of human islet amyloid polypeptide (amylin) with internal hydrophobic cavity and external wrapping morphology reveal the polymorphic nature of amyloid fibrils

The misfolding and self-assembly of human islet amyloid polypeptide (hIAPP or amylin) into amyloid fibrils is pathologically linked to type II diabetes. The polymorphic nature of both hIAPP oligomers and fibrils has been implicated for the molecular origin of hIAPP toxicity to islet β-cells, but little is known about the polymorphic structure and dynamics of these hIAPP oligomers/fibrils at the atomic level. Here, we model the polymorphism of full length hIAPP(1-37) oligomers based on experimental data from solid-state NMR, mass per length, and electron microscopy using all-atom molecular dynamics simulation with explicit solvent. As an alternative to steric zipper structures mostly presented in the 2-fold symmetrical fibrils, the most striking structural feature of our proposed hIAPP oligomers is the presence of 3-fold symmetry along the fibril growth axis, in which three β-sheet-layers wind around a hydrophobic core with different periodicities. These 3-fold triangular hIAPP structures dramatically differ in the details of the β-layer assembly and core-forming sequence at the cross section, but all display a high structural stability with favorable layer-to-layer interactions. The 3-fold hIAPP structures can also serve as templates to present triple-stranded helical fibrils via peptide elongation, with different widths from 8.7 to 9.9 nm, twists from 2.8° to 11.8°, and pitches from 14.5 to 61.1 nm, in reasonable agreement with available biophysical data. Because similar 3-fold Aβ oligomers are also observed by both NMR experiments and our previous simulations, the 3-fold structure could be a general conformation to a broad range of amyloid oligomers and fibrils. Most importantly, unlike the conventional stacking sandwich model, the proposed wrapping-cord structures can readily accommodate more than three β-layers via a two dimension conformation search by rotating and translating the β-layers to adopt different favorable packings, which can greatly enrich the polymorphism of amyloid oligomers and fibrils.     

Beta2-microglobulin amyloid fragment organization and morphology and its comparison to Abeta suggests that amyloid aggregation pathways are sequence specific

Beta2-microglobulin (beta2-m) can form dialysis-related amyloid deposits. The structure of a fragment of beta2-m (K3, Ser20-Lys41) in the oligomeric state has recently been solved. We modeled equilibrium structures of K3 oligomers with different organizations (single and double layers) and morphologies (linear-like and annular-like) for the wild-type and mutants using all-atom molecular dynamics (MD) simulations. We focused on the sheet-to-sheet association force, which is the key in the amyloid organization and morphology. For the linear-like morphology, we observed two stable organizations: (i) single-layered parallel-stranded beta-sheets and (ii) double-layered parallel-stranded antiparallel beta-sheets stacked perpendicular to the fibril axis through the hydrophobic N-terminal-N-terminal (NN) interface. No stable annular structures were observed. The structural instability of the annular morphology was mainly attributed to electrostatic repulsion of three negatively charged residues (Asp15, Glu17, and Asp19) projecting from the same beta-strand surface. Linear-like and annular-like double-layered oligomers with the NN interface are energetically more favorable than other oligomers with C-terminal-C-terminal (CC) or C-terminal-N-terminal (CN) interfaces, emphasizing the importance of hydrophobic interactions and side-chain packing in stabilizing these oligomers. Moreover, only linear-like structures, rather than annular structures, with parallel beta-strands and antiparallel beta-sheet arrangements are possible intermediate states for the K3 beta2-m amyloid fibrils in solution. Comparing the beta2-m fragment with Abeta indicates that while both adopt similar beta-strand-turn-beta-strand motifs, the final amyloid structures can be dramatically different in size, structure, and morphology due to differences in side-chain packing arrangements, intermolecular driving forces, sequence composition, and residue positions, suggesting that the mechanism leading to distinct morphologies and the aggregation pathways is sequence specific.     

Numerical simulation of aldolase tetramer stability

A theoretical study of aldolase tetramer stability, conducted by finite difference Poisson-Boltzmann (FDPB) and modified Tanford-Kirkwood (MTK) techniques using the atomic coordinates of human aldolase, is described. A method for calculating the interaction energy between subunits is proposed. An analysis of the contribution of different energy terms to the stability and oligomeric equilibria (monomer ⇔ dimer ⇔ tetramer) of aldolase is made. It is shown that the loss of solvation energy and electrostatic interactions at very high and low pH-s destabilise the oligomers. These energy terms are compensated over a wide pH range by the stabilization energy due to hydrophobic interactions. It is shown that the aldolase tetramer is energetically more preferable than other oligomers in the pH range from 5 to 11. Subunit-subunit interactions within the tetramer suggest one dimeric form to be the most stable of the possible sub-parts. For this reason the tetramer can be thought of as a “dimer of dimers”. A comparison between our theoretical results and available experimental data shows that the dissociation of the aldolase tetramer below pH 3–4 cooperatively leads to acid denaturation. A second dissociation is predicted to occur at high pH (>12) in addition to the well known acidic dissociation. The analysis suggests that a mutation of His20 or Arg257 to a neutral residue could decrease the pH of the acidic dissociation by approximately 1 pH unit. Received: 16 February 1998 / Revised version: 8 April 1998 / Accepted: 19 April 1998     

Theoretical studies of the structure and molecular dynamics of a peptide crystal

Energy minimizations and molecular dynamics simulations have been performed on the cyclic peptide cyclo-(Ala-Pro-D-Phe)2 in both the isolated and crystal states. The results of these calculations have been analyzed, both to investigate our ability to reproduce experimental data (structure and vibrational and NMR spectra) and to investigate the effects of environment on the energy, structure, and dynamics of peptides. Comparison of the minimized and time-averaged crystal systems with the experimental peptide structure shows that the calculations have closely reproduced the experimental structure. Molecular dynamics of the isolated molecule has led to a new conformation, which is approximately equal to 8.5 kcal/mol more stable than the conformation that exists in the crystal, the latter conformation being stabilized by intermolecular (packing) forces. This illustrates the considerable effect that environment can have on the conformation of peptides. The crystal environment has also been shown to significantly reduce the dynamic conformational fluctuations seen for the isolated molecule. The behavior of the peptide during the isolated simulation also supports the experimental NMR observation of a symmetric structure that differs from the asymmetric, instantaneous structures which characterize the molecule during the dynamics. Calculations of vibrational frequencies of the peptide in the crystal and isolated states show the expected shifts in bond-stretching frequencies due to intermolecular interactions. Finally, we have calculated NMR coupling constants from the dynamics simulation of the isolated peptide and have compared these with the experimental values. This has led to a possible reinterpretation of the experimental data.     

New structures help the modeling of toxic amyloidbeta ion channels

The mechanism of amyloid toxicity is poorly understood and there are two schools of thought in this hotly debated field: the first favors membrane destabilization by intermediate-to-large amyloid oligomers, with consequent thinning and non-specific ion leakage; the second favors ion-specific permeable channels lined by small amyloid oligomers. Published results currently support both mechanisms. However, the amyloidbeta (Abeta) peptide has recently been shown to form a U-shaped 'beta-strand-turn-beta-strand' structure. This structure and the available physiological data present a challenge for computational biology--to provide candidate models consistent with the experimental data. Modeling based on small Abeta oligomers containing extramembranous N-termini predicts channels with shapes and dimensions consistent with experimentally derived channel structures. These results support the hypothesis that small Abeta oligomers can form ion channels. Molecular dynamics modeling can provide blueprints of 3D structural conformations for many other amyloids whose membrane association is key to their toxicity.     

Comparison of native and non‐native ubiquitin oligomers reveals analogous structures and reactivities

Ubiquitin (Ub) chains regulate a wide range of biological processes, and Ub chain connectivity is a critical determinant of the many regulatory roles that this post‐translational modification plays in cells. To understand how distinct Ub chains orchestrate different biochemical events, we and other investigators have developed enzymatic and non‐enzymatic methods to synthesize Ub chains of well‐defined length and connectivity. A number of chemical approaches have been used to generate Ub oligomers connected by non‐native linkages; however, few studies have examined the extent to which non‐native linkages recapitulate the structural and functional properties associated with native isopeptide bonds. Here, we compare the structure and function of Ub dimers bearing native and non‐native linkages. Using small‐angle X‐ray scattering (SAXS) analysis, we show that scattering profiles for the two types of dimers are similar. Moreover, using an experimental structural library and atomistic simulations to fit the experimental SAXS profiles, we find that the two types of Ub dimers can be matched to analogous structures. An important application of non‐native Ub oligomers is to probe the activity and selectivity of deubiquitinases. Through steady‐state kinetic analyses, we demonstrate that different families of deubiquitinases hydrolyze native and non‐native isopeptide linkages with comparable efficiency and selectivity. Considering the significant challenges associated with building topologically diverse native Ub chains, our results illustrate that chains harboring non‐native linkages can serve as surrogate substrates for explorations of Ub function.     

Chemokine oligomerization and interactions with receptors and glycosaminoglycans: the role of structural dynamics in function

The first chemokine structure, that of IL-8/CXCL8, was determined in 1990. Since then, many chemokine structures have emerged. To the initial disappointment of structural biologists, the tertiary structures of these small proteins were found to be highly conserved. However, they have since proven to be much more interesting and diverse than originally expected. Somewhat like lego blocks, many chemokines oligomerize and there is significant diversity in their oligomeric forms and propensity to oligomerize. Chemokines not only interact with receptors where different oligomeric forms can induce different signaling responses, they also interact with glycosaminoglycans which can stabilize oligomers and other structures that would not otherwise form in solution. Although chemokine monomers and dimers yielded quickly to structure determination, structural information about larger chemokine oligomers, chemokines receptors, and complexes of chemokines with glycosaminoglycans and receptors has been more difficult to obtain, but recent breakthroughs suggest that this information will be forthcoming, especially with receptor structures. Equally important and challenging, will be efforts to correlate the structural information with function.     

Structural properties of trimers and tetramers of ribonuclease A

Ribonuclease A aggregates (dimers, trimers, tetramers, pentamers) can be obtained by lyophilization from 40% acetic acid solutions. Each aggregate forms two conformational isomers distinguishable by different basic net charge. The crystal structure of the two dimers has recently been determined; the structure of the higher oligomers is unknown. The results of the study of the two trimeric and tetrameric conformers can be summarized as follows: (1) RNase A trimers and tetramers form by a 3D domain-swapping mechanism. N-terminal and C-terminal types of domain swapping could coexist; (2) the secondary structures of the trimeric and tetrameric conformers do not show significant differences if compared with the secondary structure of monomeric RNase A or its two dimers; (3) a different exposure of tyrosine residues indicates that in the aggregates they have different microenvironments; (4) the two trimeric and tetrameric conformers show different susceptibility to digestion by subtilisin; (5) dimers, trimers, and tetramers of RNase A show unwinding activity on double-helical poly(dA-dT) x poly(dA-dT), that increases as a function of the size of the oligomers; (6) the less basic conformers are more stable than the more basic ones, and a low concentration in solution of trimers and tetramers favors their stability, which is definitely increased by the interaction of the aggregates with poly(dA-dT) x poly(dA-dT); (7) the products of thermal dissociation of the two trimers indicate that their structures could be remarkably different. The dissociation products of the two tetramers allow the proposal of two models for their putative structures.     

Conformational Stability of Fibrillar Amyloid-Beta Oligomers via Protofilament Pair Formation – A Systematic Computational Study

Amyloid- (A) oligomers play a crucial role in Alzheimer’s disease due to their neurotoxic aggregation properties. Fibrillar A oligomerization can lead to protofilaments and protofilament pairs via oligomer elongation and oligomer association, respectively. Small fibrillar oligomers adopt the protofilament topology, whereas fibrils contain at least protofilament pairs. To date, the underlying growth mechanism from oligomers to the mature fibril still remains to be elucidated. Here, we performed all-atom molecular dynamics simulations in explicit solvent on single layer-like protofilaments and fibril-like protofilament pairs of different size ranging from the tetramer to the 48-mer. We found that the initial U-shaped topology per monomer is maintained over time in all oligomers. The observed deviations of protofilaments from the starting structure increase significantly with size due to the twisting of the in-register parallel -sheets. This twist causes long protofilaments to be unstable and leads to a breakage. Protofilament pairs, which are stabilized by a hydrophobic interface, exhibit more fibril-like properties such as the overall structure and the twist angle. Thus, they can act as stable conformational templates for further fibril growth. Key properties like the twist angle, shape complementarity, and energetics show a size-dependent behavior so that small oligomers favor the protofilament topology, whereas large oligomers favor the protofilament pair topology. The region for this conformational transition is at the size of approximately twelve A monomers. From that, we propose the following growth mechanism from A oligomers to fibrils: (1) elongation of short protofilaments; (2) breakage of large protofilaments; (3) formation of short protofilament pairs; and (4) elongation of protofilament pairs.     

Study of the interaction of polynucleotide chains with oligomers by means of chromatography. II. Effect of magnesium ions and noncomplementary bases on the stability of complexes

The interaction of the oligonucleotides ApA, ApApA, ApApC and ApApU with poly(U) and (Ip)₅I and (Ip)₆ with poly(C) has been studied by means of equilibrium gelfiltration through Sephadex.From sorption isotherms the free energies, energies and entropies of complexing have been computed for different concentrations of magnesium ions in the medium.The stoichiometric ratio of polymers to oligomers has been measured and found equal to 2 in the case of ApApA and ApApC. This shows that the cytidylic acid residue is included in the ternary complex. But in the case of ApApU the noncomplementary base is partly squeezed out of the complex.The stacking free energy of neighbouring oligomers has been found to be in the range 1000–3000 $\text{cal}\text{mole}$ depending on the conditions.The stoichimetric ratio has been found to be 1 in the case of poly(C): oligo(I), the stacking energy is equal to 1.2 $\text{kcal}\text{mole}$. The effect of magnesium is somewhat different in the case of double and triple helices and probably reflects the formation of coordination compounds with the nitrogen bases of nucleotides.     

Modeling Amyloid-Beta as Homogeneous Dodecamers and in Complex with Cellular Prion Protein

Soluble amyloid beta (Aβ) peptide has been linked to the pathology of Alzheimer’s disease. A variety of soluble oligomers have been observed to be toxic, ranging from dimers to protofibrils. No tertiary structure has been identified as a single biologically relevant form, though many models are comprised of highly ordered β-sheets. Evidence exists for much less ordered toxic oligomers. The mechanism of toxicity remains highly debated and probably involves multiple pathways. Interaction of Aβ oligomers with the N-terminus of the cellular form of the prion protein (PrP^c) has recently been proposed. The intrinsically disordered nature of this protein and the highly polymorphic nature of Aβ oligomers make structural resolution of the complex exceptionally challenging. In this study, molecular dynamics simulations are performed for dodecameric assemblies of Aβ comprised of monomers having a single, short antiparallel β-hairpin at the C-terminus. The resulting models, devoid of any intermolecular hydrogen bonds, are shown to correlate well with experimental data and are found to be quite stable within the hydrophobic core, whereas the α-helical N-termini transform to a random coil state. This indicates that highly ordered assemblies are not required for stability and less ordered oligomers are a viable component in the population of soluble oligomers. In addition, a tentative model is proposed for the association of Aβ dimers with a double deletion mutant of the intrinsically disordered N-terminus of PrP^c. This may be useful as a conceptual working model for the binding of higher order oligomers and in the design of further experiments.     

Prediction of 1H NMR chemical shifts of DNA oligomers

A set of parameters, devised for the prediction of 1H NMR chemical shifts of heterobase and anomeric protons in the high temperature (greater than 70 degrees C) spectra of RNA oligomers has been found to be applicable to the corresponding DNA oligomers. Fifteen examples of DNA oligomers that have had high temperature spectra recorded and assigned show a mean absolute difference between predicted and assigned shifts of 0.045 ppm. The parameters for uridine H-5 are applied to the calculation of thymidine methyl group shifts and give excellent agreement with experimental assigned shifts. The RNA parameter set is a practical means of assigning heterobase and anomeric protons in DNA oligomers. A programme using the RNA parameter set has been written which enables the sequence of short DNA oligomers to be predicted from their 1H NMR spectra.     

Structural polymorphism of human islet amyloid polypeptide (hIAPP) oligomers highlights the importance of interfacial residue interactions

A 37-residue of human islet amyloid polypeptide (hIAPP or amylin) is a main component of amyloid plaques found in the pancreas of ～90% of type II diabetes patients. It is reported that hIAPP oligomers, rather than mature fibrils, are major toxic species responsible for pancreatic islet β-cell dysfunction and even cell death, but molecular structures of these oligomers remain elusive. In this work, on the basis of recent solid-state NMR and mass-per-length (MPL) data, we model a series of hIAPP oligomers with different β-layers (one, two, and three layers), symmetries (symmetry and asymmetry), and associated interfaces using molecular dynamics simulations. Three distinct interfaces formed by C-terminal β-sheet and C-terminal β-sheet (CC), N-terminal β-sheet and N-terminal β-sheet (NN), and C-terminal β-sheet and N-terminal β-sheet (CN) are identified to drive multiple cross-β-layers laterally associated together to form different amyloid organizations via different intermolecular interactions, in which the CC interface is dominated by polar interactions, the NN interface is dominated by hydrophobic interactions, and the CN interface is dominated by mixed polar and hydrophobic interactions. Overall, the structural stability of the proposed hIAPP oligomers is a result of delicate balance between maximization of favorable peptide-peptide interactions at the interfaces and optimization of solvation energy with globular structure. Different hIAPP oligomeric models indicate a general and intrinsic nature of amyloid polymorphism, driven by different interfacial side-chain interactions. The proposed models are compatible with recent experimental data in overall size, cross-section area, and molecular weight. A general hIAPP aggregation mechanism is proposed on the basis of our simulated models and experimental data.     

Prediction of atomic structure from sequence for double helical DNA oligomers

DNA can adopt different conformations depending on the base sequence, solvent, electrolyte composition and concentration, pH, temperature, and interaction with proteins. Here we present a model for calculating the three-dimensional atomic structure of double-stranded DNA oligomers. A theoretical energy function is used for calculating the interactions within the base steps and an empirical backbone function is used to restrict the conformational space accessible to the bases and to account for the conformational coupling of neighboring steps in a sequence. Conformational searching on large structures or a large number of structures is possible, because each base step can be described by just two primary degrees of freedom (slide and shift). A genetic algorithm is used to search for low-energy structures in slide-shift space, and this allows very rapid optimization of DNA oligomers. The other base step parameters have been previously optimized for all possible slide-shift sequence combinations, and a heuristic algorithm is used to add the atomic details of the backbone conformation in the final step of the calculation. The structures obtained by this method are very similar to the corresponding X-ray crystal structures observed experimentally. The average RMSD is 2.24 Angstroms for a set of 20 oligomer structures. For 15 of these sequences, the X-ray crystal structure is the global energy minimum. The other 5 are bistable sequences that have B-form global energy minima but crystallize as A-DNA.     

Vaccination with soluble Aβ oligomers generates toxicity-neutralizing antibodies

In recent studies of transgenic models of Alzheimer's disease (AD), it has been reported that antibodies to aged beta amyloid peptide 1-42 (Abeta(1-42)) solutions (mixtures of Abeta monomers, oligomers and amyloid fibrils) cause conspicuous reduction of amyloid plaques and neurological improvement. In some cases, however, neurological improvement has been independent of obvious plaque reduction, and it has been suggested that immunization might neutralize soluble, non-fibrillar forms of Abeta. It is now known that Abeta toxicity resides not only in fibrils, but also in soluble protofibrils and oligomers. The current study has investigated the immune response to low doses of Abeta(1-42) oligomers and the characteristics of the antibodies they induce. Rabbits that were injected with Abeta(1-42) solutions containing only monomers and oligomers produced antibodies that preferentially bound to assembled forms of Abeta in immunoblots and in physiological solutions. The antibodies have proven useful for assays that can detect inhibitors of oligomer formation, for immunofluorescence localization of cell-attached oligomers to receptor-like puncta, and for immunoblots that show the presence of SDS-stable oligomers in Alzheimer's brain tissue. The antibodies, moreover, were found to neutralize the toxicity of soluble oligomers in cell culture. Results support the hypothesis that immunizations of transgenic mice derive therapeutic benefit from the immuno-neutralization of soluble Abeta-derived toxins. Analogous immuno-neutralization of oligomers in humans may be a key in AD vaccines.