Comparative docking and molecular dynamics studies of molnupiravir (EIDD-2801): implications for novel mechanisms of action on influenza and SARS-CoV-2 protein targets.

Molnupiravir (EIDD-2801) (MLN) is an oral antiviral drug for COVID-19 treatment, being integrated into viral RNA through RNA-dependent RNA polymerase (RdRp). Upon ingestion, MLN is transformed into two active metabolites: ß-d-N4-hydroxycytidine (NHC) (EIDD-1931) in the host plasma, and EIDD-1931-triphosphate (MTP) within the host cells. However, recent studies provide increasing evidence of MLN's interactions with off-target proteins beyond the viral genome, suggesting that the complete mechanisms of action of MLN remain unclear. The aim of this study was therefore to investigate the molecular interactions of MLN in the form of NHC and MTP with the non-RNA structural components of avian influenza (hemagglutinin, neuraminidase) and SARS-CoV-2 (spike glycoprotein, Mpro, and RdRp) viruses and to elucidate whether these two metabolites possess the ability to form stable complexes with these major viral components. Molecular docking of NHC and MTP was performed using AutoDock 4.2.6 and the obtained protein-drug complexes were submitted to 200-ns molecular dynamics simulations in triplicate with subsequent free energy calculations using GROMACS. Docking scores, molecular dynamics and MM/GBSA results showed that MTP was tightly bound within the active site of SARS-CoV-2 RdRp and remained highly stable throughout the 200-ns simulations. Besides, it was also shown that NHC and MTP formed moderately-to-highly stable molecular complexes with off-target receptors hemagglutinin, neuraminidase and Mpro, but rather weak interactions with spike glycoprotein. Our computational findings suggest that NHC and MTP may directly inhibit these receptors, and propose that additional studies on the off-target effects of MLN, i.e. real-time protein binding assays, should be performed.Communicated by Ramaswamy H. Sarma.

A Family of Bisnaphthyl C2-Symmetric and Asymmetric Clefts: Synthesis, Solid-State Structure, and Calculation of the Interplanar Angle.

The synthesis of a new family of naphthalenoid C2-symmetric clefts has been realized through a four-step synthetic sequence giving three C2-symmetric clefts and a rare nonsymmetric example. Subsequently, stereoselective reduction of the carbonyl groups at C-8 and C-16 then provides cleft molecules with hydrogen bonding potential. Using single-crystal X-ray and computational analysis, the cleft angle of the dione has been determined.

Molecular dynamics simulations of structural and dynamical aspects of DNA hydration water.

Water is a remarkable liquid, both because of it is intriguing but also because of its importance. Water plays a key role on the structure and function of biological molecules, but on the other hand also the structure and dynamics of water are deeply influenced by its interactions with biological molecules, specially at low temperatures, where water's anomalies are enhanced. Here we present extensive molecular dynamics simulations of water hydrating a oligonucleotide down to very low temperatures (supercooled water), comparing four water models and analyzing the water structure and dynamics in different domains: water in the minor groove, water in the major groove and bulk water. We found that the water in the grooves is slowed down by the interactions with the nucleic acid and a hints of a dynamic transition regarding translational and orientational dynamics were found, specially for the water models TIP4P/2005 and TIP4P-Ew, which also showed the closest agreement with available experimental data. The behavior of water in such extreme conditions is relevant for the study of cryopreservation of biological tissues.

In silico analysis of the interactions of certain flavonoids with the receptor-binding domain of 2019 novel coronavirus and cellular proteases and their pharmacokinetic properties.

