Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “membrane potential”

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 55 records · Page 3

Plastic recognition and electrogenic uniport translocation of 1 st -, 2 nd -, and 3 rd -row transition and post-transition metals by primary-active transmembrane P 1B-2 -type ATPase pumps

Transmembrane P 1B -type ATPase pumps catalyze the extrusion of transition metal ions across cellular lipid membranes to maintain essential cellular metal homeostasis and detoxify toxic metals. Zn(II)-pumps of the P 1B-2 -type subclass, in addition to Zn 2+ , select diverse metals (Pb 2+ , Cd 2+ and Hg 2+ ) at their transmembrane binding site and feature promiscuous metal-dependent ATP hydrolysis in the presence of these metals. Yet, a comprehensive understanding of the transport of these metals, their relative translocation rates, and transport mechanism remain elusive. We developed a platform for the characterization of primary-active Zn(II)-pumps in proteoliposomes to study metal selectivity, translocation events and transport mechanism in real-time, employing a “multi-probe” approach with fluorescent sensors responsive to diverse stimuli (metals, pH and membrane potential). Together with atomic-resolution investigation of cargo selection by X-ray absorption spectroscopy (XAS), we demonstrate that Zn(II)-pumps are electrogenic uniporters that preserve the transport mechanism with 1 st -, 2 nd - and 3 rd -row transition metal substrates. Promiscuous coordination plasticity, guarantees diverse, yet defined, cargo selectivity coupled to their translocation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identification of a secretory heme‐binding protein from Nocardia seriolae involved in cell apoptosis

Abstract According to the whole‐genome bioinformatics analysis, a heme‐binding protein from Nocardia seriolae (HBP) was found. HBP was predicted to be a bacterial secretory protein, located at mitochondrial membrane in eukaryotic cells and have a similar protein structure with the heme‐binding protein of Mycobacterium tuberculosis , Rv0203. In this study, HBP was found to be a secretory protein and co‐localized with mitochondria in FHM cells. Quantitative analysis of mitochondrial membrane potential value, caspase‐3 activity, and transcription level of apoptosis‐related genes suggested that overexpression of HBP protein can induce cell apoptosis. In conclusion, HBP was a secretory protein which may target to mitochondria and involve in cell apoptosis in host cells. This research will promote the function study of HBP and deepen the comprehension of the virulence factors and pathogenic mechanisms of N. seriolae .

Wen, Yiming↗

Cation and Anion Channelrhodopsins: Sequence Motifs and Taxonomic Distribution

Cation and anion channelrhodopsins (CCRs and ACRs, respectively) primarily from two algal species, Chlamydomonas reinhardtii and Guillardia theta, have become widely used as optogenetic tools to control cell membrane potential with light. We mined algal and other protist polynucleotide sequencing projects and metagenomic samples to identify 75 channelrhodopsin homologs from four channelrhodopsin families, including one revealed in dinoflagellates in this study. We carried out electrophysiological analysis of 33 natural channelrhodopsin variants from different phylogenetic lineages and 10 metagenomic homologs in search of sequence determinants of ion selectivity, photocurrent desensitization, and spectral tuning in channelrhodopsins. Our results show that association of a reduced number of glutamates near the conductance path with anion selectivity depends on a wider protein context, because prasinophyte homologs with a glutamate pattern identical to that in cryptophyte ACRs are cation selective. Desensitization is also broadly context dependent, as in one branch of stramenopile ACRs and their metagenomic homologs, its extent roughly correlates with phylogenetic relationship of their sequences. Regarding spectral tuning, we identified two prasinophyte CCRs with red-shifted spectra to 585 nm. They exhibit a third residue pattern in their retinal-binding pockets distinctly different from those of the only two types of red-shifted channelrhodopsins known (i.e., the CCR Chrimson and RubyACRs). In cryptophyte ACRs we identified three specific residue positions in the retinal-binding pocket that define the wavelength of their spectral maxima. Lastly, we found that dinoflagellate rhodopsins with a TCP motif in the third transmembrane helix and a metagenomic homolog exhibit channel activity. Channelrhodopsins are widely used in neuroscience and cardiology as research tools and are considered prospective therapeutics, but their natural diversity and mechanisms remain poorly characterized. Genomic and metagenomic sequencing projects are producing an ever-increasing wealth of data, whereas biophysical characterization of the encoded proteins lags behind. In this study, we used manual and automated patch clamp recording of representative members of four channelrhodopsin families, including a family in dinoflagellates that we report in this study. Our results contribute to a better understanding of molecular determinants of ionic selectivity, photocurrent desensitization, and spectral tuning in channelrhodopsins.

