Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Interactions”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quantum fragmentation in the extended quantum breakdown model

We introduce a one-dimensional (1D) extended quantum breakdown model comprising a fermionic and a spin degree of freedom per site, and featuring a spatially asymmetric breakdown-type interaction between the fermions and spins. Furthermore, our model resembles the breakdown process of particles incident into a cloud chamber with nonzero quantum amplitudes of both exciting and not exciting the local vapor atoms. We analytically show that, in the absence of any magnetic field for the spins, the model exhibits Hilbert space fragmentation within each symmetry sector into exponentially many Krylov subspaces and hence displays nonthermal dynamics. Here, we demonstrate that the fragmentation naturally occurs in an entangled basis and thus provides an example of “quantum fragmentation.” Besides establishing the nature of fragmentation analytically, we also study the long-time behavior of the entanglement entropy and its deviation from the expected Page value as a probe of ergodicity in the system. Upon introducing a magnetic field for the spins, most of the Krylov subspaces merge and the model becomes chaotic. Finally, we study the effects of strong randomness on the system and observe behavior similar to that of many-body localized systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular Modes Elucidate the Nuclear Magnetic Resonance Relaxation of Viscous Fluids

The Bloembergen, Purcell, and Pound (BPP) theory of nuclear magnetic resonance (NMR) relaxation in fluids dating back to 1948 continues to be the linchpin in interpreting NMR relaxation data in applications ranging from characterizing fluids in porous media to medical imaging (MRI). The BPP theory is founded on assuming molecules are hard spheres with 1 H– 1 H dipole pairs reorienting randomly; assumptions that are severe in light of modern understanding of liquids. Nevertheless, it is intriguing to this day that the BPP theory was consistent with the original experimental data for glycerol, a hydrogen-bonding molecular fluid for which the hard-sphere-rigid-dipole assumption is inapplicable. To better understand this incongruity, atomistic molecular simulations are used to compute 1 H NMR T 1 relaxation dispersion (i.e., frequency dependence) in two contrasting cases: glycerol, and a (non hydrogen-bonding) viscosity standard. At high viscosities, simulations predict distinct functional forms of T 1 for glycerol compared to the viscosity standard, in agreement with modern measurements, yet both in contrast to BPP theory. The cause of these departures from BPP theory is elucidated, without assuming any relaxation models and without any free parameters, by decomposing the simulated T 1 response into dynamic molecular modes for both intramolecular and intermolecular interactions. The decomposition into dynamic molecular modes provides an alternative framework to understand the physics of NMR relaxation for viscous fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2015

Interphase Chromosome Conformation and Chromatin-Chromatin Interactions in Human Epithelial Cells Cultured Under Different Gravity Conditions

Microgravity has been shown to alter global gene expression patterns and protein levels both in cultured cells and animal models. It has been suggested that the packaging of chromatin fibers in the interphase nucleus is closely related to genome function, and the changes in transcriptional activity are tightly correlated with changes in chromatin folding. This study explores the changes of chromatin conformation and chromatin-chromatin interactions in the simulated microgravity environment, and investigates their correlation to the expression of genes located at different regions of the chromosome. To investigate the folding of chromatin in interphase under various culture conditions, human epithelial cells, fibroblasts, and lymphocytes were fixed in the G1 phase. Interphase chromosomes were hybridized with a multicolor banding in situ hybridization (mBAND) probe for chromosome 3 which distinguishes six regions of the chromosome as separate colors. After images were captured with a laser scanning confocal microscope, the 3-dimensional structure of interphase chromosome 3 was reconstructed at multi-mega base pair scale. In order to determine the effects of microgravity on chromosome conformation and orientation, measures such as distance between homologous pairs, relative orientation of chromosome arms about a shared midpoint, and orientation of arms within individual chromosomes were all considered as potentially impacted by simulated microgravity conditions. The studies revealed non-random folding of chromatin in interphase, and suggested an association of interphase chromatin folding with radiation-induced chromosome aberration hotspots. Interestingly, the distributions of genes with expression changes over chromosome 3 in cells cultured under microgravity environment are apparently clustered on specific loci and chromosomes. This data provides important insights into how mammalian cells respond to microgravity at molecular level.

