Hydration of Aromatic Heterocycles as an Adversary of π-Stacking

Unveiling the Role of Water on π-π Stacking Through Microwave Spectroscopy of (Thiophene)2-(Water)1-2 Clusters

There has been an ongoing debate about whether water enhances or hinders π-π stacking, a phenomenon crucial in various biological and chemical systems. In this study, the influence of water on π-π stacking is investigated by microwave spectroscopic observation of gas-phase hydrated clusters of thiophene dimers. Two isomers of (C4H4S)2-H2O and two isomers of (C4H4S)2-(H2O)2 have been unambiguously identified. These identifications are supported by quantum chemistry calculations and isotopic measurements. In each of these conformations, water molecules are situated between aromatic pairs, forming distinctive interactions. Water molecules engage with thiophene molecules either as hydrogen bond donors through OH···π interactions or as hydrogen bond acceptors through CH···O interactions. The energy decomposition analysis indicates that the bonding pattern of water molecules significantly affects the π···π interactions between aromatic rings. These findings offer valuable structural insights into the role of water in shaping π-π stacking.

Human Endocrine-Disrupting Effects of Phthalate Esters through Adverse Outcome Pathways: A Comprehensive Mechanism Analysis

Phthalate esters (PAEs) are widely exposed in the environment as plasticizers in plastics, and they have been found to cause significant environmental and health hazards, especially in terms of endocrine disruption in humans. In order to investigate the processes underlying the endocrine disruption effects of PAEs, three machine learning techniques were used in this study to build an adverse outcome pathway (AOP) for those effects on people. According to the results of the three machine learning techniques, the random forest and XGBoost models performed well in terms of prediction. Subsequently, sensitivity analysis was conducted to identify the initial events, key events, and key features influencing the endocrine disruption effects of PAEs on humans. Key features, such as Mol.Wt, Q+, QH+, ELUMO, minHCsats, MEDC-33, and EG, were found to be closely related to the molecular structure. Therefore, a 3D-QSAR model for PAEs was constructed, and, based on the three-dimensional potential energy surface information, it was discovered that the hydrophobic, steric, and electrostatic fields of PAEs significantly influence their endocrine disruption effects on humans. Lastly, an analysis of the contributions of amino acid residues and binding energy (BE) was performed, identifying and confirming that hydrogen bonding, hydrophobic interactions, and van der Waals forces are important factors affecting the AOP of PAEs' molecular endocrine disruption effects. This study defined and constructed a comprehensive AOP for the endocrine disruption effects of PAEs on humans and developed a method based on theoretical simulation to characterize the AOP, providing theoretical guidance for studying the mechanisms of toxicity caused by other pollutants.

Grid inhomogeneous solvation theory for cross-solvation in rigid solvents

Grid Inhomogeneous Solvation Theory (GIST) has proven useful to calculate localized thermodynamic properties of water around a solute. Numerous studies have leveraged this information to enhance structure-based binding predictions. We have recently extended GIST toward chloroform as a solvent to allow the prediction of passive membrane permeability. Here, we further generalize the GIST algorithm toward all solvents that can be modeled as rigid molecules. This restriction is inherent to the method and is already present in the inhomogeneous solvation theory. Here, we show that our approach can be applied to various solvent molecules by comparing the results of GIST simulations with thermodynamic integration (TI) calculations and experimental results. Additionally, we analyze and compare a matrix consisting of 100 entries of ten different solvent molecules solvated within each other. We find that the GIST results are highly correlated with TI calculations as well as experiments. For some solvents, we find Pearson correlations of up to 0.99 to the true entropy, while others are affected by the first-order approximation more strongly. The enthalpy-entropy splitting provided by GIST allows us to extend a recently published approach, which estimates higher order entropies by a linear scaling of the first-order entropy, to solvents other than water. Furthermore, we investigate the convergence of GIST in different solvents. We conclude that our extension to GIST reliably calculates localized thermodynamic properties for different solvents and thereby significantly extends the applicability of this widely used method.

