Search NASA⌕ Search

SEARCH · Search NASA

Results for “grain protein”

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 55 records · Page 3

MARTINI-Compatible Coarse-Grained Model for the Mesoscale Simulation of Peptoids

Peptoids (poly-N-substituted glycines) are a class of synthetic polymers that are regioisomers of peptides (poly-C-substituted glycines), in which the point of side-chain connectivity is shifted from the backbone C to the N atom. Peptoids have found diverse applications as peptidomimetic drugs, protein mimetic polymers, surfactants, and catalysts. Computational modeling is valuable in the understanding and design of peptoid-based nanomaterials. In this work, we report the bottom-up parameterization of coarse-grained peptoid force fields based on the MARTINI peptide force field against all-atom peptoid simulation data. Our parameterization pipeline iteratively refits coarse-grained bonded interactions using iterative Boltzmann inversion and nonbonded interactions by matching the potential of mean force for chain extension. We assure good sampling of the amide bond cis/trans isomerizations in the all-atom simulation data using parallel bias metadynamics. We develop coarse-grained models for two representative peptoids—polysarcosine (poly(N-methyl glycine)) and poly(N-((4-bromophenyl)ethyl)glycine)—and show their structural and thermodynamic properties to be in excellent accord with all-atom calculations but up to 25-fold more efficient and compatible with MARTINI force fields. Here, this work establishes a new rigorously parameterized coarse-grained peptoid force field for the understanding and design of peptoid nanomaterials at length and time scales inaccessible to all-atom calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Asymmetric oligomerization state and sequence patterning can tune multiphase condensate miscibility

Endogenous biomolecular condensates, composed of a multitude of proteins and RNAs, can organize into multiphasic structures with compositionally distinct phases. This multiphasic organization is generally understood to be critical for facilitating their proper biological function. However, the biophysical principles driving multiphase formation are not completely understood. Here we use in vivo condensate reconstitution experiments and coarse-grained molecular simulations to investigate how oligomerization and sequence interactions modulate multiphase organization in biomolecular condensates. We demonstrate that increasing the oligomerization state of an intrinsically disordered protein results in enhanced immiscibility and multiphase formation. Interestingly, we find that oligomerization tunes the miscibility of intrinsically disordered proteins in an asymmetric manner, with the effect being more pronounced when the intrinsically disordered protein, exhibiting stronger homotypic interactions, is oligomerized. Our findings suggest that oligomerization is a flexible biophysical mechanism that cells can exploit to tune the internal organization of biomolecular condensates and their associated biological functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

No evidence for transient transformation via pollen magnetofection in several monocot species

The development of rapid and efficient transformation methods for many plant species remains an obstacle in both the basic and applied plant sciences. A novel method described by Zhao et al. (2017) used magnetic nanoparticles to deliver DNA into pollen grains of several dicot species, and one monocot (lily), to achieve transformation (“pollen magnetofection”). Using the published protocol, extensive trials by two independent research groups showed no indication of transient transformation success with pollen from two monocots, maize and sorghum. To further address the feasibility of magnetofection, lily pollen was used for side-by-side trials of magnetofection with a proven methodology for transient transformation, biolistics. Using a Green Fluorescent Protein reporter plasmid, transformation efficiency with the biolistic approach averaged 0.7% over three trials. However, the same plasmid produced no recognizable transformants via magnetofection, despite screening >3500 individual pollen grains. We conclude that pollen magnetofection is not effective for transient transformation of pollen for at least three species of monocots, and suggest that efforts to replicate the magnetofection protocol in dicot species would be useful to fully assess its potential. ARISING FROM Zhao et al. Nature Plants https://doi.org/10.1038/s41477-017-0063-z (2017)

Vejlupkova, Zuzana↗

Ensuring thermodynamic consistency with invertible coarse-graining

