Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “biological processing”

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 577 records · Page 32

Thermochemical Properties of a Water-Soluble Zr(IV)-Flavonolate Complex for Fluoride Sensing

Detection of fluoride in water samples is critical for monitoring many chemical, environmental, and biological processes, however, these measurements often rely on costly instrumentation. Here, this situation has motivated the design of simple to use, fieldable fluoride sensors. Previously, a fluorometric sensor based on (EDTA)­Zr­(H 2 O) 2 (EDTA = ethylenediamine-N,N,N′,N′-tetraacetate) and 3-hydroxyflavone (FlvH) was reported for fluoride detection in ethanol-H 2 O mixtures. Herein, we report the synthesis and characterization of the previously proposed Zr­(IV)-flavonolate complex as the dimethylammonium salt, [(EDTA)­Zr­(Flv)]­[NMe 2 H 2 ] (Flv = 3-hydroxyflavonolate), and detail its utility for aqueous fluoride sensing. This species exhibits an intense blue fluorescence upon excitation with near-UV light (λ ex = 390 nm; λ em = 460 nm) that is quenched in the presence of NaF (100 mM acetate buffer, pH = 5) with good sensitivity (LOD = 20.9 nM; [F – ] from 10 –8 to 10 –4 M) under neat aqueous conditions. The equilibrium constant (K eq = 1.26 ± 0.00184 × 10 8 M –2 ) and free energy (Δ$G^°_{\textrm{rxn}}$ = –11.0 kcal mol –1 ) were determined under the conditions used for fluoride sensing, indicating a highly favorable reaction with F – in water. Finally, we developed paper-based sensing strips for qualitative fluoride sensing by immobilizing the sensor in a deposited gel enabling good selectivity and reduced interference due to slower ion diffusion.

aqueous chemistry↗

Paired Simulations and Experimental Investigations into the Calcium-Dependent Conformation of Albumin

Serum albumin is the most abundant protein in blood plasma, and it is involved in multiple biological processes. Serum albumin has recently been adapted for improving biomaterial integration with bone tissue, and studies have shown the importance of this protein in bone repair and regeneration. However, the mechanism of action is not yet clear. In stark contrast, other studies have demonstrated that albumin blocks cell adhesion to surfaces, which is seen as a limitation to its bone healing role. These apparent contradictions suggest that the conformation of albumin facilitates its bioactivity, leading to enhanced bone repair. Serum albumin is known to play a major role in maintaining the calcium ion concentration in blood plasma. Due to the prevalence of calcium at bone repair and regeneration sites, it has been hypothesized that calcium binding to serum albumin triggers a conformational change, leading to bioactivity. In the current study, molecular modeling approaches including molecular docking, atomic molecular dynamics (MD) simulation, and coarse-grained MD simulation were used to test this hypothesis by investigating the conformational changes induced in bovine serum albumin by interaction with calcium ions. The computational results were qualitatively validated with experimental Fourier-transform infrared spectroscopy analysis. We find that free calcium ions in solution transiently bind with the three major loops in albumin, triggering a conformational change where N-terminal and C-terminal domains separate from each other in a partial unfolding process. The separation distance between these domains was found to correlate with the calcium ion concentration. The experimental data support the simulation results showing that albumin has enhanced conformational heterogeneity upon exposure to intermediate levels of calcium, without any significant secondary structure changes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Defining the Specificity of Carbohydrate–Protein Interactions by Quantifying Functional Group Contributions

