Biomolecular Simulations of Halogen Bonds with a GROMOS Force Field

Delivery of rapamycin by biomimetic peptide nanoparticles targeting oxidized low-density lipoprotein in atherosclerotic plaques

Drug delivery systems based on biomimetic peptide nanoparticles are steadily gaining prominence in the treatment of diverse medical conditions. This study focused on the development of peptides that depend on ligand-receptor interactions to load rapamycin (RAPA). Furthermore, a multifunctional peptide was engineered to target oxidized low-density lipoprotein (oxLDL) within atherosclerotic plaques, facilitating the localized delivery of RAPA. The interactions between peptides and RAPA/oxLDL were analyzed by simulations and experimental approaches. Results show that the main amino acid residues on the mammalian target of rapamycin that bind to RAPA are constructed as peptides (P1 and P2), which have specific interactions with RAPA and can effectively improve the loading efficiency of RAPA. The encapsulation and drug loading efficiencies of P1/P2 were 68.0/47.9% and 48.3/36.5%, respectively. In addition, the interaction force of the multifunctional peptide (P3) and oxLDL surpassed that of their interaction with human umbilical vein endothelial cells by a factor of 3.6, conclusively establishing the specific targeting of oxLDL by these nanoparticles. The encapsulation and drug loading efficiencies of P3 for RAPA were determined to be 60.2% and 41.5%. P3 can effectively load RAPA and target oxLDL within the plaque, suggesting that P3 has potential as a therapeutic agent for atherosclerotic disease.

Development and Comprehensive Benchmark of a High-Quality AMBER-Consistent Small Molecule Force Field with Broad Chemical Space Coverage for Molecular Modeling and Free Energy Calculation

Biomolecular simulations have become an essential tool in contemporary drug discovery, and molecular mechanics force fields (FFs) constitute its cornerstone. Developing a high quality and broad coverage general FF is a significant undertaking that requires substantial expert knowledge and computing resources, which is beyond the scope of general practitioners. Existing FFs originate from only a limited number of groups and organizations, and they either suffer from limited numbers of training sets, lower than desired quality because of oversimplified representations, or are costly for the molecular modeling community to access. To address these issues, in this work, we developed an AMBER-consistent small molecule FF with extensive chemical space coverage, and we provide Open Access parameters for the entire modeling community. To validate our FF, we carried out benchmarks of quantum mechanics (QM)/molecular mechanics conformer comparison and free energy perturbation calculations on several benchmark data sets. Our FF achieves a higher level of performance at reproducing QM energies and geometries than two popular open-source FFs, OpenFF2 and GAFF2. In relative binding free energy calculations for 31 protein-ligand data sets, comprising 1079 pairs of ligands, the new FF achieves an overall root-mean-square error of 1.19 kcal/mol for ΔΔG and 0.92 kcal/mol for ΔG on a subset of 463 ligands without bespoke fitting to the data sets. The results are on par with those of the leading commercial series of OPLS FFs.

Assessment of Halogen Off-Center Point-Charge Models Using Explicit Solvent Simulations

Compounds containing halogens can form halogen bonds (XBs) with biological targets such as proteins and membranes due to their anisotropic electrostatic potential. To accurately describe this anisotropy, off-center point-charge (EP) models are commonly used in force field methods, allowing the description of XBs at the molecular mechanics and molecular dynamics level. Various EP implementations have been documented in the literature, and despite being efficient in reproducing protein-ligand geometries and sampling of XBs, it is unclear how well these EP models predict experimental properties such as hydration free energies (ΔGhyd), which are often used to validate force field performance. In this work, we report the first assessment of three EP models using alchemical free energy calculations to predict ΔGhyd values. We show that describing the halogen anisotropy using some EP models can lead to a slight improvement in the prediction of the ΔGhyd when compared with the models without EP, especially for the chlorinated compounds; however, this improvement is not related to the establishment of XBs but is most likely due to the improvement of the sampling of hydrogen bonds. We also highlight the importance of the choice of the EP model, especially for the iodinated molecules, since a slight tendency to improve the prediction is observed for compounds with a larger σ-hole but significantly worse results were obtained for compounds that are weaker XB donors.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Importance of the Positively Charged σ-Hole in Crystal Engineering of Halogenated Polypeptoids