Coronavirus Disease 2019 (COVID-19) has infected more than thirty five million people worldwide and caused nearly 1 million deaths as of October 2020. The microorganism causing COVID-19 was named as Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2 or 2019-nCoV). The aim of this study was to investigate the interactions of twenty-three phytochemicals belonging to different flavonoid subgroups with the receptor binding domain (RBD) of the spike glycoprotein of 2019-nCoV, and cellular proteases [transmembrane serine protease 2 (TMPRSS2), cathepsin B and L (CatB/L)]. The compounds interacted more strongly with CatB and CatL than with the other proteins. Van der Waals and hydrogen bonds played an important role in the receptor-ligand interactions. As a result of RBCI (relative binding capacity index) analysis conducted to rank flavonoids in terms of their interactions with the target proteins, (-)-epicatechin gallate interacted strongly with all the proteins studied. The results obtained from molecular dynamics and molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) methods also supported this data. According to Lipinski's rule of five, (-)-epicatechin gallate showed drug-likeness properties. Although this molecule is not capable of crossing the blood-brain barrier (BBB), it was concluded that (-)-epicatechin gallate can be evaluated as a candidate molecule in drug development studies against 2019-nCoV since it was not the substrate of P-gp (P-glycoprotein), did not inhibit any of the cytochrome Ps, and did not show AMES toxicity or hepatotoxicity on eukaryotic cells.

Understanding the molecular interaction of SARS-CoV-2 spike mutants with ACE2 (angiotensin converting enzyme 2).

Covid-19 is a viral disease caused by the virus SARS-CoV-2 that spread worldwide and caused more than 4.3 million deaths. Moreover, SARS-CoV-2 still continues to evolve, and specifically the E484K, N501Y, and South Africa triple (K417N + E484K + N501Y) spike protein mutants remain as the 'escape' phenotypes. The aim of this study was to compare the interaction between the receptor binding domain (RBD) of the E484K, N501Y and South Africa triple spike variants and ACE2 with the interaction between wild-type spike RBD-ACE2 and to show whether the obtained binding affinities and conformations corraborate clinical findings. The structures of the RBDs of the E484K, N501Y and South Africa triple variants were generated with DS Studio v16 and energetically minimized using the CHARMM22 force field. Protein-protein dockings were performed in the HADDOCK server and the obtained wild-type and mutant spike-ACE2 complexes were submitted to 200-ns molecular dynamics simulations with subsequent free energy calculations using GROMACS. Based on docking binding affinities and free energy calculations the E484K, N501Y and triple mutant variants were found to interact stronger with the ACE2 than the wild-type spike. Interestingly, molecular dynamics and MM-PBSA results showed that E484K and spike triple mutant complexes were more stable than the N501Y one. Moreover, the E484K and South Africa triple mutants triggered greater conformational changes in the spike glycoprotein than N501Y. The E484K variant alone, or the combination of K417N + E484K + N501Y mutations induce significant conformational transitions in the spike glycoprotein, while increasing the spike-ACE2 binding affinity.Communicated by Ramaswamy H. Sarma.

Determination of the interaction between the receptor binding domain of 2019-nCoV spike protein, TMPRSS2, cathepsin B and cathepsin L, and glycosidic and aglycon forms of some flavonols.

The novel coronavirus (COVID-19, SARS-CoV-2) is a rapidly spreading disease with a high mortality. In this research, the interactions between specific flavonols and the 2019-nCoV receptor binding domain (RBD), transmembrane protease, serine 2 (TMPRSS2), and cathepsins (CatB and CatL) were analyzed. According to the relative binding capacity index (RBCI) calculated based on the free energy of binding and calculated inhibition constants, it was determined that robinin (ROB) and gossypetin (GOS) were the most effective flavonols on all targets. While the binding free energy of ROB with the spike glycoprotein RBD, TMPRSS2, CatB, and CatL were -5.02, -7.57, -10.10, and -6.11 kcal/mol, the values for GOS were -4.67, -5.24, -8.31, and -6.76, respectively. Furthermore, both compounds maintained their stability for at least 170 ns on respective targets in molecular dynamics simulations. The molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations also corroborated these data. Considering Lipinski's rule of five, ROB and GOS exhibited 3 (MW>500, N or O>10, NH or OH>5), and 1 (NH or OH>5) violations, respectively. Neither ROB nor GOS showed AMES toxicity or hepatotoxicity. The LD50 of these compounds in rats were 2.482 and 2.527 mol/kg, respectively. Therefore, we conclude that these compounds could be considered as alternative therapeutic agents in the treatment of COVID-19. However, the possible inhibitory effects of these compounds on cytochromes (CYPs) should be verified by in vitro or in vivo tests and their adverse effects on cellular energy metabolism should be minimized by performing molecular modifications if necessary.