International Space Station↗
NASA NTRS2005

An Algorithm for the Transport of Anisotropic Neutrons

One major obstacle to human space exploration is the possible limitations imposed by the adverse effect of long-term exposure to the space environment. Even before human spaceflight began, the potentially brief exposure of astronauts to the very intense random solar particle events (SPE) were of great concern. A new challenge appears in deep space exploration from exposure to the low-intensity heavy-ion flux of the galactic cosmic rays (GCR) since the missions are of long duration and the accumulated GCR exposures can be high. Because cancer induction rates increase behind low to rather large thicknesses of aluminum shielding, according to available biological data on mammalian exposures to GCR like ions, the shield requirements for a Mars mission are prohibitively expensive in terms of mission launch costs. Therefore, a critical issue in the Human Exploration and Development of Space enterprise is cost effective mitigation of risk associated with ionizing radiation exposure. In order to estimate astronaut risk to GCR exposure and associated cancer risks and health hazards, it is necessary to do shield material studies. To determine an optimum radiation shield material it is necessary to understand nuclear interaction processes such as fragmentation and secondary particle production which is a function of energy dependent cross sections. This requires knowledge of material transmission characteristics either through laboratory testing or improved theoretical modeling. Here ion beam transport theory is of importance in that testing of materials in the laboratory environment generated by particle accelerators is a necessary step in materials development and evaluation for space use. The approximations used in solving the Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are a major emphasis of the present work.

Tweed, J.↗
NASA NTRS2001

Constitutive Modeling of Nanotube-Reinforced Polymer Composite Systems

In this study, a technique has been proposed for developing constitutive models for polymer composite systems reinforced with single-walled carbon nanotubes (SWNT). Since the polymer molecules are on the same size scale as the nanotubes, the interaction at the polymer/nanotube interface is highly dependent on the local molecular structure and bonding. At these small length scales, the lattice structures of the nanotube and polymer chains cannot be considered continuous, and the bulk mechanical properties of the SWNT/polymer composites can no longer be determined through traditional micromechanical approaches that are formulated using continuum mechanics. It is proposed herein that the nanotube, the local polymer near the nanotube, and the nanotube/polymer interface can be modeled as an effective continuum fiber using an equivalent-continuum modeling method. The effective fiber retains the local molecular structure and bonding information and serves as a means for incorporating micromechanical analyses for the prediction of bulk mechanical properties of SWNT/polymer composites with various nanotube sizes and orientations. As an example, the proposed approach is used for the constitutive modeling of two SWNT/polyethylene composite systems, one with continuous and aligned SWNT and the other with discontinuous and randomly aligned nanotubes.

Odegard, Gregory M.↗
NASA NTRS2004

Constitutive Modeling of Nanotube-Reinforced Polymer Composite Systems

In this study, a technique has been proposed for developing constitutive models for polymer composite systems reinforced with single-walled carbon nanotubes (SWNT). Since the polymer molecules are on the same size scale as the nanotubes, the interaction at the polymer/nanotube interface is highly dependent on the local molecular structure and bonding. At these small length scales, the lattice structures of the nanotube and polymer chains cannot be considered continuous, and the bulk mechanical properties of the SWNT/polymer composites can no longer be determined through traditional micromechanical approaches that are formulated using continuum mechanics. It is proposed herein that the nanotube, the local polymer near the nanotube, and the nanotube/polymer interface can be modeled as an effective continuum fiber using an equivalent-continuum modeling method. The effective fiber retains the local molecular structure and bonding information and serves as a means for incorporating micromechanical analyses for the prediction of bulk mechanical properties of SWNT/polymer composites with various nanotube sizes and orientations. As an example, the proposed approach is used for the constitutive modeling of two SWNT/polyethylene composite systems, one with continuous and aligned SWNT and the other with discontinuous and randomly aligned nanotubes.