59 BASIC BIOLOGICAL SCIENCES↗

Engineered Methylobacterium extorquens grows well on methoxylated aromatics due to its formaldehyde metabolism and stress response

ABSTRACT Lignin is a vast yet underutilized source of renewable energy. The microbial valorization of lignin is challenging due to the toxicity of its degradation intermediates, particularly formaldehyde. In this study, we engineeredMethylobacterium extorquensPA1 to metabolize lignin-derived methoxylated aromatics, vanillate (VA) and protocatechuate (PCA), by introducing thevanandpcagene clusters. Compared toPseudomonas putida,M. extorquensPA1 exhibited better formaldehyde detoxification, enabling robust growth on VA without accumulation of formaldehyde. Genetic analyses confirmed that formaldehyde oxidation and stress response systems, rather than C 1 assimilation, were important for VA metabolism. Additionally, VA and PCA were found to disrupt membrane potential, contributing to their inherent toxicity. Our findings establishM. extorquensPA1 as a promising chassis for lignin valorization and provide a framework for engineering formaldehyde-resistant microbial platforms. IMPORTANCE In developing biotechnological solutions for a circular economy, it is critical to valorize all parts of renewable resources, such as lignocellulose from vegetative components of agricultural crops and from bioenergy feedstocks. After chemical breakdown, the aromatics arising from lignin present significant challenges for use due to their toxicity. Here, we address one component of this challenge—the methoxy groups that get released as formaldehyde—and show that existing biotechnological platform organisms with strong formaldehyde metabolism, such asMethylobacterium extorquens, can be transformed into highly capable utilizers of methoxylated aromatics.

Microbiology↗

Bypassing Fast Time Scales of the Hodgkin-Huxley Neuron Model via a Thresholded Hard Reset

We propose a modification to the Hodgkin-Huxley model to reduce the numerical stiffness of the equations by introducing an explicit voltage threshold. When this threshold is crossed, the voltage and the gating variables are reset to constant values. It is found that, for all of the current stimuli considered, the proposed model accurately reproduces the behavior of the baseline Hodgkin-Huxley model while bypassing the fast dynamics of spiking. Specifically, the model accurately reproduces the spike times and, between spikes, the time courses of the membrane potential and gating variables.

97 MATHEMATICS AND COMPUTING↗

Deletion of Carboxypeptidase E in β-Cells Disrupts Proinsulin Processing but Does Not Lead to Spontaneous Development of Diabetes in Mice

Carboxypeptidase E (CPE) facilitates the conversion of prohormones into mature hormones and is highly expressed in multiple neuroendocrine tissues. Carriers of CPE mutations have elevated plasma proinsulin and develop severe obesity and hyperglycemia. We aimed to determine whether loss of Cpe in pancreatic β-cells disrupts proinsulin processing and accelerates development of diabetes and obesity in mice. Pancreatic β-cell–specific Cpe knockout mice (βCpeKO; Cpe fl/fl x Ins1 Cre/+ ) lack mature insulin granules and have elevated proinsulin in plasma; however, glucose-and KCl-stimulated insulin secretion in βCpeKO islets remained intact. High-fat diet–fed βCpeKO mice showed weight gain and glucose tolerance comparable with those of Wt littermates. Notably, β-cell area was increased in chow-fed βCpeKO mice and β-cell replication was elevated in βCpeKO islets. Transcriptomic analysis of βCpeKO β-cells revealed elevated glycolysis and Hif1α-target gene expression. On high glucose challenge, β-cells from βCpeKO mice showed reduced mitochondrial membrane potential, increased reactive oxygen species, reduced MafA, and elevated Aldh1a3 transcript levels. Following multiple low-dose streptozotocin injections, βCpeKO mice had accelerated development of hyperglycemia with reduced β-cell insulin and Glut2 expression. In conclusion, these findings suggest that Cpe and proper proinsulin processing are critical in maintaining β-cell function during the development of hyperglycemia.

