Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “well casing”

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

Mechanistic modeling of CO 2 well leakage in a generic abandoned well through a bridge plug cement-casing gap

Both known and unmapped plugged and abandoned wells are potential leakage pathways for CO 2 from geologic carbon sequestration (GCS) sites. Although many abandoned wells have cement bridge plugs installed to prevent leakage, the seal between the cement and the inner casing wall is subject to failure. In this study, we carried out detailed T2Well simulations of cases of sudden non-Darcy flow of CO 2 and brine leakage up the gap between a cement plug and the inner steel casing wall that becomes a fully connected flow path during the post-injection period. The goal of our study was two-fold: (1) to understand the dynamics, rates, and the characteristic temporal signals associated with the onset of leakage through various gap-aperture sizes, and (2) to suggest potential monitoring strategies based on the findings. Simulation results show that the leakage of CO 2 and brine upward is transient with interesting phase interference behavior. Time-dependent oscillatory flows with varying pressure, temperature, and flow rates of CO 2 and brine show strong dependence on gap aperture. Phase-change and decompression lead to very low temperatures at the top of the well for gap apertures larger than 4 mm suggesting that remote thermal monitoring at the ground surface may be an effective way of monitoring even if well locations are not known a priori. Pressure in the well is also indicative of CO 2 leakage. The temporal patterns of changing temperature and pressure may be useful diagnostic signals for leakage detection. Finally, these transient leakage signals may provide information on the cause of leakage and/or characteristics of the flow path that could inform effective remediation design and execution approaches.

58 GEOSCIENCES↗

Hot Springs and Geysers: Exploring Historical and Modern Impacts of Geothermal Energy Production on Associated Natural Surface Systems and Standardizing Management Practices

Surface thermal features, most notably hot springs and geysers are increasingly being recognized for their importance to ecosystems, indigenous cultures, and in some cases agriculture, recreation, and tourism. Geothermal project development poses a potential risk to these natural features but current regulatory requirements for assessing and managing these risks during exploration, permitting and monitoring are somewhat inconsistent and unpredictable across different geothermal fields. This has resulted in uncertainty and increases in exploration risk for geothermal energy developers that have led to costly project delays, cancellations, or hesitation to commit. Varying regulatory requirements may also influence public perception, fostering confusion, distrust and ultimately opposition to geothermal projects, further contributing to project delays or cancellations. At a time when there is an increasing urgency for reliable baseload clean energy, geothermal is a net-zero, renewable solution that additionally provides access to more equitable and environmentally just clean power. Continued integration of geothermal energy into the national energy roadmap can be facilitated through consistent and predictable permitting, providing regulators the framework they need, developers a clear path forward, and transparency that the public deserves. This project, currently in its beginning phases, seeks to address this important issue by providing a technical basis from which to build a preliminary protocol for assessing and managing potential impacts from new or existing geothermal energy projects to surface thermal features and their associated ecosystems. Development of this preliminary protocol will be informed by (1) a literature review of well-documented case studies in the western U.S. and New Zealand to understand the range of conditions that exemplify geothermal-surface thermal systems; (2) development of generic illustrative conceptual-numerical models to quantify, understand, and predict the first-order controls (e.g., pressure and permeability) on surface flows; and (3) additional independent and scientifically rigorous evaluations of geothermal-surface thermal system case studies from the Basin and Range Province that incorporate publicly available data as well as data provided by industry through data-sharing agreements. Learning from the successes of the process used to develop the Induced Seismicity Management Protocol (ISMP), we ultimately aim to use these initial efforts as a springboard for establishing a surface thermal feature management working group that will work collaboratively to finalize the protocol as well as co-create recommended best practices for implementation. We envision that the working group will primarily be composed of representatives from regulatory entities, government agencies, Tribes, academia, national laboratories, and industry, and will include early and regular engagement with community organizations and environmental groups. This will help ensure broad acceptance and implementation of the protocol, which will facilitate a more consistent, predictable, and standardized regulatory process, and help to ensure that geothermal energy continues to provide a reliable source of clean energy, and a pathway to achieving greater energy equity in the U.S.

Best Practices↗

The use of idealised experiments in testing a new convective parametrization: Performance of CoMorph-A

