Parameterization of a drug molecule with a halogen σ-hole particle using ffTK: Implementation, testing, and comparison

Broadening access to small-molecule parameterization with the force field toolkit

Access to accurate force-field parameters for small molecules is crucial for computational studies of their interactions with proteins. Although a number of general force fields for small molecules exist, e.g., CGenFF, GAFF, and OPLS, they do not cover all common chemical groups and their combinations. The Force Field Toolkit (ffTK) provides a comprehensive graphical interface that streamlines the development of classical parameters for small molecules directly from quantum mechanical (QM) calculations, allowing for force-field generation for almost any chemical group and validation of the fit relative to the target data. ffTK relies on supported external software for the QM calculations, but it can generate the necessary QM input files and parse and analyze the QM output. In previous ffTK versions, support for Gaussian and ORCA QM packages was implemented. Here, we add support for Psi4, an open-source QM package free for all users, thereby broadening user access to ffTK. We also compare the parameter sets obtained with the new ffTK version using Gaussian, ORCA, and Psi4 for three molecules: pyrrolidine, n-propylammonium cation, and chlorobenzene. Despite minor differences between the resulting parameter sets for each compound, most prominently in the dihedral and improper terms, we show that conformational distributions sampled in molecular dynamics simulations using these parameter sets are quite comparable.

Development of CHARMM compatible force field parameters and molecular dynamics simulations for the pesticide flupyradifurone

Flupyradifurone (FLU) is a novel butenolide insecticide with partial agonist activity for insect nicotinic acetylcholine receptors. Its safety for non-target organisms has been questioned in the literature, despite initial claims of its harmlessness. Detailed understanding of its toxicity and related molecular mechanisms remain under discussion. Thus, in this work, an optimized set of CHARMM compatible parameters for FLU is presented. CHARMM General Force Field program was used as a starting point while the non-bonded and bonded parameters were adjusted and optimized to reproduce MP2/6-31G(d) accuracy level results. For the validity assessment of these parameters, infrared spectrum, water-octanol partition coefficient, and normal modes were computed and compared to experimental values found in the literature. Several MD simulations of FLU in water and FLU in complex with an acetylcholine-binding protein were performed to estimate the ability of the optimized parameters to correctly describe its torsional space and reproduce observed crystallographic trends respectively.

VMD as a Platform for Interactive Small Molecule Preparation and Visualization in Quantum and Classical Simulations

Modeling and simulation of small molecules such as drugs and biological cofactors have been both a major focus of computational chemistry for decades and a growing need among computational biophysicists who seek to investigate the interaction of different types of ligands with biomolecules. Of particular interest in this regard are quantum mechanical (QM) calculations that are used to more accurately describe such small molecules, which can be of heterogeneous structures and chemistry, either in purely QM calculations or in hybrid QM/molecular mechanics (MM) simulations. QM programs are also used to develop MM force field parameters for small molecules to be used along with established force fields for biomolecules in classical simulations. With this growing need in mind, here we report a set of software tools developed and closely integrated within the broadly used molecular visualization/analysis program, VMD, that allow the user to construct, modify, and parametrize small molecules and prepare them for QM, hybrid QM/MM, or classical simulations. The tools also provide interactive analysis and visualization capabilities in an easy-to-use and integrated environment. In this paper, we briefly report on these tools and their major features and capabilities, along with examples of how they can facilitate molecular research in computational biophysics that might be otherwise prohibitively complex.

Diclofenac and other non-steroidal anti-inflammatory drugs (NSAIDs) are competitive antagonists of the human P2X3 receptor

Introduction: The P2X3 receptor (P2X3R), an ATP-gated non-selective cation channel of the P2X receptor family, is expressed in sensory neurons and involved in nociception. P2X3R inhibition was shown to reduce chronic and neuropathic pain. In a previous screening of 2000 approved drugs, natural products, and bioactive substances, various non-steroidal anti-inflammatory drugs (NSAIDs) were found to inhibit P2X3R-mediated currents.

Methods: To investigate whether the inhibition of P2X receptors contributes to the analgesic effect of NSAIDs, we characterized the potency and selectivity of various NSAIDs at P2X3R and other P2XR subtypes using two-electrode voltage clamp electrophysiology.

