Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Biomolecular dynamics”

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

Performance assessment of ensembles of in situ workflows under resource constraints

Summary Scientific breakthroughs in biomolecular methods and improvements in hardware technology have shifted from a long‐running simulation to a large set of shorter simulations running simultaneously, called an ensemble. In an ensemble, simulations are usually coupled with analyses of data produced by the simulations. In situ methods can be used to analyze large volumes of data generated by scientific simulations at runtime (i.e., simulations and analyses are performed concurrently). In this work, we study the execution of ensemble‐based simulations paired with in situ analyses using in‐memory staging methods. Using an ensemble of molecular dynamics in situ workflows with multiple simulations and analyses, we first show that collecting traditional metrics such as makespan, instructions per cycle, memory usage, or cache miss ratio is not sufficient to characterize complex behaviors of ensembles. We propose a method to evaluate the performance of ensembles of workflows that captures multiple resource usage aspects: resource efficiency, resource allocation, and resource provisioning. Experimental results demonstrate that the proposed method can effectively distinguish the performance of different component placements in an ensemble with up to 32 ensemble members. By evaluating different co‐location scenarios, our proposed performance indicators demonstrate benefits of co‐locating simulation and coupled analyses within a compute node.

Do, Tu Mai Anh↗

Data-driven prediction of α IIb β 3 integrin activation paths using manifold learning and deep generative modeling

The integrin heterodimer is a transmembrane protein critical for driving cellular process and is a therapeutic target in the treatment of multiple diseases linked to its malfunction. Activation of integrin involves conformational transitions between bent and extended states. Some of the conformations that are intermediate between bent and extended states of the heterodimer have been experimentally characterized, but the full activation pathways remain unresolved both experimentally due to their transient nature and computationally due to the challenges in simulating rare barrier crossing events in these large molecular systems. An understanding of the activation pathways can provide new fundamental understanding of the biophysical processes associated with the dynamic interconversions between bent and extended states and unveil new putative therapeutic targets. In this work, we apply nonlinear manifold learning to coarse-grained molecular dynamics simulations of bent, extended, and two intermediate states of αI I b β3 integrin to learn a low-dimensional embedding of the configurational phase space. We then train deep generative models to learn an inverse mapping between the low-dimensional embedding and high-dimensional molecular space and use these models to interpolate the molecular configurations constituting the activation pathways between the experimentally characterized states. Furthermore, this work furnishes plausible predictions of integrin activation pathways and reports a generic and transferable multi-scale technique to predict transition pathways for biomolecular systems.

97 MATHEMATICS AND COMPUTING↗

Inexact iterative numerical linear algebra for neural network-based spectral estimation and rare-event prediction

Understanding dynamics in complex systems is challenging because there are many degrees of freedom, and those that are most important for describing events of interest are often not obvious. The leading eigenfunctions of the transition operator are useful for visualization, and they can provide an efficient basis for computing statistics, such as the likelihood and average time of events (predictions). Here, we develop inexact iterative linear algebra methods for computing these eigenfunctions (spectral estimation) and making predictions from a dataset of short trajectories sampled at finite intervals. We demonstrate the methods on a low-dimensional model that facilitates visualization and a high-dimensional model of a biomolecular system. Implications for the prediction problem in reinforcement learning are discussed.

Chemistry↗

Leakless end-to-end transport of small molecules through micron-length DNA nanochannels

Designed and engineered protein and DNA nanopores can be used to sense and characterize single molecules and control transmembrane transport of molecular species. However, designed biomolecular pores are less than 100 nm in length and are used primarily for transport across lipid membranes. Nanochannels that span longer distances could be used as conduits for molecules between nonadjacent compartments or cells. We design micrometer-long, 7-nm-diameter DNA nanochannels that small molecules can traverse according to the laws of continuum diffusion. Binding DNA origami caps to channel ends eliminates transport and demonstrates that molecules diffuse from one channel end to the other rather than permeating through channel walls. These micrometer-length nanochannels can also grow, form interconnects, and interface with living cells. This work thus shows how to construct multifunctional, dynamic agents that control molecular transport, opening ways of studying intercellular signaling and modulating molecular transport between synthetic and living cells.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular dynamics simulations of the secondary-binding site in disaccharide-modified glycopeptide antibiotics

