OneOPES, a Combined Enhanced Sampling Method to Rule Them All

Leveraging Cryptic Ligand Envelopes through Enhanced Molecular Simulations

Protein-bound ligands can adopt a range of different conformations, collectively defining a ligand envelope that has proven to be crucial for the design of potent and selective drugs. Yet, the cryptic nature of this ligand envelope makes it difficult to visualize, characterize, and ultimately exploit for drug design. Using enhanced molecular dynamics simulations, here, we provide a general framework to reconstruct the cryptic ligand envelope that is dynamically accessible by protein-bound small molecules in solution. We apply this approach to quantify hidden conformational heterogeneity in structurally complex ligands including the marine natural product plitidepsin. The computed conformational heterogeneity expands the small-molecule footprint beyond that typically observed in experiments, also revealing key thermodynamic and kinetic properties of single ligand-target interactions. The model agrees quantitatively with solution NMR, X-ray crystallography, and biochemical measurements, showcasing a versatile strategy to integrate receptor-bound ligand conformational ensembles in molecular design.

Enhanced Sampling with Suboptimal Collective Variables: Reconciling Accuracy and Convergence Speed

We introduce an enhanced sampling algorithm to obtain converged free energy landscapes of molecular rare events, even when the collective variable (CV) used for biasing is not optimal. Our approach samples a time-dependent target distribution by combining the on-the-fly probability enhanced sampling and its exploratory variant, OPES Explore. This promotes more transitions between the relevant metastable states and accelerates the convergence speed of the free energy estimate. We demonstrate the successful application of this combined algorithm on the two-dimensional Wolfe-Quapp potential, millisecond time-scale ligand-receptor binding in the trypsin-benzamidine complex, and folding-unfolding transitions in chignolin mini-protein. Our proposed algorithm can compute accurate free energies at an affordable computational cost and is robust in terms of the choice of CVs, making it particularly promising for the simulation of complex biomolecular systems.

Binding Modes and Water-Mediation of Polyelectrolyte Adsorption to a Neutral CaCO3 Surface

Aqueous polyelectrolytes are effective mineralization inhibitors due to their ability to template onto crystal surfaces and chelate ions in solution. These additives have been shown to alter the morphology of calcium carbonate crystals, making them promising candidates for biological and industrial applications. However, while key to designing more effective mineralization inhibitors, the molecular mechanisms governing the interactions between polyelectrolytes and crystal surfaces remain poorly understood. In this study, we investigate the adsorption of poly(acrylic acid) (PAA) on the dominant calcite ( 101̅ 4 ) cleavage plane using all-atom molecular dynamics simulations. Although the calcite slab is electrostatically neutral, its charge distribution induces a strong electrostatic potential in an aqueous solution, which leads to significant water structuring at the interface. We observe a very favorable adsorption affinity of the polyelectrolyte chain to the surface, yet the structure of the interfacial water is not significantly affected. Direct interactions between the monomers on the polyelectrolyte and the calcite surface are infrequent, despite variations in chain length, charge density of the polyelectrolyte, and solution conditions. Intriguingly, the polyelectrolyte interaction with the calcite surface is dominantly mediated through bridging hydrogen bond interactions. As the polyelectrolyte adsorbs to the surface, the chain conformation adapts to the interfacial water structure by increasing polyelectrolyte-water contacts and integrates into pre-existing hydrogen bond networks. We found that water-mediated interactions are more dominant than direct interactions between the polyelectrolyte and the surface. This suggests an alternative pathway to the widely accepted notion that entropic effects due to water reorganization are the primary driving force. These results suggest that the polyelectrolyte binding affinity can be tuned by altering the polymer chain interactions with the interfacial water structure in addition to the surface itself.

Host-Guest Binding Free Energies à la Carte: An Automated OneOPES Protocol

Estimating absolute binding free energies from molecular simulations is a key step in computer-aided drug design pipelines, but the agreement between computational results and experiments is still very inconsistent. Both the accuracy of the computational model and the quality of the statistical sampling contribute to this discrepancy, yet disentangling the two remains a challenge. In this study, we present an automated protocol based on OneOPES, an enhanced sampling method that exploits replica exchange and can accelerate several collective variables to address the sampling problem. We apply this protocol to 37 host-guest systems. The simplicity of setting up the simulations and producing well-converged binding free energy estimates without the need to optimize simulation parameters provides a reliable solution to the sampling problem. This, in turn, allows for a systematic force field comparison and ranking according to the correlation between simulations and experiments, which can inform the selection of an appropriate model. The protocol can be readily adapted to test more force field combinations and study more complex protein-ligand systems, where the choice of an appropriate physical model is often based on heuristic considerations rather than systematic optimization.

Investigating Ligand-Mediated Conformational Dynamics of Pre-miR21: A Machine-Learning-Aided Enhanced Sampling Study

