Search NASA⌕ Search

SEARCH · Search NASA

Results for “native state simulation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Electrostatic Influence on Photoisomerization in Bacteriorhodopsin and Halorhodopsin

Bacteriorhodopsin (bR) and halorhodopsin (hR) are mutual membrane proteins that transport ions across the cell membrane in halobacteria. Their ion transport function is triggered by photoactivated isomerization of the retinal protonated Schiff base (RPSB) chromophore. In spite of their similar structures, bR and hR exhibit widely differing RPSB isomerization rates and quantum yields (with bR being both faster and more efficient than hR). Previous simulations of photoisomerization in bR and hR using ab initio multiple spawning (AIMS) with QM/MM have successfully reproduced the experimentally observed ordering of quantum yields and isomerization rates, but the origin of these differences remains elusive. In this work, we discover the role of electrostatic interactions in the protein pocket surrounding RPSB. We probe the influence of protein electrostatics by modifying the charge of the complex counterion in bR/hR to be more/less negative than the native state. We find that such modifications lead to bR-like behavior in hR and vice versa. This demonstrates the crucial role of electrostatic interactions in controlling the outcome of RPSB photoisomerization.

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

DTK C/Fortran Interface Development for NEAMS FSI Simulations

This report documents the development of DataTransferKit (DTK) C and Fortran interfaces for fluid-structure-interaction (FSI) simulations in NEAMS. In these simulations, the codes Nek5000 and Diablo are being coupled within the SHARP framework to study flow-induced vibration (FIV) in reactor steam generators. We will review the current Nek5000/Diablo coupling algorithm in SHARP and the current state of the solution transfer scheme used in this implementation. We will then present existing DTK algorithms which may be used instead to provide an improvement in both flexibility and scalability of the current SHARP implementation. We will show how these can be used within the current FSI scheme using a new set of interfaces to the algorithms developed by this work. These new interfaces currently expose the mesh-free solution transfer algorithms in DTK, a C++ library, and are written in C and Fortran to enable coupling of both Nek5000 and Diablo in their native Fortran language. They have been compiled and tested on Cooley, the test-bed machine for Mira at ALCF.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2010

Thermal-induced conformational changes in the product release area drive the enzymatic activity of xylanases 10B: Crystal structure, conformational stability and functional characterization of the xylanase 10B from Thermotoga petrophila RKU-1

Research highlights: {yields} The hyperthermostable xylanase 10B from Thermotoga petrophila RKU-1 produces exclusively xylobiose at the optimum temperature. {yields} Circular dichroism spectroscopy suggests a coupling effect of temperature-induced structural changes with its enzymatic behavior. {yields} Crystallographic and molecular dynamics studies indicate that conformational changes in the product release area modulate the enzyme action mode. -- Abstract: Endo-xylanases play a key role in the depolymerization of xylan and recently, they have attracted much attention owing to their potential applications on biofuels and paper industries. In this work, we have investigated the molecular basis for the action mode of xylanases 10B at high temperatures using biochemical, biophysical and crystallographic methods. The crystal structure of xylanase 10B from hyperthermophilic bacterium Thermotoga petrophila RKU-1 (TpXyl10B) has been solved in the native state and in complex with xylobiose. The complex crystal structure showed a classical binding mode shared among other xylanases, which encompasses the -1 and -2 subsites. Interestingly, TpXyl10B displayed a temperature-dependent action mode producing xylobiose and xylotriose at 20 {sup o}C, and exclusively xylobiose at 90 {sup o}C as assessed by capillary zone electrophoresis. Moreover, circular dichroism spectroscopy suggested a coupling effect of temperature-induced structural changes with this particular enzymatic behavior. Molecular dynamics simulations supported the CD analysis suggesting that an open conformational state adopted by the catalytic loop (Trp297-Lys326) provokes significant modifications in the product release area (+1,+2 and +3 subsites), which drives the enzymatic activity to the specific release of xylobiose at high temperatures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Metallic support films reduce optical heating in cryogenic correlative light and electron tomography