Protein–carbohydrate interactions are significant in a wide range of biological processes, disruption of which has been implicated in many different diseases. The capability of glycan-binding proteins (GBPs) to specifically bind to the corresponding glycans allows GBPs to be utilized in glycan biomarker detection or conversely to serve as targets for therapeutic intervention. However, understanding the structural origins of GBP specificity has proven to be challenging due to their typically low binding affinities (mM) and their potential to display broad or complex specificities. Here we perform molecular dynamics (MD) simulations and post-MD energy analyses with the Poisson–Boltzmann and generalized Born solvent models (MM-PB/GBSA) of the Erythrina cristagalli lectin (ECL) with its known ligands, and with new cocrystal structures reported herein. While each MM-PB/GBSA parametrization resulted in different estimates of the desolvation free energy, general trends emerged that permit us to define GBP binding preferences in terms of ligand substructure specificity. Additionally, we have further decomposed the theoretical interaction energies into contributions made between chemically relevant functional groups. Based on these contributions, the functional groups in each ligand can be assembled into a pharmacophore comprised of groups that are either critical for binding, or enhance binding, or are noninteracting. It is revealed that the pharmacophore for ECL consists of the galactopyranose (Gal) ring atoms along with C6 and the O3 and O4 hydroxyl groups. This approach provides a convenient method for identifying and quantifying the glycan pharmacophore and provides a novel method for interpreting glycan specificity that is independent of residue-level glycan nomenclature. As a result, a pharmacophore approach to defining specificity is readily transferable to molecular design software and, therefore, may be particularly useful in designing therapeutics (glycomimetics) that target GBPs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

All-Atom Biomolecular Simulation in the Exascale Era

Exascale supercomputers have opened the door to dynamic simulations, facilitated by AI/ML techniques, that model biomolecular motions over unprecedented length and time scales. This new capability holds the potential to revolutionize our understanding of fundamental biological processes. Herein we report on some of the major advances that were discussed at a recent CECAM workshop in Pisa, Italy, on the topic with a primary focus on atomic-level simulations. First, we highlight examples of current large-scale biomolecular simulations and the future possibilities enabled by crossing the exascale threshold. Next, we discuss challenges to be overcome in optimizing the usage of these powerful resources. Finally, we close by listing several grand challenge problems that could be investigated with this new computer architecture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Deep Learning-Driven Sampling Technique to Explore the Phase Space of an RNA Stem-Loop

The folding and unfolding of RNA stem-loops are critical biological processes; however, their computational studies are often hampered by the ruggedness of their folding landscape, necessitating long simulation times at the atomistic scale. Here, we adapted DeepDriveMD (DDMD), an advanced deep learning-driven sampling technique originally developed for protein folding, to address the challenges of RNA stem-loop folding. Although tempering- and order parameter-based techniques are commonly used for similar rare-event problems, the computational costs or the need for a priori knowledge about the system often present a challenge in their effective use. DDMD overcomes these challenges by adaptively learning from an ensemble of running MD simulations using generic contact maps as the raw input. DeepDriveMD enables on-the-fly learning of a low-dimensional latent representation and guides the simulation toward the undersampled regions while optimizing the resources to explore the relevant parts of the phase space. We showed that DDMD estimates the free energy landscape of the RNA stem-loop reasonably well at room temperature. Our simulation framework runs at a constant temperature without external biasing potential, hence preserving the information on transition rates, with a computational cost much lower than that of the simulations performed with external biasing potentials. Here, we also introduced a reweighting strategy for obtaining unbiased free energy surfaces and presented a qualitative analysis of the latent space. This analysis showed that the latent space captures the relevant slow degrees of freedom for the RNA folding problem of interest. Finally, throughout the manuscript, we outlined how different parameters are selected and optimized to adapt DDMD for this system. We believe this compendium of decision-making processes will help new users adapt this technique for the rare-event sampling problems of their interest.

Gupta, Ayush↗

COCOMO2: A Coarse-Grained Model for Interacting Folded and Disordered Proteins

