Free Energies of Hydration from a Generalized Born Model and an All-Atom Force Field

Real-time structural dynamics of the ultrafast solvation process around photo-excited aqueous halides

This work investigates and describes the structural dynamics taking place following charge-transfer-to-solvent photo-abstraction of electrons from I- and Br- ions in aqueous solution following single- and 2-photon excitation at 202 nm and 400 nm, respectively. A Time-Resolved X-ray Solution Scattering (TR-XSS) approach with direct sensitivity to the structure of the surrounding solvent as the water molecules adopt a new equilibrium configuration following the electron-abstraction process is utilized to investigate the structural dynamics of solvent shell expansion and restructuring in real-time. The structural sensitivity of the scattering data enables a quantitative evaluation of competing models for the interaction between the nascent neutral species and surrounding water molecules. Taking the I0-O distance as the reaction coordinate, we find that the structural reorganization is delayed by 0.1 ps with respect to the photoexcitation and completes on a time scale of 0.5-1 ps. On longer time scales we determine from the evolution of the TR-XSS difference signal that I0: e- recombination takes place on two distinct time scales of ∼20 ps and 100 s of picoseconds. These dynamics are well captured by a simple model of diffusive evolution of the initial photo-abstracted electron population where the charge-transfer-to-solvent process gives rise to a broad distribution of electron ejection distances, a significant fraction of which are in the close vicinity of the nascent halogen atoms and recombine on short time scales.

Interaction of Norsecurinine-Type Oligomeric Alkaloids with α-Tubulin: A Molecular Docking Study

The medicinal plant Securinega virosa (Roxb ex. Willd) Baill., also known as Flueggea virosa (Roxb. ex Willd.) Royle, is commonly used in traditional medicine in Africa and Asia for the management of diverse pathologies, such as parasite infections, diabetes, and gastrointestinal diseases. Numerous alkaloids have been isolated from the twigs and leaves of the plant, notably a variety of oligomeric indolizidine alkaloids derived from the monomers securinine and norsecurinine which both display anticancer properties. The recent discovery that securinine can bind to tubulin and inhibit microtubule assembly prompted us to investigate the potential binding of two series of alkaloids, fluevirosines A-H and fluevirosinine A-J, with the tubulin dimer by means of molecular modeling. These natural products are rare high-order alkaloids with tri-, tetra-, and pentameric norsecurinine motifs. Despite their large size (up to 2500 Å3), these alkaloids can bind easily to the large drug-binding cavity (about 4800 Å3) on α-tubulin facing the β-tubulin unit. The molecular docking analysis suggests that these hydrophobic macro-alkaloids can form stable complexes with α/β-tubulin. The tubulin-binding capacity varies depending on the alkaloid size and structure. Structure-binding relationships are discussed. The docking analysis identifies the trimer fluevirosine D, tetramer fluevirosinine D, and pentamer fluevirosinine H as the most interesting tubulin ligands in the series. This study is the first to propose a molecular target for these atypical oligomeric Securinega alkaloids.

A molecular docking exploration of the large extracellular loop of tetraspanin CD81 with small molecules

Tetraspanin CD81 is a transmembrane protein used as a co-receptor by different viruses and implicated in some cancer and inflammatory diseases. The design of therapeutic small molecules targeting CD81 lags behind monoclonal antibodies and peptides but different synthetic and natural products binding to CD81 have been identified. We have investigated the interaction between synthetic compounds and CD81, considering both the cholesterol-bound full-length receptor and a truncated protein corresponding to the large extracellular loop (LEL) of the tetraspanin. They represent the closed and open conformations of the protein, respectively. Stable complexes were characterized with bi-aryl compounds (notably the quinolinone-benzothiazole 6) and atypical molecules bearing a 1-amino-boraadamantane scaffold well adapted to interact with CD81 (5a-d). In each case, the mode of binding to CD81 was analyzed, the binding sites identified and the molecular contacts determined. The narrow intra-LEL binding site of CD81 can accommodate the elongated bi-aryl 6 but not a series of isosteric compounds with a bis(bicyclic) scaffold. The bora-adamantane derivatives appeared to bind well to CD81, but essentially to the external surface of the protein loop. The binding selectivity of the compounds was assessed comparing binding to the LEL of tetraspanins CD81, CD9 and Tspan15. A net preference for CD81 over CD9 was evidenced, but the LEL of Tspan15 also provided a suitable binding site for the compounds, notably for the bora-adamantane derivatives. This work provides an aid to the identification and design of tetraspanin-binding small molecules, underlining the distinct behavior of the open and closed conformation of the protein for drug binding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00203-6.

Self-Consistent Implementation of a Solvation Free Energy Framework to Predict the Salt Solubilities of Six Alkali Halides

