Machine learning-guided design of potent darunavir analogs targeting HIV-1 proteases: A computational approach for antiretroviral drug discovery

In the pursuit of novel antiretroviral therapies for human immunodeficiency virus type-1 (HIV-1) proteases (PRs), recent improvements in drug discovery have embraced machine learning (ML) techniques to guide the design process. This study employs ensemble learning models to identify crucial substructures as significant features for drug development. Using molecular docking techniques, a collection of 160 darunavir (DRV) analogs was designed based on these key substructures and subsequently screened using molecular docking techniques. Chemical structures with high fitness scores were selected, combined, and one-dimensional (1D) screening based on beyond Lipinski's rule of five (bRo5) and ADME (absorption, distribution, metabolism, and excretion) prediction implemented in the Combined Analog generator Tool (CAT) program. A total of 473 screened analogs were subjected to docking analysis through convolutional neural networks scoring function against both the wild-type (WT) and 12 major mutated PRs. DRV analogs with negative changes in binding free energy (ΔΔ G bind ) compared to DRV could be categorized into four attractive groups based on their interactions with the majority of vital PRs. The analysis of interaction profiles revealed that potent designed analogs, targeting both WT and mutant PRs, exhibited interactions with common key amino acid residues. This observation further confirms that the ML model-guided approach effectively identified the substructures that play a crucial role in potent analogs. It is expected to function as a powerful computational tool, offering valuable guidance in the identification of chemical substructures for synthesis and subsequent experimental testing.

FMO-guided design of darunavir analogs as HIV-1 protease inhibitors

The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.

Enhancing solubility and stability of piperine using β-cyclodextrin derivatives: computational and experimental investigations

Piperine (PP), a natural alkaloid found in black pepper, possesses significant bioactivities. However, its use in pharmaceutical applications is hindered by low water solubility and susceptibility to UV light degradation. To overcome these challenges, we investigated the potential of β-cyclodextrin (βCD) and its derivatives with dimethyl (DMβCD), hydroxy-propyl (HPβCD) and sulfobutyl-ether (SBEβCD) substitutions to enhance the solubility and stability of PP. This study employed computational and experimental approaches to examine the complexation between PP and βCDs. The results revealed the formation of two types of inclusion complexes: the P-form and M-form involving the insertion of piperidine moiety and the methylene-di-oxy-phenyl moiety, respectively. These complexes primarily rely on van der Waals interactions. Among the three derivatives, the PP/SBEβCD complex exhibited the highest stability followed by HPβCD, as attributed to maximum atom contacts and minimal solvent accessibility. Solubility studies confirmed the formation of inclusion complexes in a 1:1 ratio. Notably, the stability constant of the inclusion complex was approximately two-fold higher with SBEβCD and HPβCD compared to βCD. The DSC thermograms provided confirmation of the formation of the inclusion complex between the host and guest. These findings highlight the potential of βCD derivatives to effectively encapsulate PP, improving its solubility and presenting new opportunities for its pharmaceutical applications.Communicated by Ramaswamy H. Sarma.

Separation of the Chlorofluorocarbon (CFC) CCl2F2 from N2 in NaY Zeolite, in MIL-127(Fe) and in the two Carbon Nanotubes CNT (9,9) and CNT (11,11)

Four well-suited porous materials for the selective adsorption of the most prominent CFC, which is CCl2F2, from the air are carbon nanotubes CNT (9,9) and CNT (11,11), NaY zeolite, and the Metal Organic Framework MIL-125(Fe). The adsorption has been investigated through molecular simulations. Simulation results and theoretical considerations show that reasons for the extraordinarily high selectivity in all four cases were found to be the differences in the enthalpy of adsorption for the various adsorbed gases rather than steric reasons. The four adsorbate-adsorbent systems have been examined at different temperatures, pressures, and concentration ratios in the mixture. Among them, the carbon nanotube CNT (11,11) exhibited the highest selectivity, reaching up to 104.

Computational model for lipid binding regions in phospholipase (Ves a 1) from Vespa venom

The Thai banded tiger wasp (Vespa affinis) is a dangerous vespid species found in Southeast Asia, and its stings often result in fatalities due to the presence of lethal phospholipase A[Formula: see text], known as Vespapase or Ves a 1. Developing anti-venoms for Ves a 1 using chemical drugs, such as chemical drug guide, remains a challenging task. In this study, we screened 2056 drugs against the opening conformation of the venom using the ZINC 15 and e-Drug 3D databases. The binding free energy of the top five drug candidates complexed with Ves a 1 was calculated using 300-ns-MD trajectories. Our results revealed that voxilaprevir had a higher binding free energy at the catalytic sites than other drug candidates. Furthermore, the MD simulation results indicated that voxilaprevir formed stable conformations within the catalytic pocket. Consequently, voxilaprevir could act as a potent inhibitor, opening up avenues for the development of more effective anti-venom therapeutics for Ves a 1.

In Vitro and In Silico Study on the Molecular Encapsulation of α-Tocopherol in a Large-Ring Cyclodextrin