Biomolecular interactions are essential in many biological processes, including complex formation and phase separation processes. Coarse-grained computational models are especially valuable for studying such processes via simulation. Here, we present COCOMO2, an updated residue-based coarse-grained model that extends its applicability from intrinsically disordered peptides to folded proteins. This is accomplished with the introduction of a surface exposure scaling factor, which adjusts interaction strengths based on solvent accessibility, to enable the more realistic modeling of interactions involving folded domains without additional computational costs. COCOMO2 was parametrized directly with solubility and phase separation data to improve its performance on predicting concentration-dependent phase separation for a broader range of biomolecular systems compared to the original version. COCOMO2 enables new applications including the study of condensates that involve IDPs together with folded domains and the study of complex assembly processes. COCOMO2 also provides an expanded foundation for the development of multiscale approaches for modeling biomolecular interactions that span from residue-level to atomistic resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correction to “COCOMO2: A Coarse-Grained Model for Interacting Folded and Disordered Proteins”

Biomolecular interactions are essential in many biological processes, including complex formation and phase separation processes. Coarse-grained computational models are especially valuable for studying such processes via simulation. Here, we present COCOMO2, an updated residue-based coarse-grained model that extends its applicability from intrinsically disordered peptides to folded proteins. This is accomplished with the introduction of a surface exposure scaling factor, which adjusts interaction strengths based on solvent accessibility, to enable the more realistic modeling of interactions involving folded domains without additional computational costs. COCOMO2 was parametrized directly with solubility and phase separation data to improve its performance on predicting concentration-dependent phase separation for a broader range of biomolecular systems compared to the original version. COCOMO2 enables new applications including the study of condensates that involve IDPs together with folded domains and the study of complex assembly processes. COCOMO2 also provides an expanded foundation for the development of multiscale approaches for modeling biomolecular interactions that span from residue-level to atomistic resolution.

Molecular interactions↗

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗

Marcus Theory of Ion-Pairing

We present a theory for ion pair dissociation and association, motivated by the concepts of the Marcus theory of electron transfer. Despite the extensive research on ion-pairing in many chemical and biological processes, much can be learned from the exploration of collective reaction coordinates. To this end, we explore two reaction coordinates, ion pair distance and coordination number. The study of the correlation between these reaction coordinates provides a new insight into the mechanism and kinetics of ion pair dissociation and association in water. The potential of mean force on these 2D-surfaces computed from molecular dynamics simulations of different monovalent ion pairs reveal a Marcus-like mechanism for ion-pairing: Water molecules rearrange forming an activated coordination state prior to ion pair dissociation or association, followed by relaxation of the coordination state due to further water rearrangement. Like Marcus theory, we find the existence of an inverted region where the transition rates are slower with increasing exergonicity. This study provides a new perspective for the future investigations of ion-pairing and transport. SR, CJM, and GKS were supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. MDB was supported by MS3 (Materials Synthesis and Simulation Across Scales) Initiative, a Laboratory Directed Research and Development Program at Pacific Northwest National Laboratory (PNNL). The research was performed using PNNL Institutional Computing. PNNL is a multi-program national laboratory operated by Battelle for the U.S. Department of Energy.

Roy, Santanu↗

Electrostatic Control of Polymorphism in Charged Amphiphile Assemblies

Stimuli-induced structural transformations of molecular assemblies in aqueous solutions are integral to nanotechnological applications and biological processes. In particular, pH responsive amphiphiles as well as proteins with various degrees of ionization can reconfigure in response to pH variations. Here, we use in situ small and wide-angle X-ray scattering (SAXS/WAXS), transmission electron microscopy (TEM), and Monte Carlo simulations to show how charge regulation via pH induces morphological changes in the assembly of a positively charged peptide amphiphile (PA). Monte Carlo simulations and pH titration measurements reveal that ionic correlations in the PA assemblies shift the ionizable amine pK ~ 8 from pK ~ 10 in the lysine headgroup. SAXS and TEM show that with increasing pH, the assembly undergoes spherical micelle to cylindrical nanofiber to planar bilayer transitions. SAXS/WAXS reveal that the bilayer leaflets are interdigitated with the tilted PA lipid tails crystallized on a rectangular lattice. The details of the molecular packing in the membrane result from interplay between steric and van der Waals interactions. We speculate that this packing motif is a general feature of bilayers comprised of amphiphilic lipids with large ionic headgroups. Altogether, our studies correlate the molecular charge and the morphology for a pH-responsive PA system and provide insights into the Å-scale molecular packing in such assemblies.