To assess the salt solubilities of six alkali halides in aqueous systems, we proposed a thermodynamic cycle and an efficient molecular modeling methodology. The Gibbs free energy changes for vaporization, dissociation, and dissolution were calculated using the experimental data of ionic thermodynamic properties obtained from the NBS tables. Additionally, the Marcus' and Tissandier's solvation free energy data for Li+, Na+, K+, Cl-, and Br- ions were compared with the conventional solvation free energies by substituting into our self-consistent thermodynamic cycle. Furthermore, Tissandier's absolute solvation free energy data were used as the training set to refit the Lennard-Jones parameters of OPLS-AA force field for ions. To predict salt solubilities, an assumption of a pseudo-solvent was proposed to characterize the coupling work of a solute with its environment from infinitely diluted to saturated solutions, indicating that the Gibbs energy change of solvation process is a function of ionic strength. Following the self-consistency of the cycle, the newly derived formulas were used to determine the salt solubilities by interpolating the intersection of Gibbs free energy of dissolution and the zero free energy line. The refined ion parameters can also predict the structure and thermodynamic properties of aqueous electrolyte solutions, such as densities, pair correlation functions, hydration numbers, mean activity coefficients, vapor pressures, and the radial dependences of the net charge at 298.15 K and 1 bar. Our method can be used to characterize the solid-liquid equilibria of ions or charged particles in aqueous systems. Furthermore, for highly concentrated strong electrolyte systems, it is essential to introduce accurate water models and polarizable force fields.

N-glycosylation reinforces interaction of immune checkpoint TIM-3 with a small molecule ligand

N-glycosylation of eukaryotic proteins plays roles in protein folding, trafficking, and signal transduction. The biological influence of the process is well understood, whereas the pharmacological impact of protein N-glycosylation is not well under discerned. The role of N-glycosylation on drug binding to protein has been rarely studied. We have modeled the influence of a bi-antennary N-glycan introduced at position N78 on the immune checkpoint TIM-3 (T cell immunoglobulin domain and mucin domain-containing molecule 3) on the interaction with a selective drug antagonist. The bulky N-glycan introduced at the consensus sequence Asn-Val-Thr has no influence on drug binding when the glycan adopts an extended conformation. But in a folded conformation, the glycan can interact directly with the triazoloquinazolinone derivative so as to further stabilize the drug-TIM-3 complex. The non-fucosylated glycan at position N78 markedly consolidates the drug interaction, via an additional H-bond interaction with the α3-mannose residue. It provides a gain of empirical potential energy of interaction (ΔE) of about 30 %. The presence of a more rigid fucosylated N-glycan is a little less favorable, with a gain of ΔE of about 20 %. The folded N-glycan appears to protect the ligand bound to the protein cavity, with the tricyclic core of the heterocyclic molecule sandwiched between two indole rings of tryptophan residues. Similar results were obtained when using a biantennary disialyl N-glycan with a bisecting GlcNAc residue and a tetra-antennary N-glycan. The molecular models illustrate the drug-stabilizing capacity of a bulky N-glycan positioned at a validated glycosylation site (N78 corresponding to N100 for the full-length protein). The modeling approach is useful to delineate further the role of the N-glycan of the immune checkpoint TIM-3 in interaction with small molecule ligands, and to guide the design of more potent compounds. The approach is transposable to other proteins to better comprehend the influence of N-glycans on drug-receptor interactions.

Interaction of Camptothecin Anticancer Drugs with Ribosomal Proteins L15 and L11: A Molecular Docking Study

The antitumor drug topotecan (TPT) is a potent inhibitor of topoisomerase I, triggering DNA breaks lethal for proliferating cancer cells. The mechanism is common to camptothecins SN38 (the active metabolite of irinotecan) and belotecan (BLT). Recently, TPT was shown to bind the ribosomal protein L15, inducing an antitumor immune activation independent of topoisomerase I. We have modeled the interaction of four camptothecins with RPL15 derived from the 80S human ribosome. Two potential drug-binding sites were identified at Ile135 and Phe129. SN38 can form robust RPL15 complexes at both sites, whereas BLT essentially gave stable complexes with site Ile135. The empirical energy of interaction (ΔE) for SN38 binding to RPL15 is similar to that determined for TPT binding to the topoisomerase I-DNA complex. Molecular models with the ribosomal protein L11 sensitive to topoisomerase inhibitors show that SN38 can form a robust complex at a single site (Cys25), much more stable than those with TPT and BLT. The main camptothecin structural elements implicated in the ribosomal protein interaction are the lactone moiety, the aromatic system and the 10-hydroxyl group. The study provides guidance to the design of modulators of ribosomal proteins L11 and L15, both considered anticancer targets.

Molecular Docking of Cryptoconcatones to α-Tubulin and Related Pironetin Analogues

Cryptoconcatones A-L represent a series of 12 dihydropyrone derivatives isolated from the evergreen tree Cryptocarya concinna Hance, which is well distributed in southeast Asia. The lead compound in the series, cryptoconcatone L, has revealed antiproliferative activity against cultured cancer cells but its mechanism of action remains unknown. Based on a structural analogy with the anticancer natural product pironetin, which is well known for binding covalently to α-tubulin and for functioning as a microtubule polymerization inhibitor, we investigated the interaction of cryptoconcatones with tubulin dimers using molecular docking. The α-tubulin binding capacity of each compound was quantified (through calculation of the empirical energy of interaction ΔE) and structure-binding relationships were delineated. Two compounds were found to interact with α-tubulin much more potently than pironetin: cryptoconcatones F and L. In both cases, the facile formation of a covalent bond with Cys316 was evidenced, as observed with the parent compound pironetin. A few other pironetin analogues were investigated, including spicigerolide, which is an analogue of another known α-tubulin binder. Altogether, this study points to the identification of a series of 5,6-dihydro-α-pyrones as α-tubulin-binding agents. The study contributes to a better understanding of the mechanism of action of cryptoconcatones and should help the design of analogues targeting the pironetin site of α-tubulin.

