Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPLEXITY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Vibrational analysis of methyl cation—Rare gas atom complexes: CH + 3 —Rg (Rg = He, Ne, Ar, Kr)

The vibrational spectra of simple CH + 3 —Rg (Rg = He, Ne, Ar, Kr) complexes have been studied by vibrational configuration interaction theory relying on multidimensional potential energy surfaces (PESs) obtained from explicitly correlated coupled cluster calculations, CCSD(T)-F12a. In agreement with experimental results, the series of rare gas atoms leads to rather unsystematic results and indicates huge zero point vibrational energy effects for the helium complex. In order to study these sensitive complexes more consistently, we also introduce configuration averaged vibrational self-consistent field theory, which is a generalization of standard vibrational self-consistent field theory to several configurations. The vibrational spectra of the complexes are compared to that of the methyl cation, for which corrections due to scalar-relativistic effects, high-order coupled-cluster terms, e.g., quadruple excitations, and core-valence correlation have explicitly been accounted for. The occurrence of tunneling splittings for the vibrational ground-state of CH + 3 —He has been investigated on the basis of semiclassical instanton theory. Furthermore these calculations and a direct comparison of the energy profiles along the intrinsic reaction coordinates with that of the hydronium cation, H 3 O + , suggest that tunneling effects for vibrationally excited states should be very small.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Two single-reference approaches to singlet biradicaloid problems: Complex, restricted orbitals and approximate spin-projection combined with regularized orbital-optimized Møller-Plesset perturbation theory

We present a comprehensive study of two single-reference approaches to singlet biradicaloids. These two approaches are based on the recently developed regularized orbital-optimized Møller-Plesset method (κ-OOMP2). The first approach is to combine Yamaguchi's approximate projection (AP) scheme and κ-OOMP2 with unrestricted (U) orbitals (κ-UOOMP2). By capturing only essential symmetry breaking, κ-UOOMP2 can serve as a suitable basis for AP. The second approach is κ-OOMP2 with complex, restricted (cR) orbitals (κ-cROOMP2). Although its applicability is more limited due to the comparative rarity of cR solutions, κ-cROOMP2 offers a simple framework for describing singlet biradicaloids with complex polarization while removing artificial spatial symmetry breaking. We compare the scope of these two methods with numerical studies. We show that AP+κ-UOOMP2 and κ-cROOMP2 can perform similarly well in the TS12 set, a dataset that includes 12 data points for triplet-singlet gaps of several atoms and diatomic molecules with a triplet ground state. This was also found to be true for the barrier height of a reaction involving attack on a cysteine ion by a singlet oxygen molecule. However, we also demonstrate that in highly symmetric systems like C 30 (D 5h ), κ-cROOMP2 is more suitable as it conserves spatial symmetry. Finally, we present an organic biradicaloid that does not have a κ-cROOMP2 solution in which case only AP+κ-UOOMP2 is applicable. Here, we recommend κ-cROOMP2 whenever complex polarization is essential and AP+κ-UOOMP2 for biradicaloids without essential complex polarization but with essential spin-polarization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Exploring the Light-induced Dynamics in Solvated Metallogrid Complexes with Femtosecond Pulses Across the Electromagnetic Spectrum

Oligonuclear complexes of d 4 -d 7 transition metal ion centers that undergo spin-switching have long been developed for their practical role in molecular electronics. Recently, they also have appeared as promising photochemical reactants demonstrating improved stability. However, the lack of knowledge about their photophysical properties in the solution phase compared to their mononuclear congeners is currently hampering their inclusion into advanced light-driven reactions. In the present study, the ultrafast photoinduced dynamics in a solvated [2x2] iron(II) metallogrid complex are characterized by combining transient optical, infrared absorption and X-ray emission spectroscopic measurements on the femtosecond timescale. The analysis is supported by DFT calculations. The photocycle can be described in terms of intra-site transitions, where the Fe II in the low-spin (LS) state are independently photoexcited. The Franck-Condon state decays via the formation of a vibrationally-hot high-spin (HS) state that displays coherent behavior within a few picoseconds and thermalizes within tens of picosecond to yield a metastable HS state living for several hundreds of nanoseconds. Additionally, systematic comparison with the closely related mononuclear complex [Fe(terpy) 2 ] 2+ reveals that nuclearity has a profound impact on the photoinduced dynamics. Finally, this work provides guidelines for expanding the integration of oligonuclear complexes into new photoconversion schemes that may be triggered by ultrafast spin-switching.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Spin-orbit coupling in molecular complexes beyond van der Waals regime: Key factors for further splitting of 2 P 3/2 ground state