36 MATERIALS SCIENCE↗

Reaction Rate Theory in Coordination Number Space: An Application to Ion Solvation

Understanding reaction mechanisms in many chemical and biological processes require application of rare event theories. In these theories, an effective choice of a reaction coordinate to describe a reaction pathway is essential. To this end, we study ion solvation in water using molecular dynamics simulations and explore the utility of coordination number (n = number of water molecules in the first solvation shell) as the reaction coordinate. Here we compute the potential of mean force (W(n)) using umbrella sampling, predicting multiple metastable n-states for both cations and anions. We find with increasing ionic size, these states become more stable and structured for cations when compared to anions. We have extended transition state theory (TST) to calculate transition rates between n-states. TST overestimates the rate constant due to solvent-induced barrier recrossings that are not accounted for. We correct the TST rates by calculating transmission coefficients using the reactive flux method. This approach enables a new way of understanding rare events involving coordination complexes. We gratefully acknowledge Liem Dang and Panos Stinis for useful discussion. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. SR, CJM, and GKS were supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. MDB was supported by MS3 (Materials Synthesis and Simulation Across Scales) Initiative, a Laboratory Directed Research and Development Program at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy.

Roy, Santanu↗

Monovalent and Divalent Cations at the a-Al2O3(0001)/Water Interface: How Cation Identity Affects Interfacial Ordering and Vibrational Dynamics

Mineral oxide-water interfaces are important for a wide range of industrial, geochemical, and biological processes. The reactivity of these interfaces is strongly impacted by the presence of ions. Thus, it is critical to understand how ions alter the interfacial environment. This can be achieved by measuring the changes in the structure and vibrational dynamics of interfacial water induced by the presence of ions in close vicinity to the mineral surface. The a-Al2O3(0001) surface represents a flexible platform to study the effect of ions on interfacial aqueous environments at positive, neutral and negative surface charge. Using vibrational sum frequency generation (vSFG) in the frequency and time domain, we investigate how monovalent and divalent cations affect the hydrogen bonding environment of the first few layers of interfacial water next to a-Al2O3(0001). Our results indicate that monovalent cations, such as Li+, Na+, K+, and Cs+, appear to have lower binding affinities for the interface compared to Ca2+, Sr2+, and Ba2+. This leads to an interfacial region that is structured in a cation valence dependent manner. Addition of divalent cations at the negatively charged interface (pH 10) increases the spectral intensity in the 3400 cm-1 region compared to neat pH 10 H2O, in contrast to monovalent cations that only attenuate the vSFG signal. Time resolved vSFG measurements reveal that the O-H vibrational lifetime (T1) of interfacial species at pH 10 in the presence of NaCl and BaCl2 remains similar, but restructuring of the surface seen in steady state vSFG is manifested in the degree to which strongly hydrogen bonded species recover to their original populations post excitation. By tracking the accumulation of ions at the interface via the vSFG response, we can characterize the unique surface arrangements of interfacial water molecules induced by a range of monovalent and divalent cations at the a-Al2O3(0001)/water interface.

Piontek, Stefan↗

Observation of Double Excitations in the Resonant Inelastic X-ray Scattering of Nitric Oxide