Oritavancin is a semisynthetic glycopeptide antibiotic used to treat severe infections by multidrug-resistant Gram-positive pathogens. Oritavancin is known to be a thousand times more potent than vancomycin against Gram-positive bacteria due to the additional interactions with bacterial peptidoglycan (PG) facilitated by a secondary-binding site. The presence of this secondary-binding site is evident in desleucyl-oritavancin, an Edman degradation product of oritavancin, still retaining its potency against Gram-positive bacteria, whereas desleucyl-vancomycin is devoid of any antimicrobial activities. Herein, using explicit solvent molecular dynamics (MD) simulations, steered MD simulations, and umbrella sampling, we show evidence of a secondary-binding site mediated by the disaccharide-modified hydrophobic sidechain of oritavancin interactions with the pentaglycyl-bridge segment of the PG. The interactions were characterized through comparison to the interaction of PG with chloroeremomycin, vancomycin, and the desleucyl analogs of the glycopeptides. Our results show that the enhanced binding of oritavancin to PG over the binding of the other complexes studied is due to an increase in the hydrophobic effect, electrostatic and van der Waals interactions, and not the average number of hydrogen bonds. Our ranking of the binding interactions of the biomolecular complexes directly correlates with the order based on their experimental minimum inhibitory concentrations. The results of our simulations provide insight into the modification of glycopeptides to increase their antimicrobial activities or the design of novel antibiotics against pathogenic Gram-positive bacteria.

60 APPLIED LIFE SCIENCES↗

3D diffractive imaging of nanoparticle ensembles using an x-ray laser

Single particle imaging at x-ray free electron lasers (XFELs) has the potential to determine the structure and dynamics of single biomolecules at room temperature. Two major hurdles have prevented this potential from being reached, namely, the collection of sufficient high-quality diffraction patterns and robust computational purification to overcome structural heterogeneity. We report the breaking of both of these barriers using gold nanoparticle test samples, recording around 10 million diffraction patterns at the European XFEL and structurally and orientationally sorting the patterns to obtain better than 3-nm-resolution 3D reconstructions for each of four samples. With these new developments, integrating advancements in x-ray sources, fast-framing detectors, efficient sample delivery, and data analysis algorithms, we illuminate the path towards sub-nanometer biomolecular imaging. The methods developed here can also be extended to characterize ensembles that are inherently diverse to obtain their full structural landscape.

47 OTHER INSTRUMENTATION↗

Sub‐5 Ångstrom Porosity Tuning in Calixarene‐Derived Porous Liquids via Supramolecular Complexation Construction

Abstract Sub‐Ångstrom‐level porosity engineering, which is appealing in gas separations, has been demonstrated in solid carbon, polymer, and framework materials but rarely achieved in the liquid phase. In this work, a gas molecular sieving effect in the liquid phase at sub‐5 Ångstrom scale is created via sophisticated porosity tuning in calixarene‐derived porous liquids (PLs). Type II PLs are constructed via supramolecular complexation between the sodium salts of calixarene derivatives and crown ether solvents. The chemical structure variation and assembly behavior of the porous host upon PL construction are monitored by spectroscopy‐, X‐ray‐, and neutron‐scattering techniques. The presence of permanent porosity in calixarene‐derived PLs is verified by pressure swing gas uptake, altered CO 2 physisorption behavior, and molecular simulations. Sub‐5 Ångstrom porosity tuning within the PL phase is achieved by introducing bulky substituted groups on the benzene ring of the calixarene host, which then greatly affects the dynamic motion and transport behavior of CO 2 molecules and the Xe uptake performance. The approach being demonstrated in this work represents a promising pathway to tune and leverage the porosity effect for enhanced gas uptake capacity and selectivity in liquid sorbents.