Super-resolved cryogenic correlative light and electron tomography is an emerging method that provides both the single-molecule sensitivity and specificity of fluorescence imaging, and the molecular scale resolution and detailed cellular context of tomography, all in vitrified cells preserved in their native hydrated state. Technical hurdles that limit these correlative experiments need to be overcome for the full potential of this approach to be realized. Chief among these is sample heating due to optical excitation which leads to devitrification, a phase transition from amorphous to crystalline ice. Here, in this work, we show that much of this heating is due to the material properties of the support film of the electron microscopy grid, specifically the absorptivity and thermal conductivity. We demonstrate through experiment and simulation that the properties of the standard holey carbon electron microscopy grid lead to substantial heating under optical excitation. In order to avoid devitrification, optical excitation intensities must be kept orders of magnitude lower than the intensities commonly employed in room temperature super-resolution experiments. We further show that the use of metallic films, either holey gold grids, or custom made holey silver grids, alleviate much of this heating. For example, the holey silver grids permit 20× the optical intensities used on the standard holey carbon grids. Super-resolution correlative experiments conducted on holey silver grids under these increased optical excitation intensities have a corresponding increase in the rate of single-molecule fluorescence localizations. This results in an increased density of localizations and improved correlative imaging without deleterious effects from sample heating.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Thermal unfolding of barstar and the properties of interfacial water around the unfolded forms

Identification of the intermediates along the folding-unfolding pathways and probing their interactions with surrounding solvent are two important but relatively unexplored issues in protein folding. In this work, we have carried out atomistic molecular dynamics simulations to study the thermal unfolding of barstar in aqueous solution from its folded native form at two different temperatures (400 K and 450 K). The calculations at 400 K reveal partial unfolding of two α-helices (helix-1 and helix-2) and their interconnecting loop. At 450 K, on the other hand, the entire protein attains an expanded flexible conformation due to disruption of a large fraction of tertiary contacts and breaking of almost all the secondary structures. These two disordered structures obtained at such high temperatures are then studied around room temperature to probe their influence on the properties of surrounding solvent. It is found that though the unfolding of the protein in general leads to increasingly hydrated interface, but new structural motifs with locally dehydrated interface may also form during the structural transition. Additionally, independent of the conformational state of the protein, its influence on surrounding solvent has been found to be restricted to the first hydration layer.

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

Matrix Approach to Accelerate Spin-Up of CLM5

Numerical models have been developed to investigate and understand responses of biogeochemical cycle to global changes. Steady state, when a system is in dynamic equilibrium, is generally required to initialize these model simulations. However, the spin-up process that is used to achieve steady state pose a great burden to computational resources, limiting the efficiency of global modeling analysis on biogeochemical cycles. This study introduces a new Semi-Analytical Spin-Up (SASU) to tackle this grand challenge. We applied SASU to Community Land Model version 5 and examined its computational efficiency and accuracy. At the Brazil site, SASU is computationally 7 times more efficient than (or saved up to 86% computational cost in comparison with) the traditional native dynamics (ND) spin-up to reach the same steady state. Globally, SASU is computationally 8 times more efficient than the accelerated decomposition spin-up and 50 times more efficient than ND. In summary, SASU achieves the highest computational efficiency for spin-up on site and globally in comparison with other spin-up methods. It is generalizable to wide biogeochemical models and thus makes computationally costly studies (e.g., parameter perturbation ensemble analysis and data assimilation) possible for a better understanding of biogeochemical cycle under climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

On the analysis and comparison of conformer-specific essential dynamics upon ligand binding to a protein