Here, we report a joint spectroscopic and theoretical study probing spin-orbit coupling (SOC) in a variety of molecular complexes between an iodine atom and a ligand (L) with L ranging from Ar, HF to formic/acetic acids, and glycine/N-methylated glycine derivatives. Cryogenic photoelectron spectroscopy of L·I - (L=HCOOH, CH 3 COOH) reveals three distinct peaks, identified as three SOC states, denoted as X(1/2), A(3/2), and B(l/2) for the corresponding neutrals. The X and A separation ΔE XA is measured to be 0.10 eV for both, whereas the X and B gap ΔE XB is 0.98 and 0.97 eV for formic and acetic acid, respectively. These new ΔE XA values are compared with the previously reported values for the molecular complexes L·I· with L=Ar, HF, glycine, and A-methylated glycines. All together these complexes encompass a diversity of intermolecular interactions, from van der Waals to weak and strong hydrogen bonding. While the ΔE XB remains similar, the ΔE XA is shown to be extremely sensitive to the type of ligands and interactions, spanning from 5 meV to 150 meV. High-level relativistic quantum calculations including explicit SOC formulism nicely reproduce all experimental SOC splitting. A direct correlation between the magnitude of ΔE XA with the intermolecular interaction strength or bond distance of the neutral complexes—the stronger interaction (shorter bond length), the greater splitting, is established.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Amplified spontaneous emission from europium-based molecular complexes coupled to photonic crystal cavities

Rare-earth ion-based materials bear many remarkable optical properties that render them highly appealing for lighting and quantum-related applications. However, their small oscillator strength and weak emission often pose limitations. Here, in this work, we synthesize and couple Eu(III)-based molecular complexes to nanobeam photonic crystals supporting air modes. A reasonable spatial overlap between the molecular complexes and cavity modes leads to an average spontaneous emission coupling efficiency of 0.19. Our pump power-dependent photoluminescence measurements evidence amplified spontaneous emission from the molecular complexes with an amplification threshold as low as 4.4 W/cm 2 , likely benefiting from the efficient coupling. These findings suggest that integrating rare-earth ion-based molecular complexes with photonic structures could be a viable approach for regulating their emission characteristics for particular applications.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The TWD40-2 protein and the AP2 complex cooperate in the clathrin-mediated endocytosis of cellulose synthase to regulate cellulose biosynthesis

Significance Our study advances the understanding of how plants make cellulose, a critical constituent of the plant cell wall and the most abundant organic polymer in the terrestrial biosphere. Cellulose is synthesized at the plasma membrane by cellulose synthase (CESA) complexes. Here we highlight the importance of regulated CESA endocytosis in cellulose biosynthesis in Arabidopsis . Furthermore, we show that the TWD40-2 protein, a putative member of the TPLATE complex, functions in clathrin-mediated endocytosis (CME) in ways that are both distinct from and cooperative with the functionally conserved CME adaptor protein 2 (AP2) complex. Because the TPLATE complex of plants is not conserved in animals, our characterization of TWD40-2 offers the opportunity to provide insight into the evolution of eukaryotic trafficking machinery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Structure of lipoprotein lipase in complex with GPIHBP1