Crystalline nanosheets formed by amphiphilic block copolypeptoids with halogenated phenyl side chains were imaged at the atomic-scale using cryogenic transmission electron microscopy (cryo-TEM). In general, the polypeptoid molecules adopt V-shaped configurations in the crystalline state, and adjacent molecules can pack with one another in either parallel or antiparallel arrangements, depending on the chemical composition. The halogen bond, which can have characteristic energies ranging from 1 to 5 kcal/mol, is commensurate with the parallel configuration. However, cryo-TEM images show that chains in the halogenated crystals were in the antiparallel configuration. Molecular dynamics (MD) simulations show that positively charged σ-holes, which are characteristic of halogen atoms covalently bonded to carbon atoms, play an important role in determining crystal geometry. Parallel and antiparallel configurations exhibited similar stability in simulations when standard force fields that only account for the electronegativity of halogen atoms were used. However, including the σ-hole in the simulations resulted in a destabilization of the parallel configuration. This combination of imaging and simulation, which has played an important role in structural biology, has the potential to improve our understanding of factors that govern noncovalent interactions in synthetic materials.

Transmembrane Anion Transport Mediated by Halogen Bonds: Using Off-Center Charges

Synthetic anion transporters are promising therapeutic agents designed to emulate the specialized role of certain transmembrane proteins that maintain the ion concentration in cells. In the last few years, besides hydrogen bonds and ion pairs, halogen bonds have also been explored to promote the association between the synthetic molecule and the anion and their subsequent transport. This interaction is due to an anisotropic charge distribution on the halogen, and therefore, modeling halogen bonds is not a trivial task using classical force field methods that typically rely on point-charge models.Herein, a computational protocol capable of dealing with halogen bonds is presented. This protocol takes advantage of the addition of an off-center particle during the charge fitting procedure, and the resulting set of charges can be used along with the classical force field parameters from GAFF or GROMOS 54A7.

Optimized Halogen Atomic Radii for PBSA Calculations Using Off-Center Point Charges

In force-field methods, the usage of off-center point charges, also called extra points (EPs), is a common strategy to tackle the anisotropy of the electrostatic potential of covalently bonded halogens (X), thus allowing the description of halogen bonds (XBs) at the molecular mechanics/molecular dynamics (MM/MD) level. Diverse EP implementations exist in the literature differing on the charge sets and/or the X-EP distances. Poisson-Boltzmann and surface area (PBSA) calculations can be used to obtain solvation free energies (ΔG_solv) of small molecules, often to compute binding free energies (ΔG_bind) at the MM-PBSA level. This method depends, among other parameters, on the empirical assignment of atomic radii (PB radii). Given the multiplicity of off-center point-charge models and the lack of specific PB radii for halogens compatible with such implementations, in this work, we assessed the performance of PBSA calculations for the estimation of ΔG_solv values in water (ΔG_hyd), also conducting an optimization of the halogen PB radii (Cl, Br, and I) for each EP model. We not only expand the usage of EP models in the scope of the general AMBER force field (GAFF) but also provide the first optimized halogen PB radii in the context of the CHARMM general force field (CGenFF), thus contributing to improving the description of halogenated compounds in PBSA calculations.

Exploring the mechanism of the potent allosteric inhibitor compound2 on SHP2 WT and SHP2F285S by molecular dynamics study

Abnormal activation of Ras/MAPK signaling pathway could trigger excessive cell division. Src-homology 2 (SH2) domain-containing protein tyrosine phosphatase (SHP2) could promote Ras/MAPK activation by integrating growth factor signals. Thus, SHP2 inhibitors had become a hot topic in the treatment of cancer. SHP2^F285S, mutation in SHP2, was detected in leukemia variants. The compound 2 (3-benzyl-8-chloro-2-hydroxy-4H-benzo[4,5]thiazolo[3,2-a]pyrimidin-4-one) had been reported that it was a potent allosteric inhibitor of both SHP2 wild type (SHP2^WT) and the F285S mutant (SHP2^F285S). However, the mechanism of inhibition remained to be further discovered. Herein, molecular docking and molecular dynamic (MD) simulation were performed to explain the inhibition mechanism of compound 2 on SHP2^WT and SHP2^F285S. Overall, the molecular docking analysis revealed that compound 2 maintained the "close" structure of SHP2 protein probably by locking the C-SH2 and PTP domain. Next, post-analysis demonstrated that compound 2 might make TYR66-GLU76 of D'E-loop in N-SH2 and GLU258-LYS266 of B'C-loop, HIS458-ARG465 of P-loop, VAL497-THR507 of Q-loop in PTP domain regions tightly connect and much easier maintain "self-inhibited" conformation of SHP2^{F285S-compound2} than that of SHP2^WT-compound2. Importantly, GLU76 of D'E-loop could play a vital role in inhibition of SHP2^WT-compound2 and SHP2^{F285S-compound2}. This work provided a reliable clue to develop novel inhibitors for leukemia related to SHP2^F285S.

