Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “community clustering”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Distilling Knowledge from Ensembles of Cluster-Constrained-Attention Multiple-Instance Learners for Whole Slide Image Classification

The peculiar nature of whole slide imaging (WSI), digitizing conventional glass slides to obtain multiple high resolution images which capture microscopic details of a patient’s histopathological features, has garnered increased interest from the computer vision research community over the last two decades. Given the unique computational space and time complexity inherent to gigapixel-size whole slide image data, researchers have proposed novel machine learning algorithms to aid in the performance of diagnostic tasks in clinical pathology. One effective algorithm represents a Whole slide image as a bag of smaller image patches, which can be represented as low-dimension image patch embeddings. Weakly supervised deep-learning methods, such as cluster-constrained-attention multiple instance learning (CLAM), have shown promising results when combined with image patch embeddings. While traditional ensemble classifiers yield improved task performance, such methods come with a steep cost in model complexity. Through knowledge distillation, it is possible to retain some performance improvements from an ensemble, while minimizing costs to model complexity. In this work, we implement a weakly supervised ensemble using clustering-constrained-attention multiple-instance learners (CLAM), which uses attention and instance-level clustering to identify task salient regions and feature extraction in whole slides. By applying logit-based and attention-based knowledge distillation, we show it is possible to retain some performance improvements resulting from the ensemble at zero cost to model complexity.

Alamudun, Folami↗
DOE OSTI2025

dCache project status and update

The dCache project delivers an open-source, massively scalable, distributed storage system deployed internationally to satisfy today’s scientists’ ever-demanding storage requirements. Its multifaceted approach supports different use cases with the same storage, from high throughput data ingest, data sharing over wide area networks, efficient access from HPC clusters, and longterm data persistence on tertiary storage. Even though dCache was initially developed for HEP experiments, today, it is used by various scientific communities, including astrophysics, biomed, and life science, each with their specific requirements. To match the needs of these new communities and keep up with the scaling demands of existing experiments, dCache is permanently evolving. With this contribution, we would like to highlight the recent developments in dCache regarding integration with CERN Tape Archive (CTA), advanced metadata handling, token-based authorization support, bulk API for QoS transitions, REST API to control interaction with the tape system, and future development directions.

Mkrtchyan, Tigran [DESY]↗
DOE OSTI2022

Resilient Node Cluster Analysis Tool (ReNCAT) v. 2.1

SAND2021-6533 O The Resilient Node Cluster Analysis Tool (ReNCAT) identifies portions of a power distribution network that can be isolated into microgrids in order to supply services to a community during disruptive events. Candidate microgrid configurations are evaluated using a burden metric, a measure of the relative effort community members must expend to acquire services required for their well-being. Results present multiple candidate configurations, each with a different balance of cost and burden. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Melander, Darryl↗
DOE OSTI2025

dCache: The Storage System of Choice for Data-Intensive Applications

The ever-increasing volumes of data produced by modern scientific facilities like EuXFEL and LHC put significant stress on data management infrastructure operated by laboratories and research centers. The challenges to be addressed span the entire data life cycle, from ingest and efficient data analysis to long-term preservation, typically involving large tape libraries. dCache, a storage system developed in collaboration between the Deutsches Elektronen-Synchrotron (DESY), Fermi National Accelerator Laboratory, and Nordic e-Infrastructure Collaboration (NeIC), is designed to manage a large number of disk servers and to facilitate transparent data migration to and from archival storage. Its multifaceted approach offers a unified method to support a variety of scientific use cases with the same storage infrastructure, including high-throughput data ingest, data sharing over wide area networks, efficient access from HPC clusters, and long-term data preservation on tertiary storage. Initially developed for high energy physics (HEP) experiments, dCache is now used by various scientific communities, including astrophysics, biomedical research, and life sciences, each having specific requirements. This paper presents architecture, deployment strategies, performance and scalability enhancements, and recent advancements in dCache addressing the needs of scientific communities. Finally, we touch on the development and release process, ensuring the software’s high quality.