The native state of a protein consists of an equilibrium of conformational states on an energy landscape rather than existing as a single static state. The co-existence of conformers with different ligand-affinities in a dynamical equilibrium is the basis for the conformational selection model for ligand binding. In this context, the development of theoretical methods that allow us to analyze not only the structural changes but also changes in the fluctuation patterns between conformers will contribute to elucidate the differential properties acquired upon ligand binding. Molecular dynamics simulations can provide the required information to explore these features. Its use in combination with subsequent essential dynamics analysis allows separating large concerted conformational rearrangements from irrelevant fluctuations. We present a novel procedure to define the size and composition of essential dynamics subspaces associated with ligand-bound and ligand-free conformations. These definitions allow us to compare essential dynamics subspaces between different conformers. Our procedure attempts to emphasize the main similarities and differences between the different essential dynamics in an unbiased way. Essential dynamics subspaces associated to conformational transitions can also be analyzed. As a test case, we study the glutaminase interacting protein (GIP), composed of a single PDZ domain. Both GIP ligand-free state and glutaminase L peptide-bound states are analyzed. Our findings concerning the relative changes in the flexibility pattern upon binding are in good agreement with experimental Nuclear Magnetic Resonance data.

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

Single-Molecule Detection of the Encounter and Productive Electron Transfer Complexes of a Photosynthetic Reaction Center

Small, diffusible redox proteins play an essential role in electron transfer (ET) in respiration and photosynthesis, sustaining life on Earth by shuttling electrons between membrane-bound complexes via finely tuned and reversible interactions. Ensemble kinetic studies show transient ET complexes form in two distinct stages: an “encounter” complex largely mediated by electrostatic interactions, which subsequently, through subtle reorganization of the binding interface, forms a “productive” ET complex stabilized by additional hydrophobic interactions around the redox-active cofactors. Here, using single-molecule force spectroscopy (SMFS) we dissected the transient ET complexes formed between the photosynthetic reaction center-light harvesting complex 1 (RC-LH1) of Rhodobacter sphaeroides and its native electron donor cytochrome c 2 (cyt c 2 ). Importantly, SMFS resolves the distribution of interaction forces into low (~150 pN) and high (~330 pN) components, with the former more susceptible to salt concentration and to alteration of key charged residues on the RC. Thus, the low force component is suggested to reflect the contribution of electrostatic interactions in forming the initial encounter complex, whereas the high force component reflects the additional stabilization provided by hydrophobic interactions to the productive ET complex. Employing molecular dynamics simulations, we resolve five intermediate states that comprise the encounter, productive ET and leaving complexes, predicting a weak interaction between cyt c 2 and the LH1 ring near the RC-L subunit that could lie along the exit path for oxidized cyt c 2 . The multimodal nature of the interactions of ET complexes captured here may have wider implications for ET in all domains of life.

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

Activation dynamics of a water-soluble human mu-opioid receptor

The mu-opioid receptor (MOR), a class A G protein-coupled receptor mediates opioid analgesia and remains a central target for pain therapeutics. While crystal structures of MOR exist, they provide limited insight into the receptor’s dynamic conformational landscape underlying function. Here, we engineered a thermostable water-soluble MOR variant (wsMOR) that retains native-like ligand-binding and activation dynamics. This variant enables high-yield production and detailed solution-phase structural studies that are challenging with membrane-embedded MOR, providing a valuable tool for studying receptor activation and aqueous-phase drug screening. Using a combined computational and experimental approach, we performed long-timescale all-atom molecular dynamics simulations together with neutron scattering and single-molecule FRET, revealing a structurally stable receptor with a diverse ensemble of conformations at different temporal resolutions. In the ligand-free state, wsMOR displayed high conformational flexibility, which decreased upon agonist binding, particularly in transmembrane helix 6, a hallmark of G protein-coupled receptor activation. Positive allosteric modulation and G protein binding further stabilized active-like states. These findings highlight wsMOR’s conformational plasticity across picosecond to millisecond timescales and provide a foundation for structure-guided development of next-generation opioid ligands with improved efficacy and safety.

E, Agyemang [University of Tennessee Knoxville]↗
DOE OSTI2019

First-Principles Characterization of the Elusive I Fluorescent State and the Structural Evolution of Retinal Protonated Schiff Base in Bacteriorhodopsin