Understanding noncovalent bonds and their controlling forces

The fundamental underpinnings of noncovalent bonds are presented, focusing on the σ-hole interactions that are closely related to the H-bond. Different means of assessing their strength and the factors that control it are discussed. The establishment of a noncovalent bond is monitored as the two subunits are brought together, allowing the electrostatic, charge redistribution, and other effects to slowly take hold. Methods are discussed that permit prediction as to which site an approaching nucleophile will be drawn, and the maximum number of bonds around a central atom in its normal or hypervalent states is assessed. The manner in which a pair of anions can be held together despite an overall Coulombic repulsion is explained. The possibility that first-row atoms can participate in such bonds is discussed, along with the introduction of a tetrel analog of the dihydrogen bond.

Automated parameterization of quantum-mechanically derived force-fields including explicit sigma holes: A pathway to energetic and structural features of halogen bonds in gas and condensed phase

In classical molecular dynamics, general purpose atomistic force-fields (FFs) often deliver inaccurate results when dealing with halogen bonds (XBs), notwithstanding their crucial role in many fields of science, ranging from material design to drug development. Given the large dimensions of the systems of interest, it would be therefore desirable to increase the FF accuracy maintaining the simplicity of the standard Lennard-Jones (LJ) plus point charge description to avoid an excessive computational cost. A simple yet effective strategy consists in introducing a number of virtual sites able to mimic the so-called "explicit σ-hole." In this work, we present an automated FF parameterization strategy based on a global optimization of both LJ and charge parameters with respect to accurate quantum mechanical data, purposely computed for the system under investigation. As a test case, we report on two homologue series, characterized either by weak or strong XBs, namely, the di-halogenated methanes and the mono-, di-, and tri-substituted acetonitriles, taking into consideration Cl, Br, and I substituents. The resulting quantum mechanically derived FFs are validated for each compound in the gas and in the condensed phase by comparing them to general purpose and specific FFs without virtual sites and to highly accurate reference quantum mechanical data. The results strongly support the adoption of the specific FFs with virtual sites, which overcome the other investigated models in representing both gas phase energetics and the structural patterns of the liquid phase structure related to the presence of XBs.

Could Quantum Mechanical Properties Be Reflected on Classical Molecular Dynamics? The Case of Halogenated Organic Compounds of Biological Interest

Essential to understanding life, the biomolecular phenomena have been an important subject in science, therefore a necessary path to be covered to make progress in human knowledge. To fully comprehend these processes, the non-covalent interactions are the key. In this review, we discuss how specific protein-ligand interactions can be efficiently described by low computational cost methods, such as Molecular Mechanics (MM). We have taken as example the case of the halogen bonds (XB). Albeit generally weaker than the hydrogen bonds (HB), the XBs play a key role to drug design, enhancing the affinity and selectivity toward the biological target. Along with the attraction between two electronegative atoms in XBs explained by the σ-hole model, important orbital interactions, as well as relief of Pauli repulsion take place. Nonetheless, such electronic effects can be only well-described by accurate quantum chemical methods that have strong limitations dealing with supramolecular systems due to their high computational cost. To go beyond the poor description of XBs by MM methods, reparametrizing the force-fields equations can be a way to keep the balance between accuracy and computational cost. Thus, we have shown the steps to be considered when parametrizing force-fields to achieve reliable results of complex non-covalent interactions at MM level for In Silico drug design methods.

Tackling Halogenated Species with PBSA: Effect of Emulating the σ-Hole