Lipoprotein lipase (LPL) plays a central role in triglyceride (TG) metabolism. By catalyzing the hydrolysis of TGs present in TG-rich lipoproteins (TRLs), LPL facilitates TG utilization and regulates circulating TG and TRL concentrations. Until very recently, structural information for LPL was limited to homology models, presumably due to the propensity of LPL to unfold and aggregate. By coexpressing LPL with a soluble variant of its accessory protein glycosylphosphatidylinositol-anchored high-density lipoprotein binding protein 1 (GPIHBP1) and with its chaperone protein lipase maturation factor 1 (LMF1), we obtained a stable and homogenous LPL/GPIHBP1 complex that was suitable for structure determination. We report here X-ray crystal structures of human LPL in complex with human GPIHBP1 at 2.5–3.0 Å resolution, including a structure with a novel inhibitor bound to LPL. Binding of the inhibitor resulted in ordering of the LPL lid and lipid-binding regions and thus enabled determination of the first crystal structure of LPL that includes these important regions of the protein. It was assumed for many years that LPL was only active as a homodimer. The structures and additional biochemical data reported here are consistent with a new report that LPL, in complex with GPIHBP1, can be active as a monomeric 1:1 complex. The crystal structures illuminate the structural basis for LPL-mediated TRL lipolysis as well as LPL stabilization and transport by GPIHBP1.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Targeted polyelectrolyte complex micelles treat vascular complications in vivo

Vascular disease is a leading cause of morbidity and mortality in the United States and globally. Pathological vascular remodeling, such as atherosclerosis and stenosis, largely develop at arterial sites of curvature, branching, and bifurcation, where disturbed blood flow activates vascular endothelium. Current pharmacological treatments of vascular complications principally target systemic risk factors. Improvements are needed. We previously devised a targeted polyelectrolyte complex micelle to deliver therapeutic nucleotides to inflamed endothelium in vitro by displaying the peptide VHPKQHR targeting vascular cell adhesion molecule 1 (VCAM-1) on the periphery of the micelle. This paper explores whether this targeted nanomedicine strategy effectively treats vascular complications in vivo. Disturbed flow-induced microRNA-92a (miR-92a) has been linked to endothelial dysfunction. We have engineered a transgenic line (miR-92a EC-TG /Apoe –/– ) establishing that selective miR-92a overexpression in adult vascular endothelium causally promotes atherosclerosis in Apoe –/– mice. We tested the therapeutic effectiveness of the VCAM-1–targeting polyelectrolyte complex micelles to deliver miR-92a inhibitors and treat pathological vascular remodeling in vivo. VCAM-1–targeting micelles preferentially delivered miRNA inhibitors to inflamed endothelial cells in vitro and in vivo. The therapeutic effectiveness of anti–miR-92a therapy in treating atherosclerosis and stenosis in Apoe –/– mice is markedly enhanced by the VCAM-1–targeting polyelectrolyte complex micelles. These results demonstrate a proof of concept to devise polyelectrolyte complex micelle-based targeted nanomedicine approaches treating vascular complications in vivo.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Structure of the Wnt–Frizzled–LRP6 initiation complex reveals the basis for coreceptor discrimination

Wnt morphogens are critical for embryonic development and tissue regeneration. Canonical Wnts form ternary receptor complexes composed of tissue-specific Frizzled (Fzd) receptors together with the shared LRP5/6 coreceptors to initiate β-catenin signaling. The cryo-EM structure of a ternary initiation complex of an affinity-matured XWnt8–Frizzled8–LRP6 complex elucidates the basis of coreceptor discrimination by canonical Wnts by means of their N termini and linker domains that engage the LRP6 E1E2 domain funnels. Chimeric Wnts bearing modular linker “grafts” were able to transfer LRP6 domain specificity between different Wnts and enable non-canonical Wnt5a to signal through the canonical pathway. Synthetic peptides comprising the linker domain serve as Wnt-specific antagonists. The structure of the ternary complex provides a topological blueprint for the orientation and proximity of Frizzled and LRP6 within the Wnt cell surface signalosome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Large-scale conformational rearrangement of the α5-helix of Gα subunits in complex with the guanine nucleotide exchange factor Ric8A

