The Open Force Field Initiative: Open Software and Open Science for Molecular Modeling.

Force fields are a key component of physics-based molecular modeling, describing the energies and forces in a molecular system as a function of the positions of the atoms and molecules involved. Here, we provide a review and scientific status report on the work of the Open Force Field (OpenFF) Initiative, which focuses on the science, infrastructure and data required to build the next generation of biomolecular force fields. We introduce the OpenFF Initiative and the related OpenFF Consortium, describe its approach to force field development and software, and discuss accomplishments to date as well as future plans. OpenFF releases both software and data under open and permissive licensing agreements to enable rapid application, validation, extension, and modification of its force fields and software tools. We discuss lessons learned to date in this new approach to force field development. We also highlight ways that other force field researchers can get involved, as well as some recent successes of outside researchers taking advantage of OpenFF tools and data.

Structure and fragmentation of doubly ionized HNCS.

Double ionization spectra of isothiocyanic acid (HNCS) have been measured using multi-electron and multi-ion coincidence techniques combined with high-level theoretical calculations. The adiabatic double ionization energy of HNCS is found at 27.1 ± 0.1 eV and is associated with the formation of the X 3A″ ground state of HNCS2+. The characteristics of different dissociation channels are examined and compared to the results of electronic structure calculations obtained by systematically elongating the three bonds H-NCS, HN-CS, and HNC-S. For instance, the adiabatic double ionization energy of the NCS fragment is deduced to be 30.95 ± 0.5 eV. In addition, the C+ and NS+ dissociation channels are of particular interest, possibly indicating the involvement of a structural rearrangement process upon doubly ionizing HNCS.

Closing the Chasm: Understanding and Addressing the Anesthesia Workforce Supply and Demand Imbalance.

The imbalance in anesthesia workforce supply and demand has been exacerbated post-COVID due to a surge in demand for anesthesia care, especially in non-operating room anesthetizing sites, at a faster rate than the increase in anesthesia clinicians. The consequences of this imbalance or labor shortage compromise healthcare facilities, adversely affect the cost of care, worsen anesthesia workforce burnout, disrupt procedural and surgical schedules, and threaten academic missions and the ability to educate future anesthesiologists. In developing possible solutions, one must examine emerging trends that are affecting the anesthesia workforce, new technologies that will transform anesthesia care and the workforce, and financial considerations, including governmental payment policies. Possible practice solutions to this imbalance will require both short- and long-term multifactorial approaches that include increasing training positions and retention policies, improving capacity through innovations, leveraging technology, and addressing financial constraints.

Myosin-5 varies its step length to carry cargo straight along the irregular F-actin track.

Molecular motors employ chemical energy to generate unidirectional mechanical output against a track while navigating a chaotic cellular environment, potential disorder on the track, and against Brownian motion. Nevertheless, decades of nanometer-precise optical studies suggest that myosin-5a, one of the prototypical molecular motors, takes uniform steps spanning 13 subunits (36 nm) along its F-actin track. Here, we use high-resolution interferometric scattering microscopy to reveal that myosin takes strides spanning 22 to 34 actin subunits, despite walking straight along the helical actin filament. We show that cumulative angular disorder in F-actin accounts for the observed proportion of each stride length, akin to crossing a river on variably spaced stepping stones. Electron microscopy revealed the structure of the stepping molecule. Our results indicate that both motor and track are soft materials that can adapt to function in complex cellular conditions.

Comparing the quality of ChatGPT- and physician-generated responses to patients' dermatology questions in the electronic medical record.

BACKGROUND: ChatGPT is a free artificial intelligence (AI)-based natural language processing tool that generates complex responses to inputs from users. OBJECTIVES: To determine whether ChatGPT is able to generate high-quality responses to patient-submitted questions in the patient portal. METHODS: Patient-submitted questions and the corresponding responses from their dermatology physician were extracted from the electronic medical record for analysis. The questions were input into ChatGPT (version 3.5) and the outputs extracted for analysis, with manual removal of verbiage pertaining to ChatGPT's inability to provide medical advice. Ten blinded reviewers (seven physicians and three nonphysicians) rated and selected their preference in terms of 'overall quality', 'readability', 'accuracy', 'thoroughness' and 'level of empathy' of the physician- and ChatGPT-generated responses. RESULTS: Thirty-one messages and responses were analysed. Physician-generated responses were vastly preferred over the ChatGPT -responses by the physician and nonphysician reviewers and received significantly higher ratings for 'readability' and 'level of empathy'. CONCLUSIONS: The results of this study suggest that physician-generated responses to patients' portal messages are still preferred over ChatGPT, but generative AI tools may be helpful in generating the first drafts of responses and providing information on education resources for patients.