DCache↗
DOE OSTI2020

TeraChem Cloud: A High-Performance Computing Service for Scalable Distributed GPU-Accelerated Electronic Structure Calculations

The encapsulation and commoditization of electronic structure arise naturally as interoperability, and the use of nontraditional compute resources (e.g., new hardware accelerators, cloud computing) remains important for the computational chemistry community. Here, we present TERACHEM CLOUD, a high-performance computing service (HPCS) that offers on-demand electronic structure calculations on both traditional HPC clusters and cloud-based hardware. The framework is designed using off-the-shelf web technologies and containerization to be extremely scalable and portable. Within the HPCS model, users can quickly develop new methods and algorithms in an interactive environment on their laptop while allowing TERACHEM CLOUD to distribute ab initio calculations across all available resources. This approach greatly increases the accessibility of hardware accelerators such as graphics processing units (GPUs) and flexibility for the development of new methods as additional electronic structure packages are integrated into the framework as alternative backends. Cost-performance analysis indicates that traditional nodes are the most cost-effective long-term solution, but commercial cloud providers offer cutting-edge hardware with competitive rates for short-term large-scale calculations. We demonstrate the power of the TERACHEM CLOUD framework by carrying out several showcase calculations, including the generation of 300,000 density functional theory energy and gradient evaluations on medium-sized organic molecules and reproducing 300 fs of nonadiabatic dynamics on the B800-B850 antenna complex in LH2, with the latter demonstration using over 50 Tesla V100 GPUs in a commercial cloud environment in 8 h for approximately $1250.

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

Nanocluster Active Sites Formed on Heterogeneous Thermal Catalysts and Electrocatalysts by Operando Reactive Environments

Here, in this Viewpoint, we summarize our recent studies in this endeavor, ending with a forward-looking perspective into the remaining challenges and open questions for consideration by the catalysis community. Through a computationally efficient framework of systematically investigating metal atom ejection from and migration onto metal surfaces, leading to the formation of metal atom clusters, we show that the surface of a metal catalyst is dynamic and much less rigid than previously thought. We suggest that the operando formed metal atom clusters are novel active sites, which can dominate the activity of steps or kinks (let aside terraces) from where their constituent atoms were ejected. Through more readily responding to reaction conditions by forming catalytically active clusters, solid metal catalysts approach not only homogeneous catalyst motifs, thus establishing a bridge between homogeneous and heterogeneous catalysis, but also liquid metals which have been suggested to be better than traditional heterogeneous catalysts.

adsorption↗
DOE OSTI2020

In Vivo Molecular Insights into Syntrophic Geobacter Aggregates

Direct interspecies electron transfer (DIET) has been considered as a novel and highly efficient strategy in both natural anaerobic environments and artificial microbial fuel cells. A syntrophic model consisting of Geobacter metallireducens and Geobacter sulfurreducens was studied in this work. We conducted in vivo molecular mapping of the outer surface of the syntrophic community as the interface of nutrients and energy exchange. System for Analysis at the Liquid Vacuum Interface combined with time-of-flight secondary ion mass spectrometry was employed to capture the molecular distribution of syntrophic Geobacter communities in the living and hydrated state. Principal component analysis with selected peaks revealed that syntrophic Geobacter aggregates were well differentiated from other control samples, including syntrophic planktonic cells, pure cultured planktonic cells, and single population biofilms. Our in vivo imaging indicated that a unique molecular surface was formed. Specifically, aromatic amino acids, phosphatidylethanolamine components, and large water clusters were identified as key components that favored the DIET of syntrophic Geobacter aggregates. Moreover, the molecular changes in depths of the Geobacter aggregates were captured using dynamic depth profiling. Our findings shed new light on the interface components supporting electron transfer in syntrophic communities based on in vivo molecular imaging.

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

dCache: Inter-disciplinary storage system