CoMorph is a new mass-flux convection parametrization under development at the Met Office designed for use within the Unified Model and its successor model, LFRic. Use of a three-dimensional idealised model enables controlled tests of the performance of the scheme across different regimes. This includes the interaction between the physical parametrizations and the resolved dynamics, allowing study of the emergent organisation of convection on the resolved scale. A selection of well-known cases is revisited here, with the purpose of documenting the extent to which CoMorph captures a range of important, but challenging, behaviour such as the diurnal cycle and sensitivity to tropospheric moisture. Simulations using CoMorph-A, a new physics package, that has been demonstrated to perform well at numerical weather prediction (NWP) and climate scales, are compared against the current global atmosphere configuration and high-resolution results. In addition to an entirely new convection scheme, the package of changes includes significant changes to the cloud, microphysics, and boundary-layer parametrizations. Recognising that CoMorph-A is the first version of a scheme that will continue to be substantially developed and to obtain good performance, compromises in tuning have had to be made. These idealised tests therefore show what works well in this configuration, and what areas will require further work. As such, it is quite a demanding testbed and could be viewed as some of the equipment required for a “convective playground”.

54 ENVIRONMENTAL SCIENCES↗

Effect of Split-Injection Strategies on Engine Performance and Emissions under Cold-Start Operation

The recently concluded partnership for advancing combustion engines (PACE) was a US Department of Energy consortium involving multiple national laboratories focused on addressing key efficiency and emission barriers in light-duty engines. Generation of detailed experimental data and modeling capabilities to understand and predict cold-start behavior was a major pillar in this program. Cold-start, as defined by the time between first engine crank and three-way catalyst light-off, is responsible for a large percentage of NOx, unburned hydrocarbon, and particulate matter emissions in light-duty engines. Minimizing emissions during cold-start is a trade-off between achieving faster three-way catalyst light-off, and engine out emissions during that period. In this study, engine performance, emissions, and catalyst warmup potential were monitored while the engine was operated using a single direct injection (baseline case) as well as a two-way-equal-split direct injection strategy. These injection strategies were analyzed at a range of cold-start-operation relevant retarded spark timings of up to 25 degrees after top dead center of firing (dATDCf). A stoichiometric 2-bar NIMEP steady-state condition was used for all cases to simulate cold-start operation. Significant improvement in engine stability was observed with the two-way-split injection strategy at the retarded spark timings allowing for up to 2.5x increase in exhaust heat rate when engine operation is stability constrained. Similar fuel-loss-to-oil trends with exhaust heat rate were observed for both single and two-way-split injection strategies. However, the two-way split injection was observed to produce higher NOx emissions per unit exhaust heat rate. A single data point run with three-way-split direct injection at a very retarded spark-timing of 30 dATDCf pointed to further improvements in engine stability and reduction in fuel-loss-to-oil as compared to single injection strategy. Engine stability decreased as spark timing was initially retarded with a single injection but was observed to plateau and stabilize beyond spark timing of 10 dATDCf. Finally, for the two-way-split-injection strategy, retarding the start of injection (SOI) timing of the second injection led to a decrease in engine stability as well as an increase in soot emissions.

42 ENGINEERING↗

Distributed electric field sensing using fibre optics in borehole environments

In the past decade, rapid advances in distributed optical fibre sensing technologies have made it possible to record various geophysical data (e.g. strain, temperature and pressure) continuously in both time and space along the fibre, providing an unprecedented quantity and spatial density of data compared to traditional geophysical measurements as well as reducing data acquisition cost. To date, no distributed fibre-based electromagnetic field sensing system has been implemented although electromagnetic sensing could have a broad range of applications to geophysical imaging and monitoring in borehole environments. The goal of this paper is to provide a theoretical feasibility study regarding the design and use of an electromagnetic sensing optical fibre for geophysical applications. First, we present the sensitivity analysis of a ‘hypothetical’ optical fibre coated with polyvinylidene fluoride, a polymer that provides relatively high piezoelectric properties, yet unlike ceramics, is flexible. Next, using a two-dimensional electromagnetic modelling algorithm, we simulate the earth electric-field-to-fibre-strain transfer function and estimate the theoretical sensitivity of the optical fibre to electric fields. Given the state-of-the-art distributed acoustic sensing strain sensitivities in the picometres strain range, our numerical modelling analysis suggests that a perfectly coupled polyvinylidene fluoride–coated optical fibre can measure electric field values in the mV/m to V/m amplitude range. We then apply a cylindrically symmetric modelling algorithm to simulate numerical models demonstrating the applicability of such a fibre in an oilfield environment. Scenarios investigated employ an electric field source and suggest that the measurements can be used to distinguish the oil versus water ratio with a fibre mounted inside a producing steel cased oil well as well as distinguishing between brine and hydrocarbon filled reservoir zones with a fibre located outside of the casing.