Clinical and Histopathologic Spectrum of Toxic Erythema of Chemotherapy: A Series of 56 Cases From a Single Institution.

INTRODUCTION: Although many individual cases and small series of toxic erythema of chemotherapy (TEC) have been described, the full spectrum of findings is not well understood. OBJECTIVE: To provide a comprehensive review of the clinical and histopathologic features of TEC with an emphasis on novel histopathologic findings. METHODS: We searched our electronic medical record for "toxic erythema of chemotherapy" or "neutrophilic eccrine hidradenitis." Fifty-six cases meeting clinical and histopathologic criteria were identified. The electronic medical record and accompanying hematoxylin and eosin-stained slides were retrospectively reviewed. RESULTS: The clinical findings were heterogeneous but included classic presentations such as intertriginous eruptions (34%) and acral erythema (25%). The most common histopathologic features were apoptotic keratinocytes (95%), basal vacuolar change (91%), and epithelial dysmaturation (79%). Eccrine squamous syringometaplasia was seen in over half of the cases (33/56; 59%), whereas neutrophilic eccrine hidradenitis was uncommon (16%). Interestingly, many cases showed prominent interstitial histiocytes (55%). Other novel findings included irregular orthohyperkeratosis (23%), irregular epidermal hyperplasia (14%), and acantholysis (9%). LIMITATIONS: As a retrospective study, it is subject to information bias. CONCLUSION: This is the largest reported series of TEC. In addition to confirming previously reported features, we identify novel histopathologic findings to add to the spectrum of TEC.

A transferable double exponential potential for condensed phase simulations of small molecules.

The Lennard-Jones potential is the most widely-used function for the description of non-bonded interactions in transferable force fields for the condensed phase. This is not because it has an optimal functional form, but rather it is a legacy resulting from when computational expense was a major consideration and this potential was particularly convenient numerically. At present, it persists because the effort that would be required to re-write molecular modelling software and train new force fields has, until now, been prohibitive. Here, we present Smirnoff-plugins as a flexible framework to extend the Open Force Field software stack to allow custom force field functional forms. We deploy Smirnoff-plugins with the automated Open Force Field infrastructure to train a transferable, small molecule force field based on the recently-proposed double exponential functional form, on over 1000 experimental condensed phase properties. Extensive testing of the resulting force field shows improvements in transfer free energies, with acceptable conformational energetics, run times and convergence properties compared to state-of-the-art Lennard-Jones based force fields.

Myosin-5 varies its steps along the irregular F-actin track.

Molecular motors employ chemical energy to generate unidirectional mechanical output against a track. By contrast to the majority of macroscopic machines, they need to navigate a chaotic cellular environment, potential disorder in the track and Brownian motion. Nevertheless, decades of nanometer-precise optical studies suggest that myosin-5a, one of the prototypical molecular motors, takes uniform steps spanning 13 subunits (36 nm) along its F-actin track. Here, we use high-resolution interferometric scattering (iSCAT) microscopy to reveal that myosin takes strides spanning 22 to 34 actin subunits, despite walking straight along the helical actin filament. We show that cumulative angular disorder in F-actin accounts for the observed proportion of each stride length, akin to crossing a river on variably-spaced stepping stones. Electron microscopy revealed the structure of the stepping molecule. Our results indicate that both motor and track are soft materials that can adapt to function in complex cellular conditions.

Comparison of Receptor-Ligand Restraint Schemes for Alchemical Absolute Binding Free Energy Calculations.

Alchemical absolute binding free energy calculations are of increasing interest in drug discovery. These calculations require restraints between the receptor and ligand to restrict their relative positions and, optionally, orientations. Boresch restraints are commonly used, but they must be carefully selected in order to sufficiently restrain the ligand and to avoid inherent instabilities. Applying multiple distance restraints between anchor points in the receptor and ligand provides an alternative framework without inherent instabilities which may provide convergence benefits by more strongly restricting the relative movements of the receptor and ligand. However, there is no simple method to calculate the free energy of releasing these restraints due to the coupling of the internal and external degrees of freedom of the receptor and ligand. Here, a method to rigorously calculate free energies of binding with multiple distance restraints by imposing intramolecular restraints on the anchor points is proposed. Absolute binding free energies for the human macrophage migration inhibitory factor/MIF180, system obtained using a variety of Boresch restraints and rigorous and nonrigorous implementations of multiple distance restraints are compared. It is shown that several multiple distance restraint schemes produce estimates in good agreement with Boresch restraints. In contrast, calculations without orientational restraints produce erroneously favorable free energies of binding by up to approximately 4 kcal mol-1. These approaches offer new options for the deployment of alchemical absolute binding free energy calculations.