Coarse-grained models are a core computational tool in theoretical chemistry and biophysics. A judicious choice of a coarse-grained model can yield physical insights by isolating the essential degrees of freedom that dictate the thermodynamic properties of a complex, condensed-phase system. The reduced complexity of the model typically leads to lower computational costs and more efficient sampling compared with atomistic models. Designing “good” coarse-grained models is an art. Generally, the mapping from fine-grained configurations to coarse-grained configurations itself is not optimized in any way; instead, the energy function associated with the mapped configurations is. In this work, we explore the consequences of optimizing the coarse-grained representation alongside its potential energy function. We use a graph machine learning framework to embed atomic configurations into a low-dimensional space to produce efficient representations of the original molecular system. Because the representation we obtain is no longer directly interpretable as a real-space representation of the atomic coordinates, we also introduce an inversion process and an associated thermodynamic consistency relation that allows us to rigorously sample fine-grained configurations conditioned on the coarse-grained sampling. We illustrate that this technique is robust, recovering the first two moments of the distribution of several observables in proteins such as chignolin and alanine dipeptide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Coarse-Grained Model of Affinity Maturation Indicates the Importance of B-Cell Receptor Avidity in Epitope Subdominance

The elicitation of broadly neutralizing antibodies (bnAbs) is a major goal in the design of vaccines against rapidly-mutating viruses. In the case of influenza, many bnAbs that target conserved epitopes on the stem of the hemagglutinin protein (HA) have been discovered. However, these antibodies are rare, are not boosted well upon reinfection, and often have low neutralization potency, compared to strain-specific antibodies directed to the HA head. Different hypotheses have been proposed to explain this phenomenon. We use a coarse-grained computational model of the germinal center reaction to investigate how B-cell receptor binding valency affects the growth and affinity maturation of competing B-cells. We find that receptors that are unable to bind antigen bivalently, and also those that do not bind antigen cooperatively, have significantly slower rates of growth, memory B-cell production, and, under certain conditions, rates of affinity maturation. The corresponding B-cells are predicted to be outcompeted by B-cells that bind bivalently and cooperatively. We use the model to explore strategies for a universal influenza vaccine, e.g., how to boost the concentrations of the slower growing cross-reactive antibodies directed to the stem. The results suggest that, upon natural reinfections subsequent to vaccination, the protectiveness of such vaccines would erode, possibly requiring regular boosts. Collectively, our results strongly support the importance of bivalent antibody binding in immunodominance, and suggest guidelines for developing a universal influenza vaccine.

60 APPLIED LIFE SCIENCES↗

Challenges of conventional iterative all-atom and coarse-grained multiscale molecular dynamics

In this work, we evaluate the biomolecular dynamics behaviors when conventionally iterating between all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations over multiple cycles. We implemented the workflow to iterate between AA and CG in OpenMM, namely the iterative multiscale MD (iMMD) simulation workflow. In particular, we aim to identify practical applications for iterating between AA and CG simulations in a conventional manner without any constraints or model modifications. We evaluate the iMMD workflow on four representative systems, spanning folding of two soluble proteins and protein-protein as well as protein-lipid interactions of two membrane proteins. We observe that iteration between AA and CG representations could help the soluble proteins exit undesirable metastable states to fold, resulting from random protein structural distortions due to cycling. Consequently, the most reliable use of iterative AA and CG simulations appears to be to accelerating complex lipid mixing for membrane-bound protein systems rather than sampling protein conformational space. Our work explores the practical usages and limitations for iterative AA and CG simulations using readily available AA and CG force fields. The evaluated iMMD workflow in OpenMM is made available at https://github.com/lanl/iMMD.

59 BASIC BIOLOGICAL SCIENCES↗

Assembly and Analysis of Cell-Scale Membrane Envelopes

The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms-a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hexagonal Lattices of HIV Capsid Proteins Explored by Simulations Based on a Thermodynamically Consistent Model

HIV capsid proteins (CAs) may self-assemble into a variety of shapes under in vivo and in vitro conditions. Here, we employed simulations based on a residue-level coarse-grained (CG) model with full conformational flexibility to investigate hexagonal lattices, which are the underlying structural pattern for CA aggregations. Facilitated by enhanced sampling simulations to rigorously calculate CA dimerization and polymerization affinities, we calibrated our model to reproduce the experimentally measured affinities. Here using the calibrated model, we performed unbiased simulations on several large systems consisting of 1512 CA subunits, allowing reversible binding and unbinding of the CAs in a thermodynamically consistent manner. In one simulation, a preassembled hexagonal CA sheet developed spontaneous curvatures reminiscent of those observed in experiments, and the edges of the sheet exhibited local curvatures larger than those of the interior. In other simulations starting with randomly distributed CAs at different concentrations, existing CA assemblies grew by binding free capsomeres to the edges and by merging with other assemblies. At high CA concentrations, rapid establishment of predominant aggregates was followed by much slower adjustments toward more regular hexagonal lattices, with increasing numbers of intact CA hexamers and pentamers being formed. Our approach of adapting a general CG model to specific systems by using experimental binding data represents a practical and effective strategy for simulating and elucidating intricate protein aggregations.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling Interactions within and between Peptide Amphiphile Supramolecular Filaments