The role of osmolytes in the temperature-triggered conformational transition of poly(N-vinylcaprolactam): an experimental and computational study.

Biomedical industries are widely exploring the use of thermo-responsive polymers (TRPs) in the advanced development of drug delivery and in many other pharmaceutical applications. There is a great need to investigate the use of less toxic and more (bio-)compatible TRPs employing several additives, which could modify the conformational transition behavior of TRPs in aqueous solution. To move forward in this aspect, we have chosen the less toxic bio-based polymer poly(N-vinylcaprolactam) (PVCL) and three different methylamine-based osmolytes, trimethylamine N-oxide (TMAO), betaine and sarcosine, in order to investigate their particular interactions with the polymer segments in PVCL and therefore the corresponding changes in the thermo-responsive conformational behavior. Several biophysical techniques, UV-visible spectroscopy, fluorescence spectroscopy, dynamic light scattering (DLS) and laser Raman spectroscopy, as well as classical computer simulation methods such as molecular dynamics are employed in the current work. All the studied methylamines are found to favor the hydrophobic collapse of the polymer thus stabilizing the globular state of PVCL. Sarcosine is observed to cause the maximum decrease in lower critical solution temperature (LCST) of PVCL followed by TMAO and then betaine. The differences observed in the LCST values of PVCL in the presence of these molecules can be attributed to the different polymer-osmolyte interactions. The less sterically hindered N atom in the case of sarcosine causes a significant difference in the phase transition temperature values of PVCL compared to betaine and TMAO, where the nitrogen atom is buried by three methyl groups attached to it.

Growth Mechanism of Gold Nanorods: the Effect of Tip-Surface Curvature As Revealed by Molecular Dynamics Simulations.

An understanding of the anisotropic growth mechanism of gold nanorods (AuNRs) during colloidal synthesis is critical for controlling the nanocrystal size and shape and thus has implications in tuning the properties for applications in a wide range of research and technology fields. In order to investigate the role of the cetyltrimethylammonium bromide (CTAB) coating in the anisotropic growth mechanism of AuNRs, we used molecular dynamics (MD) simulations and built a computational model that considered explicitly the effect of the curvature of the gold surface on CTAB adsorption and therefore differentiated between the CTAB arrangements on flat and curved surfaces, representing the lateral and tip facets of growing AuNRs, respectively. We verified that on a curved surface, a lower CTAB coverage density and larger intermicellar channels are generated compared to those on a flat surface. Using umbrella sampling simulations, we measured the free energy profile and verified that the environment around a curved surface corresponds to an easier migration from the solution to the gold surface for the [AuBr2]- species than does a flat surface. Long unbiased molecular dynamics simulations also corroborated the umbrella sampling results. Therefore, the [AuBr2]- diffusion through the environment of the tips is much more favorable than that in the case of lateral facets. This shows that the surface curvature is an essential component of the anisotropic growth mechanism.

Single molecule force spectroscopy of a streptomycin-binding RNA aptamer: An out-of-equilibrium molecular dynamics study.

Here, we investigate the unfolding behavior of a streptomycin-binding ribonucleic acid (RNA) aptamer under application of force in shear geometry. Using Langevin out-of-equilibrium simulations to emulate the single-molecule force spectroscopy (SMFS) experiment, we were able to understand the hierarchical unfolding process that occurs in the RNA molecule under application of stretching force and the influence of streptomycin modifying this unfolding. Subsequently, the application of the Jarzynski equality to the force profiles obtained in the pulling simulations shows that the free energies for individual systems and the difference of unfolding free energy upon streptomycin binding to the RNA free aptamer are in fair agreement with the experimental values, obtained through SMFS by Nick et al. [J. Phys. Chem. B 120, 6479 (2016)].