58 GEOSCIENCES↗

Tritium Breeding Ratio Evaluation of Solid Breeder Concepts for the FESS-FNSF

This paper presents a parametric study of the Fusion Energy System Studies-Fusion Nuclear Science Facility’s (FNSF’s) tritium breeding performance for several solid breeder concepts, neutron multiplying materials, and blanket materials, assuming volume fractions based on the most recent FNSF design as a realistically representative fusion facility. In this study, we initially surveyed the tritium breeding ratio (TBR) of several solid breeder concepts by employing a simplified but efficient one-dimensional (1-D) infinite cylinder reduced-order model (ROM). Parametric studies were performed with the ROMs for the full range of breeder-to-multiplier ratios to identify the optimum mixture compositions for each breeder type that would lead to a maximum TBR. These optimized breeder-multiplier combinations were then homogenized with FNSF blanket component materials to estimate their impacts on the TBR. Subsequently, as a validation step for the optimal designs, TBR calculations were performed using a more realistic modified 1-D ROM with inner and outer breeding regions, as well as with a fully detailed 22.5-deg three-dimensional (3-D) sector of the FNSF to assess the impact of geometry details on the TBR. The differences between the two 1-D models were negligible, while the ROMs were able to correctly predict trends and identify the maximum and minimum TBR cases, as well as show consistent biases relative to the results produced by the full 3-D, 22.5-deg sector for specific breeder/multiplier combinations. Solid breeder concepts such as Li 2 O, Li 4 SiO 4 , and Li 8 ZrO 6 outperformed all others in this study in terms of TBR performance when combined with all the neutron multiplier materials selected. Here, an underlying goal of this study was to develop and improve rapid and reliable ROMs to aid designers during parametric optimizations of highly complex and computationally expensive fusion models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Identification of potential CO 2 leakage pathways and mechanisms in oil reservoirs using fault tree analysis

Geological storage of CO 2 technologies has become an important and effective way to reduce the greenhouse gas emissions, especially when it is combined with CO 2 enhanced oil recovery (EOR), which can not only trap CO 2 but also enhance oil recovery. However, the risk of CO 2 leakage has always been a prominent issue. In this paper, the mechanisms and pathways of CO 2 leakage during geological storage in oil reservoirs were analyzed using fault tree analysis (FTR). Besides, monitoring technologies were discussed and deployed in a CO 2 EOR demonstration project. The analysis results showed that the sealing failures of oil producer and CO2 injector wells, like well cement failure and casing failure, are the main reasons for the CO 2 leakage, which has been observed in the oil field monitoring project. The monitoring results indicated that there is no large-scale CO 2 leakage, while relatively high and abnormal CO 2 concentration in soil gas near some wellbores are observed, which indicates there is some leakage of CO 2 through incomplete cement ring and well casing string. FTR results provide guidelines for monitoring and preventing of CO 2 leakage during geological storage in oil reservoirs. Finally, the near-surface monitoring methods, especially the soil gas monitoring technologies, can effectively detect the leakage of CO 2 , and are a proper method for CO 2 leakage monitoring.

58 GEOSCIENCES↗

New Covariant Density Functionals of Nuclear Matter for Compact Star Simulations

We generate three families of extended covariant density functionals of nuclear matter that have varying slope of symmetry energy and skewness at nuclear saturation density, but otherwise share the same basic parameters (symmetry energy, compressibility, saturation parameters, etc.) with the standard DDME2, DD2, and MPE functionals. Tables of the parameters of these new density functionals are given, which can be straightforwardly used in DDME2, DD2, and MPE parameterization-based codes. Furthermore, we provide tables of a large number of equations of state (81 for each family) that can be used in astrophysical simulations to assess the impact of variations of the not-well-known slope of symmetry energy and skewness of nuclear systems on the astrophysics of compact objects. We also provide tables of computed integral parameters (mass, radius, and tidal deformability) that can be used, e.g., for modeling gravitational waveforms. Finally, for the extended DDME2-based parameterization, we implement a first-order phase transition to quark matter to obtain a family of equations of state that accommodates a phase transition to quark matter. Analogous tables of the equations of state and integral parameters are provided for this case as well.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of the economic implications of varied pressure drawdown strategies generated using a real-time, rapid predictive, multi-fidelity model for unconventional oil and gas wells