Resistance to inhibitors of cholinesterase 8A (Ric8A) protein is an important G protein–coupled receptor (GPCR)-independent regulator of G protein α-subunits (Gα), acting as a guanine nucleotide exchange factor (GEF) and a chaperone. Insights into the complex between Ric8A and Gα hold the key to understanding the mechanisms underlying noncanonical activation of G-protein signaling as well as the folding of nascent Gα proteins. In this work, we examined the structure of the complex of Ric8A with minimized Gαi (miniGα i ) in solution by small-angle X-ray scattering (SAXS) and exploited the scattering profile in modeling of the Ric8A/miniGα i complex by steered molecular dynamics (SMD) simulations. A small set of models of the complex featured minimal clash scores, excellent agreement with the experimental SAXS data, and a large-scale rearrangement of the signal-transducing α5-helix of Gα away from its β-sheet core. The resulting interface involved the Gα α5-helix bound to the concave surface of Ric8A and the Gα β-sheet that wraps around the C-terminal part of the Ric8A armadillo domain, leading to a severe disruption of the GDP-binding site. Further modeling of the flexible C-terminal tail of Ric8A indicated that it interacts with the effector surface of Gα. This smaller interface may enable the Ric8A-bound Gα to interact with GTP. The two-interface interaction with Gα described here distinguishes Ric8A from GPCRs and non-GPCR regulators of G-protein signaling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

Structural and Thermodynamic Basis for Weak Interactions between Dihydrolipoamide Dehydrogenase and Subunit-binding Domain of the Branched-chain [alpha]-Ketoacid Dehydrogenase Complex

The purified mammalian branched-chain {alpha}-ketoacid dehydrogenase complex (BCKDC), which catalyzes the oxidative decarboxylation of branched-chain {alpha}-keto acids, is essentially devoid of the constituent dihydrolipoamide dehydrogenase component (E3). The absence of E3 is associated with the low affinity of the subunit-binding domain of human BCKDC (hSBDb) for hE3. In this work, sequence alignments of hSBDb with the E3-binding domain (E3BD) of the mammalian pyruvate dehydrogenase complex show that hSBDb has an arginine at position 118, where E3BD features an asparagine. Substitution of Arg-118 with an asparagine increases the binding affinity of the R118N hSBDb variant (designated hSBDb*) for hE3 by nearly 2 orders of magnitude. The enthalpy of the binding reaction changes from endothermic with the wild-type hSBDb to exothermic with the hSBDb* variant. This higher affinity interaction allowed the determination of the crystal structure of the hE3/hSBDb* complex to 2.4-{angstrom} resolution. The structure showed that the presence of Arg-118 poses a unique, possibly steric and/or electrostatic incompatibility that could impede E3 interactions with the wild-type hSBDb. Compared with the E3/E3BD structure, the hE3/hSBDb* structure has a smaller interfacial area. Solution NMR data corroborated the interactions of hE3 with Arg-118 and Asn-118 in wild-type hSBDb and mutant hSBDb*, respectively. The NMR results also showed that the interface between hSBDb and hE3 does not change significantly from hSBDb to hSBDb*. Taken together, our results represent a starting point for explaining the long standing enigma that the E2b core of the BCKDC binds E3 far more weakly relative to other {alpha}-ketoacid dehydrogenase complexes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

Crystal Structure of Mitochondrial Fission Complex Reveals Scaffolding Function for Mitochondrial Division 1 (Mdv1) Coiled Coil

The mitochondrial fission machinery is best understood in the yeast Saccharomyces cerevisiae, where Fis1, Mdv1, and Dnm1 are essential components. Fis1 is a mitochondrial outer membrane protein that recruits the dynamin-related GTPase Dnm1 during the fission process. This recruitment occurs via Mdv1, which binds both Fis1 and Dnm1 and therefore functions as a molecular adaptor linking the two molecules. Mdv1 has a modular structure, consisting of an N-terminal extension that binds Fis1, a central coiled coil for dimerization, and a C-terminal WD40 repeat region that binds Dnm1. We have solved the crystal structure of a dimeric Mdv1-Fis1 complex that contains both the N-terminal extension and coiled-coil regions of Mdv1. Consistent with previous studies, Mdv1 binds Fis1 through a U-shaped helix-loop-helix motif, and dimerization of the Mdv1- Fis1 complex is mediated by the antiparallel coiled coil of Mdv1. However, the complex is surprisingly compact and rigid due to two additional contacts mediated by the surface of the Mdv1 coiled coil. The coiled coil packs against both Fis1 and the second helix of the Mdv1 helix-loop-helix motif. Mutational analyses showed that these contacts are important for mitochondrial fission activity. These results indicate that, in addition to dimerization, the unusually long Mdv1 coiled coil serves a scaffolding function to stabilize the Mdv1-Fis1 complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Integrative Structure–Function Mapping of the Nucleoporin Nup133 Suggests a Conserved Mechanism for Membrane Anchoring of the Nuclear Pore Complex