59 BASIC BIOLOGICAL SCIENCES↗

Scaling and Benchmarking an Evolutionary Algorithm for Constructing Biophysical Neuronal Models

Single neuron models are fundamental for computational modeling of the brain's neuronal networks, and understanding how ion channel dynamics mediate neural function. A challenge in defining such models is determining biophysically realistic channel distributions. Here, we present an efficient, highly parallel evolutionary algorithm for developing such models, named NeuroGPU-EA. NeuroGPU-EA uses CPUs and GPUs concurrently to simulate and evaluate neuron membrane potentials with respect to multiple stimuli. We demonstrate a logarithmic cost for scaling the stimuli used in the fitting procedure. NeuroGPU-EA outperforms the typically used CPU based evolutionary algorithm by a factor of 10 on a series of scaling benchmarks. We report observed performance bottlenecks and propose mitigation strategies. Finally, we also discuss the potential of this method for efficient simulation and evaluation of electrophysiological waveforms.

59 BASIC BIOLOGICAL SCIENCES↗

Energy coupling and stoichiometry of Zn 2+ /H + antiport by the prokaryotic cation diffusion facilitator YiiP

YiiP from Shewanella oneidensis is a prokaryotic Zn 2+ /H + antiporter that serves as a model for the Cation Diffusion Facilitator (CDF) superfamily, members of which are generally responsible for homeostasis of transition metal ions. Previous studies of YiiP as well as related CDF transporters have established a homodimeric architecture and the presence of three distinct Zn 2+ binding sites named A, B, and C. In this study, we use cryo-EM, microscale thermophoresis and molecular dynamics simulations to address the structural and functional roles of individual sites as well as the interplay between Zn 2+ binding and protonation. Structural studies indicate that site C in the cytoplasmic domain is primarily responsible for stabilizing the dimer and that site B at the cytoplasmic membrane surface controls the structural transition from an inward facing conformation to an occluded conformation. Binding data show that intramembrane site A, which is directly responsible for transport, has a dramatic pH dependence consistent with coupling to the proton motive force. A comprehensive thermodynamic model encompassing Zn 2+ binding and protonation states of individual residues indicates a transport stoichiometry of 1 Zn 2+ to 2–3 H + depending on the external pH. This stoichiometry would be favorable in a physiological context, allowing the cell to use the proton gradient as well as the membrane potential to drive the export of Zn 2+ .

59 BASIC BIOLOGICAL SCIENCES↗

Evidence for long-term potentiation in phospholipid membranes

Biological supramolecular assemblies, such as phospholipid bilayer membranes, have been used to demonstrate signal processing via short-term synaptic plasticity (STP) in the form of paired pulse facilitation and depression, emulating the brain’s efficiency and flexible cognitive capabilities. However, STP memory in lipid bilayers is volatile and cannot be stored or accessed over relevant periods of time, a key requirement for learning. Using droplet interface bilayers (DIBs) composed of lipids, water and hexadecane, and an electrical stimulation training protocol featuring repetitive sinusoidal voltage cycling, we show that DIBs displaying memcapacitive properties can also exhibit persistent synaptic plasticity in the form of long-term potentiation (LTP) associated with capacitive energy storage in the phospholipid bilayer. The time scales for the physical changes associated with the LTP range between minutes and hours, and are substantially longer than previous STP studies, where stored energy dissipated after only a few seconds. STP behavior is the result of reversible changes in bilayer area and thickness. On the other hand, LTP is the result of additional molecular and structural changes to the zwitterionic lipid headgroups and the dielectric properties of the lipid bilayer that result from the buildup of an increasingly asymmetric charge distribution at the bilayer interfaces.

59 BASIC BIOLOGICAL SCIENCES↗

Ion Separations Based on Spontaneously Arising Streaming Potentials in Rotating Isoporous Membranes