The dCache project provides open-source software deployed internationally to satisfy ever more demanding storage requirements. Its multifaceted approach provides an integrated way of supporting different use-cases with the same storage, from high throughput data ingest, data sharing over wide area networks, efficient access from HPC clusters and long term data persistence on a tertiary storage. Though it was originally developed for the HEP experiments, today it is used by various scientific communities, including astrophysics, biomed, life science, which have their specific requirements. In this paper we describe some of the new requirements as well as demonstrate how dCache developers are addressing them.

Mkrtchyan, Tigran↗
DOE OSTI2022

Theoretical Perspective on Operando Spectroscopy of Fluxional Nanocatalysts

Improvements in operando spectroscopy have enabled the catalysis community to investigate the dynamic nature of catalysts under operating conditions with increasing detail. Still, the highly dynamic nature of some catalysts, such as fluxional supported subnano clusters, presents a formidable challenge even for the most state-of-the-art techniques. Here, the reason is that such fluxional catalytic interfaces contain a variety of thermally accessible states. Operando spectroscopies used in catalysis generally fall into two categories: ensemble-based techniques, which provide spectra containing the signals of the entire ensemble of states of the catalyst and are not necessarily dominated by the most active species, and localized techniques, which provide atomistic-level information about the dynamics of active sites in a very small area, which might not include the most active species. Combining many different kinds of techniques can provide detailed insight; however, we propose that effective utilization of specific computational techniques and approaches within the fluxionality paradigm can fill the gap and enable atomistic characterization of the most relevant catalytic sites.

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

US Department of Energy, Office of Science High Performance Computing Facility Operational Assessment 2019 Oak Ridge Leadership Computing Facility

Oak Ridge National Laboratory's (ORNL's) Leadership Computing Facility (OLCF) continues to surpass its operational target goals: supporting users; delivering fast, reliable computational ecosystems; creating innovative solutions for high performance computing (HPC) needs; and managing risks, safety, and security associated with operating some of the most powerful computers in the world. The results can be seen in the cutting-edge science conducted by users and the praise from the research community. Calendar year (CY) 2019 was a big year as OLCF staff ran five world-class resources (the leadershipclass computers Titan and Summit, the large analysis cluster called Eos, and the massive parallel filesystems called Atlas and Alpine)) and also began power and cooling upgrades for a 2021 exascale system called Frontier. While continuing exceptional operation of Titan, Eos, and Rhea, the OLCF released the Summit supercomputer for production on January 1, 2019. Summit debuted as the most capable and efficient system in its class and has been recognized as the most powerful system in the world for its performance on both the high performance linpack (HPL) and conjugate gradient (HPCG) benchmark applications since June 2018 according to TOP500. Summit represents the culmination of a multiyear effort between the OLCF, IBM, NVIDIA, and Mellanox to deliver a system that is unmatched for modeling, simulation, data analysis, and learning. To hit the ground running with science-ready applications on day one, application teams worked closely with the OLCF through the Center for Accelerated Application Readiness (CAAR) program for years in advance of the Summit deployment. CY 2019 was filled with outstanding results and accomplishments: a very high rating from users on overall satisfaction for the sixth year in a row; a tremendous amount of core-hours delivered to researchers from two leadership-class systems; and success in delivering on the allocation split of roughly 60%, 30%, and 10% of core-hours offered for the Innovative and Novel Computational Impact on Theory and Experiment (INCITE), Advanced Scientific Computing Research Leadership Computing Challenge (ALCC), and Director's Discretionary (DD) programs, respectively (see Operational Performance section). These accomplishments, coupled with the high utilization rates (overall and capability usage), represent the fulfillment of the promise of both leadership-class machines: efficient facilitation of leadership-class computational applications. Table ES.1 presents a summary of the 2019 OLCF metric targets and the associated results. More information can be found in the Operational Performance section for each OLCF resource. The scientific accomplishments of OLCF users are a strong indication of long-term operational success, with publications this year in such notable journals and publications as Nature, Nature Physics, Nature Plants, Physical Review X, Journal of the American Physical Society, Cell, Nano Letters, and Trends in Biotechnology. Crucial domain-specific discoveries facilitated by resources at the OLCF are described in the High Performance Computing Facility Operational Assessment 2019 Oak Ridge Leadership Computing Facility (OAR) Strategic Results section. For example, researchers used Summit to pinpoint and understand the production of proteins from genetic information, including mutations and the functional expression of disease (Section 8.2).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Clonality, local population structure, and gametophyte sex ratios in cryptic species of the Sphagnum magellanicum complex

