Characterization of the In Vivo Deuteration of Native Phospholipids by Mass Spectrometry Yields Guidelines for Their Regiospecific Customization.

Customization of deuterated biomolecules is vital for many advanced biological experiments including neutron scattering. However, because it is challenging to control the proportion and regiospecificity of deuterium incorporation in live systems, often only two or three synthetic lipids are mixed together to form simplistic model membranes. This limits the applicability and biological accuracy of the results generated with these synthetic membranes. Despite some limited prior examination of deuterating Escherichia coli lipids in vivo, this approach has not been widely implemented. Here, an extensive mass spectrometry-based profiling of E. coli phospholipid deuteration states with several different growth media was performed, and a computational method to describe deuterium distributions with a one-number summary is introduced. The deuteration states of 36 lipid species were quantitatively profiled in 15 different growth conditions, and tandem mass spectrometry was used to reveal deuterium localization. Regressions were employed to enable the prediction of lipid deuteration for untested conditions. Small-angle neutron scattering was performed on select deuterated lipid samples, which validated the deuteration states calculated from the mass spectral data. Based on these experiments, guidelines for the design of specifically deuterated phospholipids are described. This unlocks even greater capabilities from neutron-based techniques, enabling experiments that were formerly impossible.

Ectopic Production of 3,4-Dihydroxybenzoate in Planta Affects Cellulose Structure and Organization.

Lignocellulosic biomass is a highly sustainable and largely carbon dioxide neutral feedstock for the production of biofuels and advanced biomaterials. Although thermochemical pretreatment is typically used to increase the efficiency of cell wall deconstruction, genetic engineering of the major plant cell wall polymers, especially lignin, has shown promise as an alternative approach to reduce biomass recalcitrance. Poplar trees with reduced lignin content and altered composition were previously developed by overexpressing bacterial 3-dehydroshikimate dehydratase (QsuB) enzyme to divert carbon flux from the shikimate pathway. In this work, three transgenic poplar lines with increasing QsuB expression levels and different lignin contents were studied using small-angle neutron scattering (SANS) and wide-angle X-ray scattering (WAXS). SANS showed that although the cellulose microfibril cross-sectional dimension remained unchanged, the ordered organization of the microfibrils progressively decreased with increased QsuB expression. This was correlated with decreasing total lignin content in the QsuB lines. WAXS showed that the crystallite dimensions of cellulose microfibrils transverse to the growth direction were not affected by the QsuB expression, but the crystallite dimensions parallel to the growth direction were decreased by ∼20%. Cellulose crystallinity was also decreased with increased QsuB expression, which could be related to high levels of 3,4-dihydroxybenzoate, the product of QsuB expression, disrupting microfibril crystallization. In addition, the cellulose microfibril orientation angle showed a bimodal distribution at higher QsuB expression levels. Overall, this study provides new structural insights into the impact of ectopic synthesis of small-molecule metabolites on cellulose organization and structure that can be used for future efforts aimed at reducing biomass recalcitrance.

Synthesis of Perdeuterated Alkyl Amines/Amides with Pt/C as Catalyst under Mild Conditions.

A convenient method for the synthesis of perdeuterated alkyl amides/amines is disclosed. Perdeuterated acetyl amides can be achieved by a hydrogen-deuterium (H/D) exchange protocol with Pt/C as a catalyst and D2O as a deuterium source under mild conditions. After removal or reduction of the acetyl group, this protocol can provide perdeuterated primary, secondary, and tertiary amines, which are difficult to achieve via other methods.

Disordered Domain Shifts the Conformational Ensemble of the Folded Regulatory Domain of the Multidomain Oncoprotein c-Src.

c-Src kinase is a multidomain non-receptor tyrosine kinase that aberrantly phosphorylates several signaling proteins in cancers. Although the structural properties of the regulatory domains (SH3-SH2) and the catalytic kinase domain have been extensively characterized, there is less knowledge about the N-terminal disordered region (SH4UD) and its interactions with the other c-Src domains. Here, we used domain-selective isotopic labeling combined with the small-angle neutron scattering contrast matching technique to study SH4UD interactions with SH3-SH2. Our results show that in the presence of SH4UD, the radius of gyration (Rg) of SH3-SH2 increases, indicating that it has a more extended conformation. Hamiltonian replica exchange molecular dynamics simulations provide a detailed molecular description of the structural changes in SH4UD-SH3-SH2 and show that the regulatory loops of SH3 undergo significant conformational changes in the presence of SH4UD, while SH2 remains largely unchanged. Overall, this study highlights how a disordered region can drive a folded region of a multidomain protein to become flexible, which may be important for allosteric interactions with binding partners. This may help in the design of therapeutic interventions that target the regulatory domains of this important family of kinases.