Interaction of the renin inhibitor aliskiren with the SARS-CoV-2 main protease: a molecular docking study

The renin protein is an upstream enzymatic regulator of the renin-aldosterone-angiotensin system (RAAS) essential for the maintenance of blood pressure. The angiotensin-converting enzyme-2 (ACE2) is a major component of the RAAS and a cell surface receptor exploited by the SARS-CoV-2 virus to enter host cells. A recent molecular modeling study has revealed that the direct renin peptide inhibitor remikiren can bind to the catalytic site of SARS-CoV-2 main protease (M^pro). By analogy, we postulated that the non-peptidic drug aliskiren, a more potent renin inhibitor than remikiren and a drug routinely used to treat hypertension, may also be able to interact with M^pro. An in silico comparison of the binding of the two compounds to M^pro indicates that aliskiren (ΔE = -75.9 kcal/mol) can form stable complexes with the main viral protease, binding to the active site, as remikiren (ΔE = -83.2 kcal/mol). The comparison with a panoply of 30 references compounds (mainly antiviral drugs) indicated that remikiren is a potent M^pro binder comparable to drugs like glecaprevir and pibrentasvir (ΔE = -96.5 kcal/mol). The energy of interaction (ΔE) of aliskiren with M^pro is about 10% lower than with remikiren, comparable to that calculated with drugs like velpatasvir and sofosbuvir. A model is proposed to define the drug binding site, with the best binders (including remikiren) penetrating deeply into the site, whereas the less potent binders (including aliskiren) interact more superficially with the protein.Communicated by Ramaswamy H. Sarma.

Direct Observation of Contact Ion-Pair Formation in La3+ Methanol Solution

An approach combining molecular dynamics (MD) simulations and X-ray absorption spectroscopy (XAS) has been used to carry out a comparative study about the solvation properties of dilute La(NO₃)₃ solutions in water and methanol, with the aim of elucidating the still elusive coordination of the La³⁺ ion in the latter medium. The comparison between these two systems enlightened a different behavior of the nitrate counterions in the two environments: while in water the La(NO₃)₃ salt is fully dissociated and the La³⁺ ion is coordinated by water molecules only, the nitrate anions are able to enter the metal first solvation shell to form inner-sphere complexes in methanol solution. The speciation of the formed complexes showed that the 10-fold coordination is preferential in methanol solution, where the nitrate anions coordinate the La³⁺ cations in a monodentate fashion and the methanol molecules complete the solvation shell to form an overall bicapped square antiprism geometry. This is at variance with the aqueous solution where a more balanced situation is observed between the 9- and 10-fold coordination. An experimental confirmation of the MD results was obtained by La K-edge XAS measurements carried out on 0.1 M La(NO₃)₃ solutions in the two solvents, showing the distinct presence of the nitrate counterions in the La³⁺ ion first solvation sphere of the methanol solution. The analysis of the extended X-ray absorption fine structure (EXAFS) part of the absorption spectrum collected on the methanol solution was carried out starting from the MD results and confirmed the structural arrangement observed by the simulations.

The formation of contact ion pairs between the La³⁺ ions and the nitrate anions has been demonstrated in diluted methanol solution using a combined approach using Molecular Dynamics simulations and X-ray absorption spectroscpy.

Binding of Vialinin A and p-Terphenyl Derivatives to Ubiquitin-Specific Protease 4 (USP4): A Molecular Docking Study

The para-terphenyl derivative vialinin A (Vi-A), isolated from Thelephora fungi, has been characterized as a potent inhibitor of the ubiquitin-specific protease 4 (USP4). Blockade of USP4 contributes to the anti-inflammatory and anticancer properties of the natural product. We have investigated the interaction of Vi-A with USP4 by molecular modeling, to locate the binding site (around residue V98 within the domain in USP segment) and to identify the binding process and interaction contacts. From this model, a series of 32 p-terphenyl compounds were tested as potential USP4 binders, mainly in the vialinin, terrestrin and telephantin series. We identified 11 compounds presenting a satisfactory USP4 binding capacity, including two fungal products, vialinin B and aurantiotinin A, with a more favorable empirical energy of USP4 interaction (ΔE) than the reference product Vi-A. The rare p-terphenyl aurantiotinin A, isolated from the basidiomycete T. aurantiotincta, emerged as a remarkable USP4 binder. Structure-binding relationships have been identified and discussed, to guide the future design of USP4 inhibitors based on the p-terphenyl skeleton. The docking study should help the identification of other protease inhibitors from fungus.

Interaction of panduratin A and derivatives with the SARS-CoV-2 main protease (mpro): a molecular docking study