The nitrogen K-edge resonant inelastic X-ray scattering (RIXS) map of nitric oxide (NO) has been measured and simulated to provide a detailed analysis of the observed features. High-resolution experimental RIXS maps were collected using an in situ gas flow cell and a high-transmission soft X-ray spectrometer. Accurate descriptions of the ground, excited, and core-excited states are based upon restricted active space self-consistent-field calculations using second order multiconfigurational perturbation theory. The nitrogen K-edge RIXS map of NO shows a range of features that can be assigned to intermediate states arising from 1s → π* and 1s → Rydberg excitations; additional bands are attributed to doubly excited intermediate states comprising 1s → π* and π → π* excitations. These results provide a detailed picture of RIXS for an open-shell molecule and an extensive description of the core-excited electronic structure of NO, an important molecule in many chemical and biological processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair P D1 –P D2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state P D1 – –P D2 + stabilizes energetically below the excitonic state. Importantly, this underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic Signature of Metal-hydroxo and Peroxo Species in K-edge X-ray Absorption Spectra

Metal-dioxygen species are important intermediates formed during dioxygen activations by metalloenzymes in various biological processes, by catalysts in fuel cells, and prior to O 2 evolution by photosystem II. In this work, we focus on manganese-porphyrin complexes using tetramesitylporphyrin ligand (TMP) to explore changes in Mn K-edge X-ray absorption spectroscopy (XAS) associated with the formation of Mn-hydroxide and Mn–O 2 peroxide species. With limited spectroscopic characterization of these compounds, Mn K β X-ray emission spectroscopy (XES), XAS, density functional theory (DFT), and time-dependent DFT (TD-DFT) analysis will enhance our understanding of their complex electronic structure. We show that the shape of the pre-edge in the K-edge Mn X-ray absorption near-edge structure (XANES) can serve as a spectroscopic signature of the Mn III -peroxo formation and thus can be used to track the presence of the side-on peroxide as an intermediate in time-resolved or in situ experiments. In conclusion, our results will help to further summarize the spectroscopic fingerprints for peroxo and hydroxo species, addressing the challenge of identifying the reactive metal species in catalytic reactions.

Extended X-ray absorption fine structure↗

Direct Molecular Evidence of Proton Transfer and Mass Dynamics at the Electrode–Electrolyte Interface

Proton transfer has been widely regarded as a key step in many electrochemical and biological processes. However, direct molecular evidence has long been lacking. In this work, we chose the electrochemical oxidation of acetaminophen (APAP) as a model system and utilized in situ liquid time-of-flight secondary ion mass spectroscopy (ToF-SIMS) to molecularly examine proton solvation and transfer in this process. In addition, we successfully captured and identified the transient radical intermediate, providing solid molecular evidence to resolve an important debate in electron transfer-proton transfer oxidation mechanism of APAP. Moreover, the potential-dependent behaviors of both inert ions and electroactive species during the dynamic potential scanning were chemically monitored in real time and the mass diffusion mechanism regarding the electroactive and non-electroactive species was revealed under polarized conditions. The results are well consistent with our computer simulation. The observations in this work greatly improved the comprehensive understanding of proton transfer and mass dynamics occurring at the electrode-electrolyte interface involved in complex electrochemical processes.

Wang, Jun-Gang↗

Stimulated X-ray Raman and Absorption Spectroscopy of Iron–Sulfur Dimers

Iron–sulfur complexes play an important role in biological processes such as metabolic electron transport. A detailed understanding of the mechanism of long-range electron transfer requires knowledge of the electronic structure of the complexes, which has traditionally been challenging to obtain, either by theory or by experiment, but the situation has begun to change with advances in quantum chemical methods and intense free electron laser light sources. We compute the spectra for stimulated X-ray Raman spectroscopy (SXRS) and absorption spectroscopy of homovalent and mixed-valence [2Fe–2S] complexes, using the ab initio density matrix renormalization group algorithm. The simulated spectra show clear signatures of the theoretically predicted dense low-lying excited states within the d–d manifold. Finally, the difference in spectral intensity between the absorption-active and Raman-active states provides a potential mechanism to selectively excite states by a proper tuning of the excitation pump, to access the electronic dynamics within this manifold.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