Many peptides are able to self-assemble into one-dimensional (1D) nanostructures, such as cylindrical fibers or ribbons of variable widths, but the relationship between the morphology of 1D objects and their molecular structure is not well understood. Here, we use coarse-grained molecular dynamics (CG-MD) simulations to study the nanostructures formed by self-assembly of different peptide amphiphiles (PAs). The results show that ribbons are hierarchical superstructures formed by laterally assembled cylindrical fibers. Simulations starting from bilayer structures demonstrate the formation of filaments, whereas other simulations starting from filaments indicate varying degrees of interaction among them depending on chemical structure. These interactions are verified by observations using atomic force microscopy of the various systems. The interfilament interactions are predicted to be strongest in supramolecular assemblies that display hydrophilic groups on their surfaces, while those with hydrophobic ones are predicted to interact more weakly as confirmed by viscosity measurements. The simulations also suggest that peptide amphiphiles with hydrophobic termini bend to reduce their interfacial energy with water, which may explain why these systems do not collapse into superstructures of bundled filaments. The simulations suggest that future experiments will need to address mechanistic questions about the self-assembly of these systems into hierarchical structures, namely, the preformation of interactive filaments vs equilibration of large assemblies into superstructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗

Identification of Novel Genomic Associations and Gene Candidates for Grain Starch Content in Sorghum

Starch accumulated in the endosperm of cereal grains as reserve energy for germination serves as a staple in human and animal nutrition. Unraveling genetic control for starch metabolism is important for breeding grains with high starch content. In this study, we used a sorghum association panel with 389 individuals and 141,557 single nucleotide polymorphisms (SNPs) to fit linear mixed models (LMM) for identifying genomic regions and potential candidate genes associated with starch content. Three associated genomic regions, one in chromosome (chr) 1 and two novel associations in chr-8, were identified using combination of LMM and Bayesian sparse LMM. All significant SNPs were located within protein coding genes, with SNPs ∼ 52 Mb of chr-8 encoding a Casperian strip membrane protein (CASP)-like protein (Sobic.008G111500) and a heat shock protein (HSP) 90 (Sobic.008G111600) that were highly expressed in reproductive tissues including within the embryo and endosperm. The HSP90 is a potential hub gene with gene network of 75 high-confidence first interactors that is enriched for five biochemical pathways including protein processing. The first interactors of HSP90 also showed high transcript abundance in reproductive tissues. The candidates of this study are likely involved in intricate metabolic pathways and represent candidate gene targets for source-sink activities and drought and heat stress tolerance during grain filling.

genomics↗

Cooperative multivalent receptor binding promotes exposure of the SARS-CoV-2 fusion machinery core

The molecular events that permit the spike glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to bind and enter cells are important to understand for both fundamental and therapeutic reasons. Spike proteins consist of S1 and S2 domains, which recognize angiotensin-converting enzyme 2 (ACE2) receptors and contain the viral fusion machinery, respectively. Ostensibly, the binding of spike trimers to ACE2 receptors promotes dissociation of the S1 domains and exposure of the fusion machinery, although the molecular details of this process have yet to be observed. We report the development of bottom-up coarse-grained (CG) models consistent with cryo-electron tomography data, and the use of CG molecular dynamics simulations to investigate viral binding and S2 core exposure. We show that spike trimers cooperatively bind to multiple ACE2 dimers at virion-cell interfaces in a manner distinct from binding between soluble proteins, which processively induces S1 dissociation. We also simulate possible variant behavior using perturbed CG models, and find that ACE2-induced S1 dissociation is primarily sensitive to conformational state populations and the extent of S1/S2 cleavage, rather than ACE2 binding affinity. These simulations reveal an important concerted interaction between spike trimers and ACE2 dimers that primes the virus for membrane fusion and entry.

59 BASIC BIOLOGICAL SCIENCES↗