Odegard, Gregory M.↗
NASA NTRS2009

Measurements of the Flowfield Interaction Between Tandem Cylinders

This paper presents the most recent measurements from an ongoing investigation of the unsteady wake interference between a pair of circular cylinders in tandem. The purpose of this investigation is to help build an in-depth experimental database for this canonical flow configuration that embodies the effects of component interaction in landing gear noise. This new set of measurements augments the previous database at the primary Reynolds number (based on tunnel speed and cylinder diameter) of 1.66 105 in four important respects. First, better circumferential resolution of surface pressure fluctuations is obtained via cylinder "clocking". Second, higher resolution particle image velocimetry measurements of the shear layer separating from the cylinders are achieved. Third, the effects of simultaneous boundary layer trips along both the front and rear cylinders, versus front cylinder alone in the previous measurements, are studied. Lastly, on-surface and off-surface characteristics of unsteady flow near the "critical" cylinder spacing, wherein the flow switches intermittently between two states that are characteristic of lower and higher spacings, are examined. This critical spacing occurs in the middle of a relatively sudden change in the drag of either cylinder and is characterized by a loud intermittent noise and a flow behavior that randomly transitions between shear layer attachment to the rear cylinder and constant shedding and rollup in front of it. Analysis of this bistable flow state reveals much larger spanwise correlation lengths of surface pressure fluctuations than those at larger and smaller values of the cylinder spacing.

Neuhart, Dan H.↗
NASA NTRS2012

Processes that Drove the Transition from Chemistry to Biology: Concepts and Evidence

Two properties are particularly germane to the transition from chemistry to biology. One is the emergence of complex molecules (polymers) capable of performing non-trivial functions, such as catalysis, energy transduction or transport across cell walls. The other is the ability of several functions to work in concert to provide reproductive advantage to systems hosting these functions. Biological systems exhibit these properties at remarkable levels of efficiency and accuracy in a way that appears effortless. However, dissection of these properties reveals great complexities that are involved. This opens a question: how a simple, ancestral system could have acquired the required properties? Other questions follow. What are the chances that a functional polymer emerges at random? What is the minimum structural complexity of a polymer to carry out a function at a reasonable level of efficiency? Can we identify concrete, protobiologically plausible mechanisms that yield advantageous coupling between different functions? These and similar questions are at the core of the main topic of this session: how soulless chemistry became life? Clearly, we do not have complete answers to any of these questions. However, in recent years a number of new and sometimes unexpected clues have been brought to light. Of particular interest are proteins because they are the main functional polymers in contemporary cells. The emergence of protein functions is a puzzle. It is widely accepted that a well ]defined, compact structure (fold) is a prerequisite for function. It is equally widely accepted that compact folds are rare among random amino acid polymers. Then, how did protein functionality start? According to one hypothesis well folded were preceded by their poorly folded, yet still functional ancestors. Only recently, however, experimental evidence supporting this hypothesis has been presented. In particular, a small enzyme capable of ligating two RNA fragments with the rate of 106 above background was evolved in vitro. This enzyme does not look like any contemporary protein. It is very flexible and its structure is kept together just by a single salt bridge between a charged residue and a coordinating zinc. A similar picture emerges from studies of simple transmembrane channels that mimic those in ancestral cells. Again, they are extremely flexible and do not form a conventional pore. Yet, they efficiently mediate ion transport. Studies on simple proteins that are on-going in several laboratories hold promise of revealing crucial links between chemical and biological catalysis and other ubiquitous cell functions. Interaction between composition, growth and division of protobiologically relevant vesicles and metabolic reactions that they encapsulate is an example of coupling between simple functions that promotes reproduction and evolution. Recent studies have demonstrated possible mechanisms by which vesicles might have evolved their composition from fatty acids to phospholipids, thus facilitating a number of new cellular functions. Conversely, it has been also demonstrated that an encapsulated metabolism might drive vesicle division. These are, again, examples of processes that might have driven the transition from chemistry to biology.

Pohorille, Andrew↗
DOE OSTI2019