α-tocopherol is the physiologically most active form of vitamin E, with numerous biological activities, such as significant antioxidant activity, anticancer capabilities, and anti-aging properties. However, its low water solubility has limited its potential use in the food, cosmetic, and pharmaceutical industries. One possible strategy for addressing this issue is the use of a supramolecular complex with large-ring cyclodextrins (LR-CDs). In this study, the phase solubility of the CD26/α-tocopherol complex was investigated to assess the possible ratios between host and guest in the solution phase. Next, the host-guest association of the CD26/α-tocopherol complex at different ratios of 1:2, 1:4, 1:6, 2:1, 4:1, and 6:1 was studied by all-atom molecular dynamics (MD) simulations. At 1:2 ratio, two α-tocopherol units interact spontaneously with CD26, forming an inclusion complex, as supported by the experimental data. In the 2:1 ratio, a single α-tocopherol unit was encapsulated by two CD26 molecules. In comparison, increasing the number of α-tocopherol or CD26 molecules above two led to self-aggregation and consequently limited the solubility of α-tocopherol. The computational and experimental results indicate that a 1:2 ratio could be the most suitable stoichiometry to use in the CD26/α-tocopherol complex to improve α-tocopherol solubility and stability in inclusion complex formation.

Exceptionally high selectivity in the separation of light hydrocarbons by adsorption on MIL-127(Fe) and on a (9,9) carbon nanotube

Porous solids with channel sizes that are not much above the size of small hydrocarbons can yield extremely large adsorption selectivity. Our Grand Canonical Monte-Carlo simulations indicate exceptionally high selectivity for the separation of methane, ethane and propane from natural gas. At 250 K the C₃H₈/CH₄ separation on MIL-127 at low pressure has a selectivity of more than 1000 and the C₃H₈/CH₄ separation on CNT (9,9) is even above 2000. This is due to the strong molecule lattice interaction in narrow channels which leads to large enthalpies of adsorption. The Arrhenius law for the Henry coefficients is analysed in order to show that the effect is due to this enthalpy rather than to steric reasons.

Identification of repurposing therapeutics toward SARS-CoV-2 main protease by virtual screening

SARS-CoV-2 causes the current global pandemic coronavirus disease 2019. Widely-available effective drugs could be a critical factor in halting the pandemic. The main protease (3CL^pro) plays a vital role in viral replication; therefore, it is of great interest to find inhibitors for this enzyme. We applied the combination of virtual screening based on molecular docking derived from the crystal structure of the peptidomimetic inhibitors (N3, 13b, and 11a), and experimental verification revealed FDA-approved drugs that could inhibit the 3CL^pro of SARS-CoV-2. Three drugs were selected using the binding energy criteria and subsequently performed the 3CL^pro inhibition by enzyme-based assay. In addition, six common drugs were also chosen to study the 3CL^pro inhibition. Among these compounds, lapatinib showed high efficiency of 3CL^pro inhibition (IC₅₀ value of 35 ± 1 μM and K_i of 23 ± 1 μM). The binding behavior of lapatinib against 3CL^pro was elucidated by molecular dynamics simulations. This drug could well bind with 3CL^pro residues in the five subsites S1’, S1, S2, S3, and S4. Moreover, lapatinib’s key chemical pharmacophore features toward SAR-CoV-2 3CL^pro shared important HBD and HBA with potent peptidomimetic inhibitors. The rational design of lapatinib was subsequently carried out using the obtained results. Our discovery provides an effective repurposed drug and its newly designed analogs to inhibit SARS-CoV-2 3CL^pro.

Halogenated Baicalein as a Promising Antiviral Agent toward SARS-CoV-2 Main Protease

The coronavirus disease pandemic is a constant reminder that global citizens are in imminent danger of exposure to emerging infectious diseases. Therefore, developing a technique for inhibitor discovery is essential for effective drug design. Herein, we proposed fragment molecular orbital (FMO)-based virtual screening to predict the molecular binding energy of potential severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease inhibitors. The integration of quantum mechanical approaches and trajectory analysis from a microsecond molecular dynamics simulation was used to identify potential inhibitors. We identified brominated baicalein as a potent inhibitor of the SARS-CoV-2 main protease and confirmed its inhibitory activity in an in vitro assay. Brominated baicalein did not demonstrate significant toxicity in either in vitro or in vivo studies. The pair interaction energy from FMO-RIMP2/PCM and inhibitory constants based on the protease enzyme assay suggested that the brominated baicalein could be further developed into novel SARS-CoV-2 protease inhibitors.

In Silico Elucidation of Potent Inhibitors and Rational Drug Design against SARS-CoV-2 Papain-like Protease