Development and Benchmarking of Open Force Field 2.0.0: The Sage Small Molecule Force Field.

We introduce the Open Force Field (OpenFF) 2.0.0 small molecule force field for drug-like molecules, code-named Sage, which builds upon our previous iteration, Parsley. OpenFF force fields are based on direct chemical perception, which generalizes easily to highly diverse sets of chemistries based on substructure queries. Like the previous OpenFF iterations, the Sage generation of OpenFF force fields was validated in protein-ligand simulations to be compatible with AMBER biopolymer force fields. In this work, we detail the methodology used to develop this force field, as well as the innovations and improvements introduced since the release of Parsley 1.0.0. One particularly significant feature of Sage is a set of improved Lennard-Jones (LJ) parameters retrained against condensed phase mixture data, the first refit of LJ parameters in the OpenFF small molecule force field line. Sage also includes valence parameters refit to a larger database of quantum chemical calculations than previous versions, as well as improvements in how this fitting is performed. Force field benchmarks show improvements in general metrics of performance against quantum chemistry reference data such as root-mean-square deviations (RMSD) of optimized conformer geometries, torsion fingerprint deviations (TFD), and improved relative conformer energetics (ΔΔE). We present a variety of benchmarks for these metrics against our previous force fields as well as in some cases other small molecule force fields. Sage also demonstrates improved performance in estimating physical properties, including comparison against experimental data from various thermodynamic databases for small molecule properties such as ΔHmix, ρ(x), ΔGsolv, and ΔGtrans. Additionally, we benchmarked against protein-ligand binding free energies (ΔGbind), where Sage yields results statistically similar to previous force fields. All the data is made publicly available along with complete details on how to reproduce the training results at https://github.com/openforcefield/openff-sage.

What is the Optimal Dipole Moment for Nonpolarizable Models of Liquids?

In classical nonpolarizable models, electrostatic interactions are usually described by assigning fixed partial charges to interaction sites. Despite the multitude of methods and theories proposed over the years for partial charge assignment, a fundamental question remains─what is the correct degree of polarization that a fixed-charge model should possess to provide the best balance of interactions (including induction effects) and yield the best description of the potential energy surface of a liquid phase? We address this question by approaching it from two separate and independent viewpoints: the QUantum mechanical BEspoke (QUBE) approach, which assigns bespoke force field parameters for individual molecules from ab initio calculations with minimal empirical fitting, and the Polarization-Consistent Approach (PolCA) force field, based on empirical fitting of force field parameters with an emphasis on transferability by rigorously accounting for polarization effects in the parameterization process. We show that the two approaches yield consistent answers to the above question, namely, that the dipole moment of the model should be approximately halfway between those of the gas and the liquid phase. Crucially, however, the reference liquid-phase dipole needs to be estimated using methods that explicitly consider both mean-field and local contributions to polarization. In particular, continuum dielectric models are inadequate for this purpose because they cannot account for local effects and therefore significantly underestimate the degree of polarization of the molecule. These observations have profound consequences for the development, validation, and testing of nonpolarizable models.

An open-source molecular builder and free energy preparation workflow.

Automated free energy calculations for the prediction of binding free energies of congeneric series of ligands to a protein target are growing in popularity, but building reliable initial binding poses for the ligands is challenging. Here, we introduce the open-source FEgrow workflow for building user-defined congeneric series of ligands in protein binding pockets for input to free energy calculations. For a given ligand core and receptor structure, FEgrow enumerates and optimises the bioactive conformations of the grown functional group(s), making use of hybrid machine learning/molecular mechanics potential energy functions where possible. Low energy structures are optionally scored using the gnina convolutional neural network scoring function, and output for more rigorous protein-ligand binding free energy predictions. We illustrate use of the workflow by building and scoring binding poses for ten congeneric series of ligands bound to targets from a standard, high quality dataset of protein-ligand complexes. Furthermore, we build a set of 13 inhibitors of the SARS-CoV-2 main protease from the literature, and use free energy calculations to retrospectively compute their relative binding free energies. FEgrow is freely available at https://github.com/cole-group/FEgrow, along with a tutorial.

Open Force Field BespokeFit: Automating Bespoke Torsion Parametrization at Scale.