The nuclear pore complex (NPC) is the sole passageway for the transport of macromolecules across the nuclear envelope. Nup133, a major component in the essential Y-shaped Nup84 complex, is a large scaffold protein of the NPC's outer ring structure. Here, we describe an integrative modeling approach that produces atomic models for multiple states of Saccharomyces cerevisiae (Sc) Nup133, based on the crystal structures of the sequence segments and their homologs, including the related Vanderwaltozyma polyspora (Vp) Nup133 residues 55 to 502 (VpNup133 55–502 ) determined in this study, small angle X-ray scattering profiles for 18 constructs of ScNup133 and one construct of VpNup133, and 23 negative-stain electron microscopy class averages of ScNup1332–1157. Using our integrative approach, we then computed a multi-state structural model of the full-length ScNup133 and validated it with mutational studies and 45 chemical cross-links determined via mass spectrometry. Finally, the model of ScNup133 allowed us to annotate a potential ArfGAP1 lipid packing sensor (ALPS) motif in Sc and VpNup133 and discuss its potential significance in the context of the whole NPC; we suggest that ALPS motifs are scattered throughout the NPC's scaffold in all eukaryotes and play a major role in the assembly and membrane anchoring of the NPC in the nuclear envelope. Our results are consistent with a common evolutionary origin of Nup133 with membrane coating complexes (the protocoatomer hypothesis); the presence of the ALPS motifs in coatomer-like nucleoporins suggests an ancestral mechanism for membrane recognition present in early membrane coating complexes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Can Third-Body Stabilization of Bimolecular Collision Complexes in Cold Molecular Clouds Happen?

For more than half a century, networks of radiative association, dissociative recombination, and bimolecular reactions have been postulated to drive the low-temperature chemistry of cold molecular clouds. Third-body stabilizations of collision complexes have been assumed to be 'irrelevant' due to short lifetime of such complexes. Here, we conduct crossed molecular beam studies of ground state atomic silicon with diacetylene in combination with electronic structure calculations and microcanonical kinetics models operating under cold molecular cloud conditions. Our combined experimental, electronic structure, and microcanonical kinetics modelling investigations provide compelling evidence that three-body collisions of molecular hydrogen with long-lived reaction intermediates accessed through intersystem crossing are prevalent deep inside molecular clouds. This concept might be exportable to reactions involving polycyclic aromatic hydrocarbons thus affording a versatile machinery to complex organics via third-body stabilizations of bimolecular collision complexes deep inside cold molecular clouds.

astrochemistry↗
DOE OSTI2013

TWO NOVEL PARAMETERS TO EVALUATE THE GLOBAL COMPLEXITY OF THE SUN'S MAGNETIC FIELD AND TRACK THE SOLAR CYCLE

Since the unusually prolonged and weak solar minimum between solar cycles 23 and 24 (2008-2010), the sunspot number is smaller and the overall morphology of the Sun's magnetic field is more complicated (i.e., less of a dipole component and more of a tilted current sheet) compared with the same minimum and ascending phases of the previous cycle. Nearly 13 yr after the last solar maximum ({approx}2000), the monthly sunspot number is currently only at half the highest value of the past cycle's maximum, whereas the polar magnetic field of the Sun is reversing (north pole first). These circumstances make it timely to consider alternatives to the sunspot number for tracking the Sun's magnetic cycle and measuring its complexity. In this study, we introduce two novel parameters, the standard deviation (SD) of the latitude of the heliospheric current sheet (HCS) and the integrated slope (SL) of the HCS, to evaluate the complexity of the Sun's magnetic field and track the solar cycle. SD and SL are obtained from the magnetic synoptic maps calculated by a potential field source surface model. We find that SD and SL are sensitive to the complexity of the HCS: (1) they have low values when the HCS is flat at solar minimum, and high values when the HCS is highly tilted at solar maximum; (2) they respond to the topology of the HCS differently, as a higher SD value indicates that a larger part of the HCS extends to higher latitude, while a higher SL value implies that the HCS is wavier; (3) they are good indicators of magnetically anomalous cycles. Based on the comparison between SD and SL with the normalized sunspot number in the most recent four solar cycles, we find that in 2011 the solar magnetic field had attained a similar complexity as compared to the previous maxima. In addition, in the ascending phase of cycle 24, SD and SL in the northern hemisphere were on the average much greater than in the southern hemisphere, indicating a more tilted and wavier HCS in the north than the south, associated with the early reversal of the polar magnetic field in the north relative to the south.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Non-analyticity in holographic complexity near critical points