Panduratin A (Pa-A) is a prenylated cyclohexenyl chalcone isolated from the rhizomes of the medicinal and culinary plant Boesenbergia rotunda (L.) Mansf., commonly called fingerroots. Both an ethanolic plant extract and Pa-A have shown a marked antiviral activity against the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), responsible for the COVID-19 pandemic disease. Pa-A functions as a protease inhibitor inhibiting infection of human cells by the virus. We have modeled the interaction of Pa-A, and 26 panduratin analogues with the main protease (M^pro) of SARS-CoV-2 using molecular docking. The natural product 4-hydroxypanduratin showed a higher M^pro binding capacity than Pa-A and isopanduratin A. The interaction with M^Pro of all known panduratin derivatives (Pa-A to Pa-Y) have been compared, together with more than 60 reference products. Three compounds emerged as potential robust M^Pro binders: Pa-R, Pa-V, Pa-S, with a binding capacity significantly higher than 4-OH-Pa-A and Pa-A. The empirical energy of interaction (ΔE) calculated with the best compound in the panduratin series, Pa-R bound to M^pro, surpassed that measured with the top reference protease inhibitors such a ruprintrivir, lufotrelvir, and glecaprevir. Structure-binding relationships are discussed. Compounds with a flavanone moiety (PA-R/S) are the best binders, better than those with a chromene unit (Pa-F/G). The extended molecules (such as Pa-V) exhibit good M^pro binding, but the dimeric compound Pa-Y is too long and protrudes outside the binding cavity. The work provides novel ideas to guide the design of new molecules interacting with M^pro.Communicated by Ramaswamy H. Sarma.

Molecular docking study of GSK-3β interaction with nomilin, kihadanin B, and related limonoids and triterpenes with a furyl-δ-lactone core

Glycogen synthase kinase-3β (GSK-3β) is a target enzyme considered for the treatment of multiple human diseases, from neurodegenerative pathologies to viral infections and cancers. Numerous inhibitors of GSK-3β have been discovered but thus far only a few have reached clinical trials and only one drug, tideglusib (1), has been registered. Natural products targeting GSK-3β have been identified, including the two anticancer limonoids obacunone (5) and gedunin (4), both presenting a furyl-δ-lactone core. To help identifying novel GSK-3β ligands, we have performed a molecular docking study with 15 complementary natural products bearing a furyl-δ-lactone unit (such as limonin (6) and kihadanins A (8) and B (9)) or a closely related structure (such as cedrelone (10) and nimbolide (11)). The formation of GSK-3β-binding complexes for those natural products was compared to reference GSK-3β ATP-competitive inhibitors LY2090314 (3) and AR-A014418 (2). Our in silico analysis led to the identification of two new GSK-3β-binding natural products: kihadanin B (9) and nomilin (7). The latter surpassed the reference compounds in terms of calculated empirical energy of interaction (ΔE). Nomilin (7) can possibly bind to the active site of GSK-3β, notably via the furyl-δ-lactone core and its 1-acetyl group, implicated in the protein interaction. Compound structure-binding relationships are discussed. The study should help the discovery of novel natural products targeting GSK-3β.

Spectrophotometric studies of charge-transfer complexes formed with ions N,N'-alkyldiyl-bis(pyridinium) derivatives and iodide

The charge-transfer complexes (CTC) between electron donor iodide and pyridinium dimers and monomers as acceptors have been studied spectrophotometrically in acetonitrile, DMSO and water. The dimers were: N,N'-alkyldiyl-bis(4-cyanopyridinium) and N,N'-alkyldiyl-bis(2-bromopyridinium), and the monomers were: N-alkyl-4-cyanopyridinium and N-alkyl-2-bromopyridinium, bridged by n methylene spacers. The formation constant (K_CTC), molar absorptivity (ε_CTC), fluorescence-quenching constant (K_SV) were determined. The results indicate that the stoichiometry of the CTC for dimers and monomers is 1:1 (equimolar ratio), and its formation is favored for dimers more than for monomers, especially dimers with short methylene spacers forming a "sandwiched type-complex" in which the iodide is close to the two pyridinium rings. Solvents with low polarity favored the complex, which was destroyed by the presence of water. The experimental studies were complemented with theoretical studies with quantum mechanics (QM) calculations using density functional theory (DFT) and molecular mechanics (MM) simulations with Molecular Dynamics for identify the most stable conformers in acetonitrile solution. The electronic excitations were calculated with sequential QM/MM hybrid method, showing a charge-transfer wavelength in agreement with the UV-Vis absorption spectra obtained experimentally. It confirms the "sandwiched type-complex" conformers favoring the interaction of the iodide with one pyridinium rings and simultaneously forming one, or two, hydrogen bonds with the alkyl chain and additionally, for the case of dimers, with the other pyridinium ring.

Determining the charge distribution and the direction of bond cleavage with femtosecond anisotropic x-ray liquidography

Energy, structure, and charge are fundamental quantities characterizing a molecule. Whereas the energy flow and structure change in chemical reactions are experimentally characterized, determining the atomic charges of a molecule in solution has been elusive, even for a triatomic molecule such as triiodide ion, I3-. Moreover, it remains to be answered how the charge distribution is coupled to the molecular geometry; which I-I bond, if two I-I bonds are unequal, dissociates depending on the electronic state. Here, femtosecond anisotropic x-ray solution scattering allows us to provide the following answers in addition to the overall rich structural dynamics. The analysis unravels that the negative charge of I3- is highly localized on the terminal iodine atom forming the longer bond with the central iodine atom, and the shorter I-I bond dissociates in the excited state, whereas the longer one in the ground state. We anticipate that this work may open a new avenue for studying the atomic charge distribution of molecules in solution and taking advantage of orientational information in anisotropic scattering data for solution-phase structural dynamics.

Ligand-Assisted Charge-Transfer Mechanism: The Case of CdSe/Cysteine/MoS2 Heterostructures