Proton Beam Emittance Growth due to Centroid Fluctuations in Strongly Disrupted Electron Bunches in the Linac-Ring eRHIC

In the energy recovery linac (ERL) based eRHIC design, the electron beam, accelerated in a superconducting ERL, collides with the proton beam circulating in the RHIC ring ones, undergoing a very strong beam-beam interaction with the protons, which warrants careful examination. Following the analysis of the transverse disruption and linear mismatch effects in the electron beam caused by beam-beam collision discussed in [2], we adopt the beam parameters listed in Table I for the eRHIC linac-ring design. After a review of the electron disruption process, we perform long term strong-strong beam-beam simulations with a modification of K. Ohmi’s strong-strong beam-beam code BBSS, with the modification consisting of the regeneration of the electron beam after every collision with the proton beam. With the parameters listed in Table I, a kink instability is observed in the proton beam. By reducing the proton bunch length, thus decreasing the electron disruption parameter, the kink instability is observed to disappear for a proton bunch length σ z,p < 10 cm, or electron disruption parameter d e < 1.5. The effect of random offsets in the electron bunch centroids on the long term stability of the proton beam is studied next, with a proton bunch length σ z,p = 6.3 cm (d e = 0.9), thus below the kink instability threshold. The random offsets of the electron bunch centroids are assumed to be Gaussian distributed, with standard deviation offs, and "delta-correlated" in time (white noise). Studies to take into account finite time-correlation effects and random offsets generated by non-Gaussian random processes will be considered in the future. The effect from fluctuations in the electron bunch population has been studied analytically and will be discussed numerically in a companion paper.

36 MATERIALS SCIENCE↗
NASA NTRS2001

The Use of CASES-97 Observations to Assess and Parameterize the Impact of Land-Surface Heterogeneity on Area-Average Surface Heat Fluxes for Large-Scale Coupled Atmosphere-Hydrology Models

To understand the effects of land-surface heterogeneity and the interactions between the land-surface and the planetary boundary layer at different scales, we develop a multiscale data set. This data set, based on the Cooperative Atmosphere-Surface Exchange Study (CASES97) observations, includes atmospheric, surface, and sub-surface observations obtained from a dense observation network covering a large region on the order of 100 km. We use this data set to drive three land-surface models (LSMs) to generate multi-scale (with three resolutions of 1, 5, and 10 kilometers) gridded surface heat flux maps for the CASES area. Upon validating these flux maps with measurements from surface station and aircraft, we utilize them to investigate several approaches for estimating the area-integrated surface heat flux for the CASES97 domain of 71x74 square kilometers, which is crucial for land surface model development/validation and area water and energy budget studies. This research is aimed at understanding the relative contribution of random turbulence versus organized mesoscale circulations to the area-integrated surface flux at the scale of 100 kilometers, and identifying the most important effective parameters for characterizing the subgrid-scale variability for large-scale atmosphere-hydrology models.

Chen, Fei↗
DOE OSTI2018

An atomistic investigation of the interaction of dislocations with Guinier-Preston zones in Al-Cu alloys

The interaction between edge dislocations and Guinier-Preston zones in an Al-Cu alloy was analyzed by means of atomistic simulations. The different thermodynamic functions that determine the features of these obstacles for the dislocation glide were computed using molecular statics, molecular dynamics and the nudged elastic band method. It was found that Guinier-Preston zones are sheared by dislocations and the rate at which dislocations overcome the precipitate is controlled by the activation energy, Δ U , in agreement with the postulates of the harmonic transition state theory. Moreover, the entropic contribution to the Helmholtz activation free energy was in the range 1.3–1.8 k b , which can be associated with the typical vibrational entropy of solids. Finally, an estimation of the initial shear flow stress as a function of temperature was carried out from the thermodynamic data provided by the atomistic simulations. Furthermore, comparison with experimental results showed that the effect of the random precipitate distribution and of the dislocation character and dislocation/precipitation orientation has to be taken into account in the simulations to better reproduce experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shear-induced lift force on spheres in a viscous linear shear flow at finite volume fractions

