SEARCH · Engineering Papers
Results for “Molecular conformation”
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.
MOF-enabled confinement and related effects for chemical catalyst presentation and utilization
A defining characteristic of nearly all catalytically functional MOFs is uniform, molecular-scale porosity. MOF pores, linkers and nodes that define them, help regulate reactant and product transport, catalyst siting, catalyst accessibility, catalyst stability, catalyst activity, co-catalyst proximity, composition of the chemical environment at and beyond the catalytic active site, chemical intermediate and transition-state conformations, thermodynamic affinity of molecular guests for MOF interior sites, framework charge and density of charge-compensating ions, pore hydrophobicity/hydrophilicity, pore and channel rigidity vs. flexibility, and other features and properties. Collectively and individually, these properties help define overall catalyst functional behaviour. Here, this review focuses on how porous, catalyst-containing MOFs capitalize on molecular-scale confinement, containment, isolation, environment modulation, energy delivery, and mobility to accomplish desired chemical transformations with potentially superior selectivity or other efficacy, especially in comparison to catalysts in homogeneous solution environments.
Conformational Isomerization of Imide Anions Governs Solvation and Transport in Water-in-Salt Electrolytes
The behavior of highly concentrated electrolytes departs radically from the dilute-solution theory, yet the molecular origin of this transformation remains unresolved. Here, we identify the conformational isomerization of molecular ions as a decisive, previously unrecognized control parameter governing structure and transport in crowded aqueous electrolytes. Across a series of fluorosulfonimide anions, we show that increasing concentration drives a collective shift from extended transoid to compact cisoid conformers, revealed by small-angle X-ray scattering, vibrational spectroscopy, pulsed-field gradient NMR, and molecular dynamics simulations. This conformational transition triggers a collapse of the hydrogen-bonded water network and the emergence of densely packed ionic domains with confined water, producing a qualitative change in Li+ transport from solvent-mediated diffusion to network-confined hopping. Anion size and asymmetry systematically tune the onset of this transition, demonstrating that molecular geometry dictates mesoscale organization and dynamics in the ion-rich regime. Our results establish ion conformation, not merely composition or coordination, as a fundamental thermodynamic variable in concentrated solutions, providing a chemical framework that unifies solvation structure and transport in water-in-salt electrolytes and suggesting new principles for designing dense ionic media.
Designer molecules of the synaptic organizer MDGA1 reveal 3D conformational control of biological function
MDGAs (MAM domain-containing glycosylphosphatidylinositol anchors) are synaptic cell surface molecules that regulate the formation of trans-synaptic bridges between neurexins (NRXNs) and neuroligins (NLGNs), which promote synaptic development. Mutations in MDGAs are implicated in various neuropsychiatric diseases. MDGAs bind NLGNs in cis on the postsynaptic membrane and physically block NLGNs from binding to NRXNs. In crystal structures, the six immunoglobulin (Ig) and single fibronectin III domains of MDGA1 reveal a striking compact, triangular shape, both alone and in complex with NLGNs. Whether this unusual domain arrangement is required for biological function or other arrangements occur with different functional outcomes is unknown. Here, we show that WT MDGA1 can adopt both compact and extended 3D conformations that bind NLGN2. Designer mutants targeting strategic molecular elbows in MDGA1 alter the distribution of 3D conformations while leaving the binding affinity between soluble ectodomains of MDGA1 and NLGN2 intact. In contrast, in a cellular context, these mutants result in unique combinations of functional consequences, including altered binding to NLGN2, decreased capacity to conceal NLGN2 from NRXN1β, and/or suppressed NLGN2-mediated inhibitory presynaptic differentiation, despite the mutations being located far from the MDGA1-NLGN2 interaction site. Thus, the 3D conformation of the entire MDGA1 ectodomain appears critical for its function, and its NLGN-binding site on Ig1-Ig2 is not independent of the rest of the molecule. As a result, global 3D conformational changes to the MDGA1 ectodomain via strategic elbows may form a molecular mechanism to regulate MDGA1 action within the synaptic cleft.
Conformational heterogeneity of the BTK PHTH domain drives multiple regulatory states
Full-length Bruton’s tyrosine kinase (BTK) has been refractory to structural analysis. The nearest full-length structure of BTK to date consists of the autoinhibited SH3–SH2–kinase core. Precisely how the BTK N-terminal domains (the Pleckstrin homology/Tec homology [PHTH] domain and proline-rich regions [PRR] contain linker) contribute to BTK regulation remains unclear. We have produced crystals of full-length BTK for the first time but despite efforts to stabilize the autoinhibited state, the diffraction data still reveal only the SH3–SH2–kinase core with no electron density visible for the PHTH–PRR segment. Cryo-electron microscopy (cryoEM) data of full-length BTK, on the other hand, provide the first view of the PHTH domain within full-length BTK. CryoEM reconstructions support conformational heterogeneity in the PHTH–PRR region wherein the globular PHTH domain adopts a range of states arrayed around the autoinhibited SH3–SH2–kinase core. On the way to activation, disassembly of the SH3–SH2–kinase core opens a new autoinhibitory site on the kinase domain for PHTH domain binding that is ultimately released upon interaction of PHTH with phosphatidylinositol (3,4,5)-trisphosphate. Membrane-induced dimerization activates BTK and we present here a crystal structure of an activation loop swapped BTK kinase domain dimer that likely represents the conformational state leading to trans-autophosphorylation. Together, these data provide the first structural elucidation of full-length BTK and allow a deeper understanding of allosteric control over the BTK kinase domain during distinct stages of activation.
Conformational diversity of 1-phenylpiperidin-4-one in the gas phase
Conformational preferences of substituted, saturated cyclic ketones may deviate from those of their parent analogues. 1-Phenylpiperidin-4-one was found to exist in several forms differing by not only position (axial or equatorial) of the aryl-substituent but also by a configuration of the piperidine cycle itself: chair-Eq, chair-Ax and twist-Tw. Here, the ratios of the co-existing species contributions in gas phase predicted by various quantum chemical calculations were confirmed by synchronous gas-phase electron diffraction and mass spectrometric experiments carried out for the vapor at 337 K demonstrating a best fit with the B3LYP-D3/cc-pVTZ combination: Eq:Ax:Tw=55(13):22(9):23(10) vs. 40:35:25.
Robust and scalable uncertainty estimation with conformal prediction for machine-learned interatomic potentials
Abstract Uncertainty quantification (UQ) is important to machine learning (ML) force fields to assess the level of confidence during prediction, as ML models are not inherently physical and can therefore yield catastrophically incorrect predictions. Established a-posteriori UQ methods, including ensemble methods, the dropout method, the delta method, and various heuristic distance metrics, have limitations such as being computationally challenging for large models due to model re-training. In addition, the uncertainty estimates are often not rigorously calibrated. In this work, we propose combining the distribution-free UQ method, known as conformal prediction (CP), with the distances in the neural network’s latent space to estimate the uncertainty of energies predicted by neural network force fields. We evaluate this method (CP+latent) along with other UQ methods on two essential aspects, calibration, and sharpness, and find this method to be both calibrated and sharp under the assumption of independent and identically-distributed (i.i.d.) data. We show that the method is relatively insensitive to hyperparameters selected, and test the limitations of the method when the i.i.d. assumption is violated. Finally, we demonstrate that this method can be readily applied to trained neural network force fields with traditional and graph neural network architectures to obtain estimates of uncertainty with low computational costs on a training dataset of 1 million images to showcase its scalability and portability. Incorporating the CP method with latent distances offers a calibrated, sharp and efficient strategy to estimate the uncertainty of neural network force fields. In addition, the CP approach can also function as a promising strategy for calibrating uncertainty estimated by other approaches.
Comparison of Chain Conformation of Poly(vinyl alcohol) in Solutions and Melts from Quantum Chemistry Based Molecular Dynamics Simulations
Confirmations of 2,4-dihydroxypentane (DHP), a model molecule for poly(vinyl alcohol), have been studied by quantum chemistry (QC) calculations and molecular dynamics (MD) simulations. QC calculations at the 6-311G MP2 level show the meso tt conformer to be lowest in energy followed by the racemic tg, due to intramolecular hydrogen bond between the hydroxy groups. The Dreiding force field has been modified to reproduce the QC conformer energies for DHP. MD simulations using this force field have been carried out for DHP molecules in the gas phase, melt, and CHCl3 and water solutions. Extensive intramolecular hydrogen bonding is observed for the gas phase and CHCl3 solution, but not for the melt or aqueous solution, Such a condensed phase effect due to intermolecular interactions results in a drastic change in chain conformations, in agreement with experiments.
AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2
Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic’s evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (M pro ), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. Here, this work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic “warhead” to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of M pro with micromolar affinities (K I of 5.27 μM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.
Deconstructing Chirality: Probing Local and Nonlocal Effects in Azobenzene Derivatives with X-ray Circular Dichroism
Resolving molecular chirality on the atomic scale remains a critical challenge in chemistry. Conventional Optical Circular Dichroism spectroscopy often overlooks subtle and localized structural features. Here, we computationally investigate site-specific X-ray circular dichroism (XCD) across a series of trans-azobenzene derivatives to deconstruct and interpret chiroptical signals at the atomic level. Our modeling reveals that XCD is capable of distinguishing dichroic contributions arising from both a local chiral center and a global molecular twist, revealing their intricate interplay and potential for constructive or destructive interference. We show that sterically induced global distortions can dominate the XCD signal in some cases, even suppressing the response from the chiral center itself. This insight suggests a new molecular design principle for tuning the chiroptical activity, which we extend by proposing strategies to achieve unidirectional photoisomerization through steric gearing. Altogether, this work establishes a quantitative framework for engineering chiroptical responses, laying the foundation for the design of functional chiral systems utilizing principles of unidirectional molecular motor-like conformational dynamics.
Direct Observation of the Pressure-Induced Structural Variation in Gold Nanoclusters and the Correlated Optical Response
The ability to gradually modify the atomic structures of nanomaterials and directly identify such structural variation is important in nanoscience research. Here, we present the first example of a high-pressure single-crystal X-ray diffraction analysis of atomically precise metal nanoclusters. The pressure-dependent, subangstrom structural evolution of an ultrasmall gold nanoparticle, Au 25 S 18 , has been directly identified. We found that a 0.1 Å decrease of the Au–Au bond length could induce a blue-shift of 30 nm in the photoluminescence spectra of gold nanoclusters. From theoretical calculations, the origins of the blue-shift and enhanced photoluminescence under pressure are investigated, which are ascribed to molecular orbital symmetry and conformational locking, respectively. Here, the combination of the high-pressure in situ X-ray results with both theoretical and experimental optical spectra provides a direct and generalizable avenue to unveil the underlying structure–property relations for nanoclusters and nanoparticles which cannot be obtained through traditional physical chemistry measurements.
Tailoring Hierarchical Structure and Rare Earth Affinity of Compositionally Identical Polymers via Sequence Control
Macromolecule sequence, structure, and function are inherently intertwined. While well-established relationships exist in proteins, they are more challenging to define for synthetic polymer nanoparticles due to their molecular weight, sequence, and conformational dispersities. Furthermore, to explore the impact of sequence on nanoparticle structure, we synthesized a set of 16 compositionally identical, sequence-controlled polymers with distinct monomer patterning of dimethyl acrylamide and a bioinspired, structure-driving di(phenylalanine) acrylamide (FF). Sequence control was achieved through multiblock polymerizations, yielding unique ensembles of polymer sequences which were simulated by kinetic Monte Carlo simulations. Systematic analysis of the global (tertiary- and quaternary-like) structure in this amphiphilic copolymer series revealed the effect of multiple sequence descriptors: the number of domains, the hydropathy of terminal domains, and the patchiness (density) of FF within a domain, each of which impacted both chain collapse and the distribution of single- and multichain assemblies. Furthermore, both the conformational freedom of chain segments and local-scale, β-sheet-like interactions were sensitive to the patchiness of FF. To connect sequence, structure, and target function, we evaluated an additional series of nine sequence-controlled copolymers as sequestrants for rare earth elements (REEs) by incorporating a functional acrylic acid monomer into select polymer scaffolds. We identified key sequence variables that influence the binding affinity, capacity, and selectivity of the polymers for REEs. Collectively, these results highlight the potential of and boundaries of sequence control via multiblock polymerizations to drive primary sequence ensembles hierarchical structures, and ultimately the functionality of compositionally identical polymeric materials.
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.
Characterization of the structural forces governing the reversibility of the thermal unfolding of the human acidic fibroblast growth factor
Human acidic fibroblast growth factor (hFGF1) is an all beta-sheet protein that is involved in the regulation of key cellular processes including cell proliferation and wound healing. hFGF1 is known to aggregate when subjected to thermal unfolding. In this study, we investigate the equilibrium unfolding of hFGF1 using a wide array of biophysical and biochemical techniques. Systematic analyses of the thermal and chemical denaturation data on hFGF1 variants (Q54P, K126N, R136E, K126N/R136E, Q54P/K126N, Q54P/R136E, and Q54P/K126N/R136E) indicate that nullification of charges in the heparin-binding pocket can significantly increase the stability of wtFGF1. Triple variant (Q54P/K126N/R136E) was found to be the most stable of all the hFGF1 variants studied. With the exception of triple variant, thermal unfolding of wtFGF1 and the other variants is irreversible. Thermally unfolded triple variant refolds completely to its biologically native conformation. Microsecond-level molecular dynamic simulations reveal that a network of hydrogen bonds and salt bridges linked to Q54P, K126N, and R136E mutations, are responsible for the high stability and reversibility of thermal unfolding of the triple variant. In our opinion, the findings of the study provide valuable clues for the rational design of a stable hFGF1 variant that exhibits potent wound healing properties.
Unique Structural Features Relate to Evolutionary Adaptation of Cytochrome P450 in the Abyssal Zone
Cytochromes P450 (CYPs) form one of the largest enzyme superfamilies, with similar structural folds yet biological functions varying from synthesis of physiologically essential compounds to metabolism of myriad xenobiotics. Sterol 14α-demethylases (CYP51s) represent a very special P450 family, regarded as a possible evolutionary progenitor for all currently existing P450s. In metazoans CYP51 is critical for the biosynthesis of sterols including cholesterol. Here we determined the crystal structures of ligand-free CYP51s from the abyssal fish Coryphaenoides armatus and human-. Comparative sequence–structure–function analysis revealed specific structural elements that imply elevated conformational flexibility, uncovering a molecular basis for faster catalytic rates, lower substrate selectivity, and intrinsic resistance to inhibition. In addition, the C. armatus structure displayed a large-scale repositioning of structural segments that, in vivo, are immersed in the endoplasmic reticulum membrane and border the substrate entrance (the FG arm, >20 Å, and the β4 hairpin, >15 Å). The structural distinction of C. armatus CYP51, which is the first structurally characterized deep sea P450, suggests stronger involvement of the membrane environment in regulation of the enzyme function. We interpret this as a co-adaptation of the membrane protein structure with membrane lipid composition during evolutionary incursion to life in the deep sea.
Replication Data for: Deconstructing Chirality: Probing Local and Non-local Effects in Azobenzene Derivatives with X-ray Circular Dichroism
Resolving molecular chirality at the atomic scale remains a critical challenge in chemistry. Conventional Optical Circular Dichroism spectroscopy often overlooks subtle and localized structural features. Here, we computationally investigate site-specific X-ray Circular Dichroism (XCD) across a series of trans-azobenzene derivatives to deconstruct and interpret chiroptical signals at the atomic level. Our modeling reveals that XCD is capable of distinguishing dichroic contributions arising from both a local chiral center and global molecular twist, revealing their intricate interplay and potential for constructive or destructive interference. We show that sterically-induced global distortions can dominate the XCD signal in some cases, even suppressing the response from the chiral center itself. This insight suggests a new molecular design principle for tuning chiroptical activity, which we extend by proposing strategies to achieve unidirectional photoisomerization through steric gearing. Altogether, this work establishes a quantitative framework for engineering chiroptical responses, laying the foundation for the design of functional chiral systems utilizing principles of unidirectional molecular motor-like conformational dynamics.
Conformational and oligomeric states of SPOP from small-angle X-ray scattering and molecular dynamics simulations
Speckle-type POZ protein (SPOP) is a substrate adaptor in the ubiquitin proteasome system, and plays important roles in cell-cycle control, development, and cancer pathogenesis. SPOP forms linear higher-order oligomers following an isodesmic self-association model. Oligomerization is essential for SPOP’s multivalent interactions with substrates, which facilitate phase separation and localization to biomolecular condensates. Structural characterization of SPOP in its oligomeric state and in solution is, however, challenging due to the inherent conformational and compositional heterogeneity of the oligomeric species. Here, we develop an approach to simultaneously and self-consistently characterize the conformational ensemble and the distribution of oligomeric states of SPOP by combining small-angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations. We build initial conformational ensembles of SPOP oligomers using coarse-grained molecular dynamics simulations, and use a Bayesian/maximum entropy approach to refine the ensembles, along with the distribution of oligomeric states, against a concentration series of SAXS experiments. Our results suggest that SPOP oligomers behave as rigid, helical structures in solution, and that a flexible linker region allows SPOP’s substrate-binding domains to extend away from the core of the oligomers. Additionally, our results are in good agreement with previous characterization of the isodesmic self-association of SPOP. In the future, the approach presented here can be extended to other systems to simultaneously characterize structural heterogeneity and self-assembly.
Identification of structural transitions in bacterial fatty acid binding proteins that permit ligand entry and exit at membranes
Fatty acid (FA) transfer proteins extract FA from membranes and sequester them to facilitate their movement through the cytosol. Detailed structural information is available for these soluble protein–FA complexes, but the structure of the protein conformation responsible for FA exchange at the membrane is unknown. Staphylococcus aureus FakB1 is a prototypical bacterial FA transfer protein that binds palmitate within a narrow, buried tunnel. Here, we define the conformational change from a “closed” FakB1 state to an “open” state that associates with the membrane and provides a path for entry and egress of the FA. Using NMR spectroscopy, we identified a conformationally flexible dynamic region in FakB1, and X-ray crystallography of FakB1 mutants captured the conformation of the open state. In addition, molecular dynamics simulations show that the new amphipathic α-helix formed in the open state inserts below the phosphate plane of the bilayer to create a diffusion channel for the hydrophobic FA tail to access the hydrocarbon core and place the carboxyl group at the phosphate layer. The membrane binding and catalytic properties of site-directed mutants were consistent with the proposed membrane docked structure predicted by our molecular dynamics simulations. Finally, the structure of the bilayer-associated conformation of FakB1 has local similarities with mammalian FA binding proteins and provides a conceptual framework for how these proteins interact with the membrane to create a diffusion channel from the FA location in the bilayer to the protein interior.