SERCA plays a crucial role in the toxicity of a betulinic acid derivative with potential antimalarial activity.

Malaria is one of the most significant infectious diseases that affect poor populations in tropical areas throughout the world. Plants have been shown to be a good source for the development of new antimalarial chemotherapeutic agents, as shown for the discovery of quinine and artemisinin derivatives. Our research group has been working with semisynthetic triterpene derivatives that show potential antimalarial activity toward different strains of Plasmodium falciparum by specifically modulating calcium pathways in the parasite. Promising results were obtained for nanomolar concentrations of the semisynthetic betulinic acid derivative LAFIS13 against the P. falciparum 3D7 strain in vitro, with a selectivity index of 18 compared to a mammalian cell line. Continuing these studies, we present here in vitro and in vivo toxicological evaluations of this compound, followed by docking studies with PfATP6, a sarco/endoplasmic reticulum Ca+2-ATPase (SERCA) protein. LAFIS13 showed an LD50 between 300 and 50 mg/kg, and the acute administration of 50 mg/kg (i.p.) had no negative effects on hematological, biochemical and histopathological parameters. Based on the results of the in vitro assays, LAFIS13 not exerted significant effects on coagulation parameters of human peripheral blood, but a hemolytic activity was verified at higher concentrations. According to the molecular docking study, the PfATP6 protein may be a target for LAFIS13, which corroborates its previously reported modulatory effects on calcium homeostasis in the parasite. Notably, LAFIS13 showed a higher selectivity for the mammalian SERCA protein than for PfATP6, thus impairing the selectivity between parasite and host. In summary, the direct interaction with calcium pumps and the hemolytic potential of the compound proved to be plausible mechanism of LAFIS13 toxicity.

Molecular dynamics study of the LCST transition in aqueous poly(N-n-propylacrylamide).

The breadth of technological applications of smart polymers relies on the possibility of tuning their molecular structure to respond to external stimuli. In this context, N-substituted acrylamide-based polymers are widely studied thermoresponsive polymers. Poly(N-n-propylacrylamide) (PNnPAm), which is a structural isomer of the poly(N-isopropylacrylamide) (PNIPAm) exhibits however, a lower phase transition in aqueous solution. In this work, we use all-atom molecular dynamics simulations of PNnPAm in aqueous solutions to study, from a microscopic point-of-view, the influence of chain size and concentration on the LCST of PNnPAm. Our analysis shows that the collapse of a single oligomer of PNnPAm upon heating is dependent on the chain length and corresponds to a complex interplay between hydration and intermolecular interactions. Analysis of systems with multiple chains shows an aggregation of PNnPAm chains above the LCST.

Conjugation with polyamines enhances the antitumor activity of naphthoquinones against human glioblastoma cells.

Glioblastoma multiform (GBM) is the most common and devastating type of primary brain tumor, being considered the deadliest of human cancers. In this context, extensive efforts have been undertaken to develop new drugs that exhibit both antiproliferation and antimetastasis effects on GBM. 1,4-Naphthoquinone (1,4-NQ) scaffold has been found in compounds able to inhibit important biological targets associated with cancer, which includes DNA topoisomerase, Hsp90 and monoamine oxidase. Among potential antineoplastic 1,4-NQs is the plant-derived lapachol (2-hydroxy-3-prenyl-1,4-naphthoquinone) that was found to be active against the Walker-256 carcinoma and Yoshida sarcoma. In the present study, we examined the effect of polyamine (PA)-conjugated derivatives of lapachol, nor-lapachol and lawsone on the growth and invasion of the human GBM cells. The conjugation with PA (a spermidine analog) resulted in dose-dependent and time-dependent increase of cytotoxicity of the 1,4-NQs. In addition, in-vitro inhibition of GBM cell invasion by lapachol was increased upon PA conjugation. Previous biochemical experiments indicated that these PA-1,4-NQs are capable of inhibiting DNA human topoisomerase II-α (topo2α), a major enzyme involved in maintaining DNA topology. Herein, we applied molecular docking to investigate the binding of PA-1,4-NQs to the ATPase site of topo2α. The most active molecules preferentially bind at the ATP-binding site of topo2α, which is energetically favored by the conjugation with PA. Taken together, these findings suggested that the PA-1,4-NQ conjugates might represent potential molecules in the development of new drugs in chemotherapy for malignant brain tumors.