Li, Errui [Department of Chemistry University of T↗

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The development and applications of multidimensional biomolecular spectroscopy illustrated by photosynthetic light harvesting

The parallel and synergistic developments of atomic resolution structural information, new spectroscopic methods, their underpinning formalism, and the application of sophisticated theoretical methods have led to a step function change in our understanding of photosynthetic light harvesting, the process by which photosynthetic organisms collect solar energy and supply it to their reaction centers to initiate the chemistry of photosynthesis. The new spectroscopic methods, in particular multidimensional spectroscopies, have enabled a transition from recording rates of processes to focusing on mechanism. We discuss two ultrafast spectroscopies – two-dimensional electronic spectroscopy and two-dimensional electronic-vibrational spectroscopy – and illustrate their development through the lens of photosynthetic light harvesting. Both spectroscopies provide enhanced spectral resolution and, in different ways, reveal pathways of energy flow and coherent oscillations which relate to the quantum mechanical mixing of, for example, electronic excitations (excitons) and nuclear motions. The new types of information present in these spectra provoked the application of sophisticated quantum dynamical theories to describe the temporal evolution of the spectra and provide new questions for experimental investigation. While multidimensional spectroscopies have applications in many other areas of science, we feel that the investigation of photosynthetic light harvesting has had the largest influence on the development of spectroscopic and theoretical methods for the study of quantum dynamics in biology, hence the focus of this review. We conclude with key questions for the next decade of this review.

59 BASIC BIOLOGICAL SCIENCES↗

Current limitations of solid-state NMR in carbohydrate and cell wall research

High-resolution investigation of cell wall materials has emerged as an important application of biomolecular solid-state NMR (ssNMR). Multidimensional correlation experiments have become a standard method for obtaining sufficient spectral resolution to determine the polymorphic structure of carbohydrates and address biochemical questions regarding the supramolecular organization of cell walls. Here, using plant cellulose and matrix polysaccharides as examples, we will review how the multifaceted complexity of polysaccharide structure is impeding the resonance assignment process and assess the available biochemical and spectroscopic approaches that could circumvent this barrier. We will emphasize the ineffectiveness of the current methods in reconciling the ever-growing dataset and deriving structural information. We will evaluate the protocols for achieving efficient and homogeneous hyperpolarization across the cell wall material using magic-angle spinning dynamic nuclear polarization (MAS-DNP). Critical questions regarding the line-broadening effects of cell wall molecules at cryogenic temperature and by paramagnetic biradicals will be considered. Finally, the MAS-DNP method will be placed into a broader context with other structural characterization techniques, such as cryo-electron microscopy, to advance ssNMR research in carbohydrate and cell wall biomaterials.

09 BIOMASS FUELS↗

Deuteron off-resonance rotating frame relaxation for the characterization of slow motions in rotating and static solid-state proteins

Here, we demonstrate the feasibility of deuterium solid-state NMR off-resonance rotating frame relaxation measurements for studies of slow motions in biomolecular solids. The pulse sequence, which includes adiabatic pulses for magnetization alignment, is illustrated for static and magic-angle spinning conditions away from rotary resonances. We apply the measurements for three systems with selective deuterium labels at methyl groups: a) a model compound, Fluorenylmethyloxycarbonyl methionine-D 3 amino acid, for which the principles of the measurements and corresponding motional modeling based on rotameric interconversions are demonstrated; b) amyloid-β 1-40 fibrils labeled at a single alanine methyl group located in the disordered N-terminal domain. This system has been extensively studied in prior work and here serves as a test of the method for complex biological systems. The essential features of the dynamics consist of large-scale rearrangements of the disordered N-terminal domain and the conformational exchange between the free and bound forms of the domain, the latter one due to transient interactions with the structured core of the fibrils. and c) a 15-residue helical peptide which belongs to the predicted α-helical domain near the N-terminus of apolipoprotein B. The peptide is solvated with triolein and incorporates a selectively labeled leucine methyl groups. The method permits model refinement, indicating rotameric interconversions with a distribution of rate constants.

59 BASIC BIOLOGICAL SCIENCES↗

Performance efficient macromolecular mechanics via sub-nanometer shape based coarse graining