Highly selective ion separations are vital for producing pure salts, and membrane-based separations are promising alternatives to conventional ion-separation techniques. Our previous work demonstrated that simple pressure-driven flow through negatively charged isoporous membranes can separate Li+ and K+ with selectivities as high as 70 in dilute solutions. The separation mechanism relies on spontaneously arising streaming potentials that induce electromigration, which opposes advection and separates cations based on differences in their electrophoretic mobilities. Although the separation technique is simple, this work shows that high selectivities are possible only with careful consideration of experimental conditions including transmembrane pressure, solution ionic strength, the K+/Li+ ratio in the feed, and the extent of concentration polarization. Separations conducted with a rotating membrane show Li+/K+ selectivities as high as 150 with a 1000 rpm membrane rotation rate, but the selectivity decreases to 1.3 at 95 rpm. These results demonstrate the benefits and necessity of quantitative control of concentration polarization in highly selective separations. Increases in solution ionic strength or the K+/Li+ feed ratio can also decrease selectivities more than an order of magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Membranes for Lithium Recovery From Conventional and Unconventional Sources

Lithium has been deemed a critical mineral of national importance that finds uses in a wide range of applications, and its demand has been rising significantly in recent years. The urgency of meeting this demand requires lithium extraction from various aqueous sources such as continental brines, geothermal brines, seawater, produced water, and battery waste. While direct lithium extraction (DLE) technologies such as adsorption, ion exchange, and solvent extraction have emerged as possible solutions, membrane technologies are also being investigated for various sources and at different stages of the recovery process. Here, we analyze the application of membranes for pretreatment of lithium source waters, bring management, lithium/magnesium separation, lithium/sodium separation, and lithium hydroxide conversion, and evaluate performance metrics for critical lithium separations from the literature. We explore the potential of membranes at every stage of the recovery process and describe their current status and future prospects. We describe hypothetical process trains with integrated membrane technologies for each source type and address their feasibility and challenges. The potential energy and water impacts of membrane-integrated and conventional DLE processes are also critically considered alongside performance and selectivity metrics, and this is illustrated using examples and calculated from published technical reports. This paper thus provides a comprehensive overview of the application of membranes along every stage of the lithium recovery process, emphasizing the versatility and potential of membrane technologies for critical mineral recovery.

36 MATERIALS SCIENCE↗

A Pro‐Angiogenic Immunoprotective Membrane for Cell Therapies

Abstract Immunoisolation strategies that rely on porous membranes play an important role in cell transplantation therapies to protect cells from the host's immune system. These membranes must possess immunoprotective properties while facilitating the transport of nutrients and cell products to maintain the functional integrity of encapsulated cells. An easy and scalable process is described to fabricate a dual function porous polymeric membrane that shields cells against immune cell attack and promotes vascularization to address the nutritional and oxygen requirements of transplanted cells. The fabrication process results in a membrane cross‐section with a gradient of nanopores to micropores that support cell immunoisolation and interfacial vascularization requirements, respectively. The membranes demonstrate excellent cell compatibility and effectively prevent T cell transmigration without compromising glucose diffusion and oxygen permeability. In a murine subcutaneous implantation model, membranes are stable for 60 days and exhibit significantly reduced fibrous capsules, with enhanced vascularization near the membrane. These porous polymeric membranes can potentially be used as pro‐angiogenic immunoprotective membranes for cell transplantation applications where maximizing cell viability and function is of critical importance.

Engineering↗

Calcium Sulfate and Calcium Carbonate Scaling of Thin-Film Composite Polyamide Reverse Osmosis Membranes with Surface-Tethered Polyacrylic Acid Chains