Structure and mobility of water confined in AlPO4-54 nanotubes.

We performed molecular dynamics simulations of water confined within AlPO4-54 nanotubes. AlPO4-54 is an artificial material made of AlO4 and of PO4 in tetrahedra arranged in a periodic structure forming pores of approximately 1.3 nm in diameter. This makes AlPO4-54 an excellent candidate for practical applications, such as for water filtration and desalination. In this work, the structural and dynamical properties of the confined water are analyzed for various temperatures and water loadings. We find that the water structure is controlled by the heterogeneity of the nanopore surface with the water molecules located preferentially next to the surface of oxygens of AlPO4-54; consequently, at very low densities, water forms helicoidal structures in string-like arrangements.

Ru-Catalyzed Estragole Isomerization under Homogeneous and Ionic Liquid Biphasic Conditions.

The isomerization of estragole to trans-anethole is an important reaction and is industrially performed using an excess of NaOH or KOH in ethanol at high temperatures with very low selectivity. Simple Ru-based transition-metal complexes, under homogeneous, ionic liquid (IL)-supported (biphasic) and "solventless" conditions, can be used for this reaction. The selectivity of this reaction is more sensitive to the solvent/support used than the ligands associated with the metal catalyst. Thus, under the optimized reaction conditions, 100% conversion can be achieved in the estragole isomerization, using as little as 4 × 10-3 mol % (40 ppm) of [RuHCl(CO)(PPh3)3] in toluene, reflecting a total turnover number (TON) of 25 000 and turnover frequencies (TOFs) of up to 500 min-1 at 80 °C. Using a dimeric Ru precursor, [RuCl(µ-Cl)(η3:η3-C10H16)]2, in ethanol associated with P(OEt)3, a TON of 10 000 and a TOF of 125 min-1 are obtained with 100% conversion and 99% selectivity. These two Ru catalytic systems can be transposed to biphasic IL systems by using ionic-tagged P-ligands such as 1-(3-(diphenylphosphanyl)propyl)-2,3-dimethylimidazolium bis(trifluoromethanesulfonyl)imide immobilized in 1-(3-hydroxypropyl)-2,3-dimethylimidazolium bis(trifluoromethanesulfonyl) imide with up to 99% selectivity and almost complete estragole conversion. However, the reaction is much slower than that performed under solventless or homogeneous conditions. The use of ionic-tagged ligands significantly reduces the Ru leaching to the organic phase, compared to that in reactions performed under homogeneous conditions, where the catalytic system loses catalytic performance after the second recycling. Detailed kinetic investigations of the reaction catalyzed by [RuHCl(CO)(PPh3)3] indicate that a simplified kinetic model (a monomolecular reversible first-order reaction) is adequate for fitting the homogeneous reaction at 80 °C and under biphasic conditions. However, the kinetics of the reaction are better described if all of the elementary steps are taken into consideration, especially at 40 °C.

Organocatalytic Imidazolium Ionic Liquids H/D Exchange Catalysts.

Simple 1,2,3-trialkylimidazolium cation associated with basic anions, such as hydrogen carbonate, prolinate, and imidazolate, is an active catalyst for the H/D exchange reaction of various substrates using CDCl3 as D source, without the addition of any extra bases or metal. High deuterium incorporation (up to 49%) in acidic C-H bonds of ketone and alkyne substrates (pKa from 18.7 to 28.8) was found at room temperature. The reaction proceeds through the fast and reversible deuteration of the 2-methyl H of the imidazolium cation followed by D transfer to the substrate. The IL acts as a neutral base catalyst in which the contact ion pair is maintained in the course of the reaction. The basic active site is due to the presence of a remote basic site in the anion namely, OH of bicarbonate, NH of prolinate, and activated water in the imidazolate anion. Detailed kinetic experiments demonstrate that the reaction is first order on the substrate and pseudozero order relative to the ionic liquid, due to the fast reversible reaction involving the deuteration of the ionic liquid by the solvent.