Several studies have shown a significant increase in drag on a distribution of solid spherical particles within a fluid with increasing particle volume fraction. As a result, many empirical drag laws accounting for the dependence on the Reynolds number and volume fraction can be found in the literature. This study investigates the possibility of a similar effect of the particle volume fraction on the mean hydrodynamic lift force on randomly distributed spherical particles in a linear shear flow. Particle-resolved direct numerical simulations are performed to evaluate the mean lift force, and the results are compared with the case of an isolated particle in a linear shear flow for the same Reynolds number and shear rate. The mean lift force acting on the particles appears to remain nearly the same as that on an isolated particle. However, due to the influence of neighboring particles, there is a substantial force variation in transverse directions on each individual particle, whose magnitude is comparable to the mean drag force. The distribution of drag force in a linear shear flow is shown to be nearly the same as in a uniform flow at the same volume fraction and Reynolds number. A simple stochastic model based on a Gaussian distribution is presented for the lift force variation, and its performance is compared to the prediction of the deterministic pairwise interaction extended point-particle model.

42 ENGINEERING↗
DOE OSTI2022

Neural network reconstruction of the dense matter equation of state from neutron star observables

The Equation of State (EoS) of strongly interacting cold and hot ultra-dense QCD matter remains a major challenge in the field of nuclear astrophysics. With the advancements in measurements of neutron star masses, radii, and tidal deformabilities, from electromagnetic and gravitational wave observations, neutron stars play an important role in constraining the ultra-dense QCD matter EoS. Here, in this work, we present a novel method that exploits deep learning techniques to reconstruct the neutron star EoS from mass-radius (M-R) observations. We employ neural networks (NNs) to represent the EoS in a model-independent way, within the range of ~1-7 times the nuclear saturation density. The unsupervised Automatic Differentiation (AD) framework is implemented to optimize the EoS, so as to yield through TOV equations, an M-R curve that best fits the observations. We demonstrate that this method works by rebuilding the EoS on mock data, i.e., mass-radius pairs derived from a randomly generated polytropic EoS. The reconstructed EoS fits the mock data with reasonable accuracy, using just 11 mock M-R pairs observations, close to the current number of actual observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Simulating the filamentation of smoothed laser beams with three-dimensional nonlinear dynamics

In a plasma, the ponderomotive force of an inhomogeneous electromagnetic field expels plasma from regions of high intensity. When a laser propagates through a plasma, this force creates density wells that subsequently modify the index of refraction. The beam refracts and focuses into these wells and may filament. In extreme cases, the laser beam will spray due to increasing angular divergence of the beam. The threshold for ponderomotive self-focusing is well established for isolated laser hotspots or speckles. Here, we define a practical threshold for characterizing the filamentation of thousands of speckles that are found in the focal plane of high-power laser beams spatially smoothed with random phase plates as used at high energy and power laser facilities studying inertial confinement fusion. This threshold is tested against three-dimensional simulations of speckled laser light propagating through plasma. Four metrics are applied to assess filamentation: the fraction of power above five times the average intensity, an effective f-number, the mean-squared perpendicular wavenumber, and the fraction of rarefied density with deviation from the initial density exceeding | δn/ n| = 0.1. The speckled beams studied are generated by random phase plates, both with and without additional polarization smoothing, in a parameter regime of relevance to indirect drive experiments. While filamentation has been discussed extensively in the literature, we believe this to be the first published simulation study with three-dimensional nonlinear hydrodynamics that addresses the onset threshold of ponderomotive filamentation and establishes the lengths and time scales necessary to reach a statistical steady state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Functional traits driving species role in the structure of terrestrial vertebrate scavenger networks