Deconstruction of biomass enabled by local demixing of cosolvents at cellulose and lignin surfaces.

A particularly promising approach to deconstructing and fractionating lignocellulosic biomass to produce green renewable fuels and high-value chemicals pretreats the biomass with organic solvents in aqueous solution. Here, neutron scattering and molecular-dynamics simulations reveal the temperature-dependent morphological changes in poplar wood biomass during tetrahydrofuran (THF):water pretreatment and provide a mechanism by which the solvent components drive efficient biomass breakdown. Whereas lignin dissociates over a wide temperature range (>25 °C) cellulose disruption occurs only above 150 °C. Neutron scattering with contrast variation provides direct evidence for the formation of THF-rich nanoclusters (Rg ∼ 0.5 nm) on the nonpolar cellulose surfaces and on hydrophobic lignin, and equivalent water-rich nanoclusters on polar cellulose surfaces. The disassembly of the amphiphilic biomass is thus enabled through the local demixing of highly functional cosolvents, THF and water, which preferentially solvate specific biomass surfaces so as to match the local solute polarity. A multiscale description of the efficiency of THF:water pretreatment is provided: matching polarity at the atomic scale prevents lignin aggregation and disrupts cellulose, leading to improvements in deconstruction at the macroscopic scale.

Assessment of antimicrobial activity of melittin encapsulated in bicontinuous microemulsions prepared using renewable oils.

The objective of this study is to demonstrate that melittin, a well-studied antimicrobial peptide (AMP), can be solubilized in an active form in bicontinuous microemulsions (BMEs) that employ biocompatible oils. The systems investigated consisted of Winsor-III and -IV BME phases composed of Water/Aerosol-OT (AOT)/Polysorbate 85/isopropyl myristate and a Winsor-IV BME employing Polysorbate 80 and limonene. We found that melittin resided in an α-helix-rich configuration and was in an apolar environment for the AOT/Polysorbate 85 Winsor-III system, suggesting that melittin interacted with the surfactant monolayer and was in an active conformation. An apolar environment was also detected for melittin in the two Winsor-IV systems, but to a lesser extent than the Winsor-III system. Small-angle X-ray scattering analysis indicated that melittin at a concentration of 1.0 g/Laq in the aqueous subphase of the Winsor-IV systems led to the greatest impact on the BME structure (e.g., decrease of quasi-periodic repeat distance and correlation length and induction of interfacial fluidity). The antimicrobial activity of the Polysorbate 80 Winsor-IV system was evaluated against several bacteria prominent in chronic wounds and surgical site infections (SSIs). Melittin-free BMEs inhibited the growth of all tested bacteria due to its oil, limonene, while the inclusion of 1.0 g/Laq of melittin in the BMEs enhanced the activity against several bacteria. A further increase of melittin concentration in the BMEs had no further enhancement. These results demonstrate the potential utility of BMEs as a delivery platform for AMPs and other hydrophilic and lipophilic drugs to inhibit antibiotic-resistant microorganisms in chronic wounds and SSIs.

Evidence for Plant-Conserved Region Mediated Trimeric CESAs in Plant Cellulose Synthase Complexes.

Higher plants synthesize cellulose using membrane-bound, six-lobed cellulose synthase complexes, each lobe containing trimeric cellulose synthases (CESAs). Although molecular biology reports support heteromeric trimers composed of different isoforms, a homomeric trimer was reported for in vitro studies of the catalytic domain of CESA1 of Arabidopsis (AtCESA1CatD) and confirmed in cryoEM structures of full-length CESA8 and CESA7 of poplar and cotton, respectively. In both structures, a small portion of the plant-conserved region (P-CR) forms the only contacts between catalytic domains of the monomers. We report inter-subunit lysine-crosslinks that localize to the small P-CR, negative-stain EM structure, and modeling data for homotrimers of AtCESA1CatD. Molecular dynamics simulations for AtCESA1CatD trimers based on the CESA8 cryoEM structure were stable and dependent upon a small set of residue contacts. The results suggest that homomeric CESA trimers may be important for the synthesis of primary and secondary cell walls and identify key residues for future mutagenic studies.