Global public health has been a critical problem by the sudden increase of the COVID-19 outbreak. The papain-like protease (PLpro) of SARS-CoV-2 is a key promising target for antiviral drug development since it plays a pivotal role in viral replication and innate immunity. Here, we employed the all-atom molecular dynamics (MD) simulations and binding free energy calculations based on MM-PB(GB)SA and SIE methods to elucidate and compare the binding behaviors of five inhibitors derived from peptidomimetic inhibitors (VIR250 and VIR251) and naphthalene-based inhibitors (GRL-0617, compound 3, and compound Y96) against SARS-CoV-2 PLpro. The obtained results revealed that all inhibitors interacting within the PLpro active site are mostly driven by vdW interactions, and the hydrogen bond formation in residues G163 and G271 with peptidomimetics and the Q269 residue with naphthalene-based inhibitors was essential for stabilizing the protein-ligand complexes. Among the five studied inhibitors, VIR250 exhibited the most binding efficiency with SARS-CoV-2 PLpro, and thus, it was chosen for the rational drug design. Based on the computationally designed ligand-protein complexes, the replacement of aromatic rings including heteroatoms (e.g., thiazolopyridine) at the P2 and P4 sites could help to improve the inhibitor-binding efficiency. Furthermore, the hydrophobic interactions with residues at P1-P3 sites can be increased by enlarging the nonpolar moieties (e.g., ethene) at the N-terminal of VIR250. We expect that the structural data obtained will contribute to the development of new PLpro inhibitors with more inhibitory potency for COVID-19.

Adsorption and the Chemical Reaction N2O4 ↔ 2NO2 in the Presence of N2 in a Gas Phase Connected with a Carbon Nanotube

The paper shows, by molecular simulations, that a CNT (9,9) carbon nanotube allows very efficient separation of nitrogen oxides (NO _x ) from N₂, that has in good approximation properties of the complete air mixture. Gibbs ensemble Monte Carlo simulations are used to describe the adsorption. The permanent chemical reaction between N₂O₄ and NO₂, which occurs simultaneously to adsorption, is treated by the reactive Monte Carlo simulation. A very high selectivity has been found. For a low pressure and at T = 298 K, an adsorption/reaction selectivity between NO _x and N₂ can reach values up to 3 × 10³.

Structural Basis of Specific Glucoimidazole and Mannoimidazole Binding by Os3BGlu7

β-Glucosidases and β-mannosidases hydrolyze substrates that differ only in the epimer of the nonreducing terminal sugar moiety, but most such enzymes show a strong preference for one activity or the other. Rice Os3BGlu7 and Os7BGlu26 β-glycosidases show a less strong preference, but Os3BGlu7 and Os7BGlu26 prefer glucosides and mannosides, respectively. Previous studies of crystal structures with glucoimidazole (GIm) and mannoimidazole (MIm) complexes and metadynamic simulations suggested that Os7BGlu26 hydrolyzes mannosides via the B_2,5 transition state (TS) conformation preferred for mannosides and glucosides via their preferred ⁴H₃/⁴E TS conformation. However, MIm is weakly bound by both enzymes. In the present study, we found that MIm was not bound in the active site of crystallized Os3BGlu7, but GIm was tightly bound in the -1 subsite in a ⁴H₃/⁴E conformation via hydrogen bonds with the surrounding residues. One-microsecond molecular dynamics simulations showed that GIm was stably bound in the Os3BGlu7 active site and the glycone-binding site with little distortion. In contrast, MIm initialized in the B_2,5 conformation rapidly relaxed to a E₃/⁴H₃ conformation and moved out into a position in the entrance of the active site, where it bound more stably despite making fewer interactions. The lack of MIm binding in the glycone site in protein crystals and simulations implies that the energy required to distort MIm to the B_2,5 conformation for optimal active site residue interactions is sufficient to offset the energy of those interactions in Os3BGlu7. This balance between distortion and binding energy may also provide a rationale for glucosidase versus mannosidase specificity in plant β-glycosidases.

Combined Adsorption and Reaction in the Ternary Mixture N2, N2O4, NO2 on MIL-127 Examined by Computer Simulations

A high selectivity of NO _x over N₂ (simulating air) is found in silico when studying the adsorption of the ternary mixture N₂O₄/NO₂/N₂ on the metal-organic framework MIL-127(Fe) by molecular simulations under consideration of the recombination reaction N₂O₄ ↔ 2NO₂. The number of N atoms in nitrogen oxides NO _x and that in N₂ is used to define a selectivity of the combined adsorption and chemical recombination that can reach values of about 1000.

Delivery of Alpha-Mangostin Using Cyclodextrins through a Biological Membrane: Molecular Dynamics Simulation

α-Mangostin (MGS) exhibits various pharmacological activities, including antioxidant, anticancer, antibacterial, and anti-inflammatory properties. However, its low water solubility is the major obstacle for its use in pharmaceutical applications. To increase the water solubility of MGS, complex formation with beta-cyclodextrins (βCDs), particularly with the native βCD and/or its derivative 2,6-dimethyl-β-CD (DMβCD) is a promising technique. Although there have been several reports on the adsorption of βCDs on the lipid bilayer, the release of the MGS/βCDs inclusion complex through the biological membrane remains unclear. In this present study, the release the MGS from the two different βCDs (βCD and DMβCD) across the lipid bilayer was investigated. Firstly, the adsorption of the free MGS, free βCDs, and inclusion complex formation was studied by conventional molecular dynamics simulation. The MGS in complex with those two βCDs was able to spontaneously release free MGS into the inner membrane. However, both MGS and DMβCD molecules potentially permeated into the deeper region of the interior membrane, whereas βCD only adsorbed at the outer membrane surface. The interaction between secondary rim of βCD and the 1-palmitoeyl-2-oleoyl-glycero-3-phosphocholine (POPC) phosphate groups showed the highest number of hydrogen bonds (up to 14) corresponding to the favorable location of βCD on the POPC membrane. Additionally, the findings suggested that electrostatic energy was the main driving force for βCD adsorption on the POPC membrane, while van der Waals interactions played a predominant role in DMβCD adsorption. The release profile of MGS from the βCDs pocket across the lipid bilayer exhibited two energy minima along the reaction coordinate associated with the permeation of the MGS molecule into the deeper region of the POPC membrane.