Abstract Species assemblages often have a non‐random nested organization, which in vertebrate scavenger (carrion‐consuming) assemblages is thought to be driven by facilitation in competitive environments. However, not all scavenger species play the same role in maintaining assemblage structure, as some species are obligate scavengers (i.e., vultures) and others are facultative, scavenging opportunistically. We used a database with 177 vertebrate scavenger species from 53 assemblages in 22 countries across five continents to identify which functional traits of scavenger species are key to maintaining the scavenging network structure. We used network analyses to relate ten traits hypothesized to affect assemblage structure with the “role” of each species in the scavenging assemblage in which it appeared. We characterized the role of a species in terms of both the proportion of monitored carcasses on which that species scavenged, or scavenging breadth (i.e., the species “normalized degree”), and the role of that species in the nested structure of the assemblage (i.e., the species “paired nested degree”), therefore identifying possible facilitative interactions among species. We found that species with high olfactory acuity, social foragers, and obligate scavengers had the widest scavenging breadth. We also found that social foragers had a large paired nested degree in scavenger assemblages, probably because their presence is easier to detect by other species to signal carcass occurrence. Our study highlights differences in the functional roles of scavenger species and can be used to identify key species for targeted conservation to maintain the ecological function of scavenger assemblages.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Diaminopurine in Nonenzymatic RNA Template Copying

In the RNA World before the emergence of an RNA polymerase, nonenzymatic template copying would have been essential for the transmission of genetic information. However, the products of chemical copying with the canonical nucleotides (A, U, C, and G) are heavily biased toward the incorporation of G and C, which form a more stable base pair than A and U. We therefore asked whether replacing adenine (A) with diaminopurine (D) might lead to more efficient and less biased nonenzymatic template copying by making a stronger version of the A:U pair. As expected, primer extension substrates containing D bound to U in the template more tightly than substrates containing A. However, primer extension with D exhibited elevated reaction rates on a C template, leading to concerns about fidelity. Our crystallographic studies revealed the nature of the D:C mismatch by showing that D can form a wobble-type base pair with C. We then asked whether competition with G would decrease the mismatched primer extension. We performed nonenzymatic primer extension with all four activated nucleotides on randomized RNA templates containing all four letters and used deep sequencing to analyze the products. We found that the DUCG genetic system exhibited a more even product distribution and a lower mismatch frequency than the canonical AUCG system. Furthermore, primer extension is greatly reduced following all mismatches, including the D:C mismatch. Our study suggests that D deserves further attention for its possible role in the RNA World and as a potentially useful component of artificial nonenzymatic RNA replication systems.

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

Distance‐based reconstruction of protein quaternary structures from inter‐chain contacts

Abstract Predicting the quaternary structure of protein complex is an important problem. Inter‐chain residue‐residue contact prediction can provide useful information to guide the ab initio reconstruction of quaternary structures. However, few methods have been developed to build quaternary structures from predicted inter‐chain contacts. Here, we develop the first method based on gradient descent optimization (GD) to build quaternary structures of protein dimers utilizing inter‐chain contacts as distance restraints. We evaluate GD on several datasets of homodimers and heterodimers using true/predicted contacts and monomer structures as input. GD consistently performs better than both simulated annealing and Markov Chain Monte Carlo simulation. Starting from an arbitrarily quaternary structure randomly initialized from the tertiary structures of protein chains and using true inter‐chain contacts as input, GD can reconstruct high‐quality structural models for homodimers and heterodimers with average TM‐score ranging from 0.92 to 0.99 and average interface root mean square distance from 0.72 Å to 1.64 Å. On a dataset of 115 homodimers, using predicted inter‐chain contacts as restraints, the average TM‐score of the structural models built by GD is 0.76. For 46% of the homodimers, high‐quality structural models with TM‐score ≥ 0.9 are reconstructed from predicted contacts. There is a strong correlation between the quality of the reconstructed models and the precision and recall of predicted contacts. Only a moderate precision or recall of inter‐chain contact prediction is needed to build good structural models for most homodimers. Moreover, GD improves the quality of quaternary structures predicted by AlphaFold2 on a Critical Assessment of Techniques for Protein Structure Prediction–Critical Assessments of Predictions of Interactions dataset.

59 BASIC BIOLOGICAL SCIENCES↗