Generation of the configurational ensemble of an intrinsically disordered protein from unbiased molecular dynamics simulation.

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampled molecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function.

Small Angle Neutron Scattering Shows Nanoscale PMMA Distribution in Transparent Wood Biocomposites.

Transparent wood biocomposites based on PMMA combine high optical transmittance with excellent mechanical properties. One hypothesis is that despite poor miscibility the polymer is distributed at the nanoscale inside the cell wall. Small-angle neutron scattering (SANS) experiments are performed to test this hypothesis, using biocomposites based on deuterated PMMA and "contrast-matched" PMMA. The wood cell wall nanostructure soaked in heavy water is quantified in terms of the correlation distance d between the center of elementary cellulose fibrils. For wood/deuterated PMMA, this distance d is very similar as for wood/heavy water (correlation peaks at q ≈ 0.1 Å-1). The peak disappears when contrast-matched PMMA is used, indeed proving nanoscale polymer distribution in the cell wall. The specific processing method used for transparent wood explains the nanocomposite nature of the wood cell wall and can serve as a nanotechnology for cell wall impregnation of polymers in large wood biocomposite structures.

Transient and stabilized complexes of Nsp7, Nsp8, and Nsp12 in SARS-CoV-2 replication.

The replication transcription complex (RTC) from the virus SARS-CoV-2 is responsible for recognizing and processing RNA for two principal purposes. The RTC copies viral RNA for propagation into new virus and for ribosomal transcription of viral proteins. To accomplish these activities, the RTC mechanism must also conform to a large number of imperatives, including RNA over DNA base recognition, basepairing, distinguishing viral and host RNA, production of mRNA that conforms to host ribosome conventions, interfacing with error checking machinery, and evading host immune responses. In addition, the RTC will discontinuously transcribe specific sections of viral RNA to amplify certain proteins over others. Central to SARS-CoV-2 viability, the RTC is therefore dynamic and sophisticated. We have conducted a systematic structural investigation of three components that make up the RTC: Nsp7, Nsp8, and Nsp12 (also known as RNA-dependent RNA polymerase). We have solved high-resolution crystal structures of the Nsp7/8 complex, providing insight into the interaction between the proteins. We have used small-angle x-ray and neutron solution scattering (SAXS and SANS) on each component individually as pairs and higher-order complexes and with and without RNA. Using size exclusion chromatography and multiangle light scattering-coupled SAXS, we defined which combination of components forms transient or stable complexes. We used contrast-matching to mask specific complex-forming components to test whether components change conformation upon complexation. Altogether, we find that individual Nsp7, Nsp8, and Nsp12 structures vary based on whether other proteins in their complex are present. Combining our crystal structure, atomic coordinates reported elsewhere, SAXS, SANS, and other biophysical techniques, we provide greater insight into the RTC assembly, mechanism, and potential avenues for disruption of the complex and its functions.

Unusual zwitterionic catalytic site of SARS-CoV-2 main protease revealed by neutron crystallography.

The main protease (3CL Mpro) from SARS-CoV-2, the etiological agent of COVID-19, is an essential enzyme for viral replication. 3CL Mpro possesses an unusual catalytic dyad composed of Cys145 and His41 residues. A critical question in the field has been what the protonation states of the ionizable residues in the substrate-binding active-site cavity are; resolving this point would help understand the catalytic details of the enzyme and inform rational drug development against this pernicious virus. Here, we present the room-temperature neutron structure of 3CL Mpro, which allowed direct determination of hydrogen atom positions and, hence, protonation states in the protease. We observe that the catalytic site natively adopts a zwitterionic reactive form in which Cys145 is in the negatively charged thiolate state and His41 is doubly protonated and positively charged, instead of the neutral unreactive state usually envisaged. The neutron structure also identified the protonation states, and thus electrical charges, of all other amino acid residues and revealed intricate hydrogen-bonding networks in the active-site cavity and at the dimer interface. The fine atomic details present in this structure were made possible by the unique scattering properties of the neutron, which is an ideal probe for locating hydrogen positions and experimentally determining protonation states at near-physiological temperature. Our observations provide critical information for structure-assisted and computational drug design, allowing precise tailoring of inhibitors to the enzyme's electrostatic environment.