Molecular linkers, such as cysteine, are used to stabilize colloidal quantum dots (QDs) and anchor them. Despite the typically large molecular HOMO/LUMO gap of linkers, they can increase the quantum yield and provide an effective charge-transfer channel. Through first-principles calculations, we investigate the ligand binding and the implications for charge transfer using a prototypical CdSe-Cysteine-MoS₂ three-way heterostructure. We find that the deprotonated ligand interacts with both sides of the heterostructure, which allows for successful self-passivation of the cysteine ligand molecule and the formation of dative bonds with a greatly reduced molecular gap compared with the gas phase. This leads to the formation of a charge-transfer state that is delocalized across the ligand and can directly assist electron transfer from the conduction band of colloidal CdSe QDs to the underlying MoS₂ substrate, which is a mechanism that could extend far beyond 0D-2D hybrid systems.

Interaction of obtusilactone B and related butanolide lactones with the barrier-to-autointegration factor 1 (BAF1). A computational study

The barrier-to-autointegration factor 1 (BAF1) protein is a DNA-binding protein implicated in nuclear envelop repair and reformation after mitosis. This nuclear protein is frequently overexpressed in cancer cells and plays a role in the occurrence and development of different tumors. It is a potential therapeutic target for gastric cancer, breast cancer and other malignancies. For this reason, BAF1 inhibitors are searched. The butanolide lactone obtusilactone B (Ob-B) has been found to inhibit VRK1-dependent phosphorylation of BAF1, upon direct binding to the nuclear protein. Taking advantage of the known crystallographic structure of BAF1, we have elaborated molecular models of Ob-B bound to BAF1 to delimit the binding site and binding configuration. The long endoolefinic alkyl side chain of Ob-B extends into a small groove on the protein surface, and the adjacent exomethylene-γ-lactone moiety occupies a pocket comprising to the Ser-4 phosphorylation site of BAF1. Twenty butanolide lactones structurally close to ObB were screened for BAF1 binding. Several natural products with BAF1-binding capacity potentially superior to Ob-B were identified, including mahubanolide, kotomolide B, epilitsenolide D2, and a few other known anticancer plant natural products. Our study provides new ideas to guide the discovery and design of BAF1 inhibitors.

•

Obtusilactone B (Ob-B) is an anticancer inhibitor of VRK1-mediated BAF1 phosphorylation.

•

Molecular models of Ob-B bound to BAF1 have been constructed and the binding site determined.

•

Screening of 20 butanolide lactones led to the identification of new potential BAF1 binders.

•

Mahubanolide, kotomolide B and epilitsenolide D2 emerge as potential BAF1 inhibitors.

Molecular docking study of britannin binding to PD-L1 and related anticancer pseudoguaianolide sesquiterpene lactones

The pseudoguaianolide-type sesquiterpene lactone (SL) britannin (BRT), found in different Inula species, has been characterized as a potent anticancer agent acting via modulation of the transcription factor NFkB and the Nrf2-Keap1 signaling pathway. In addition, a BRT-induced down-regulation of the immune checkpoint PD-L1 (programmed cell death ligand 1) expressed on cancer cells has been evidenced. Here we have performed a docking analysis of the direct binding of BRT to the PD-L1 protein, both in its monomeric and dimeric state. BRT appears to form stable complexes with PD-L1, with a preference for the dimeric form, binding at the interface of the two monomers. The calculated empirical energy of interaction (ΔE) value reaches -63.1 kcal/mol for the BRT-PD-L1 dimer complex, not far from the value calculated with the reference PD-L1 ligand BMS-202 (ΔE = -73.4 kcal/mol) under identical conditions. We also studied the potential PD-L1 dimer binding of 15 pseudoguaianolide sesquiterpene lactones analogues to BRT, including helenalin, gaillardin, bigelovin, coronopilin, and others. The docking analysis predicted that the SL chamissonolide (CHM) can also form equally stable complexes with PD-L1 dimer (ΔE = -64.8 kcal/mol). Preliminary compound structure-PD-L1 binding relationships have been delineated. This computational study supports the proposed interaction of BRT with PD-L1 and provides a guidance to the design of novel PD-L1 binders incorporating a SL-like tricyclic core unit.

On the Coordination Chemistry of the lanthanum(III) Nitrate Salt in EAN/MeOH Mixtures

A thorough structural characterization of the La(NO₃)₃ salt dissolved into several mixtures of ethyl ammonium nitrate (EAN) and methanol (MeOH) with EAN molar fraction χ_EAN ranging from 0 to 1 has been carried out by combining molecular dynamics (MD) and X-ray absorption spectroscopy (XAS). The XAS and MD results show that changes take place in the La³⁺ first solvation shell when moving from pure MeOH to pure EAN. With increasing the ionic liquid content of the mixture, the La³⁺ first-shell complex progressively loses MeOH molecules to accommodate more and more nitrate anions. Except in pure EAN, the La³⁺ ion is always able to coordinate both MeOH and nitrate anions, with a ratio between the two ligands that changes continuously in the entire concentration range. When moving from pure MeOH to pure EAN, the La³⁺ first solvation shell passes from a 10-fold bicapped square antiprism geometry where all the nitrate anions act only as monodentate ligands to a 12-coordinated icosahedral structure in pure EAN where the nitrate anions bind the La³⁺ cation both in mono- and bidentate modes. The La³⁺ solvation structure formed in the MeOH/EAN mixtures shows a great adaptability to changes in the composition, allowing the system to reach the ideal compromise among all of the different interactions that take place into it.