Why Are Lopinavir and Ritonavir Effective against the Newly Emerged Coronavirus 2019? Atomistic Insights into the Inhibitory Mechanisms

Since the emergence of a novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported from Wuhan, China, neither a specific vaccine nor an antiviral drug against SARS-CoV-2 has become available. However, a combination of two HIV-1 protease inhibitors, lopinavir and ritonavir, has been found to be effective against SARS-CoV, and both drugs could bind well to the SARS-CoV 3C-like protease (SARS-CoV 3CLpro). In this work, molecular complexation between each inhibitor and SARS-CoV-2 3CLpro was studied using all-atom molecular dynamics simulations, free energy calculations, and pair interaction energy analyses based on MM/PB(GB)SA and FMO-MP2/PCM/6-31G* methods. Both anti-HIV drugs interacted well with the residues at the active site of SARS-CoV-2 3CLpro. Ritonavir showed a somewhat higher number atomic contacts, a somewhat higher binding efficiency, and a somewhat higher number of key binding residues compared to lopinavir, which correspond with the slightly lower water accessibility at the 3CLpro active site. In addition, only ritonavir could interact with the oxyanion hole residues N142 and G143 via the formation of two hydrogen bonds. The interactions in terms of electrostatics, dispersion, and charge transfer played an important role in the drug binding. The obtained results demonstrated how repurposed anti-HIV drugs could be used to combat COVID-19.

In silico study directed towards identification of the key structural features of GyrB inhibitors targeting MTB DNA gyrase: HQSAR, CoMSIA and molecular dynamics simulations

Mycobacterium tuberculosis DNA gyrase subunit B (GyrB) has been identified as a promising target for rational drug design against fluoroquinolone drug-resistant tuberculosis. In this study, we attempted to identify the key structural feature for highly potent GyrB inhibitors through 2D-QSAR using HQSAR, 3D-QSAR using CoMSIA and molecular dynamics (MD) simulations approaches on a series of thiazole urea core derivatives. The best HQSAR and CoMSIA models based on IC₅₀ and MIC displayed the structural basis required for good activity against both GyrB enzyme and mycobacterial cell. MD simulations and binding free energy analysis using MM-GBSA and waterswap calculations revealed that the urea core of inhibitors has the strongest interaction with Asp79 via hydrogen bond interactions. In addition, cation-pi interaction and hydrophobic interactions of the R₂ substituent with Arg82 and Arg141 help to enhance the binding affinity in the GyrB ATPase binding site. Thus, the present study provides crucial structural features and a structural concept for rational design of novel DNA gyrase inhibitors with improved biological activities against both enzyme and mycobacterial cell, and with good pharmacokinetic properties and drug safety profiles.

Lymphocyte proliferation induced by high-affinity peptides for HLA-B*51:01 in Behçet's uveitis

Several proteins have been proposed as candidate auto-antigens in the pathogenesis of Behçet’s disease (BD). In this study, we aimed to confirm the cellular responses to candidate peptide autoantigens with high affinity for the HLA-B*51:01 molecule using computerized binding predictions and molecular dynamics simulations. We identified two new candidate peptides (HSP65PD, derived from heat shock protein-65, and B51PD, derived from HLA-B*51:01) with high-affinity to the HLA-B*51:01 binding pocket using the Immune Epitope Database for Major Histocompatibility Complex-I Binding Prediction and molecular dynamics simulations. The peptide-induced proliferation of lymphocytes from patients with BD, sarcoidosis, Vogt–Koyanagi–Harada disease (VKH) with panuveitis, systemic scleroderma (SSc) without uveitis, and healthy controls (HC) was investigated using the bromodeoxyuridine assay. The proliferative response of leukocytes to HSP65PD was significantly higher in BD (SI 1.92 ± 0.65) than that in sarcoidosis (SI 1.38 ± 0.46), VKH (SI 1.40 ± 0.33), SSc (SI 1.32 ± 0.31), and HC (SI 1.27 ± 0.28) (P = 0.0004, P = 0.0007, P < 0.0001, P < 0.0001, respectively, Mann-Whitney’s U-test). The proliferative response of leukocytes to B51PD was also higher in BD than that in sarcoidosis, VKH, SSc without uveitis, and HC, whereas no significant differences were observed among the five groups in response to a control peptide derived from topoisomerase 1. A significantly higher response to HPS65PD and B51PD was observed in the HLA-B*51:01-positive patients with BD than in the HLA-B*51:01-negative patients. In conclusion, two peptides that had high affinity to HLA-B*51:01 in computerized binding prediction showed significantly higher response in HLA-B*51:01-positive patients with BD, indicating the usefulness of computerized simulations for identifying autoreactive peptides to HLAs.

