Engineering Papers⌕ Search

Engineering topics

Ferguson, Andrew L.

Publications and source records attributed to Ferguson, Andrew L..

12 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Permutationally Invariant Networks for Enhanced Sampling (PINES): Discovery of Multimolecular and Solvent-Inclusive Collective Variables

The typically rugged nature of molecular free energy landscapes can frustrate efficient sampling of the thermodynamically relevant phase space due to the presence of high free energy barriers. Enhanced sampling techniques can improve phase space exploration by accelerating sampling along particular collective variables (CVs). A number of techniques exist for data-driven discovery of CVs parameterizing the important large scale motions of the system. A challenge to CV discovery is learning CVs invariant to symmetries of the molecular system, frequently rigid translation, rigid rotation, and permutational relabeling of identical particles. Of these, permutational invariance has proved a persistent challenge in frustrating the data-driven discovery of multi-molecular CVs in systems of self-assembling particles and solvent-inclusive CVs for solvated systems. In this work, we integrate Permutation Invariant Vector (PIV) featurizations with autoencoding neural networks to learn nonlinear CVs invariant to translation, rotation, and permutation, and perform interleaved rounds of CV discovery and enhanced sampling to iteratively expand sampling of configurational phase space and obtain converged CVs and free energy landscapes. Here, we demonstrate the Permutationally Invariant Network for Enhanced Sampling (PINES) approach in applications to the self-assembly of a 13-atom Argon cluster, association/dissociation of a NaCl ion pair in water, and hydrophobic collapse of a C 45 H 92 n-pentatetracontane polymer chain. We make the approach freely available as a new module within the PLUMED2 enhanced sampling libraries.

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

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↗
DOE OSTI2023

Data-Driven Discovery of Linear Molecular Probes with Optimal Selective Affinity for PFAS in Water

Approaches to tackle the wide and growing variety of highly persistent per- and polyfluoroalkyl substances (PFAS) are of pressing global need because of their detrimental human health effects, such as cancer, birth defects, and hormone imbalance. Sensitive, selective, and easy-to-use real-time sensors to monitor and detect PFAS and sorbents to extract them are critical to meeting government-mandated environmental concentrations. In this work, we combine all-atom molecular dynamics simulations, enhanced sampling, deep representational learning, and Bayesian optimization to perform high-throughput virtual screening for highly sensitive and selective molecular probes. Our molecular design space consists of 3850 linear hydrocarbon chains with varying degrees of halogenation with and without amine- and phosphine-based headgroups. By employing a data-driven search process, we efficiently explore the molecular design space to optimize the sensitivity to perfluorooctanesulfonic acid (PFOS) as a prototypical PFAS analyte and selectivity relative to a sodium dodecyl sulfate (SDS) interferent. We calculate 504 Gibbs free energies of probe-analyte and probe-interferent interactions and identify probes with PFOS association free energies of up to (-ΔG PFOS ) = 9.8 ± 0.2 kJ/mol and selectivities relative to SDS of (-ΔΔG PFOS–SDS ) = 3.1 ± 1.5 kJ/mol. A C 11 Br 23 P(CH 3 ) 2 probe containing 11 backbone brominated carbons and a tertiary phosphine headgroup possesses the most sensitive binding constant to PFOS within the defined search space of K b PFOS = 177.4 ± 12.7, and a semibrominated probe C 5 H 11 C 7 Br 14 N(CH 3 ) 2 containing 12 backbone carbons and a tertiary amine headgroup possesses the highest selectivity relative to SDS of K b PFOS /K b SDS = 4.6 ± 1.7. A retrospective analysis of our data to extract interpretable design rules reveals that the sensitivity of linear hydrogenated probes increases by approximately 1 kJ/mol per C–C bond. The addition or removal of halogen atoms and amine or phosphine headgroups produces nonmonotonic changes in both sensitivity and selectivity with changes to the sensitivity of up to 2.5 kJ/mol. Finally, this work places empirical limitations on the performance of a wide range of linear probes for PFOS detection and offers a generic strategy for high-throughput computational screening to promote selective and sensitive binding.

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

Thermodynamic Basis for the Stabilization of Helical Peptoids by Chiral Sidechains

Peptoids are a class of highly customizable biomimetic foldamers that retain properties from both proteins and polymers. It has been shown that peptoids can adopt peptide-like secondary structures through the careful selection of sidechain chemistries, but the underlying conformational landscapes that drive these assemblies at the molecular level remain poorly understood. Given the high flexibility of the peptoid backbone, it is essential that methods applied to study peptoid secondary structure formation possess the requisite sensitivity to discriminate between structurally similar yet energetically distinct microstates. In this work, a generalizable simulation scheme is used to robustly sample the complex folding landscape of various 12mer polypeptoids, resulting in a predictive model that links sidechain chemistry with preferential assembly into one of 12 accessible backbone motifs. Using a variant of the metadynamics sampling method, four peptoid dodecamers are simulated in water: sarcosine, N-(1-phenylmethyl)glycine (Npm), (S)-N-(1-phenylethyl)glycine (Nspe), and (R)-N-(1-phenylethyl)glycine (Nrpe)–to determine the underlying entropic and energetic impacts of hydrophobic and chiral peptoid sidechains on secondary structure formation. Our results indicate that the driving forces to assemble Nrpe and Nspe sequences into polyproline type-I helices in water are found to be enthalpically driven, with small benefits from an entropic gain for isomerization and steric strain due to the presence of the chiral center. The minor entropic gains from bulky chiral sidechains in Nrpe- and Nspe-containing peptoids can be explained through increased configurational entropy in the cis state. However, overall assembly into a helix is found to be overall entropically unfavorable. Furthermore, these results highlight the importance of considering the many various competing interactions in the rational design of peptoid secondary structure building blocks.

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

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH 2 -capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young’s modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. Here, we construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a “Goldilocks” tube radius at which capillary wave fluctuations in the tube wall are minimized.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hierarchical Materials from High Information Content Macromolecular Building Blocks: Construction, Dynamic Interventions, and Prediction