MicroRNAs (miRs) are short, noncoding RNA strands that regulate the activity of mRNAs by affecting the repression of protein translation, and their dysregulation has been implicated in several pathologies. miR21 in particular has been implicated in tumorigenesis and anticancer drug resistance, making it a critical target for drug design. miR21 biogenesis involves precise biochemical pathways, including the cleavage of its precursor, pre-miR21, by the enzyme Dicer. The present work investigates the conformational dynamics of pre-miR21, focusing on the role of adenine29 in switching between Dicer-binding-prone and inactive states. We also investigated the effect of L50, a cyclic peptide binder of pre-miR21 and a weak inhibitor of its processing. Using time series data and our novel collective variable-based enhanced sampling technique, OneOPES, we simulated these conformational changes and assessed the effect of L50 on the conformational plasticity of pre-miR21. Our results provide insight into peptide-induced conformational changes and pave the way for the development of a computational platform for the screening of inhibitors of pre-miR21 processing that considers RNA flexibility, a stepping stone for effective structure-based drug design, with potentially broad applications in drug discovery.

Absolute Binding Free Energies with OneOPES

The calculation of absolute binding free energies (ABFEs) for protein-ligand systems has long been a challenge. Recently, refined force fields and algorithms have improved the quality of the ABFE calculations. However, achieving the level of accuracy required to inform drug discovery efforts remains difficult. Here, we present a transferable enhanced sampling strategy to accurately calculate absolute binding free energies using OneOPES with simple geometric collective variables. We tested the strategy on two protein targets, BRD4 and Hsp90, complexed with a total of 17 chemically diverse ligands, including both molecular fragments and drug-like molecules. Our results show that OneOPES accurately predicts protein-ligand binding affinities with a mean unsigned error within 1 kcal mol-1 of experimentally determined free energies, without the need to tailor the collective variables to each system. Furthermore, our strategy effectively samples different ligand binding modes and consistently matches the experimentally determined structures regardless of the initial protein-ligand configuration. Our results suggest that the proposed OneOPES strategy can be used to inform lead optimization campaigns in drug discovery and to study protein-ligand binding and unbinding mechanisms.

Revolutionizing Molecular Design for Innovative Therapeutic Applications through Artificial Intelligence

The field of computational protein engineering has been transformed by recent advancements in machine learning, artificial intelligence, and molecular modeling, enabling the design of proteins with unprecedented precision and functionality. Computational methods now play a crucial role in enhancing the stability, activity, and specificity of proteins for diverse applications in biotechnology and medicine. Techniques such as deep learning, reinforcement learning, and transfer learning have dramatically improved protein structure prediction, optimization of binding affinities, and enzyme design. These innovations have streamlined the process of protein engineering by allowing the rapid generation of targeted libraries, reducing experimental sampling, and enabling the rational design of proteins with tailored properties. Furthermore, the integration of computational approaches with high-throughput experimental techniques has facilitated the development of multifunctional proteins and novel therapeutics. However, challenges remain in bridging the gap between computational predictions and experimental validation and in addressing ethical concerns related to AI-driven protein design. This review provides a comprehensive overview of the current state and future directions of computational methods in protein engineering, emphasizing their transformative potential in creating next-generation biologics and advancing synthetic biology.

Entropy drives the ligand recognition in G-protein-coupled receptor subtypes

Achieving ligand subtype selectivity within highly homologous subtypes of G-protein-coupled receptor (GPCR) is critical yet challenging for GPCR drug discovery, primarily due to the unclear mechanism underlying ligand subtype selectivity, which hampers the rational design of subtype-selective ligands. Herein, we disclose an unusual molecular mechanism of entropy-driven ligand recognition in cannabinoid (CB) receptor subtypes, revealed through atomic-level molecular dynamics simulations, cryoelectron microscopy structure, and mutagenesis experiments. This mechanism is attributed to the distinct conformational dynamics of the receptor's orthosteric pocket, leading to variations in ligand binding entropy and consequently, differential binding affinities, which culminate in specific ligand recognition. We experimentally validated this mechanism and leveraged it to design ligands with enhanced or ablated subtype selectivity. One such ligand demonstrated favorable pharmacokinetic properties and significant efficacy in rodent inflammatory analgesic models. More importantly, it is precisely due to the high subtype selectivity obtained based on this mechanism that this ligand does not show addictive properties in animal models. Our findings elucidate the unconventional role of entropy in CB receptor subtype selectivity and suggest a strategy for rational design of ligands to achieve entropy-driven subtype selectivity for many pharmaceutically important GPCRs.

Addressing Suboptimal Poses in Nonequilibrium Alchemical Calculations

Alchemical transformations can be used to quantitatively estimate absolute binding free energies at a reasonable computational cost. However, most of the approaches currently in use require knowledge of the correct (crystallographic) pose. In this paper, we present a combined Hamiltonian replica exchange nonequilibrium alchemical method that allows us to reliably calculate absolute binding free energies, even when starting from suboptimal initial binding poses. Performing a preliminary Hamiltonian replica exchange enhances the sampling of slow degrees of freedom of the ligand and the target, allowing the system to populate the correct binding pose when starting from an approximate docking pose. We apply the method on 6 ligands of the first bromodomain of the BRD4 bromodomain-containing protein. For each ligand, we start nonequilibrium alchemical transformations from both the crystallographic pose and the top-scoring docked pose that are often significantly different. We show that the method produces statistically equivalent binding free energies, making it a useful tool for computational drug discovery pipelines.