Atomistic mechanisms underlying the activation of the G protein-coupled sweet receptor heterodimer by sugar alcohol recognition

The human T1R2-T1R3 sweet taste receptor (STR) plays an important role in recognizing various low-molecular-weight sweet-tasting sugars and proteins, resulting in the release of intracellular heterotrimeric G protein that in turn leads to the sweet taste perception. Xylitol and sorbitol, which are naturally occurring sugar alcohols (polyols) found in many fruits and vegetables, exhibit the potential caries-reducing effect and are widely used for diabetic patients as low-calorie sweeteners. In the present study, computational tools were applied to investigate the structural details of binary complexes formed between these two polyols and the T1R2-T1R3 heterodimeric STR. Principal component analysis revealed that the Venus flytrap domain (VFD) of T1R2 monomer was adapted by the induced-fit mechanism to accommodate the focused polyols, in which α-helical residues 233-268 moved significantly closer to stabilize ligands. This finding likely suggested that these structural transformations might be the important mechanisms underlying polyols-STR recognitions. The calculated free energies also supported the VFD of T1R2 monomer as the preferential binding site for such polyols, rather than T1R3 region, in accord with the lower number of accessible water molecules in the T1R2 pocket. The E302 amino acid residue in T1R2 was found to be the important recognition residue for polyols binding through a strongly formed hydrogen bond. Additionally, the binding affinity of xylitol toward the T1R2 monomer was significantly higher than that of sorbitol, making it a sweeter tasting molecule.

Interactions of HLA-DR and Topoisomerase I Epitope Modulated Genetic Risk for Systemic Sclerosis

The association of systemic sclerosis with anti-Topoisomerase 1 antibody (ATASSc) with specific alleles of human leukocyte antigen (HLA)-DR has been observed among various ethnics. The anti-Topoisomerase 1 antibody is a common autoantibody in SSc with diffuse cutaneous scleroderma, which is one of the clinical subtypes of SSc. On the other hand, an immunodominant peptide of topoisomerase 1 (Top1) self-protein (residues 349-368) was reported to have strong association with ATASSc. In this study, molecular dynamics simulation was performed on the complexes of Top1 peptide with various HLA-DR subtypes divided into ATASSc-associated alleles (HLA-DRB1*08:02, HLA-DRB1*11:01 and HLA-DRB1*11:04), suspected allele (HLA-DRB5*01:02), and non-associated allele (HLA-DRB1*01:01). The unique interaction for each system was compared to the others in terms of dynamical behaviors, binding free energies and solvation effects. Our results showed that three HLA-DR/Top1 complexes of ATASSc association mostly exhibited high protein stability and increased binding efficiency without solvent interruption, in contrast to non-association. The suspected case (HLA-DRB5*01:02) binds Top1 as strongly as the ATASSc association case, which implied a highly possible risk for ATASSc development. This finding might support ATASSc development mechanism leading to a guideline for the treatment and avoidance of pathogens like Top1 self-peptide risk for ATASSc.

The role of conserved QXG and binding affinity of S23G & S26G receptors on avian H5, swine H1 and human H1 of influenza A virus hemagglutinin

Outbreaks of avian, human and swine influenza are a serious concern for public health. In the reproductive cycle of the influenza virus, hemagglutinin (HA) is the primary protein responsible for binding to glycan receptor sites on the host cell surface. MD simulations of avian H5, swine H1 and human H1 complexed with S23G and S26G receptors were performed to study the role of key residues on the receptor conformational behaviors, hydrogen bond formation, binding free energy and residue-wise energy contribution. The obtained results indicated that the relative energies of swH1_S23G and swH1_S26G were found to be close to each other (3.1 kcal/mol) while the relative energies of AvH5 and HuH1 were found to be significantly different (11.1 ± 6.8 and 29.0 ± 8.2 kcal/mol for AvH5 and HuH1, respectively).

Effect of D168V mutation in NS3/4A HCV protease on susceptibilities of faldaprevir and danoprevir

Hepatitis C virus (HCV) is a serious cause of liver inflammation, cirrhosis and the development of hepatocellular carcinoma. Its NS3/4A serine protease functions to cleave a specific peptide bond, which is an important step in HCV replication. Thus the NS3/4A protease has become one of the main drug-targets in the design and development of anti-HCV agents. Unfortunately, high mutation rates in HCV have been reported due to the lack of RNA proofreading activity resulting in drug resistance. Herein, all-atom molecular dynamics simulations were employed to understand and illustrate the effects of the NS3/4A D168V mutation on faldaprevir (FDV) and danoprevir (DNV) binding efficiency. The D168V mutation was shown to interrupt the hydrogen bonding network of Q80R155D168R123 embedded in the extended S2 and partial S4 subsites of the NS3 protein and as a result the R123 side chain was displaced and moved out from the binding pocket. By means of MM/PBSA and MM/GBSA binding free energy calculations, the FDV and DNV binding affinities were shown to be significantly reduced by ∼10-15 kcal mol^-1 and ∼4-9 kcal mol^-1 relative to the wild-type complexes, respectively, which somewhat agrees with the experimental resistance folds.