Dimensionality reduction via coarse grain modeling is a valuable tool in biomolecular research. For large assemblies, ultra coarse models are often knowledge-based, relying on a priori information to parameterize models thus hindering general predictive capability. Here, we present substantial advances to the shape based coarse graining (SBCG) method, which we refer to as SBCG2. SBCG2 utilizes a revitalized formulation of the topology representing network which makes high-granularity modeling possible, preserving atomistic details that maintain assembly characteristics. Further, we present a method of granularity selection based on charge density Fourier Shell Correlation and have additionally developed a refinement method to optimize, adjust and validate high-granularity models. We demonstrate our approach with the conical HIV-1 capsid and heteromultimeric cofilin-2 bound actin filaments. Our approach is available in the Visual Molecular Dynamics (VMD) software suite, and employs a CHARMM-compatible Hamiltonian that enables high-performance simulation in the GPU-resident NAMD3 molecular dynamics engine.

59 BASIC BIOLOGICAL SCIENCES↗

CryoDRGN-AI: neural ab initio reconstruction of challenging cryo-EM and cryo-ET datasets

Proteins and other biomolecules form dynamic macromolecular machines that are tightly orchestrated to move, bind, and perform chemistry. Cryo-electron microscopy (cryo-EM) and cryo-electron tomography (cryo-ET) can access the intrinsic heterogeneity of these complexes and are therefore key tools for understanding their function. However, 3D reconstruction of the collected imaging data presents a challenging computational problem, especially without any starting information, a setting termed ab initio reconstruction. Here, in this study, we introduce cryoDRGN-AI, a method leveraging an expressive neural representation and combining an exhaustive search strategy with gradient-based optimization to process challenging heterogeneous datasets. Using cryoDRGN-AI, we reveal new conformational states in large datasets, reconstruct previously unresolved motions from unfiltered datasets, and demonstrate ab initio reconstruction of biomolecular complexes from in situ data. With this expressive and scalable model for structure determination, we hope to unlock the full potential of cryo-EM and cryo-ET as a high-throughput tool for structural biology and discovery.