The structural properties of the La(NO₃)₃ salt dissolved into EAN/methanol mixtures were characterized by molecular dynamics and X-ray absorption spectroscopy. The La³⁺ solvation shell undergoes significant changes with increasing the ionic liquid content of the mixture, progressively losing methanol molecules to accommodate more and more nitrate anions. The La³⁺ solvation structure shows great adaptability to composition changes, passing from a 10-fold bicapped square antiprism geometry in pure methanol to a 12-coordinated icosahedral complex in EAN.

In silico analysis of echinocandins binding to the main proteases of coronaviruses PEDV (3CLpro) and SARS-CoV-2 (Mpro)

The porcine epidemic diarrhea virus (PEDV) and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are two highly pathogenic viruses causing tremendous damages to the swine and human populations, respectively. Vaccines are available to prevent contamination and to limit dissemination of these two coronaviruses, but efficient and widely affordable treatments are needed. Recently, four natural products targeting the 3C-like protease (3CLpro) of PEDV and inhibiting replication of the virus in vitro have been identified: tomatidine, epigallocatechin-3-gallate, buddlejasaponin IVb and pneumocandin B0. We have evaluated the interaction of these compounds with 3CLpro of PEDV and with the structurally similar main protease (Mpro) of SARS-CoV-2. The molecular docking analysis indicated that the echinocandin-type lipopeptide pneumocandin B0 can generate much more stable complexes with both proteases compared to tomatidine. The empirical energy of interaction (ΔE) calculated with pneumocandin B0 bound to Mpro is extremely high, comparable to that measured with known antiviral drugs. Pneumocandin B0 and its analogue capsofungin appeared a little less adapted to interact with 3CLpro compared to Mpro. In contrast, the antifungal drug micafungin bearing an unfused tricyclic side chain, emerges as a better ligand of 3CLpro of PEDV compared to Mpro of SARS-CoV-2, based on our calculations. Collectively, the analysis underlines the benefit of echinocandin-type antifungal drugs as potential inhibitors of PEDV and SARS-CoV-2 main proteases. These clinically important antifungal natural products deserve further studies as antiviral agents.

Binding of hydroxychloroquine and chloroquine dimers to palmitoyl-protein thioesterase 1 (PPT1) and its glycosylated forms: a computational approach

The lysosomal enzyme palmitoyl-protein thioesterase 1 (PPT1) removes thioester-linked fatty acid groups from membrane-bound proteins to facilitate their proteolysis. A lack of PPT1 (due to gene mutations) causes the progressive death of cortical neurons and is responsible for infantile neural ceroid lipofuscinosis (INCL), a severe neurodegenerative disorder in children. Conversely, PPT1 is often over-expressed in cancer, and considered as a valid target to control tumor growth. Potent and selective inhibitors of PPT1 have been designed, in particular 4-amino-7-chloro-quinoline derivatives such as hydroxychloroquine (HCQ) and the dimeric analogues Lys05 and DC661. We have modeled the interaction of these three compounds with the enzyme, taking advantage of the PPT1 crystallographic structure. The molecules can fit into the palmitate site of the protein, with the dimeric compounds forming more stable complexes than the monomer. But the molecular modeling suggests that the most favorable binding sites are located outside the active site. Two sites centered on residues Met112 and Gln144 were identified, offering suitable cavities for drug binding. According to the calculated empirical energies of interaction (ΔE), the dimer DC661 forms the most stable complex at site Met112 of palmitate-bound PPT1. N-glycosylated forms of PPT1 were elaborated. Paucimannosidic glycans (M2FA and M3F) and a bulkier tetra-antennary complex glycan were introduced at asparagine residues N197, N212 and N232. These N-glycans do not impede drug binding, thus suggesting that all glycoforms of PPT1 can be targeted with these compounds.Communicated by Ramaswamy H. Sarma.

In silico analysis of the antidiabetic terpenoid acankoreagenin binding to PPARγ

Acankoreagenin (ACK) is a lupane triterpene found in several Acanthopanax and Schefflera plant species. ACK, also known as acankoreanogenin or HLEDA, bears a major structural analogy with other lupane triterpenoids such as impressic acid (IA) and the largely used phytochemical betulinic acid (BA). These compounds display marked anti-inflammatory, anti-diabetes, and anti-cancer properties. BA can form stable complexes with the peroxisome proliferator-activated receptor gamma (PPARγ). The tridimensional structure of the BA-PPARγ complex was used to perform a molecular docking analysis of the binding of ACK and IA to the protein. The 3-hydroxyl epimers (R/S) of each natural product were also modeled to examine the role of the C3-OH stereochemistry that distinguishes BA [3(S)] from ACK and AI [3(R)]. Calculations indicate that ACK can form more stable complexes with PPARγ than BA, upon insertion of the drug into the same binding pocket. The inversion of the C3-OH stereochemistry is not an obstacle for binding and the additional carboxy group of ACK at C23 position seems to reinforce the protein interaction. The 3-hydroxyl group does not play a major role in the geometry of the protein-drug complex, which is preserved between BA and ACK. Additional structure-binding relationships are provided, through the evaluation of the PPARγ binding capacity of ACK derivatives. Binding of ACK to PPARγ would account for its marked antidiabetic effect, at least partially. ACK can be used as a platform to design new antidiabetic compounds.

Influence of Solvent on the Drug-Loading Process of Amphiphilic Nanogel Star Polymers