The conversion of light energy into work is essential to life on earth. Bacteriorhodopsin (bR), a light-activated proton pump in Archae, has served for many years as a model system for the study of this process in photoactive proteins. Upon absorption of a photon, its chromophore, the retinal protonated Schiff base (RPSB), isomerizes from its native all-trans form to a 13-cis form and pumps a proton out of the cell in a process that is coupled to eventual ATP synthesis. Despite numerous time-resolved spectroscopic studies over the years, the details of the photodynamics of bR on the excited state, particularly the characterization of the I fluorescent state, the time-resolved reaction mechanism, and the role of the counterion cluster of RPSB, remain uncertain. Here, in this work, we use ab initio multiple spawning (AIMS) with spin-restricted ensemble Kohn–Sham (REKS) theory to simulate the nonadiabatic dynamics of the ultrafast photoreaction in bR. The excited state dynamics can be partitioned into three distinct phases: (1) relaxation away from the Franck–Condon region dominated by changes in retinal bond length alternation, (2) dwell time on the excited state in the I fluorescent state featuring an untwisted, bond length inverted RPSB, and (3) rapid torsional evolution to the conical intersection after overcoming a small excited state barrier. We fully characterize the I fluorescent state and the excited state barrier that hinders direct evolution to the conical intersection following photoexcitation. We also find that photoisomerization is accompanied by weakening of the interaction between RPSB and its counterion cluster. However, in contradiction with a recent time-resolved X-ray experiment, hydrogen bond cleavage is not necessary to reproduce the observed photoisomerization dynamics.

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

Control And Optimization Modular Modeling Application For Nuclear Deployment