Susceptibility of inhibitors against 3C protease of coxsackievirus A16 and enterovirus A71 causing hand, foot and mouth disease: A molecular dynamics study

Hand foot and mouth disease (HFMD) epidemic has occurred in many countries. Coxsackievirus A16 (CV-A16) and Enterovirus A71 (EV-A71) are the main causes of HFMD. Up to now, there are no anti-HFMD drugs available. Rupintrivir, a broad-spectrum inhibitor, is a drug candidate for HFMD treatment, while other HFMD inhibitors designed from several studies have a relatively low efficiency. Therefore, in this work we aim to study the binding mechanisms of rupintrivir and a peptidic α,β-unsaturated ethyl ester (SG85) against both CV-A16 and EV-A71 3C proteases (3Cpro) using all-atoms molecular dynamics simulation. The obtained results indicate that SG85 shows a stronger binding affinity than rupintrivir against CV-A16. Both inhibitors exhibit a comparable affinity against EV-A71 3Cpro. The molecular information of the binding of the two inhibitors to the proteases will be elucidated. Thus, it is implied that these two compounds may be used as leads for further anti-HFMD drug design and development.

Role of R292K mutation in influenza H7N9 neuraminidase toward oseltamivir susceptibility: MD and MM/PB(GB)SA study

The H7N9 avian influenza virus is a novel re-assortment from at least four different strains of virus. Neuraminidase, which is a glycoprotein on the surface membrane, has been the target for drug treatment. However, some H7N9 strains that have been isolated from patient after drug treatment have a R292K mutation in neuraminidase. This substitution was found to facilitate drug resistance using protein- and virus- assays, in particular it gave a high resistance to the most commonly used drug, oseltamivir. The aim of this research is to understand the source of oseltamivir resistance using MD simulations and the MM/PB(GB)SA binding free energy approaches. Both methods can predict the reduced susceptibility of oseltamivir in good agreement to the IC ₅₀ binding energy, although MM/GBSA underestimates this prediction compared to the MM/PBSA calculation. Electrostatic interaction is the main contribution for oseltamivir binding in terms of both interaction and solvation. We found that the source of the drug resistance is a decrease in the binding interaction combined with the reduction of the dehydration penalty. The smaller K292 mutated residue has a larger binding pocket cavity compared to the wild-type resulting in the loss of drug carboxylate-K292 hydrogen bonding and an increased accessibility for water molecules around the K292 mutated residue. In addition, oseltamivir does not bind well to the R292K mutant complex as shown by the high degree of fluctuation in ligand RMSD during the simulation and the change in angular distribution of bulky side chain groups.

Protein-protein interactions between SWCNT/chitosan/EGF and EGF receptor: a model of drug delivery system

Epidermal growth factor (EGF) was used as the targeting ligand to enhance the specificity of a cancer drug delivery system (DDS) via its specific interaction with the EGF receptor (EGFR) that is overexpressed on the surface of some cancer cells. To investigate the intermolecular interaction and binding affinity between the EGF-conjugated DDS and the EGFR, 50 ns molecular dynamics simulations were performed on the complex of tethered EGFR and EGF linked to single-wall carbon nanotube (SWCNT) through a biopolymer chitosan wrapping the tube outer surface (EGFR·EGF-CS-SWCNT-Drug complex), and compared to the EGFR·EGF complex and free EGFR. The binding pattern of the EGF-CS-SWCNT-Drug complex to the EGFR was broadly comparable to that for EGF, but the binding affinity of the EGF-CS-SWCNT-Drug complex was predicted to be somewhat better than that for EGF alone. Additionally, the chitosan chain could prevent undesired interactions of SWCNT at the binding pocket region. Therefore, EGF connected to SWCNT via a chitosan linker is a seemingly good formulation for developing a smart DDS served as part of an alternative cancer therapy.

A 3D-RISM/RISM study of the oseltamivir binding efficiency with the wild-type and resistance-associated mutant forms of the viral influenza B neuraminidase

The binding affinity of oseltamivir to the influenza B neuraminidase and to its variants with three single substitutions, E119G, R152K, and D198N, is investigated by the MM/3D-RISM method. The binding affinity or the binding free energy of ligand to receptor was found to be determined by a subtle balance of two major contributions that largely cancel out each other: the ligand-receptor interactions and the dehydration free energy. The theoretical results of the binding affinity of the drug to the mutants reproduced the observed trend in the resistivity, measured by IC50 ; the high-level resistance of E119G and R152K, and the low-level resistance of D198N. For E119G and R152K, reduction of the direct drug-target interaction, especially at the mutated residue, is the main source of high-level oseltamivir resistance. This phenomenon, however, is not found in the D198N strain, which is located in the framework of the active-site.

Binding specificity of polypeptide substrates in NS2B/NS3pro serine protease of dengue virus type 2: A molecular dynamics Study