The gypsum and calcite scaling propensities of the thin-film composite polyamide (PA-TFC) reverse osmosis (RO) membrane, modified with a tethered surface layer of polyacrylic acid (PAA) chains, was evaluated and compared to the scaling of selected commercial RO membranes. The tethered PAA layer was synthesized onto a commercial polyamide membrane (i.e., base-PA) via atmospheric pressure plasma-induced graft polymerization (APPIGP). The PAA nano-structured (SNS) base-PA membrane (SNS-PAA-PA) was scaled to a lesser degree, as quantified by a lower permeate flux decline and surface imaging, relative to the tested commercial membranes (Dow SW30, Toray SWRO, and BWRO). The cleaning of gypsum-scaled membranes with D.I. water flushing achieved 100% water permeability recovery for both the SNS-PAA-PA and Dow SW30 membranes, relative to 92–98% permeability restoration for the Toray membranes. The calcium carbonate scaling of SNS-PAA-PA membranes was also lower relative to the commercial membranes, but permeability recovery after D.I. water cleaning was somewhat lower (94%) but consistent with the level of surface scale coverage. In contrast, the calcite and gypsum-scaled membrane areas of the commercial membranes post-cleaning were significantly higher than for the SNS-PAA-PA membrane but with 100% permeability recovery, suggesting the potential for membrane damage when mineral scaling is severe.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of Carbon Molecular Sieves Hollow Fiber Membranes based on Polybenzimidazole Doped with Polyprotic Acids with Superior H 2 /CO 2 Separation Properties (Final Report)

The goal of this project was to develop a highly efficient membrane-based process to capture CO 2 from coal-derived syngas with 95% CO 2 purity, achieving the cost of electricity (COE) 30% below the baseline capture approaches (i.e., Selexol process) when coupled with the advancement in other areas of the power generation facility. Our core approach is based on high-permeance hollow fiber membranes (HFMs) with superior H 2 /CO 2 separation properties at the syngas process conditions, which can then be further utilized to design membrane reactors for process intensification of H 2 production and purification and CO 2 capture. Three organizations with complementary skills collaborated to achieve the goal, including the University at Buffalo (UB), Los Alamos National Laboratory (LANL), and Trimeric Corporation (Trimeric). We formulated logical steps to bring the membrane technology from Technology Readiness Level (TRL) 3 (Experimental proof of concept) to TRL 4 (Laboratory scale validation in relevant environment). During the budget period (BP) 1, we screened various polymeric materials and identified polybenzimidazole doped with inorganic polyprotic acids as the desirable platform. The acid doping increases the H 2 /CO 2 selectivity, and the sequential carbonization increases H 2 permeability while retaining the high selectivity. By manipulating the acid type and doping level and the carbonization temperature, we demonstrated advanced carbon molecular sieving (CMS) materials with H 2 permeability of above 200 Barrer (1 Barrer = 10 -10 cm 3 (STP) cm cm -2 s -2 cmHg -1 ) and H 2 /CO 2 selectivity of above 40 at 200-300°C with simulated syngas containing CO and water vapor. For example, the PBI-(H 3 PO 4 ) 0.11 carbonized at 700 °C exhibits H 2 permeability of 200 Barrer and H 2 /CO 2 selectivity of 60 at 200 °C, and H 2 permeability of 240 Barrer and H 2 /CO 2 selectivity of 54 at 225 °C, which meets the targeted properties and surpasses Robeson’s upper bound. During the BP2, we focused on the conversion of the advanced CMS materials to stable HFMs. Membranes with H 2 permeance of 1,090 GPU (1 GPU = 10 -6 cm 3 (STP) cm -2 s -2 cmHg -1 ) and H 2 /CO 2 selectivity of 57 at 300 °C were successfully fabricated. The effects of temperature, gas compositions, pressure, and time on the separation properties were systematically investigated. Pencil modules were continuously evaluated for 219 h (dry pure gas) and 669 h (dry simulated syngas) and showed initial decline in permeance and increased H 2 /CO 2 selectivity over time, ultimately achieving a steady state stable value, indicating that the ageing phenomena in the nanoporous structures of the membranes during the test. Membrane reactors were fabricated based on the CMS membranes and evaluated for water-gas shift (WGS) reaction. The use of membranes slightly improves the conversion of the CO. However, parametric tests of the membrane reactors at various temperatures and flow rates need to be conducted, as well as the membranes with improved separation performance. We performed a sensitivity analysis on the impact of H 2 /CO 2 selectivity on the COE based on a hybrid process of a membrane unit and cryogenic unit developed by Membrane Technology and Research, Inc. (MTR). Three H 2 /CO 2 selectivity (40, 60, and 15) cases were developed and compared with the baseline capture process (Case B5B) provided by the DOE report. Increasing the membrane H 2 /CO 2 selectivity reduces COE, but the rate of the decrease of COE also diminishes. The COE values for H 2 /CO 2 selectivities of 40 and 60 were nearly the same. As the H 2 /CO 2 selectivity increases, the inert recycling decreases, leading to smaller equipment, less auxiliary power requirements, and less heating, cooling, and refrigeration duty. The refrigeration system used to liquefy the CO 2 is the most expensive piece of equipment and consumes the most electricity within the CO 2 capture process. Increasing the CO 2 concentration in the recycle stream would improve the economics of the process by reducing the refrigeration duty requirement of the unit and also allow for higher liquefaction temperatures. The high H 2 -selective membrane developed by our team may be applicable in other separation processes where lower pressure H 2 retains value. Typically, hydrogen retains its pressure when it is separated from syngas components. Residual components may be used as low-quality fuel and then vented to the atmosphere. Applications might include control of H 2 /CO ratios or mitigation of the water gas shift reaction by CO 2 recycling to the feed of a gasifier or steam methane reformer. To summarize, we have developed industrial HFMs with the best H 2 /CO 2 separation performance reported in the literature. The membranes demonstrate stability with simulated syngas and show great potential for membrane reactors for WGS reactions, lowering the cost of blue H 2 production.