Sphagnum (peatmoss) comprises a moss (Bryophyta) clade with approximately 300-500 species. The genus has unparalleled ecological importance because Sphagnum-dominated peatlands store almost a third of the terrestrial carbon pool and peatmosses engineer the formation and microtopography of peatlands. Genomic resources for Sphagnum are being actively expanded, but many aspects of their biology are still poorly known. Among these are the degree to which Sphagnum species reproduce asexually, and the relative frequencies of male and female gametophytes in these haploid-dominant plants. Here, we assess clonality and gametophyte sex ratios and test hypotheses about the local-scale distribution of clones and sexes in four North American species of the S. magellanicum complex. These four species are difficult to distinguish morphologically and are very closely related. We also assess microbial communities associated with Sphagnum host plant clones and sexes at two sites. 405 samples of the four species, representing 57 populations, were subjected to RADseq. Analyses of population structure and clonality based on the molecular data utilized both phylogenetic and phenetic approaches. Multi-locus genotypes (genets) were identified using the RADseq data. Sexes of sampled ramets were determined using a molecular approach that utilized coverage of loci on the sex chromosomes after the method was validated using a sample of plants that expressed sex phenotypically. Sex ratios were estimated for each species, and populations within species. Differences in fitness between genets was estimated as the numbers of ramets each genet comprised. Degree of clonality (numbers of genets/numbers of ramets [samples]) within species, among sites, and between gametophyte sexes were estimated. Sex ratios were estimated for each species, and populations within species. Sphagnum-associated microbial communities were assessed at two sites in relation to Sphagnum clonality and sex. All four species appear to engage in a mixture of sexual and asexual (clonal) reproduction. A single ramet represents most genets but 2-8 ramets were detected for some genets. Only one genet is represented by ramets in multiple populations; all other genets are restricted to a single population. Within populations ramets of individual genets are spatially clustered, suggesting limited dispersal even within peatlands. Sex ratios are male-biased in S. diabolicum but female-biased in the other three species, although significantly so only in S. divinum. Neither species nor males/females differ in levels of clonal propagation. At St. Regis Lake (NY) and Franklin Bog (VT), microbial community composition is strongly differentiated between the sites, but differences between species, genets, and sexes were not detected. Within S. divinum, however, female gametophytes harbored 2-3 times the number oi of microbial taxa as males. These four Sphagnum species all exhibit a similar reproductive patterns that result from a mixture of sexual and asexual reproduction. The spatial patterns of clonally replicated ramets of genets suggest that these species fall between the so-called phalanx patterns where genets abut one another but do not extensively mix, because of limited ramet fragmentation, and the guerrilla patterns where extensive genet fragmentation and dispersal results in greater mixing of different genets. Although sex ratios in bryophytes are most often female-biased, both male and female biases occur in this complex of closely related species. The association of far greater microbial diversity for female gametophytes in S. divinum, which has a female-biased sex ratio, suggests additional research to determine if levels of microbial diversity are consistently correlated with differing patterns of sex ratio biases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Interventions to Disrupt Coronavirus Disease Transmission at a University, Wisconsin, USA, August–October 2020