The development of accurate transferable force fields is key to realizing the full potential of atomistic modeling in the study of biological processes such as protein-ligand binding for drug discovery. State-of-the-art transferable force fields, such as those produced by the Open Force Field Initiative, use modern software engineering and automation techniques to yield accuracy improvements. However, force field torsion parameters, which must account for many stereoelectronic and steric effects, are considered to be less transferable than other force field parameters and are therefore often targets for bespoke parametrization. Here, we present the Open Force Field QCSubmit and BespokeFit software packages that, when combined, facilitate the fitting of torsion parameters to quantum mechanical reference data at scale. We demonstrate the use of QCSubmit for simplifying the process of creating and archiving large numbers of quantum chemical calculations, by generating a dataset of 671 torsion scans for druglike fragments. We use BespokeFit to derive individual torsion parameters for each of these molecules, thereby reducing the root-mean-square error in the potential energy surface from 1.1 kcal/mol, using the original transferable force field, to 0.4 kcal/mol using the bespoke version. Furthermore, we employ the bespoke force fields to compute the relative binding free energies of a congeneric series of inhibitors of the TYK2 protein, and demonstrate further improvements in accuracy, compared to the base force field (MUE reduced from 0.560.390.77 to 0.420.280.59 kcal/mol and R2 correlation improved from 0.720.350.87 to 0.930.840.97).

Direct observation of the molecular mechanism underlying protein polymerization.

Protein assembly is a main route to generating complexity in living systems. Revealing the relevant molecular details is challenging because of the intrinsic heterogeneity of species ranging from few to hundreds of molecules. Here, we use mass photometry to quantify and monitor the full range of actin oligomers during polymerization with single-molecule sensitivity. We find that traditional nucleation-based models cannot account for the observed distributions of actin oligomers. Instead, the key step of filament formation is a slow transition between distinct states of an actin filament mediated by cation exchange or ATP hydrolysis. The resulting model reproduces important aspects of actin polymerization, such as the critical concentration for filament formation and bulk growth behavior. Our results revise the mechanism of actin nucleation, shed light on the role and function of actin-associated proteins, and introduce a general and quantitative means to studying protein assembly at the molecular level.

Exploration and validation of force field design protocols through QM-to-MM mapping.

The scale of the parameter optimisation problem in traditional molecular mechanics force field construction means that design of a new force field is a long process, and sub-optimal choices made in the early stages can persist for many generations. We hypothesise that careful use of quantum mechanics to inform molecular mechanics parameter derivation (QM-to-MM mapping) should be used to significantly reduce the number of parameters that require fitting to experiment and increase the pace of force field development. Here, we design and train a collection of 15 new protocols for small, organic molecule force field derivation, and test their accuracy against experimental liquid properties. Our best performing model has only seven fitting parameters, yet achieves mean unsigned errors of just 0.031 g cm-3 and 0.69 kcal mol-1 in liquid densities and heats of vaporisation, compared to experiment. The software required to derive the designed force fields is freely available at https://github.com/qubekit/QUBEKit.

Perioperative Brain Health in the Older Adult: A Patient Safety Imperative.

While people 65 years of age and older represent 16% of the population in the United States, they account for >40% of surgical procedures performed each year. Maintaining brain health after anesthesia and surgery is not only important to our patients, but it is also an increasingly important patient safety imperative for the specialty of anesthesiology. Aging is a complex process that diminishes the reserve of every organ system and often results in a patient who is vulnerable to the stress of surgery. The brain is no exception, and many older patients present with preoperative cognitive impairment that is undiagnosed. As we age, a number of changes occur in the human brain, resulting in a patient who is less resilient to perioperative stress, making older adults more susceptible to the phenotypic expression of perioperative neurocognitive disorders. This review summarizes the current scientific and clinical understanding of perioperative neurocognitive disorders and recommends patient-centered, age-focused interventions that can better mitigate risk, prevent harm, and improve outcomes for our patients. Finally, it discusses the emerging topic of sleep and cognitive health and other future frontiers of scientific inquiry that might inform clinical best practices.

High-Fidelity Indirect Readout of Trapped-Ion Hyperfine Qubits.

We propose and demonstrate a protocol for high-fidelity indirect readout of trapped ion hyperfine qubits, where the state of a ^{9}Be^{+} qubit ion is mapped to a ^{25}Mg^{+} readout ion using laser-driven Raman transitions. By partitioning the ^{9}Be^{+} ground-state hyperfine manifold into two subspaces representing the two qubit states and choosing appropriate laser parameters, the protocol can be made robust to spontaneous photon scattering errors on the Raman transitions, enabling repetition for increased readout fidelity. We demonstrate combined readout and back-action errors for the two subspaces of 1.2_{-0.6}^{+1.1}×10^{-4} and 0_{-0}^{+1.9}×10^{-5} with 68% confidence while avoiding decoherence of spectator qubits due to stray resonant light that is inherent to direct fluorescence detection.