Effects of stereochemistry and copolymerization on the LCST of PNIPAm.

Poly(N-isopropylacrylamide) (PNIPAm) is a smart polymer that presents a lower critical transition temperature (LCST) of 305 K. Interestingly, this transition point falls within the range of the human body temperature, making PNIPAm a highly suitable candidate for bio-medical applications. However, it is sometimes desirable to have a rather flexible tuning of the LCST of these polymers to further increase their range of applications. In this work, we use all-atom molecular dynamics simulations to study the LCST of PNIPAm-based (co-)polymers. We study different molecular architectures where the polymer sequences are tuned either by modifying its stereochemistry or by the co-polymerization of PNIPAm with acrylamide (Am) units. Our analysis connects global polymer conformations with the microscopic intermolecular interactions. These findings suggest that the collapse of a PNIPAm chain upon heating is dependent on the hydration structure around the monomers, which is strongly dependent on the tacticity and the presence of more hydrophilic acrylamide monomers. Our results are found to be in good agreement with the existing experimental data.

Adaptive resolution simulation of oligonucleotides.

Nucleic acids are characterized by a complex hierarchical structure and a variety of interaction mechanisms with other molecules. These features suggest the need of multiscale simulation methods in order to grasp the relevant physical properties of deoxyribonucleic acid (DNA) and RNA using in silico experiments. Here we report an implementation of a dual-resolution modeling of a DNA oligonucleotide in physiological conditions; in the presented setup only the nucleotide molecule and the solvent and ions in its proximity are described at the atomistic level; in contrast, the water molecules and ions far from the DNA are represented as computationally less expensive coarse-grained particles. Through the analysis of several structural and dynamical parameters, we show that this setup reliably reproduces the physical properties of the DNA molecule as observed in reference atomistic simulations. These results represent a first step towards a realistic multiscale modeling of nucleic acids and provide a quantitatively solid ground for their simulation using dual-resolution methods.

C-IBI: Targeting cumulative coordination within an iterative protocol to derive coarse-grained models of (multi-component) complex fluids.

We present a coarse-graining strategy that we test for aqueous mixtures. The method uses pair-wise cumulative coordination as a target function within an iterative Boltzmann inversion (IBI) like protocol. We name this method coordination iterative Boltzmann inversion (C-IBI). While the underlying coarse-grained model is still structure based and, thus, preserves pair-wise solution structure, our method also reproduces solvation thermodynamics of binary and/or ternary mixtures. Additionally, we observe much faster convergence within C-IBI compared to IBI. To validate the robustness, we apply C-IBI to study test cases of solvation thermodynamics of aqueous urea and a triglycine solvation in aqueous urea.

Why does high pressure destroy co-non-solvency of PNIPAm in aqueous methanol?

It is well known that poly(N-isopropylacrylamide) (PNIPAm) exhibits an interesting, yet puzzling, phenomenon of co-non-solvency. Co-non-solvency occurs when two competing good solvents for PNIPAm, such as water and alcohol, are mixed together. As a result, the same PNIPAm collapses within intermediate mixing ratios. This complex conformational transition is driven by preferential binding of methanol with PNIPAm. Interestingly, co-non-solvency can be destroyed when applying high hydrostatic pressures. In this work, using a large scale molecular dynamics simulation employing high pressures, we propose a microscopic picture behind the suppression of the co-non-solvency phenomenon. Based on thermodynamic and structural analysis, our results suggest that the preferential binding of methanol with PNIPAm gets partially lost at high pressures, making the background fluid reasonably homogeneous for the polymer. This is consistent with the hypothesis that the co-non-solvency phenomenon is driven by preferential binding and is not based on depletion effects.