20 FOSSIL-FUELED POWER PLANTS↗

Quaternized poly(arylene ether benzonitrile) membranes for vanadium redox flow batteries

In this work, a series of quaternized poly (arylene ether benzonitrile)s (QAxCN; x = percent quaternization) were synthesized for the vanadium redox flow battery (VRFB) application. The presence of highly polar nitrile groups incorporated in the quaternized polymer chain provides desirable polymer properties (low water uptake, high dimensional stability and low VO 2+ permeability) without compromising ion conductivity. VRFB single cells using the aforementioned membranes showed excellent performance. As a result, the cell using a 15 μm-thick QA50CN membrane exhibited superior coulombic efficiency (98.2%) and comparable energy efficiency (86.5%) at 40 mA cm -2 compared to the single cell using Nafion-212 with a coulombic efficiency of 95.5% and an energy efficiency of 86.2%. The cell cycling stability and membrane oxidative stability data confirmed that the QAxCN membranes have the potential to replace Nafion membranes for energy storage in VRFB applications.

25 ENERGY STORAGE↗

Genome-wide characterization of the soybean DOMAIN OF UNKNOWN FUNCTION 679 membrane protein gene family highlights their potential involvement in growth and stress response

The DMP (DUF679 membrane proteins) family is a plant-specific gene family that encodes membrane proteins. The DMP family genes are suggested to be involved in various programmed cell death processes and gamete fusion during double fertilization in Arabidopsis. However, their functional relevance in other crops remains unknown. This study identified 14 genes from the DMP family in soybean (Glycine max) and characterized their physiochemical properties, subcellular location, gene structure, and promoter regions using bioinformatics tools. Additionally, their tissue-specific and stress-responsive expressions were analyzed using publicly available transcriptome data. Phylogenetic analysis of 198 DMPs from monocots and dicots revealed six clades, with clade-I encoding senescence-related AtDMP1/2 orthologues and clade-II including pollen-specific AtDMP8/9 orthologues. The largest clade, clade-III, predominantly included monocot DMPs, while monocot- and dicot-specific DMPs were assembled in clade-IV and clade-VI, respectively. Evolutionary analysis suggests that soybean GmDMPs underwent purifying selection during evolution. Using 68 transcriptome datasets, expression profiling revealed expression in diverse tissues and distinct responses to abiotic and biotic stresses. The genes Glyma.09G237500 and Glyma.18G098300 showed pistil-abundant expression by qPCR, suggesting they could be potential targets for female organ-mediated haploid induction. Furthermore, cis-acting regulatory elements primarily related to stress-, hormone-, and light-induced pathways regulate GmDMPs, which is consistent with their divergent expression and suggests involvement in growth and stress responses. Overall, our study provides a comprehensive report on the soybean GmDMP family and a framework for further biological functional analysis of DMP genes in soybean or other crops.

59 BASIC BIOLOGICAL SCIENCES↗

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↗