We present an all-atom molecular dynamics study of the effect of a range of organic solvents (dichloromethane, diethyl ether, toluene, methanol, dimethyl sulfoxide, and tetrahydrofuran) on the conformations of a nanogel star polymeric nanoparticle with solvophobic and solvophilic structural elements. These nanoparticles are of particular interest for drug delivery applications. As drug loading generally takes place in an organic solvent, this work serves to provide insight into the factors controlling the early steps of that process. Our work suggests that nanoparticle conformational structure is highly sensitive to the choice of solvent, providing avenues for further study as well as predictions for both computational and experimental explorations of the drug-loading process. Our findings suggest that when used in the drug-loading process, dichloromethane, tetrahydrofuran, and toluene allow for a more extensive and increased drug-loading into the interior of nanogel star polymers of the composition studied here. In contrast, methanol is more likely to support shallow or surface loading and, consequently, faster drug release rates. Finally, diethyl ether should not work in a formulation process since none of the regions of the nanogel star polymer appear to be sufficiently solvated by it.

Enhanced Monte Carlo Methods for Modeling Proteins Including Computation of Absolute Free Energies of Binding

The generation of a complete ensemble of geometrical configurations is required to obtain reliable estimations of absolute binding free energies by alchemical free energy methods. Molecular dynamics (MD) is the most popular sampling method, but the representation of large biomolecular systems may be incomplete owing to energetic barriers that impede efficient sampling of the configurational space. Monte Carlo (MC) methods can possibly overcome this issue by adapting the attempted movement sizes to facilitate transitions between alternative local-energy minima. In this study, we present an MC statistical mechanics algorithm to explore the protein-ligand conformational space with emphasis on the motions of the protein backbone and side chains. The parameters for each MC move type were optimized to better reproduce conformational distributions of 18 dipeptides and the well-studied T4-lysozyme L99A protein. Next, the performance of the improved MC algorithms was evaluated by computing absolute free energies of binding for L99A lysozyme with benzene and seven analogs. Results for benzene with L99A lysozyme from MD and the optimized MC protocol were found to agree within 0.6 kcal/mol, while a mean unsigned error of 1.2 kcal/mol between MC results and experiment was obtained for the seven benzene analogs. Significant advantages in computation speed are also reported with MC over MD for similar extents of configurational sampling.

Computational insight into the origin of unexpected contrast in chiral markers as revealed by STM

Internal substituents can serve the double purpose of generating stereogenic centers and (potentially) being identifiable with Scanning Tunneling Microscopy (STM) in 2D self-assembled molecular layers. We investigate computationally the origin of stark contrast variations in STM images of chirally substituted self-assembled organic films. STM images of alkyl derivatives with secondary -CH₃ and -OH groups have been simulated. Density functional theory calculations reveal bias-dependent contrast reversals in the substituent regions: a lack of local density of states in the relevant energy regime results in 'dark spots' in the simulated STM images, which turn bright upon increasing the bias voltage.

Ion Pairing and Diffusion in Magnesium Electrolytes Based on Magnesium Borohydride

One obstacle to realizing a practical, rechargeable magnesium-ion battery is the development of efficient Mg electrolytes. Electrolytes based on simple Mg(BH₄)₂ salts suffer from poor salt solubility and/or low conductivity, presumably due to strong ion pairing. Understanding the molecular-scale processes occurring in these electrolytes would aid in overcoming these performance limitations. Toward this goal, the present study examines the solvation, agglomeration, and transport properties of a family of Mg electrolytes based on the Mg(BH₄)₂ salt using classical molecular dynamics. These properties were examined across five different solvents (tetrahydrofuran and the glymes G1-G4) and at four salt concentrations ranging from the dilute limit up to 0.4 M. Significant and irreversible salt agglomeration was observed in all solvents at all nondilute Mg(BH₄)₂ concentrations. The degree of clustering observed in these divalent Mg systems is much larger than that reported for electrolytes containing monovalent cations, such as Li. The salt agglomeration rate and diffusivity of Mg²⁺ were both observed to correlate with solvent self-diffusivity: electrolytes using longer- (shorter-) chain solvents had the lowest (highest) Mg²⁺ diffusivity and agglomeration rates. Incorporation of Mg²⁺ into Mg²⁺-BH₄^- clusters significantly reduces the diffusivity of Mg²⁺ by restricting displacements to localized motion within largely immobile agglomerates. Consequently, diffusion is increasingly impeded with increasing Mg(BH₄)₂ concentration. These data are consistent with the solubility limitations observed experimentally for Mg(BH₄)₂-based electrolytes and highlight the need for strategies that minimize salt agglomeration in electrolytes containing divalent cations.

Influence of Dielectric Constant on Ionic Transport in Polyether-Based Electrolytes

We use all-atom molecular dynamics simulations to study the effect of polymer polarity, as quantified by the dielectric constant, on the transport properties of lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) doped polyethers. Our results indicate that increasing the host dielectric constant leads to a decrease in ionic cluster sizes and reduction in correlated motion of oppositely charged ions. This causes the ionic conductivity to more closely approach the Nernst-Einstein limit in which ionic conductivity is only limited by the diffusivities of Li⁺ and TFSI^-. We compare our results to recent experimental observations which demonstrate similar qualitative trends in host polarity.

Harmonic Force Constants for Molecular Mechanics Force Fields via Hessian Matrix Projection