The pathogenic dengue virus (DV) is a growing global threat, particularly in South East Asia, for which there is no specific treatment available. The virus possesses a two-component (NS2B/NS3) serine protease that cleaves the viral precursor proteins. Here, we performed molecular dynamics simulations of the NS2B/NS3 protease complexes with six peptide substrates (capsid, intNS3, 2A/2B, 4B/5, 3/4A and 2B/3 containing the proteolytic site between P(1) and P(1)' subsites) of DV type 2 to compare the specificity of the protein-substrate binding recognition. Although all substrates were in the active conformation for cleavage reaction by NS2B/NS3 protease, their binding strength was somewhat different. The simulated results of intermolecular hydrogen bonds and decomposition energies suggested that among the ten substrate residues (P(5)-P(5)') the P(1) and P(2) subsites play a major role in the binding with the focused protease. The arginine residue at these two subsites was found to be specific preferential binding at the active site with a stabilization energy of <-10 kcal mol(-1). Besides, the P(3), P(1)', P(2)' and P(4)' subsites showed a less contribution in binding interaction (<-2 kcal mol(-1)). The catalytic water was detected nearby the carbonyl oxygen of the P(1) reacting center of the capsid, intNS3, 2A/2B and 4B/5 peptides. These results led to the order of absolute binding free energy (ΔGbind) between these substrates and the NS2B/NS3 protease ranked as capsid>intNS3>2A/2B>4B/5>3/4A>2B/3 in a relative correspondence with previous experimentally derived values.

High-level QM/MM calculations support the concerted mechanism for Michael addition and covalent complex formation in thymidylate synthase

Thymidylate synthase (TS) is a promising cancer target, due to its crucial function in thymine synthesis. It performs the reductive methylation of 2'-deoxyuridine-5'-phosphate (dUMP) to thymidine-5'-phosphate (dTMP), using N-5,10-methylene-5,6,7,8-tetrahydrofolate (mTHF) as a cofactor. After the formation of the dUMP/mTHF/TS noncovalent complex, and subsequent conformational activation, this complex has been proposed to react via nucleophilic attack (Michael addition) by Cys146, followed by methylene-bridge formation to generate the ternary covalent intermediate. Herein, QM/MM (B3LYP-D/6-31+G(d)-CHARMM27) methods are used to model the formation of the ternary covalent intermediate. A two-dimensional potential energy surface reveals that the methylene-bridged intermediate is formed via a concerted mechanism, as indicated by a single transition state on the minimum energy pathway and the absence of a stable enolate intermediate. A range of different QM methods (B3LYP, MP2 and SCS-MP2, and different basis sets) are tested for the calculation of the activation energy barrier for the formation of the methylene-bridged intermediate. We test convergence of the QM/MM results with respect to size of the QM region. Inclusion of Arg166, which interacts with the nucleophilic thiolate, in the QM region is important for reliable results; the MM model apparently does not reproduce energies for distortion of the guanidinium side chain correctly. The spin component scaled-Møller-Plessett perturbation theory (SCS-MP2) approach was shown to be in best agreement (within 1.1 kcal/mol) while the results obtained with MP2 and B3LYP also yielded acceptable values (deviating by less than 3 kcal/mol) compared with the barrier derived from experiment. Our results indicate that using a dispersion-corrected DFT method, or a QM method with an accurate treatment of electron correlation, increases the agreement between the calculated and experimental activation energy barriers, compared with the semiempirical AM1 method. These calculations provide important insight into the reaction mechanism of TS and may be useful in the design of new TS inhibitors.

Cluster Formation in a Concentrated Lithium-Liquid Ammonia Solution. A Monte Carlo Study

Results of a Monte Carlo study of a lithium-liquid ammonia solution at 240 K are reported. The basic cube contained 135 Li+ and 1025 NH3 . With an experimental density of 0.554 g/cm3 a side-length of 37.89 Å resulted. The pseudopotential theory is employed, which permits the exclusion of the electrons from an explicit consideration. The structure of the solution is described by various site-site radial distribution functions. The six ammonia molecules in the first solvation shell of the lithium ion are arranged octahedrally. Clusters are formed which consist almost exclusively of two solvated Li+ which have simultaneously either one ammonia molecule or an octahedral edge or an octahedral plane in common. About one third of the ammonia molecules belong to the bulk phase.

Three-body Effects in Calcium(II)-ammonia Solutions: Molecular Dynamics Simulations

Molecular dynamics simulations have been performed with and without three-body corrections at an average temperature of 240 K using a flexible ammonia model. The system consists of one calcium ion and 215 ammonia molecules. The calcium(II)-ammonia interactions were newly developed, based on ab initio calculations with a basis set of double zeta quality. The role of three-body interactions on the structural and dynamical properties of the solution has been investigated. The presence of three-body corrections leads to the reduction of the first shell coordination number of Ca(II) in liquid ammonia from 9 to 8, the increase of the size of the solvation shell by 0.33 A and the disappearance of the sec-ond solvation shell.

Effect of silver nanoparticles on rice (Oryza sativa L. cv. KDML 105) seed germination and seedling growth