Hierarchical materials that exhibit order over multiple length scales are ubiquitous in nature. Because hierarchy gives rise to unique properties and functions, many have sought inspiration from nature when designing and fabricating hierarchical matter. More and more, however, nature’s own high-information content building blocks, proteins, peptides, and peptidomimetics, are being coopted to build hierarchy because the information that determines structure, function, and interfacial interactions can be readily encoded in these versatile macromolecules. Here, we take stock of recent progress in the rational design and characterization of hierarchical materials produced from high-information content blocks with a focus on stimuli-responsive and “smart” architectures. We also review advances in the use of computational simulations and data-driven predictions to shed light on how the side chain chemistry and conformational flexibility of macromolecular blocks drive the emergence of order and the acquisition of hierarchy and also on how ionic, solvent, and surface effects influence the outcomes of assembly. Furthermore, continued progress in the above areas will ultimately usher in an era where an understanding of designed interactions, surface effects, and solution conditions can be harnessed to achieve predictive materials synthesis across scale and drive emergent phenomena in the self-assembly and reconfiguration of high-information content building blocks.

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

Active learning of polarizable nanoparticle phase diagrams for the guided design of triggerable self-assembling superlattices

Polarizable nanoparticles are of interest in materials science because of their rich and complex phase behavior that can be used to engineer nanostructured materials with long-range crystalline order. To understand and rationally navigate the design space of polarizable nanoparticles for self-assembling highly ordered superlattices, we developed a coarse-grained computational model to describe the nanoparticle-nanoparticle interactions in implicit solvent and employ the computationally efficient image method to model many-body polarization interactions. We conducted high-throughput virtual screening over a five-dimensional particle design space spanned by temperature, particle size, particle charge, particle dielectric, and solvent dielectric using enhanced sampling molecular dynamics calculations within an active learning framework to efficiently map out the regions of thermodynamic stability of the self-assembled aggregates. We validate our predictions in comparisons against small angle x-ray scattering measurements of gold nanoparticles surface functionalized with metal chalcogenide ligands. Lastly, we use our validated phase maps to computationally design switchable nanostructured materials capable of triggered assembly and disassembly as a function of temperature and solvent dielectric with potential applications as sensors, smart windows, optoelectronic devices, and in medical diagnostics.

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

Hierarchical Self-Assembly Pathways of Peptoid Helices and Sheets

Peptoids (N-substituted glycines) are a class of tailorable synthetic peptidomic polymers. Amphiphilic diblock peptoids have been engineered to assemble 2D crystalline lattices with applications in catalysis and molecular separations. Assembly is induced in an organic solvent/water mixture by evaporating the organic phase, but the assembly pathways remain uncharacterized. We conduct all-atom molecular dynamics simulations of Nbrpe6Nc6 as a prototypical amphiphilic diblock peptoid comprising an NH 2 -capped block of six hydrophobic N-((4-bromophenyl)ethyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail. We identify a thermodynamically controlled assembly mechanism by which monomers assemble into disordered aggregates that self-order into 1D chiral helical rods then 2D achiral crystalline sheets. We support our computational predictions with experimental observations of 1D rods using small-angle X-ray scattering, circular dichroism, and atomic force microscopy and 2D crystalline sheets using X-ray diffraction and atomic force microscopy. Furthermore, this work establishes a new understanding of hierarchical peptoid assembly and principles for the design of peptoid-based nanomaterials.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Pushing and Pulling: A Dual pH Trigger Controlled by Varying the Alkyl Tail Length in Heme Coordinating Peptide Amphiphiles

Some organisms in nature that undergo anaerobic respiration utilize 1D nanoscale arrays of densely packed cytochromes containing the molecule heme. The assemblies can be mimicked with 1D nanoscale fibrils composed of peptide amphiphiles designed to coordinate heme in dense arrays. To create such materials and assemblies it is critical to understand the assembly process and what controls the various aspects of hierarchical assembly. MD simulations suggest that shorter alkyl chains on the peptide lead to more dynamic structures than the peptides with longer chains that yield kinetically trapped states. The hydration parameters manifest themselves experimentally through the observation of a dual pH trigger which controls the peptide assembly rate, the heme binding affinity, and heme organization kinetics. Great strides in understanding the relative complexity of the self-assembly process in relation to incorporating a functional moiety like heme opens up many possibilities in developing abiotic assemblies for bioelectronic devices and assemblies

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

Molecular latent space simulators

Latent space simulators learn kinetic models for atomistic simulations and generate novel trajectories at six orders of magnitude lower cost.

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

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