Biochemical and structural analyses reveal that the tumor suppressor neurofibromin (NF1) forms a high-affinity dimer.

Neurofibromin is a tumor suppressor encoded by the NF1 gene, which is mutated in Rasopathy disease neurofibromatosis type I. Defects in NF1 lead to aberrant signaling through the RAS-mitogen-activated protein kinase pathway due to disruption of the neurofibromin GTPase-activating function on RAS family small GTPases. Very little is known about the function of most of the neurofibromin protein; to date, biochemical and structural data exist only for its GAP domain and a region containing a Sec-PH motif. To better understand the role of this large protein, here we carried out a series of biochemical and biophysical experiments, including size-exclusion chromatography-multiangle light scattering (SEC-MALS), small-angle X-ray and neutron scattering, and analytical ultracentrifugation, indicating that full-length neurofibromin forms a high-affinity dimer. We observed that neurofibromin dimerization also occurs in human cells and likely has biological and clinical implications. Analysis of purified full-length and truncated neurofibromin variants by negative-stain EM revealed the overall architecture of the dimer and predicted the potential interactions that contribute to the dimer interface. We could reconstitute structures resembling high-affinity full-length dimers by mixing N- and C-terminal protein domains in vitro The reconstituted neurofibromin was capable of GTPase activation in vitro, and co-expression of the two domains in human cells effectively recapitulated the activity of full-length neurofibromin. Taken together, these results suggest how neurofibromin dimers might form and be stabilized within the cell.

Cellulose synthase interactive1- and microtubule-dependent cell wall architecture is required for acid growth in Arabidopsis hypocotyls.

Auxin-induced cell elongation relies in part on the acidification of the cell wall, a process known as acid growth that presumably triggers expansin-mediated wall loosening via altered interactions between cellulose microfibrils. Cellulose microfibrils are a major determinant for anisotropic growth and they provide the scaffold for cell wall assembly. Little is known about how acid growth depends on cell wall architecture. To explore the relationship between acid growth-mediated cell elongation and plant cell wall architecture, two mutants (jia1-1 and csi1-3) that are defective in cellulose biosynthesis and cellulose microfibril organization were analyzed. The study revealed that cell elongation is dependent on CSI1-mediated cell wall architecture but not on the overall crystalline cellulose content. We observed a correlation between loss of crossed-polylamellate walls and loss of auxin- and fusicoccin-induced cell growth in csi1-3. Furthermore, induced loss of crossed-polylamellate walls via disruption of cortical microtubules mimics the effect of csi1 in acid growth. We hypothesize that CSI1- and microtubule-dependent crossed-polylamellate walls are required for acid growth in Arabidopsis hypocotyls.

Production of deuterated biomass by cultivation of Lemna minor (duckweed) in D2O.

MAIN CONCLUSION: Common duckweed Lemna minor was cultivated in 50% D2O to produce biomass with 50-60% deuterium incorporation containing cellulose with degree of polymerization close (85%) to that of H2O-grown controls. The small aquatic plant duckweed, particularly the genus Lemna, widely used for toxicity testing, has been proposed as a potential source of biomass for conversion into biofuels as well as a platform for production of pharmaceuticals and specialty chemicals. Ability to produce deuterium-substituted duckweed can potentially extend the range of useful products as well as assist process improvement. Cultivation of these plants under deuterating conditions was previously been reported to require addition of kinetin to induce growth and was hampered by anomalies in cellular morphology and protein metabolism. Here, we report the production of biomass with 50-60% deuterium incorporation by long-term photoheterotrophic growth of common duckweed Lemna minor in 50% D2O with 0.5% glucose. L. minor grown in 50% D2O without addition of kinetin exhibited a lag phase twice that of H2O-grown controls, before start of log phase growth at 40% of control rates. Compared to continuous white fluorescent light, growth rates increased fivefold for H2O and twofold for 50% D2O when plants were illuminated at higher intensity with a metal halide lamp and a diurnal cycle of 12-h light/12-h dark. Deuterium incorporation was determined by a combination of 1H and 2H nuclear magnetic resonance (NMR) to be 40-60%. The cellulose from the deuterated plants had an average-number degree of polymerization (DPn) and polydispersity index (PDI) close to that of H2O-grown controls, while Klason lignin content was reduced. The only major gross morphological change noted was root inhibition.