University settings have demonstrated potential for coronavirus disease (COVID-19) outbreaks; they combine congregate living, substantial social activity, and a young population predisposed to mild illness. Using genomic and epidemiologic data, we describe a COVID-19 outbreak at the University of Wisconsin-Madison, Madison, Wisconsin, USA. During August-October 2020, a total of 3,485 students, including 856/6,162 students living in dormitories, tested positive. Case counts began rising during move-in week, August 25–31, 2020, then rose rapidly during September 1–11, 2020. The university initiated multiple prevention efforts, including quarantining 2 dormitories; a subsequent decline in cases was observed. Genomic surveillance of cases from Dane County, in which the university is located, did not find evidence of transmission from a large cluster of cases in the 2 quarantined dorms during the outbreak. Coordinated implementation of prevention measures can reduce COVID-19 spread in university settings and may limit spillover to the surrounding community.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Profiles of upcoming HPC Applications and their Impact on Reservation Strategies

With the expected convergence between HPC, BigData and AI, new applications with different profiles are coming to HPC infrastructures. Here, we aim at better understanding the features and needs of these applications in order to be able to run them efficiently on HPC platforms. The approach followed is bottom-up: we study thoroughly an emerging application from the neuroscience community (SLANT) to understand its behavior. Based on these observations, we derive a generic, yet simple, application model (namely, a linear sequence of stochastic jobs). We expect this model to be representative for a large set of upcoming applications that require the computational power of HPC clusters without fitting the typical behavior of large-scale traditional applications. In a second step, we show how one can manipulate this generic model in a scheduling framework. Specifically we consider the problem of making reservations (both time and memory) for an execution on an HPC platform. We derive solutions using the model of the first step of this work. We experimentally show the robustness of the model, even with very few data or with another application, to generate the model, and provide performance gains with regards to standard and more recent approaches used in the neuroscience community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Chemolithoautotroph distributions across the subsurface of a convergent margin

Subducting oceanic crusts release fluids rich in biologically relevant compounds into the overriding plate, fueling subsurface chemolithoautotrophic ecosystems. To understand the impact of subsurface geochemistry on microbial communities, we collected fluid and sediments from 14 natural springs across a ~200 km transect across the Costa Rican convergent margin and performed shotgun metagenomics. The resulting 404 metagenome-assembled genomes (MAGs) cluster into geologically distinct regions based on MAG abundance patterns: outer forearc-only (25% of total relative abundance), forearc/arc-only (38% of total relative abundance), and delocalized (37% of total relative abundance) clusters. In the outer forearc, Thermodesulfovibrionia, Candidatus Bipolaricaulia, and Firmicutes have hydrogenotrophic sulfate reduction and Wood-Ljungdahl (WL) carbon fixation pathways. In the forearc/arc, Anaerolineae, Ca. Bipolaricaulia, and Thermodesulfovibrionia have sulfur oxidation, nitrogen cycling, microaerophilic respiration, and WL, while Aquificae have aerobic sulfur oxidation and reverse tricarboxylic acid carbon fixation pathway. Transformation-based canonical correspondence analysis shows that MAG distribution corresponds to concentrations of aluminum, iron, nickel, dissolved inorganic carbon, and phosphate. While delocalized MAGs appear surface-derived, the subsurface chemolithoautotrophic, metabolic, and taxonomic landscape varies by the availability of minerals/metals and volcanically derived inorganic carbon. Furthermore, the WL pathway persists across all samples, suggesting that this versatile, energy-efficient carbon fixation pathway helps shape convergent margin subsurface ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Network Analysis of Academic Medical Center Websites in the United States