Results: We identified diclofenac as a hP2X3R and hP2X2/3R antagonist with micromolar potency (with IC50 values of 138.2 and 76.7 µM, respectively). A weaker inhibition of hP2X1R, hP2X4R, and hP2X7R by diclofenac was determined. Flufenamic acid (FFA) inhibited hP2X3R, rP2X3R, and hP2X7R (IC50 values of 221 µM, 264.1 µM, and ∼900 µM, respectively), calling into question its use as a non-selective ion channel blocker, when P2XR-mediated currents are under study. Inhibition of hP2X3R or hP2X2/3R by diclofenac could be overcome by prolonged ATP application or increasing concentrations of the agonist α,β-meATP, respectively, indicating competition of diclofenac and the agonists. Molecular dynamics simulation showed that diclofenac largely overlaps with ATP bound to the open state of the hP2X3R. Our results suggest a competitive antagonism through which diclofenac, by interacting with residues of the ATP-binding site, left flipper, and dorsal fin domains, inhibits the gating of P2X3R by conformational fixation of the left flipper and dorsal fin domains. In summary, we demonstrate the inhibition of the human P2X3 receptor by various NSAIDs. Diclofenac proved to be the most effective antagonist with a strong inhibition of hP2X3R and hP2X2/3R and a weaker inhibition of hP2X1R, hP2X4R, and hP2X7R.

Discussion: Considering their involvement in nociception, inhibition of hP2X3R and hP2X2/3R by micromolar concentrations of diclofenac, which are rarely reached in the therapeutic range, may play a minor role in analgesia compared to the high-potency cyclooxygenase inhibition but may explain the known side effect of taste disturbances caused by diclofenac.

Mechanisms by which small molecules of diverse chemotypes arrest Sec14 lipid transfer activity

Phosphatidylinositol (PtdIns) transfer proteins (PITPs) enhance the activities of PtdIns 4-OH kinases that generate signaling pools of PtdIns-4-phosphate. In that capacity, PITPs serve as key regulators of lipid signaling in eukaryotic cells. Although the PITP phospholipid exchange cycle is the engine that stimulates PtdIns 4-OH kinase activities, the underlying mechanism is not understood. Herein, we apply an integrative structural biology approach to investigate interactions of the yeast PITP Sec14 with small-molecule inhibitors (SMIs) of its phospholipid exchange cycle. Using a combination of X-ray crystallography, solution NMR spectroscopy, and atomistic MD simulations, we dissect how SMIs compete with native Sec14 phospholipid ligands and arrest phospholipid exchange. Moreover, as Sec14 PITPs represent new targets for the development of next-generation antifungal drugs, the structures of Sec14 bound to SMIs of diverse chemotypes reported in this study will provide critical information required for future structure-based design of next-generation lead compounds directed against Sec14 PITPs of virulent fungi.

On the force field optimisation of [Formula: see text]-lactam cores using the force field Toolkit

When employing molecular dynamics (MD) simulations for computer-aided drug design, the quality of the used force fields is highly important. Here we present reparametrisations of the force fields for the core molecules from 9 different [Formula: see text]-lactam classes, for which we utilized the force field Toolkit and Gaussian calculations. We focus on the parametrisation of the dihedral angles, with the goal of reproducing the optimised quantum geometry in MD simulations. Parameters taken from CGenFF turn out to be a good initial guess for the multiplicity of each dihedral angle, but the key to a successful parametrisation is found to lie in the phase shifts. Based on the optimised quantum geometry, we come up with a strategy for predicting the phase shifts prior to the dihedral potential fitting. This allows us to successfully parameterise 8 out of the 11 molecules studied here, while the remaining 3 molecules can also be parameterised with small adjustments. Our work highlights the importance of predicting the dihedral phase shifts in the ligand parametrisation protocol, and provides a simple yet valuable strategy for improving the process of parameterising force fields of drug-like molecules.

All-Atom MD Simulations of the HBV Capsid Complexed with AT130 Reveal Secondary and Tertiary Structural Changes and Mechanisms of Allostery

The hepatitis B virus (HBV) capsid is an attractive drug target, relevant to combating viral hepatitis as a major public health concern. Among small molecules known to interfere with capsid assembly, the phenylpropenamides, including AT130, represent an important antiviral paradigm based on disrupting the timing of genome packaging. Here, all-atom molecular dynamics simulations of an intact AT130-bound HBV capsid reveal that the compound increases spike flexibility and improves recovery of helical secondary structure in the spike tips. Regions of the capsid-incorporated dimer that undergo correlated motion correspond to established sub-domains that pivot around the central chassis. AT130 alters patterns of correlated motion and other essential dynamics. A new conformational state of the dimer is identified, which can lead to dramatic opening of the intradimer interface and disruption of communication within the spike tip. A novel salt bridge is also discovered, which can mediate contact between the spike tip and fulcrum even in closed conformations, revealing a mechanism of direct communication across these sub-domains. Altogether, results describe a dynamical connection between the intra- and interdimer interfaces and enable mapping of allostery traversing the entire core protein dimer.