Experience has suggested that pressure maintenance in hydraulically fractured reservoirs via lower, more sustained production drawdowns may offer improved cumulative recovery and overall resource extraction efficiency compared to more rapid drawdown approaches aimed at generating high initial production. However, given the inherent variability of oil and natural gas markets, operators pursue production strategies that maximize profitability over resource extraction efficiency. This study focuses on evaluating the implications of contrasting pressure drawdown strategies on the long-term production and resulting economics for a real, producing unconventional gas well in the Marcellus Shale of the Appalachian Basin using a techno-economic analysis approach. Our research combines elements of well-specific horizontal well design, production forecasting, equipment sizing and capital cost estimation, operating cost estimation, and revenue and tax calculations. Gas production forecast outlook scenarios were generated under varying pressure drawdowns using two approaches: 1) a novel physics-informed machine learning workflow and 2) traditional reservoir simulation. A discounted cash flow model was used to evaluate the resulting economic implications for each drawdown scenario—generating output for exploring the coupled effect of factors like the timing and volume of gas production, prevailing economic and market conditions for natural gas, and overall estimated ultimate recovery on profitability metrics such as internal rate of return and net present value. Results show that there is potential to maximize the cumulative gas produced in the specific case study well by employing a lower pressure drawdown. Conversely, the greatest profitability is achieved using rapid drawdown as signified by a small, specific subset of our outlook scenarios. On an averaging basis, we find that the combinations of highest cumulative producing and most profitable scenarios occur under lower drawdowns with long (>40 years) producing timeframes, but require higher relative gas price and lower discounting considerations. Further, the machine learning predictive outlooking capability proved effective for enabling rapid generation of a multitude of scenario forecasts. As a result, a variety of prominent example cases could be generated to strike the balance of greater productivity and economic return given their associated producing features and economic conditions when compared to similar producing scenarios—critical insight that offers improved decision support for unconventional oil and gas operations.

42 ENGINEERING↗

Gargantuan Hail in Argentina

On 8 February 2018, a supercell storm produced gargantuan (>15 cm or >6 inches in maximum dimension) hail as it moved over the heavily populated city of Villa Carlos Paz in Cordoba Province, Argentina, South America. Observations of gargantuan hail are quite rare, but the large population density here yielded numerous witnesses and social media pictures and videos from this event that document multiple large hailstones. The storm was also sampled by the newly installed operational polarimetric C-band radar in Cordoba. During the RELAMPAGO campaign, the authors interviewed local residents about their accounts of the storm, and uncovered additional social media video and photographs revealing extremely large hail at multiple locations in town. This article documents the case, including the meteorological conditions supporting the storm (with the aid of a high-resolution WRF simulation), the storm's observed radar signatures, and three noteworthy hailstones observed by residents. These hailstones include a freezer-preserved 4.48-inch (11.38-cm) maximum dimension stone that was scanned with a 3D infrared laser scanner, a 7.1-inch (18-cm) maximum dimension stone, and a hailstone photogrammetrically estimated to be between 7.4 and 9.3 inches (18.8-23.7-cm) in maximum dimension, which is close to or exceeds the world record for maximum dimension. Such a well-observed case is an important step forward in understanding environments and storms that produce gargantuan hail, and ultimately how to anticipate and detect such extreme events.

Kumjian, Matthew R.↗

Unraveling TeV halos with the Cherenkov Telescope Array