From Farm to Flight: CoverCress as a Low Carbon Intensity Cash Cover Crop for Sustainable Aviation Fuel Production. A Review of Progress Towards Commercialization

Thlaspi arvense L. (Field Pennycress; pennycress) is being converted into a winter-annual oilseed crop that confers cover crop benefits when grown throughout the 12 million-hectares U.S. Midwest. To ensure a fit with downstream market demand, conversion involves not only improvements in yield and maturity through traditional breeding, but also improvements in the composition of the oil and protein through gene editing tools. The conversion process is similar to the path taken to convert rapeseed into Canola. In the case of field pennycress, the converted product that is suitable as a rotational crop is called CoverCress™ as marketed by CoverCress Inc. or golden pennycress if marketed by others. Off-season integration of a CoverCress crop into existing corn and soybean hectares would extend the growing season on established croplands and avoid displacement of food crops or ecosystems while yielding up to 1 billion liters of seed oil annually by 2030, with the potential to grow to 8 billion liters from production in the U.S. Midwest alone. The aviation sector is committed to carbon-neutral growth and reducing emissions of its global market, which in 2019 approached 122 billion liters of consumption in the U.S. and 454 billion liters globally. The oil derived from a CoverCress crop is ideally suited as a new bioenergy feedstock for the production of drop-in Sustainable Aviation Fuel (SAF), renewable diesel, biodiesel and other value-added coproducts. Through a combination of breeding and genomics-enabled mutagenesis approaches, considerable progress has been made in genetically improving yield and other agronomic traits. With USDA-NIFA funding and continued public and private investments, improvements to CoverCress germplasm and agronomic practices suggest that field-scale production can surpass 1,680 kg ha -1 (1,500 lb ac -1 ) in the near term. At current commodity prices, economic modeling predicts this level of production can be profitable across the entire supply chain. Two-thirds of the grain value is in oil converted to fuels and chemicals, and the other one-third is in the meal used as an animal feed, industrial applications, and potential plant-based protein products. In addition to strengthening rural communities by providing income to producers and agribusinesses, cultivating a CoverCress crop potentially offers a myriad of ecosystem services. The most notable service is water quality protection through reduced nutrient leaching and reduced soil erosion. Biodiversity enhancement by supporting pollinators’ health is also a benefit. While the efforts described herein are focused on the U.S., cultivation of a CoverCress crop will likely have a broader application to regions around the world with similar agronomic and environmental conditions.

09 BIOMASS FUELS↗

Comprehensive Genome-Wide Natural Variation and Expression Analysis of Tubby-like Proteins Gene Family in Brachypodium distachyon

The Tubby-like proteins (TLPs) gene family is a group of transcription factors found in both animals and plants. In this study, we identified twelve B. distachyon TLPs, divided into six groups based on conserved domains and evolutionary relationships. We predicted cis-regulatory elements involved in light, hormone, and biotic and abiotic stresses. The expression patterns in response to light and hormones revealed that BdTLP3, 4, 7, and 14 are involved in light responses, and BdTLP1 is involved in ABA responses. Furthermore, BdTLP2, 7, 9, and 13 are expressed throughout vegetative and reproductive stages, whereas BdTLP1, 3, 5, and 14 are expressed at germinating grains and early vegetative development, and BdTLP4, 6, 8, and 10 are expressed at the early reproduction stage. The natural variation in the eleven most diverged B. distachyon lines revealed high conservation levels of BdTLP1-6 to high variation in BdTLP7-14 proteins. Based on diversifying selection, we identified amino acids in BdTLP1, 3, 8, and 13, potentially substantially affecting protein functions. This analysis provided valuable information for further functional studies to understand the regulation, pathways involved, and mechanism of BdTLPs.

Plant Sciences↗

The confluence of machine learning and multiscale simulations

Multiscale modeling has a long history of use in structural biology, as computational biologists strive to overcome the time- and length-scale limits of atomistic molecular dynamics. Contemporary machine learning techniques, such as deep learning, have promoted advances in virtually every field of science and engineering and are revitalizing the traditional notions of multiscale modeling. Deep learning has found success in various approaches for distilling information from fine-scale models, such as building surrogate models and guiding the development of coarse-grained potentials. However, perhaps its most powerful use in multiscale modeling is in defining latent spaces that enable efficient exploration of conformational space. In conclusion, this confluence of machine learning and multiscale simulation with modern high-performance computing promises a new era of discovery and innovation in structural biology.