With the advances in nanotechnology, silver nanoparticles (AgNPs) have been applied in many industries, increasing their potential exposure level in the environment, yet their environmental safety remains poorly evaluated. The possible effects of different sized AgNPs (20, 30-60, 70-120 and 150nm diameter) on jasmine rice, Oryza sativa L. cv. KDML 105, were investigated at different concentrations (0.1, 1, 10, 100 and 1000mg/L) upon seed germination and seedling growth. The results revealed that the level of seed germination and subsequent growth of those seedlings that germinated were both decreased with increasing sizes and concentrations of AgNPs. Based on the analysis of AgNPs accumulation in plant tissues, it implied that the higher uptake was found when the seeds were treated with the smaller AgNPs, 20nm diameter AgNPs, but it was trapped in the roots rather than transported to the leaves. These resulted in the less negative effects on seedling growth, when compared to the seed soaking with the larger AgNPs with 150nm diameter. The negative effects of AgNPs were supported by leaf cell deformation when rice seeds were treated with 150-nm-diameter AgNP at the concentration of 10 or 100mg/L during seed germination. These results further strengthen our understanding of environmental safety information with respect to nanomaterials.

Rational design of InhA inhibitors in the class of diphenyl ether derivatives as potential anti-tubercular agents using molecular dynamics simulations

A series of diphenyl ether derivatives were developed and showed promising potency for inhibiting InhA, an essential enoyl acyl carrier protein reductase involved in mycolic acid biosynthesis, leading to the lysis of Mycobacterium tuberculosis. To understand the structural basis of diphenyl ether derivatives for designing more potent inhibitors, molecular dynamics (MD) simulations were performed. Based on the obtained results, the dynamic behaviour in terms of flexibility, binding free energy, binding energy decomposition, conformation, and the inhibitor-enzyme interaction of diphenyl ether inhibitors were elucidated. Phe149, Tyr158, Met161, Met199, Val203 and NAD+ are the key residues for binding of diphenyl ether inhibitors in the InhA binding pocket. Our results could provide the structural concept to design new diphenyl ether inhibitors with better enzyme inhibitory activity against M. tuberculosis InhA. The present work facilitates the design of new and potentially more effective anti-tuberculosis agents.

Key binding and susceptibility of NS3/4A serine protease inhibitors against hepatitis C virus

Hepatitis C virus (HCV) causes an infectious disease that manifests itself as liver inflammation, cirrhosis, and can lead to the development of liver cancer. Its NS3/4A serine protease is a potent target for drug design and development since it is responsible for cleavage of the scissile peptide bonds in the polyprotein important for the HCV life cycle. Herein, the ligand-target interactions and the binding free energy of the four current NS3/4A inhibitors (boceprevir, telaprevir, danoprevir, and BI201335) were investigated by all-atom molecular dynamics simulations with three different initial atomic velocities. The per-residue free energy decomposition suggests that the key residues involved in inhibitor binding were residues 41-43, 57, 81, 136-139, 155-159, and 168 in the NS3 domain. The van der Waals interactions yielded the main driving force for inhibitor binding at the protease active site for the cleavage reaction. In addition, the highest number of hydrogen bonds was formed at the reactive P1 site of the four studied inhibitors. Although the hydrogen bond patterns of these inhibitors were different, their P3 site was most likely to be recognized by the A157 backbone. Both molecular mechanic (MM)/Poisson-Boltzmann surface area and MM/generalized Born surface area approaches predicted the relative binding affinities of the four inhibitors in a somewhat similar trend to their experimentally derived biological activities.

Effect of halogen substitutions on dUMP to stability of thymidylate synthase/dUMP/mTHF ternary complex using molecular dynamics simulation

The stability of the thymidylate synthase (TS)/2-deoxyuridine-5-monophosphate (dUMP)/5,10-methylene-5,6,7,8-tetrahydrofolate (mTHF) ternary complex formation and Michael addition are considered as important steps that are involved in the inhibition mechanism of the anticancer prodrug 5-fluorouracil (5-FU). Here, the effect of three different halogen substitutions on the C-5 position of the dUMP (XdUMPs = FdUMP, CldUMP, and BrdUMP), the normal substrate, on the stability of the TS/dUMP and TS/dUMP/mTHF binary and ternary complexes, respectively, was investigated via molecular dynamics simulation. The simulated results revealed that the stability of all the systems was substantially increased by mTHF binding to the catalytic pocket. In the ternary complex, a much greater stabilization of the dUMP and XdUMPs through electrostatic interactions, including charge-charge and hydrogen bond interactions, was found compared to mTHF. An additional unique hydrogen bond between the substituted fluorine of FdUMP and the hydroxyl group of the TS Y94 residue was observed in both the binary and ternary complexes. The distance between the S(-) atom of the TS C146 residue and the C6 atom of dUMP, at <4 Å in all systems, suggested that a Michael addition with the formation of a S-C6 covalent bond potentially occurred, although the hydrogen atom on C6 of dUMP is substituted by a halogen atom. The MM/PBSA binding free energy revealed the significant role of the bridging waters around the ligands in the increased binding affinity (∼10 kcal/mol) of dUMP/XdUMP, either alone or together with mTHF, toward TS. The order of the averaged binding affinity in the ternary systems was found to be CldUMP ≈ FdUMP > dUMP > BrdUMP, suggesting that CldUMP could be a potent candidate TS inhibitor, the same as FdUMP (the metabolite form of 5-FU).