Healthcare resources are published annually in repositories such as the AHA Annual Survey Database TM . However, these data repositories are created via manual surveying techniques which are cumbersome in collection and not updated as frequently as website information of the respective hospital systems represented. Also, this resource is not widely available to patients in an easy-to-use format. Network analysis techniques have the potential to create topological maps which serve to aid in pathfinding for patients in their search for healthcare services. This study explores the topological structure of forty United States academic health center websites. Network analysis is utilized to analyze and visualize 48,686 webpages. Several elements of network structure are examined including basic network properties, and centrality measures distributions. The Louvain community detection algorithm is used to examine the extent to which these techniques allow identification of healthcare resources within networks. The results indicate that websites with related healthcare services tend to form observable clusters useful in mapping key resources within a hospital system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Genome Analysis of a Verrucomicrobial Endosymbiont With a Tiny Genome Discovered in an Antarctic Lake

Organic Lake in Antarctica is a marine-derived, cold (–13°C), stratified (oxic-anoxic), hypersaline (>200 gl –1 ) system with unusual chemistry (very high levels of dimethylsulfide) that supports the growth of phylogenetically and metabolically diverse microorganisms. Symbionts are not well characterized in Antarctica. However, unicellular eukaryotes are often present in Antarctic lakes and theoretically could harbor endosymbionts. Here, we describe Candidatus Organicella extenuata, a member of the Verrucomicrobia with a highly reduced genome, recovered as a metagenome-assembled genome with genetic code 4 (UGA-to-Trp recoding) from Organic Lake. It is closely related to Candidatus Pinguicocccus supinus (163,218 bp, 205 genes), a newly described cytoplasmic endosymbiont of the freshwater ciliate Euplotes vanleeuwenhoeki (Serra et al., 2020). At 158,228 bp (encoding 194 genes), the genome of Ca. Organicella extenuata is among the smallest known bacterial genomes and similar to the genome of Ca. Pinguicoccus supinus (163,218 bp, 205 genes). Ca. Organicella extenuata retains a capacity for replication, transcription, translation, and protein-folding while lacking any capacity for the biosynthesis of amino acids or vitamins. Notably, the endosymbiont retains a capacity for fatty acid synthesis (type II) and iron–sulfur (Fe-S) cluster assembly. Metagenomic analysis of 150 new metagenomes from Organic Lake and more than 70 other Antarctic aquatic locations revealed a strong correlation in abundance between Ca. Organicella extenuata and a novel ciliate of the genus Euplotes. Like Ca. Pinguicoccus supinus, we infer that Ca. Organicella extenuata is an endosymbiont of Euplotes and hypothesize that both Ca. Organicella extenuata and Ca. Pinguicocccus supinus provide fatty acids and Fe-S clusters to their Euplotes host as the foundation of a mutualistic symbiosis. The discovery of Ca. Organicella extenuata as possessing genetic code 4 illustrates that in addition to identifying endosymbionts by sequencing known symbiotic communities and searching metagenome data using reference endosymbiont genomes, the potential exists to identify novel endosymbionts by searching for unusual coding parameters.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Assessing CESM2 Clouds and Their Response to Climate Change Using Cloud Regimes

Abstract The Community Earth System Model, version 2 (CESM2), has a very high climate sensitivity driven by strong positive cloud feedbacks. To evaluate the simulated clouds in the present climate and characterize their response with climate warming, a clustering approach is applied to three independent satellite cloud products and a set of coupled climate simulations. Using k -means clustering with a Wasserstein distance cost function, a set of typical cloud configurations is derived for the satellite cloud products. Using satellite simulator output, the model clouds are classified into the observed cloud regimes in both current and future climates. The model qualitatively reproduces the observed cloud configurations in the historical simulation using the same time period as the satellite observations, but it struggles to capture the observed heterogeneity of clouds which leads to an overestimation of the frequency of a few preferred cloud regimes. This problem is especially apparent for boundary layer clouds. Those low-level cloud regimes also account for much of the climate response in the late twenty-first century in four shared socioeconomic pathway simulations. The model reduces the frequency of occurrence of these low-cloud regimes, especially in tropical regions under large-scale subsidence, in favor of regimes that have weaker cloud radiative effects.

58 GEOSCIENCES↗
DOE OSTI2021

OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.

59 BASIC BIOLOGICAL SCIENCES↗