Abstract The region near a critical point is studied using holographic models of second-order phase transitions. In a previous paper, we argued that the quantum circuit complexity of the vacuum ( C 0 ) is the largest at the critical point. When deforming away from the critical point by a term the complexity C ( τ ) has a piece non-analytic in τ , namely C 0 − C ( τ ) ∼ | τ − τ c | ν ( d − 1 ) + a n a l y t i c . Here, as usual, ν = 1 d − Δ and ξ is the correlation length ξ ∼ | τ − τ c | − ν and there are possible logarithmic corrections to this expression. That was derived using numerical results for the Bose–Hubbard model and general scaling considerations. In this paper, we show that the same is valid in the case of holographic complexity providing evidence that the results are universal, and at the same time providing evidence for holographic computations of complexity.

Physics↗
DOE OSTI2019

Regulation of LC3 lipidation by the autophagy-specific class III phosphatidylinositol-3 kinase complex

Autophagy is a conserved eukaryotic pathway critical for cellular adaptation to changes in nutrition levels and stress. The class III phosphatidylinositol (PI)3-kinase complexes I and II (PI3KC3-C1 and -C2) are essential for autophagosome initiation and maturation, respectively, from highly curved vesicles. We used a cell-free reaction that reproduces a key autophagy initiation step, LC3 lipidation, as a biochemical readout to probe the role of autophagy-related gene (ATG)14, a PI3KC3-C1-specific subunit implicated in targeting the complex to autophagy initiation sites. We reconstituted LC3 lipidation with recombinant PI3KC3-C1, -C2, or various mutant derivatives added to extracts derived from a CRISPR/Cas9-generated ATG14-knockout cell line. Both complexes C1 and C2 require the C-terminal helix of VPS34 for activity on highly curved membranes. However, only complex C1 supports LC3 lipidation through the curvature-targeting amphipathic lipid packing sensor (ALPS) motif of ATG14. Furthermore, the ALPS motif and VPS34 catalytic activity are required for downstream recruitment of WD-repeat domain phosphoinositide-interacting protein (WIPI)2, a protein that binds phosphatidylinositol 3-phosphate and its product phosphatidylinositol 3, 5-bisphosphate, and a WIPI-binding protein, ATG2A, but do not affect membrane association of ATG3 and ATG16L1, enzymes contributing directly to LC3 lipidation. These data reveal the nuanced role of the ATG14 ALPS in membrane curvature sensing, suggesting that the ALPS has additional roles in supporting LC3 lipidation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Regional heritability advanced complex trait analysis for GPU and traditional parallel architectures

Motivation: Quantification of the contribution of genetic variation to phenotypic variation for complex traits becomes increasingly computationally demanding with increasing numbers of single-nucleotide polymorphisms and individuals. To meet the challenges in making feasible large-scale studies, we present the REgional heritability advanced complex trait analysis software. Adapted from advanced complex trait analysis (and, in turn, genome-wide complex trait analysis), it is tailored to exploit the parallelism present in modern traditional and graphics processing unit (GPU)-accelerated machines, from workstations to supercomputers. Results: We adapt the genetic relationship matrix estimation algorithm to remove limitations on memory, allowing the analysis of large datasets. Further, we build on this to develop a version of the code able to efficiently exploit GPU-accelerated systems for both the genetic relationship matrix and REstricted maximum likelihood (REML) parts of the analysis, offering substantial speedup over the traditional central processing unit version. We develop the ability to analyze multiple small regions of the genome across multiple compute nodes in parallel, following the ‘regional heritability’ approach. We demonstrate the new software using 1024 GPUs in parallel on one of the world’s fastest supercomputers. Availability: The code is freely available at http://www.epcc.ed.ac.uk/software-products

59 BASIC BIOLOGICAL SCIENCES↗