Pulsars are observed to emit bright and spatially extended gamma-ray emission at multi-TeV energies. These so-called "TeV halos" are now understood to be a nearly universal feature of middle-aged pulsars. However, many of the key physical processes that govern these systems, particularly those affecting particle diffusion, remain poorly constrained. We aim to evaluate the ability of the Cherenkov Telescope Array (CTA) to probe the physical properties of TeV halos, with a focus on the nearby and well-studied case of the Geminga pulsar. We simulate gamma-ray emission from various TeV halo models, incorporating different assumptions for the injected electron spectrum, spin-down evolution, and energy-dependent diffusion. These models are then used to forecast CTA's sensitivity to spectral and spatial differences, based on realistic mock observations and instrument response simulations. We find that CTA will be able to distinguish between a wide range of TeV halo models that are currently consistent with existing data. In particular, CTA observations can constrain the normalization, energy dependence, and spatial extent of the diffusion coefficient surrounding Geminga, as well as the spectral shape of the injected electron population.

79 ASTRONOMY AND ASTROPHYSICS↗

Wabash CarbonSAFE. Geologic Analysis of the Potosi Dolomite Reservoir Interval and Potential Confining Units, at the Wabash #1 Well, Vigo Co., Indiana, Subtask 7.2, Technical Report

The Wabash CarbonSAFE project drilled the Wabash #1 stratigraphic test well (ID# 168045) at the Wabash Valley Resources (WVR) IGCC facility in Vigo County, Indiana, to evaluate the feasibility of commercial-scale CO 2 storage near the site. The Mt. Simon Sandstone was the initial target for storage evaluation and was found to have generally poor reservoir qualities in the Wabash #1 well. Focus was then placed on the secondary target, the Potosi Dolomite. For storage and containment evaluation an extensive suite of geophysical logs was collected throughout the borehole, and whole cores were obtained from the Mt. Simon Sandstone, Eau Claire Formation, and Maquoketa Group. The Maquoketa Group is identified as a regionally extensive secondary sealing interval above the Mt. Simon Sandstone and as a primary seal for the Potosi Dolomite storage complex. Due to concerns over potential lost circulation zones in the Potosi Dolomite during drilling, no core samples or Formation Micro Imager (FMI) logs were acquired in the Potosi Dolomite. However, cased hole well testing was performed over a 20 feet (6 m) interval within the Potosi Dolomite during which a fluid swab sample was obtained. Approximately 35 miles of 2D seismic information was acquired in the project area near Wabash #1 to aid in reservoir and caprock characterization. The highly porous and permeable zones in the Potosi Dolomite and lower Oneota Dolomite observed within the Wabash #1 well are, for the purpose of this report, considered the same reservoir package and referred to as the Potosi reservoir interval. This report presents a regional depositional context for the Potosi reservoir interval, the overlying confining zones in the Knox Group, and Maquoketa Group primary seal. This context is developed through detailed information obtained from the Wabash #1 well such as lithologic data collected from cuttings and core, geophysical logging, geomechanical analysis of core samples, and well testing and fluid sampling within the Potosi reservoir interval. Information discussed in this report was compiled as a basis for the development of static and dynamic geologic models of the Potosi Dolomite storage complex. The static and dynamic modeling of CO 2 injection in the Potosi reservoir interval are discussed in a separate report (Dessenberger et al., 2022) under the Wabash CarbonSAFE project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Cements and a Modeling Tool to Calculate their Viability under Various Exploitation Conditions of HT RTES Systems

Primary cementing is the most important operation performed on a subterranean reservoir. Cement, placed in the annulus between the casing and the formations serves as a hydraulic seal preventing fluids and gas migrations, protecting steel casing from corrosion and supporting the well structure. Poor cementing jobs can be directly responsible for wells not reaching their full capacity, casing corrosion, compromised well integrity and, in the worst-case scenarios, well collapse.

15 GEOTHERMAL ENERGY↗

Bare mass effects on the reheating process after inflation

We consider the effects of a bare mass term for the inflaton, when the inflationary potential takes the form V(Φ)=λΦ k about its minimum with k ≥ 4. We concentrate on k =4, but discuss general cases as well. Further, we assume $λΦ$$^{2}_{end}$ >> $m$$^{2}_{Φ}$, where Φ end is the inflaton field value when the inflationary expansion ends. We show that the presence of a mass term (which may be present due to radiative corrections or supersymmetry breaking) can significantly alter the reheating process, as the equation of state of the inflaton condensate changes from w Φ = $\frac{1}{3}$ to w Φ = 0 when λ⁢Φ 2 drops below $m$$^{2}_{Φ}$. We show that, for a mass mΦ ≳ 3⁢λ $\frac{1}{4}$ T RH , the mass term will dominate at reheating. The value of λ is relatively model independent as it is normalized by the cosmic microwave background perturbation spectrum. For T models of inflation, this leads to m Φ ≳ T RH /250. We compute the effects on the reheating temperature for cases where reheating is due to inflaton decay (to fermions, scalars, or vectors) or to inflaton scattering (to scalars or vectors). For scattering to scalars and in the absence of a decay, there is always a residual inflaton background that acts as cold dark matter. In this case, we derive a strong upper limit to the inflaton bare mass which for T models is m Φ < 350 MeV⁢(T RH /10 10 GeV) 3/5 . We also consider the effect of the bare mass term on the fragmentation of the inflaton condensate.