Levy, Axel [Stanford Univ., CA (United States); SL↗

Out of the Dark and Into the Light: A New View of Phytochrome Photobodies

Light is a critical environmental stimulus for plants, serving as an energy source via photosynthesis and a signal for developmental programming. Plants perceive light through various light-responsive proteins, termed photoreceptors. Phytochromes are red-light photoreceptors that are highly conserved across kingdoms. In the model plant Arabidopsis thaliana , phytochrome B serves as a light and thermal sensor, mediating physiological processes such as seedling germination and establishment, hypocotyl growth, chlorophyll biogenesis, and flowering. In response to red light, phytochromes convert to a biologically active form, translocating from the cytoplasm into the nucleus and further compartmentalizes into subnuclear compartments termed photobodies. PhyB photobodies regulate phytochrome-mediated signaling and physiological outputs. However, photobody function, composition, and biogenesis remain undefined since their discovery. Based on photobody cellular dynamics and the properties of internal components, photobodies have been suggested to undergo liquid-liquid phase separation, a process by which some membraneless compartments form. Here, we explore photobodies as environmental sensors, examine the role of their protein constituents, and outline the biophysical perspective that photobodies may be undergoing liquid-liquid phase separation. Understanding the molecular, cellular, and biophysical processes that shape how plants perceive light will help in engineering improved sunlight capture and fitness of important crops.

54 ENVIRONMENTAL SCIENCES↗

Spray Mechanism of Contained-Electrospray Ionization

Analytical characteristics of contained electrospray ionization (ESI) are summarized in terms of its potential to modify the analyte solution during the stages of droplet formation to provide opportunities to generate native versus denatured biomolecular gas-phase ions, without the need for bulk-phase analyte modifications. The real-time modification of the charged microdroplets occurs in a cavity that is included in the outlet of the contained-ESI ion source. Close examination of the inside of the cavity using a high-speed camera revealed the formation of discrete droplets as well as thin liquid films in the droplets wake. When operated at 20 psi N2 pressure, the droplets were observed to move at an average speed of 8 mm/s providing ~1 s mixing time in a 10 mm cavity length. Evidence is provided for the presence of highly reactive charged droplets based on myoglobin charge state distribution, apo-myoglobin contents, and ion mobility drift time profiles under different spray conditions. Mechanistic insights for the capture of vapor-phase reagents and droplet dynamics as influenced by different operational modes are also described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear Magnetic Resonance Dynamics of LiTFSI–Pyrazole Eutectic Solvents

Deep Eutectic Solvents (DESs) have emerged as promising candidates to replace conventional organic solvents in various technological applications due to their low vapor pressure, non-flammability, and ease of preparation at low costs. In particular, Type IV DESs, which are composed of metal salts and hydrogen bond donors, are possible replacements for lithium-ion battery electrolytes. In this study, we investigate the molecular dynamics of solvents of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and pyrazole (PYR) at varying LiTFSI:PYR molar ratios (1:2, 1:3, 1:4, 1:5) using Nuclear Magnetic Resonance Dispersion (NMRD) and Pulsed Field Gradient (PFG) Nuclear Magnetic Resonance (NMR). PFG NMR reveals composition-dependent diffusion trends, while NMRD provides molecular-level insights into the longitudinal relaxation rate (R 1 = 1/T 1 ). Notably, the LiTFSI:PYR (1:2) sample shows distinct behavior across both techniques, exhibiting enhanced relaxation rates and lower self-diffusion for 1 H compared to the other nuclei ( 19 F and 7 Li), suggestive of stronger and more efficient Li + –pyrazole interactions, as confirmed by the modeling of the relaxation profiles. Our study advances understanding of ion dynamics in azole-based eutectic solvents, supporting their potential use in safer battery electrolytes.

FFC-NMR↗

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

Alchemical free energy simulations have long been utilized to predict free energy changes for binding affinity and solubility of small molecules. However, while the theoretical foundation of these methods is well established, seamlessly handling many of the practical aspects regarding the preparation of the different thermodynamic end states of complex molecular systems and the numerous processing scripts often remains a burden for successful applications. Here, we present CHARMM-GUI Free Energy Calculator (http://www.charmm-gui.org/input/fec) that provides various alchemical free energy perturbation molecular dynamics (FEP/MD) systems with input and post-processing scripts for NAMD and GENESIS. Four submodules are available: Absolute Ligand Binder (for absolute ligand binding FEP/MD), Relative Ligand Binder (for relative ligand binding FEP/MD), Absolute Ligand Solvator (for absolute ligand solvation FEP/MD), and Relative Ligand Solvator (for relative ligand solvation FEP/MD). Each module is designed to build multiple systems of a set of selected ligands at once for high-throughput FEP/MD simulations. The capability of Free Energy Calculator is illustrated by absolute and relative solvation FEP/MD of a set of ligands and absolute and relative binding FEP/MD of a set of ligands for T4-lysozyme in solution and the adenosine A 2A receptor in a membrane. The calculated free energy values are overall consistent with the experimental and published free energy results (within ~1 kcal/mol). We hope that Free Energy Calculator is useful to carry out high-throughput FEP/MD simulations in the field of biomolecular sciences and drug discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mol* Viewer: modern web app for 3D visualization and analysis of large biomolecular structures

Abstract Large biomolecular structures are being determined experimentally on a daily basis using established techniques such as crystallography and electron microscopy. In addition, emerging integrative or hybrid methods (I/HM) are producing structural models of huge macromolecular machines and assemblies, sometimes containing 100s of millions of non-hydrogen atoms. The performance requirements for visualization and analysis tools delivering these data are increasing rapidly. Significant progress in developing online, web-native three-dimensional (3D) visualization tools was previously accomplished with the introduction of the LiteMol suite and NGL Viewers. Thereafter, Mol* development was jointly initiated by PDBe and RCSB PDB to combine and build on the strengths of LiteMol (developed by PDBe) and NGL (developed by RCSB PDB). The web-native Mol* Viewer enables 3D visualization and streaming of macromolecular coordinate and experimental data, together with capabilities for displaying structure quality, functional, or biological context annotations. High-performance graphics and data management allows users to simultaneously visualise up to hundreds of (superimposed) protein structures, stream molecular dynamics simulation trajectories, render cell-level models, or display huge I/HM structures. It is the primary 3D structure viewer used by PDBe and RCSB PDB. It can be easily integrated into third-party services. Mol* Viewer is open source and freely available at https://molstar.org/.

Sehnal, David↗