A modification to the Seminario method [ Int. J. Quantum Chem. 1996 , 60 , 1271 - 1277 ] is proposed, which derives accurate harmonic bond and angle molecular mechanics force field parameters directly from the quantum mechanical Hessian matrix. The new method reduces the average error in the reproduction of quantum mechanical normal-mode frequencies of a benchmark set of 70 molecules from 12.3% using the original method, to 6.3%. The modified Seminario method is fully automated, and all parameters are computed directly from quantum mechanical data, thereby avoiding interdependency between bond and angle parameters and other components of the force field. A complete set of bond and angle force field parameters for the 20 naturally occurring amino acids is also provided for use in the future development of protein force fields.

Modeling Halogen Bonds in Ionic Liquids: A Force Field for Imidazolium and Halo-Imidazolium Derivatives

In this work, a force field for molecular dynamics and Monte Carlo simulations of ionic liquids containing imidazolium and halo-imidazolium derivatives is presented. This force field is an extension of the well-known CL&P and OPLS-AA models and was validated by comparing predicted crystalline structures for 22 ionic liquid compounds with the corresponding data deposited at the Cambridge Structural Database. The obtained results indicate that the proposed force field extension allows the reproduction of the crystal data with an absolute average deviation lower than 2.4%. Finally, it was also established that the halogen atoms covalently bound to the studied imidazolium cations are positively charged and do not exhibit a so-called σ-hole feature. For this reason, the formation of halogen bonds in the proposed force field appears naturally from the parametrized atomic point-charge distribution, without the necessity of any extra interaction sites.

1.14*CM1A-LBCC: Localized Bond-Charge Corrected CM1A Charges for Condensed-Phase Simulations

The quality of the 1.14*CM1A and 1.20*CM5 charge models was evaluated for calculations of free energies of hydration. For a set of 426 neutral molecules, 1.14*CM1A and 1.20*CM5 yield MADs of 1.26 and 1.21 kcal/mol, respectively. The 1.14*CM1A charges, which can be readily obtained for large systems, exhibit large deviations only for a subset of functional groups. The results for these cases were systematically improved using localized bond-charge corrections (LBCC) by which offsetting adjustments are made to the partial charges for atoms in specified bond types. Only 19 LBCCs were needed to yield 1.14*CM1A-LBCC charges that reduce the errors for the 426 ΔG_hyd values to only 0.61 kcal/mol. The modified charge method was also tested in computation of heats of vaporization and densities for pure organic liquids, yielding average errors of 1.40 kcal/mol and 0.024 g/cm³, similar to those for 1.14*CM1A.

Bayesian Model Averaging for Ensemble-Based Estimates of Solvation-Free Energies

This paper applies the Bayesian Model Averaging statistical ensemble technique to estimate small molecule solvation free energies. There is a wide range of methods available for predicting solvation free energies, ranging from empirical statistical models to ab initio quantum mechanical approaches. Each of these methods is based on a set of conceptual assumptions that can affect predictive accuracy and transferability. Using an iterative statistical process, we have selected and combined solvation energy estimates using an ensemble of 17 diverse methods from the fourth Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) blind prediction study to form a single, aggregated solvation energy estimate. Methods that possess minimal or redundant information are pruned from the ensemble and the evaluation process repeats until aggregate predictive performance can no longer be improved. We show that this process results in a final aggregate estimate that outperforms all individual methods by reducing estimate errors by as much as 91% to 1.2 kcal mol^-1 accuracy. This work provides a new approach for accurate solvation free energy prediction and lays the foundation for future work on aggregate models that can balance computational cost with prediction accuracy.

A polymer dataset for accelerated property prediction and design

Emerging computation- and data-driven approaches are particularly useful for rationally designing materials with targeted properties. Generally, these approaches rely on identifying structure-property relationships by learning from a dataset of sufficiently large number of relevant materials. The learned information can then be used to predict the properties of materials not already in the dataset, thus accelerating the materials design. Herein, we develop a dataset of 1,073 polymers and related materials and make it available at http://khazana.uconn.edu/. This dataset is uniformly prepared using first-principles calculations with structures obtained either from other sources or by using structure search methods. Because the immediate target of this work is to assist the design of high dielectric constant polymers, it is initially designed to include the optimized structures, atomization energies, band gaps, and dielectric constants. It will be progressively expanded by accumulating new materials and including additional properties calculated for the optimized structures provided.

Hydration Properties and Solvent Effects for All-Atom Solutes in Polarizable Coarse-Grained Water

Due to the importance of water in chemical and biological systems, a coarse-grained representation of the solvent can greatly simplify the description of the system while retaining key thermodynamic properties of the medium. A multiscale solvation model that couples all-atom solutes and polarizable Martini coarse-grained water (AAX/CGS) is developed to reproduce free energies of hydration of organic solutes. Using Monte Carlo/free energy perturbation (MC/FEP) calculations, results from multiscale and all-atom simulations are compared. Improved accuracy is obtained with the AAX/CGS approach for hydrophobic and sulfur- or halogen-containing solutes, but larger deviations are found for polar solute molecules where hydrogen bonding is featured. Furthermore, solvent effects on conformational and tautomeric equilibria of AA solutes were investigated using AA, CG, and GB/SA solvent models. It is found that the CG solvent model can reproduce well the medium effects from experiment and AA simulations; however, the GB/SA solvent model fails in some cases. A 7-30-fold reduction in computational cost is found for the present AAX/CGS multiscale simulations compared to the AA alternative.