79 ASTRONOMY AND ASTROPHYSICS↗

Draft Rules for Alloy 617 Creep-Fatigue Design Using an EPP+SMT Approach

The report describes an improved new creep-fatigue design method for structural components in high temperature nuclear service. The new method uses an integrated elastic-perfectly plastic (EPP) analysis and Simplified Model Test (SMT) approach that greatly simplifies the design process by avoiding the separate evaluation of creep and fatigue damages and eliminating the requirement of stress classification. The report includes draft design rules, presented in a format compatible with an ASME nuclear Code Case, as well as a commentary on the rules and the validation of the rules using pressurized SMT (p-SMT) test data. Areas of potential improvement have been identified for further evaluation. The report also develops EPP+SMT design charts for Alloy 617 at temperatures between 800°C and 950°C by extrapolating high strain range, short hold test data to low strain range and long hold times, typical for structural components in high temperature nuclear service. The design charts are currently limited to 10,000 cycles to avoid excessive extrapolation outside the test database. More robust alternate extrapolation procedures will be considered in future work. The report includes several verification problems comparing the new EPP+SMT design method with current ASME creep-fatigue design methods. The comparison demonstrates that the EPP+SMT method significantly reduces the over-conservatism in current methods. Finally, the report includes a fully-documented sample problem with multiple load cases detailing the application of proposed EPP+SMT design rules.

36 MATERIALS SCIENCE↗

Deep learning uncertainty quantification for clinical text classification

Machine learning algorithms are expected to work side-by-side with humans in decision-making pipelines. Thus, the ability of classifiers to make reliable decisions is of paramount importance. Deep neural networks (DNNs) represent the state-of-the-art models to address real-world classification. Although the strength of activation in DNNs is often correlated with the network’s confidence, in-depth analyses are needed to establish whether they are well calibrated. In this paper, we demonstrate the use of DNN-based classification tools to benefit cancer registries by automating information extraction of disease at diagnosis and at surgery from electronic text pathology reports from the US National Cancer Institute (NCI) Surveillance, Epidemiology, and End Results (SEER) population-based cancer registries. In particular, we introduce multiple methods for selective classification to achieve a target level of accuracy on multiple classification tasks while minimizing the rejection amount—that is, the number of electronic pathology reports for which the model’s predictions are unreliable. We evaluate the proposed methods by comparing our approach with the current in-house deep learning-based abstaining classifier. Overall, all the proposed selective classification methods effectively allow for achieving the targeted level of accuracy or higher in a trade-off analysis aimed to minimize the rejection rate. On in-distribution validation and holdout test data, with all the proposed methods, we achieve on all tasks the required target level of accuracy with a lower rejection rate than the deep abstaining classifier (DAC). Interpreting the results for the out-of-distribution test data is more complex; nevertheless, in this case as well, the rejection rate from the best among the proposed methods achieving 97% accuracy or higher is lower than the rejection rate based on the DAC. We show that although both approaches can flag those samples that should be manually reviewed and labeled by human annotators, the newly proposed methods retain a larger fraction and do so without retraining—thus offering a reduced computational cost compared with the in-house deep learning-based abstaining classifier.

59 BASIC BIOLOGICAL SCIENCES↗

Assessing models of force-dependent unbinding rates via infrequent metadynamics

Protein–ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent “bonds” being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed—slip bonds, where the unbinding rate increases, and catch bonds, where the unbinding rate decreases. Despite much theoretical work, we cannot predict for which protein–ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin–streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite the use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predicts an overall slip bond for larger forces consistent with prior experimental results although there are substantial deviations at small forces that require further investigation. Furthermore, this work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