Predicting the Relative Binding Affinity for Reversible Covalent Inhibitors by Free Energy Perturbation Calculations

Covalent inhibitors are assuming central importance in drug discovery projects, especially in this pandemic scenario. Many research groups have focused their attention on inhibiting viral proteases or human proteases such as cathepsin L (hCatL). The inhibition of these critical enzymes may impair viral replication. However, molecular modeling of covalent ligands is challenging since covalent and noncovalent ligand-bound states must be considered in the binding process. In this work, we evaluated the suitability of free energy perturbation (FEP) calculations as a tool for predicting the binding affinity of reversible covalent inhibitors of hCatL. Our strategy relies on the relative free energy calculated for both covalent and noncovalent complexes and the free energy changes have been compared with experimental data for eight nitrile-based inhibitors, including three new inhibitors of hCatL. Our results demonstrate that the covalent complex can be employed to properly rank the inhibitors. Nevertheless, a comparison of the free energy changes in both noncovalent and covalent states is valuable to interpret the effect triggered by the formation of the covalent bond on the interactions played by functional groups distant from the warhead. Overall, FEP can be employed as a powerful predictor tool in developing and understanding the activity of reversible covalent inhibitors.

Conformational Shifts of Stacked Heteroaromatics: Vacuum vs. Water Studied by Machine Learning

Stacking interactions play a crucial role in drug design, as we can find aromatic cores or scaffolds in almost any available small molecule drug. To predict optimal binding geometries and enhance stacking interactions, usually high-level quantum mechanical calculations are performed. These calculations have two major drawbacks: they are very time consuming, and solvation can only be considered using implicit solvation. Therefore, most calculations are performed in vacuum. However, recent studies have revealed a direct correlation between the desolvation penalty, vacuum stacking interactions and binding affinity, making predictions even more difficult. To overcome the drawbacks of quantum mechanical calculations, in this study we use neural networks to perform fast geometry optimizations and molecular dynamics simulations of heteroaromatics stacked with toluene in vacuum and in explicit solvation. We show that the resulting energies in vacuum are in good agreement with high-level quantum mechanical calculations. Furthermore, we show that using explicit solvation substantially influences the favored orientations of heteroaromatic rings thereby emphasizing the necessity to include solvation properties starting from the earliest phases of drug design.

Quantum Chemical Microsolvation by Automated Water Placement

We developed a quantitative approach to quantum chemical microsolvation. Key in our methodology is the automatic placement of individual solvent molecules based on the free energy solvation thermodynamics derived from molecular dynamics (MD) simulations and grid inhomogeneous solvation theory (GIST). This protocol enabled us to rigorously define the number, position, and orientation of individual solvent molecules and to determine their interaction with the solute based on physical quantities. The generated solute-solvent clusters served as an input for subsequent quantum chemical investigations. We showcased the applicability, scope, and limitations of this computational approach for a number of small molecules, including urea, 2-aminobenzothiazole, (+)-syn-benzotriborneol, benzoic acid, and helicene. Our results show excellent agreement with the available ab initio molecular dynamics data and experimental results.

Conformational Ensembles of Antibodies Determine Their Hydrophobicity

A major challenge in the development of antibody biotherapeutics is their tendency to aggregate. One root cause for aggregation is exposure of hydrophobic surface regions to the solvent. Many current techniques predict the relative aggregation propensity of antibodies via precalculated scales for the hydrophobicity or aggregation propensity of single amino acids. However, those scales cannot describe the nonadditive effects of a residue’s surrounding on its hydrophobicity. Therefore, they are inherently limited in their ability to describe the impact of subtle differences in molecular structure on the overall hydrophobicity. Here, we introduce a physics-based approach to describe hydrophobicity in terms of the hydration free energy using grid inhomogeneous solvation theory (GIST). We apply this method to assess the effects of starting structures, conformational sampling, and protonation states on the hydrophobicity of antibodies. Our results reveal that high-quality starting structures, i.e., crystal structures, are crucial for the prediction of hydrophobicity and that conformational sampling can compensate errors introduced by the starting structure. On the other hand, sampling of protonation states only leads to good results when combined with high-quality structures, whereas it can even be detrimental otherwise. We conclude by pointing out that a single static homology model may not be adequate for predicting hydrophobicity.