Hemicellulose-Cellulose Composites Reveal Differences in Cellulose Organization after Dilute Acid Pretreatment.

Model hemicellulose-cellulose composites that mimic plant cell wall polymer interactions were prepared by synthesizing deuterated bacterial cellulose in the presence of glucomannan or xyloglucan. Dilute acid pretreatment (DAP) of these materials was studied using small-angle neutron scattering, X-ray diffraction, and sum frequency generation spectroscopy. The macrofibril dimensions of the pretreated cellulose alone were smaller but with similar entanglement of macrofibrillar network as native cellulose. In addition, the crystallite size dimension along the (010) plane increased. Glucomannan-cellulose underwent similar changes to cellulose, except that the macrofibrillar network was more entangled after DAP. Conversely, in xyloglucan-cellulose the macrofibril dimensions and macrofibrillar network were relatively unchanged after pretreatment, but the cellulose Iß content was increased. Our results point to a tight interaction of xyloglucan with microfibrils while glucomannan only interacts with macrofibril surfaces. This study provides insight into roles of different hemicellulose-cellulose interactions and may help in improving pretreatment processes or engineering plants with decreased recalcitrance.

Bicontinuous microemulsions as a biomembrane mimetic system for melittin.

Antimicrobial peptides effectively kill antibiotic-resistant bacteria by forming pores in prokaryotes' biomembranes via penetration into the biomembranes' interior. Bicontinuous microemulsions, consisting of interdispersed oil and water nanodomains separated by flexible surfactant monolayers, are potentially valuable for hosting membrane-associated peptides and proteins due to their thermodynamic stability, optical transparency, low viscosity, and high interfacial area. Here, we show that bicontinuous microemulsions formed by negatively-charged surfactants are a robust biomembrane mimetic system for the antimicrobial peptide melittin. When encapsulated in bicontinuous microemulsions formed using three-phase (Winsor-III) systems, melittin's helicity increases greatly due to penetration into the surfactant monolayers, mimicking its behavior in biomembranes. But, the threshold melittin concentration required to achieve these trends is lower for the microemulsions. The extent of penetration was decreased when the interfacial fluidity of the microemulsions was increased. These results suggest the utility of bicontinuous microemulsions for isolation, purification, delivery, and host systems for antimicrobial peptides.

Gradients in Wall Mechanics and Polysaccharides along Growing Inflorescence Stems.

At early stages of Arabidopsis (Arabidopsis thaliana) flowering, the inflorescence stem undergoes rapid growth, with elongation occurring predominantly in the apical ∼4 cm of the stem. We measured the spatial gradients for elongation rate, osmotic pressure, cell wall thickness, and wall mechanical compliances and coupled these macroscopic measurements with molecular-level characterization of the polysaccharide composition, mobility, hydration, and intermolecular interactions of the inflorescence cell wall using solid-state nuclear magnetic resonance spectroscopy and small-angle neutron scattering. Force-extension curves revealed a gradient, from high to low, in the plastic and elastic compliances of cell walls along the elongation zone, but plots of growth rate versus wall compliances were strikingly nonlinear. Neutron-scattering curves showed only subtle changes in wall structure, including a slight increase in cellulose microfibril alignment along the growing stem. In contrast, solid-state nuclear magnetic resonance spectra showed substantial decreases in pectin amount, esterification, branching, hydration, and mobility in an apical-to-basal pattern, while the cellulose content increased modestly. These results suggest that pectin structural changes are connected with increases in pectin-cellulose interaction and reductions in wall compliances along the apical-to-basal gradient in growth rate. These pectin structural changes may lessen the ability of the cell wall to undergo stress relaxation and irreversible expansion (e.g. induced by expansins), thus contributing to the growth kinematics of the growing stem.