To model halogen-bond phenomena using classical force fields, an extra point (EP) of charge is frequently introduced at a given distance from the halogen (X) to emulate the σ-hole. The resulting molecular dynamics (MD) trajectories can be used in subsequent molecular mechanics (MM) combined with Poisson-Boltzmann and surface area calculations (PBSA) to estimate protein-ligand binding free energies (Δ G_bind). While EP addition improves the MM/MD description of halogen-containing systems, its effect on the calculation of solvation free energies (Δ G_solv) using the PBSA approach is yet to be assessed. As the PBSA calculations depend, among other parameters, on the empirical assignment of radii (PB radii), a problematic issue arises, since standard halogen radii are smaller than the typical X···EP distances, thus placing the EP within the solvent dielectric. Herein, we took a common literature EP parametrization scheme, which uses X···EP = R_min and RESP charges in the context of GAFF, and performed a comprehensive study on the performance of PBSA (using three different setups) in the calculation of Δ G_solv values for 142 halogenated compounds (bearing Cl, Br, or I) for which the experimental values are known. By conducting an optimization (minimizing the error against experimental values), we provide a new optimized set of halogen PB radii, for each PBSA setup, that should be used in the context of the aforementioned scenario. A simultaneous optimization of PB radii and X···EP distances shows that a wide range of distance/radius pairs can be used without significant loss of accuracy, therefore laying the basis for expanding this halogen radii optimization strategy to other force fields and EP implementations. As ligand Δ G_solv estimation is an important term in the determination of protein-ligand Δ G_bind, this work is particularly relevant in the framework of structure-based virtual screening and related computer-aided drug design routines.

Halogen bonding in halocarbon-protein complexes and computational tools for rational drug design

Introduction: Halogens have a prominent role in drug design. Often used as a mean to improve ADME properties, they are also becoming a tool in protein-ligand recognition given their ability to form a non-covalent interaction, termed halogen bond, where halogens act as electrophilic species interacting with electron-rich partners. Rational drug design of halogen-bonding lead molecules requires an accurate description of halocarbon-protein complexes by computational tools though not all methods are able to tackle this non-covalent interaction. Areas covered: The authors present a review of computational methodologies that can be used to properly describe halogen bonds in the context of protein-ligand complexes, providing also insights on how these methods can be used in the context of computer-aided drug design. Expert opinion: Although in the last few years many computational tools, ranging from fast screening methods to the more expensive QM calculations, have been developed to tackle the halogen bonding phenomenon, they are not yet standard in the literature. This will eventually change as official software distributions are including support for halogen bonding in their methods. Tackling desolvation of halogenated species seems to be a good strategy to improve the accuracy of computational methods, that will be more commonly used prior to laboratory work in the future.

Recent studies of docking and molecular dynamics simulation for liquid-phase enantioseparations

Liquid-phase enantioseparations have been fruitfully applied in several fields of science. Various applications along with technical and theoretical advancements contributed to increase significantly the knowledge in this area. Nowadays, chromatographic techniques, in particular HPLC on chiral stationary phase, are considered as mature technologies. In the last thirty years, CE has been also recognized as one of the most versatile technique for analytical scale separation of enantiomers. Despite the huge number of papers published in these fields, understanding mechanistic details of the stereoselective interaction between selector and selectand is still an open issue, in particular for high-molecular weight chiral selectors like polysaccharide derivatives. With the ever growing improvement of computer facilities, hardware and software, computational techniques have become a basic tool in enantioseparation science. In this field, molecular docking and dynamics simulations proved to be extremely adaptable to model and visualize at molecular level the spatial proximity of interacting molecules in order to predict retention, selectivity, enantiomer elution order, and profile noncovalent interaction patterns underlying the recognition process. On this basis, topics and trends in using docking and molecular dynamics as theoretical complement of experimental LC and CE chiral separations are described herein. The basic concepts of these computational strategies and seminal studies performed over time are presented, with a specific focus on literature published between 2015 and November 2018. A systematic compilation of all published literature has not been attempted.

A Simple Model for Halogen Bond Interaction Energies

Halogen bonds are prevalent in many areas of chemistry, physics, and biology. We present a statistical model for the interaction energies of halogen-bonded systems at equilibrium based on high-accuracy ab initio benchmark calculations for a range of complexes. Remarkably, the resulting model requires only two fitted parameters, X and B—one for each molecule—and optionally the equilibrium separation, R e , between them, taking the simple form E = X B / R e n . For n = 4 , it gives negligible root-mean-squared deviations of 0.14 and 0.28 kcal mol - 1 over separate fitting and validation data sets of 60 and 74 systems, respectively. The simple model is shown to outperform some of the best density functionals for non-covalent interactions, once parameters are available, at essentially zero computational cost. Additionally, we demonstrate how it can be transferred to completely new, much larger complexes and still achieve accuracy within 0.5 kcal mol - 1 . Using a principal component analysis and symmetry-adapted perturbation theory, we further show how the model can be used to predict the physical nature of a halogen bond, providing an efficient way to gain insight into the behavior of halogen-bonded systems. This means that the model can be used to highlight cases where induction or dispersion significantly affect the underlying nature of the interaction.