Solvation Thermodynamics in Different Solvents: Water-Chloroform Partition Coefficients from Grid Inhomogeneous Solvation Theory

Reliable information on partition coefficients plays a key role in drug development, as solubility decisively affects bioavailability. In a physicochemical context, the partition coefficient of a solute between two different solvents can be described as a function of solvation free energies. Hence, substantial scientific efforts have been made toward accurate predictions of solvation free energies in various solvents. The grid inhomogeneous solvation theory (GIST) facilitates the calculation of solvation free energies. In this study, we introduce an extended version of the GIST algorithm, which enables the calculation for chloroform in addition to water. Furthermore, GIST allows localization of enthalpic and entropic contributions. We test our approach by calculating partition coefficients between water and chloroform for a set of eight small molecules. We report a Pearson correlation coefficient of 0.96 between experimentally determined and calculated partition coefficients. The capability to reliably predict partition coefficients between water and chloroform and the possibility to localize their contributions allow the optimization of a compound's partition coefficient. Therefore, we presume that this methodology will be of great benefit for the efficient development of pharmaceuticals.

Macrocycle Cell Permeability Measured by Solvation Free Energies in Polar and Apolar Environments

The relation of surface polarity and conformational preferences is decisive for cell permeability and thus bioavailability of macrocyclic drugs. Here, we employ grid inhomogeneous solvation theory (GIST) to calculate solvation free energies for a series of six macrocycles in water and chloroform as a measure of passive membrane permeability. We perform accelerated molecular dynamics simulations to capture a diverse structural ensemble in water and chloroform, allowing for a direct profiling of solvent-dependent conformational preferences. Subsequent GIST calculations facilitate a quantitative measure of solvent preference in the form of a transfer free energy, calculated from the ensemble-averaged solvation free energies in water and chloroform. Hence, the proposed method considers how the conformational diversity of macrocycles in polar and apolar solvents translates into transfer free energies. Following this strategy, we find a striking correlation of 0.92 between experimentally determined cell permeabilities and calculated transfer free energies. For the studied model systems, we find that the transfer free energy exceeds the purely water-based solvation free energies as a reliable estimate of cell permeability and that conformational sampling is imperative for a physically meaningful model. We thus recommend this purely physics-based approach as a computational tool to assess cell permeabilities of macrocyclic drug candidates.

STACKED - Solvation Theory of Aromatic Complexes as Key for Estimating Drug Binding

The use of fragments to biophysically characterize a protein binding pocket and determine the strengths of certain interactions is a computationally and experimentally commonly applied approach. Almost all drug like molecules contain at least one aromatic moiety forming stacking interactions in the binding pocket. In computational drug design, the strength of stacking and the resulting optimization of the aromatic core or moiety is usually calculated using high level quantum mechanical approaches. However, as these calculations are performed in a vacuum, solvation properties are neglected. We close this gap by using Grid Inhomogeneous Solvation Theory (GIST) to describe the properties of individual heteroaromatics and complexes and thereby estimate the desolvation penalty. In our study, we investigated the solvation free energies of heteroaromatics frequently occurring in drug design projects in complex with truncated side chains of phenylalanine, tyrosine, and tryptophan. Furthermore, we investigated the properties of drug-fragments crystallized in a fragment-based lead optimization approach investigating PDE-10-A. We do not only find good correlation for the estimated desolvation penalty and the experimental binding free energy, but our calculations also allow us to predict prominent interaction sites. We highlight the importance of including the desolvation penalty of the respective heteroaromatics in stacked complexes to explain the gain or loss in affinity of potential lead compounds.