Direct Determination of Hydroxymethyl Conformations of Plant Cell Wall Cellulose Using 1H Polarization Transfer Solid-State NMR.

In contrast to the well-studied crystalline cellulose of microbial and animal origins, cellulose in plant cell walls is disordered due to its interactions with matrix polysaccharides. Plant cell wall (PCW) is an undisputed source of sustainable global energy; therefore, it is important to determine the molecular structure of PCW cellulose. The most reactive component of cellulose is the exocyclic hydroxymethyl group: when it adopts the tg conformation, it stabilizes intrachain and interchain hydrogen bonding, while gt and gg conformations destabilize the hydrogen-bonding network. So far, information about the hydroxymethyl conformation in cellulose has been exclusively obtained from 13C chemical shifts of monosaccharides and oligosaccharides, which do not reflect the environment of cellulose in plant cell walls. Here, we use solid-state Nuclear Magnetic Resonance (ssNMR) spectroscopy to measure the hydroxymethyl torsion angle of cellulose in two model plants, by detecting distance-dependent polarization transfer between H4 and H6 protons in 2D 13C-13C correlation spectra. We show that the interior crystalline portion of cellulose microfibrils in Brachypodium and Arabidopsis cell walls exhibits H4-H6 polarization transfer curves that are indicative of a tg conformation, whereas surface cellulose chains exhibit slower H4-H6 polarization transfer that is best fit to the gt conformation. Joint constraints by the H4-H6 polarization transfer curves and 13C chemical shifts indicate that it is unlikely for interior cellulose to have a significant population of the gt and gg conformation mixed with the tg conformation, while surface cellulose may adopt a small percentage of the gg conformation. These results provide new constraints to the structure and matrix interactions of cellulose in plant cell walls, and represent the first direct determination of a torsion angle in an important noncrystalline carbohydrate polymer.

Observation of a structural gradient in Winsor-III microemulsion systems.

We demonstrate here for the first time via small-angle neutron scattering (SANS) that the middle, bicontinuous microemulsion (BµE) phase of Winsor-III systems undergoes a gradual change of structure and composition in the vertical direction, contrary to the commonly held belief of uniform structure and composition. A vertical stage was deployed to enable precise alignment of a custom-designed rectangular cell containing the WIII system with respect to the neutron beam, allowing for several different vertical positions to be analyzed. For the water/AOT/CK-2,13 (two-tailed alkyl ethoxylate containing a 1,3-dioxolane linkage)/heptane Winsor-III system, the quasi-periodic repeat distance (d) and correlation length (ξ), obtained from the Teubner-Strey model applied to the SANS data, decreased and the surface area per volume of the surfactant monolayer (via Porod analysis) increased in the downward direction, trends that reflect an increase of surfactant concentration, consistent with the ultralow interfacial tension that often occurs for the lower liquid-liquid interface of many WIII systems. The water/sodium dodecyl sulfate (SDS)/1-pentanol/dodecane system shared the same trend with regard to d as observed for AOT/CK-2,13. In contrast, for SDS/pentanol, ξ increased and the amphiphilicity factor (fa) decreased in the downward direction, trends consistent with a decrease of cosurfactant (pentanol) concentration in the downward direction. Non-uniformity in the vertical direction has implications in the transport of solutes between WIII phases during the extractive purification of proteins or the removal of heavy metals and pollutants from wastewater, or the deposition of BµEs onto hydrophilic vs. hydrophobic surfaces as thin coatings.

The global pattern of trace-element distributions in ocean floor basalts.

The magmatic layers of the oceanic crust are created at constructive plate margins by partial melting of the mantle as it wells up. The chemistry of ocean floor basalts, the most accessible product of this magmatism, is studied for the insights it yields into the compositional heterogeneity of the mantle and its thermal structure. However, before eruption, parental magma compositions are modified at crustal pressures by a process that has usually been assumed to be fractional crystallization. Here we show that the global distributions of trace elements in ocean floor basalts describe a systematic pattern that cannot be explained by simple fractional crystallization alone, but is due to cycling of magma through the global ensemble of magma chambers. Variability in both major and incompatible trace-element contents about the average global pattern is due to fluctuations in the magma fluxes into and out of the chambers, and their depth, as well as to differences in the composition of the parental magmas.