59 BASIC BIOLOGICAL SCIENCES↗

Reweighting configurations generated by transferable, machine learned models for protein sidechain backmapping

Multiscale modeling requires the linking of models at different levels of detail, with the goal of gaining accelerations from lower fidelity models while recovering fine details from higher resolution models. Communication across resolutions is particularly important in modeling soft matter, where tight couplings exist between molecular-level details and mesoscale structures. While multiscale modeling of biomolecules has become a critical component in exploring their structure and self-assembly, backmapping from coarse-grained to fine-grained, or atomistic, representations presents a challenge, despite recent advances through machine learning. A major hurdle, especially for strategies utilizing machine learning, is that backmappings can only approximately recover the atomistic ensemble of interest. We demonstrate conditions for which backmapped configurations may be reweighted to exactly recover the desired atomistic ensemble. By training separate decoding models for each sidechain type, we develop an algorithm based on normalizing flows and geometric algebra attention to autoregressively propose backmapped configurations for any protein sequence. Critical for reweighting with modern protein force fields, our trained models include all hydrogen atoms in the backmapping and make probabilities associated with atomistic configurations directly accessible. We also demonstrate, however, that reweighting is extremely challenging despite state-of-the-art performance on recently developed metrics and generation of configurations with low energies in atomistic protein force fields. Through detailed analysis of configurational weights, we show that machine-learned backmappings must not only generate configurations with reasonable energies, but also correctly assign relative probabilities under the generative model. These are broadly important considerations in generative modeling of atomistic molecular configurations.

Monroe, Jacob I. [Univ. of Arkansas, Fayetteville,↗

Moving the needle: Employing deep reinforcement learning to push the boundaries of coarse-grained vaccine models

Highly mutable infectious disease pathogens (hm-IDPs) such as HIV and influenza evolve faster than the human immune system can contain them, allowing them to circumvent traditional vaccination approaches and causing over one million deaths annually. Agent-based models can be used to simulate the complex interactions that occur between immune cells and hm-IDP-like proteins (antigens) during affinity maturation—the process by which antibodies evolve. Compared to existing experimental approaches, agent-based models offer a safe, low-cost, and rapid route to study the immune response to vaccines spanning a wide range of design variables. However, the highly stochastic nature of affinity maturation and vast sequence space of hm-IDPs render brute force searches intractable for exploring all pertinent vaccine design variables and the subset of immunization protocols encompassed therein. To address this challenge, we employed deep reinforcement learning to drive a recently developed agent-based model of affinity maturation to focus sampling on immunization protocols with greater potential to improve the chosen metrics of protection, namely the broadly neutralizing antibody (bnAb) titers or fraction of bnAbs produced. Using this approach, we were able to coarse-grain a wide range of vaccine design variables and explore the relevant design space. Our work offers new testable insights into how vaccines should be formulated to maximize protective immune responses to hm-IDPs and how they can be minimally tailored to account for major sources of heterogeneity in human immune responses and various socioeconomic factors. Our results indicate that the first 3 to 5 immunizations, depending on the metric of protection, should be specially tailored to achieve a robust protective immune response, but that beyond this point further immunizations require only subtle changes in formulation to sustain a durable bnAb response.

60 APPLIED LIFE SCIENCES↗

Arbitrarily accurate, nonparametric coarse graining with Markov renewal processes and the Mori–Zwanzig formulation

Stochastic dynamics, such as molecular dynamics, are important in many scientific applications. However, summarizing and analyzing the results of such simulations is often challenging due to the high dimension in which simulations are carried out and, consequently, due to the very large amount of data that are typically generated. Coarse graining is a popular technique for addressing this problem by providing compact and expressive representations. Coarse graining, however, potentially comes at the cost of accuracy, as dynamical information is, in general, lost when projecting the problem in a lower-dimensional space. This article shows how to eliminate coarse-graining error using two key ideas. First, we represent coarse-grained dynamics as a Markov renewal process. Second, we outline a data-driven, non-parametric Mori–Zwanzig approach for computing jump times of the renewal process. Numerical tests on a small protein illustrate the method.

36 MATERIALS SCIENCE↗