The purpose of the COMMAND code is to provide a flexible, scalable tool for use in developing, integrating, and testing the technologies necessary for achieving autonomous operations of advanced nuclear reactors. The code enables users to efficiently implement custom simulations and experiments by combining key methods from different software modules. These modules are focused on: modeling and simulation tools, such as nuclear simulation tools used for high-fidelity modeling (e.g., Reactor Excursion and Leak Analysis Program [RELAP5-3D] and Monte Carlo N-Particle [MCNP]); machine learning and optimization tools (e.g., anomaly detection and data-driven modeling techniques); advanced control in its digital, high-performance, and supervisory control forms (e.g., proportional integral derivative (PID) control and model predictive control (MPC); and integration with hardware through industrial communication protocols. To ensure flexibility and scalability, COMMAND was designed to be both modular—the software “pieces” all inherit from generic building blocks and can be combined and connected to create complicated simulations—and high performing—designed for parallel processing, enabling simulations and experiments to take advantage of multi-core computers, servers, and nodes. The code is written in the Python programming language due to the language's popularity, active community, and open-source and cross-platform nature. Maintaining consistency with other simulation tools used within the nuclear energy community, users implement simulations and experiments through text input files, which define components, parameters, connections, etc., through lines of text. Given that COMMAND is written in Python, these input files are native Python scripts, and so use the standard Python structure and formatting. This also enables users to take advantage of Python's extensive package library to develop custom capabilities for their specific use cases.

Faber, Jacob [Idaho National Laboratory (INL), Ida↗
DOE OSTI2017

Collapse Precedes Folding in Denaturant-Dependent Assembly of Ubiquitin

The folding of small protein ubiquitin (Ub), which plays an indispensable role in targeting proteins for degradation and DNA damage response, is complex. A number of experiments on Ub folding have reached differing conclusions regarding the relation between collapse and folding, and whether intermediates are populated. In order to resolve these vexing issues, we elucidate the denaturant-dependent thermodynamics and kinetics of Ub folding at low and neutral pH as a function of guanidinium chloride and urea using coarse-grained molecular simulations. The changes in the fraction of the folded Ub, and the radius of gyration (R g ) as a function of the denaturant concentration, [C], are in quantitative agreement with experiments. Under conditions used in experiments, R g of the unfolded state at neutral pH changes only by ≈17% as the [GdmCl] decreases from 6 to 0 M. We predict that the extent of compaction of the unfolded state increases as temperature decreases. A two-dimensional folding landscape as a function of R g and a measure of similarity to the folded state reveals unambiguously that the native state assembly is preceded by collapse, as discovered in fast mixing experiments on several proteins. Assessment of the folding trajectories, under mildly denaturing conditions ([GdmCl] = 1.0 M or [Urea] = 1.0 M), shows that Ub folds by collision between preformed secondary structural elements involving kinetic intermediates that are primarily stabilized by long-range contacts. Here, we explain the results of small angle X-ray scattering (SAXS) experiments on Ub quantitatively, and establishes that evolved globular proteins in the unfolded ensemble are poised to collapse as the solvent conditions for the biopolymer changes from good solvent to Θ-solvent like conditions on denaturant dilution. In the process, we explain the discrepancy between SAXS and single molecule fluorescent resonant energy transfer (smFRET) experiments, which have arrived at a contradicting conclusion concerning the collapse of polypeptide chains.

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

Thermal Spectrum Molten Salt-Fueled Reactor Reference Plant Model

This report details a fully coupled neutronics thermal hydraulics reference plant model for a thermal spectrum molten salt-fueled reactor. The Multiphysics model is developed on the Nuclear Regulatory Commission's Comprehensive Reactor Analysis Bundle available on the Idaho National Laboratory high-performance computer, which natively and seamlessly couples Griffin and Sam codes of the Multiphysics Object-Oriented Simulation Environment (MOOSE) based applications. Griffin provides the reactor physics capabilities, k-eigenvalue, delayed neutron precursor concentrations, and transient solutions for stationary and flowing fuels. Sam provides the solution of mass, momentum, and energy conservation equations of the whole system with a set of system components for one-dimensional single-phase flow. The neutronics feedback model relies primarily on fuel and moderator temperatures, fuel salt density, and delayed neutron precursors distributions. This report presents our results for the coupled steady state core and an unprotected loss of flow event. Although this model is prototypical regarding capabilities in the Comprehensive Reactor Analysis Bundle, its results are consistent with experimental data of the pump transient tests that were performed at zero power of the Molten Salt Reactor Experiment (MSRE).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2013

Translating MAPGEN to ASPEN for MER

This software translates MAPGEN (Europa and APGEN) domains to ASPEN, and the resulting domain can be used to perform planning for the Mars Exploration Rover (MER). In other words, this is a conversion of two distinct planning languages (both declarative and procedural) to a third (declarative) planning language in order to solve the problem of faithful translation from mixed-domain representations into the ASPEN Modeling Language. The MAPGEN planning system is an example of a hybrid procedural/declarative system where the advantages of each are leveraged to produce an effective planner/scheduler for MER tactical planning. The adaptation of the planning system (ASPEN) was investigated, and, with some translation, much of the procedural knowledge encoding is amenable to declarative knowledge encoding. The approach was to compose translators from the core languages used for adapting MAGPEN, which consists of Europa and APGEN. Europa is a constraint- based planner/scheduler where domains are encoded using a declarative model. APGEN is also constraint-based, in that it tracks constraints on resources and states and other variables. Domains are encoded in both constraints and code snippets that execute according to a forward sweep through the plan. Europa and APGEN communicate to each other using proxy activities in APGEN that represent constraints and/or tokens in Europa. The composition of a translator from Europa to ASPEN was fairly straightforward, as ASPEN is also a declarative planning system, and the specific uses of Europa for the MER domain matched ASPEN s native encoding fairly closely. On the other hand, translating from APGEN to ASPEN was considerably more involved. On the surface, the types of activities and resources one encodes in APGEN appear to match oneto- one to the activities, state variables, and resources in ASPEN. But, when looking into the definitions of how resources are profiled and activities are expanded, one sees code snippets that access various information available during planning for the moment in time being planned to decide at the time what the appropriate profile or expansion is. APGEN is actually a forward (in time) sweeping discrete event simulator, where the model is composed of code snippets that are artfully interleaved by the engine to produce a plan/schedule. To solve this problem, representative code is simulated as a declarative series of task expansions. Predominantly, three types of procedural models were translated: loops, if statements, and code blocks. Loops and if statements were handled using controlled task expansion, and code blocks were handled using constraint networks that maintained the generation of results based on what the order of execution would be for a procedural representation. One advantage with respect to performance for MAPGEN is the use of APGEN s GUI. This GUI is written in C++ and Motif, and performs very well for large plans.

Rabideau, Gregg R.↗
DOE OSTI2019

Development of a Stochastic Infiltration Model for PPA Implementation - 19293

Performance Assessment models often require inputs representing specific hydrologic processes and properties such as infiltration, water content, groundwater flow rate, streamflow, and sediment flux. This paper details the development of a process-level infiltration model and its use in a Probabilistic Performance Assessment (PPA) model for the West Valley Site in western New York. Additionally, this paper demonstrates a method for obtaining distributions of the long-term average hydrologic processes and properties for a PA model. For this PPA, an infiltration model was developed using HELP (Hydrologic Evaluation of Landfill Performance) to evaluate the rates of infiltration, percolation, and drainage through several site soil and proposed decommissioning cover systems, as well as the volumetric water contents of each of the unsaturated layers in a modeled soil column. A particular challenge of probabilistic modeling over long periods of time is scaling parameters appropriately so that outcomes are not grossly over- or underestimated. For these models, distributions of climatic parameters such as average precipitation rate, temperature, wind speed, humidity, and solar radiation were scaled to 1000-year averages and were used to force steady-state models that predict 1000-year average infiltration rates, drainage rates, and water contents. For the HELP model, distributions were developed for the physical and hydrologic properties of each material in a given soil column or cover system, and these properties were allowed to vary stochastically along with climatic inputs. The HELP model was run 5000 times for each of seven different soil and cover systems representing various existing and proposed conditions (i.e. native soil, backfill, existing covers, and a proposed close-in-place cover). The results of the HELP model provide ranges of 1000-year average infiltration rates, drainage rates, and water contents through and within the various layers of each of the different soil and cover conditions. The model outcomes also demonstrate the effectiveness of a proposed close-in-place cover, which essentially eliminates infiltration during its operational lifetime. Following model simulation, a sensitivity analysis was performed to establish the relative importance of each stochastic input with respect to each output. The outcomes of the model simulations and the sensitivity analysis were implemented into a system-level GoldSim model via a novel sensitivity-informed nearest neighbor approach. These infiltration rates, drainage rates, and water contents are used to parameterize mixing cells within the GoldSim Contaminant Transport module, which is the primary system-level model used in the PPA. The infiltration rates are also used as boundary conditions for a stochastic groundwater flow model that is also utilized in the PPA, but not discussed in this paper. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2013

Protein Conformational Populations and Functionally Relevant Substates

Functioning proteins do not remain fixed in a unique structure, rather they sample a range of conformations facilitated by motions within the protein. Even in the native state, a protein exists as a collection of interconverting conformations driven by thermodynamic fluctuations. Motions on the fast time scale allow a protein to sample conformations in the nearby area of its conformational landscape, while motions on slower time scales give it access to conformations in distal areas of the landscape. Emerging evidence indicates that protein landscapes contain conformational substates with dynamic and structural features that support the designated function of the protein. Nuclear magnetic resonance (NMR) experiments provide information about conformational ensembles of proteins. X-ray crystallography allows researchers to identify the most populated states along the landscape, and computational simulations give atom-level information about the conformational substates of different proteins. This ability to characterize and obtain quantitative information about the conformational substates and the populations of proteins within them is allowing researchers to better understand the relationship between protein structure and dynamics and the mechanisms of protein function. In this Account, we discuss recent developments and challenges in the characterization of functionally relevant conformational populations and substates of proteins. In some enzymes, the sampling of functionally relevant conformational substates is connected to promoting the overall mechanism of catalysis. For example, the conformational landscape of the enzyme dihydrofolate reductase hasmultiple substates, which facilitate the binding and the release of the cofactor and substrate and catalyze the hydride transfer. For the enzyme cyclophilin A, computational simulations reveal that the long time scale conformational fluctuations enable the enzyme to access conformational substates that allow it to attain the transition state, therefore promoting the reaction mechanism. Long-term, this emerging view of proteins with conformational substates has broad implications for improving our understanding of enzymes, enzyme engineering, and better drug design. Researchers have already used photoactivation to modulate protein conformations as a strategy to develop a hypercatalytic enzyme. Furthermore, the alteration of the conformational substates through binding of ligands at locations other than the active site provides the basis for the design of new medicines through allosteric modulation.

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

Protein Conformational Populations and Functionally Relevant Sub-states

Functioning proteins do not remain fixed in a unique structure, but instead they sample a range of conformations facilitated by motions within the protein. Even in the native state, a protein exists as a collection of interconverting conformations driven by thermodynamic fluctuations. Motions on the fast time scale allow a protein to sample conformations in the nearby area of its conformational landscape, while motions on slower time scales give it access to conformations in distal areas of the landscape. Emerging evidence indicates that protein landscapes contain conformational substates with dynamic and structural features that support the designated function of the protein. Nuclear magnetic resonance (NMR) experiments provide information about conformational ensembles of proteins. X-ray crystallography allows researchers to identify the most populated states along the landscape, and computational simulations give atom-level information about the conformational substates of different proteins. This ability to characterize and obtain quantitative information about the conformational substates and the populations of proteins within them is allowing researchers to better understand the relationship between protein structure and dynamics and the mechanisms of protein function. In this Account, we discuss recent developments and challenges in the characterization of functionally relevant conformational populations and substates of proteins. In some enzymes, the sampling of functionally relevant conformational substates is connected to promoting the overall mechanism of catalysis. For example, the conformational landscape of the enzyme dihydrofolate reductase has multiple substates, which facilitate the binding and the release of the cofactor and substrate and catalyze the hydride transfer. For the enzyme cyclophilin A, computational simulations reveal that the long time scale conformational fluctuations enable the enzyme to access conformational substates that allow it to attain the transition state, therefore promoting the reaction mechanism. In the long term, this emerging view of proteins with conformational substates has broad implications for improving our understanding of enzymes, enzyme engineering, and better drug design. Researchers have already used photoactivation to modulate protein conformations as a strategy to develop a hypercatalytic enzyme. In addition, the alteration of the conformational substates through binding of ligands at locations other than the active site provides the basis for the design of new medicines through allosteric modulation.

Agarwal, Pratul K↗
DOE OSTI2014

Coarse-grain model for natural cellulose fibrils in explicit water

Understanding biomass structure and dynamics on multiple time and length scales is important for the development of cellulosic biofuels. To this aim, we have developed a coarse-grain (CG) model for molecular dynamics (MD) simulations of cellulose fibrils in explicit water based on target observables from fully atomistic simulations. This model examines the significance of the presence of explicit solvent and compares results with the previous, implicit solvent CG cellulose models. The present, constraint-free CG model is used to generate a series of noncrystalline fibril structures using a coupling parameter, λ, between fully crystalline and fully amorphous potentials. By exploring various structural parameters, including the root-mean-square deviation, root-mean-square fluctuations, radius of gyration, and persistence length, we find the crystalline-to-amorphous state transition takes place at λ ≈ 0.386. The persistence length of cellulose fibril in the transition region corresponds to that of native cellulose fibrils. The transition between crystalline and amorphous fibrils occurs at larger values of λ in explicit water than in the implicit case. Detailed analysis of individual energetic contribution to the transition reveals that the nonbonded interactions, in particular, that of cellulose water interaction, plays a significant role in the observed crystalline to amorphous transition of cellulose fibril. The present study thus highlights the importance of solvent presence that cannot be adequately described with the previous implicit solvent model. In conclusion, the present method provides an accurate and constraint-free approach for deriving a variety of structures of cellulose in water, with a wide range of crystallinity, suitable for incorporation into large-scale models of lignocellulosic biomass.

09 BIOMASS FUELS↗