Steered molecular dynamics approach for promising drugs for influenza A virus targeting M2 channel proteins

Steered Molecular Dynamics Simulations Study on FABP4 Inhibitors

Ordinary small molecule de novo drug design is time-consuming and expensive. Recently, computational tools were employed and proved their efficacy in accelerating the overall drug design process. Molecular dynamics (MD) simulations and a derivative of MD, steered molecular dynamics (SMD), turned out to be promising rational drug design tools. In this paper, we report the first application of SMD to evaluate the binding properties of small molecules toward FABP4, considering our recent interest in inhibiting fatty acid binding protein 4 (FABP4). FABP4 inhibitors (FABP4is) are small molecules of therapeutic interest, and ongoing clinical studies indicate that they are promising for treating cancer and other diseases such as metabolic syndrome and diabetes.

Electrostatic Interactions Explain the Higher Binding Affinity of the CR3022 Antibody for SARS-CoV-2 than the 4A8 Antibody

A structural understanding of the mechanism by which antibodies bind SARS-CoV-2 at the atomic level is highly desirable as it can tell the development of more effective antibodies to treat Covid-19. Here, we use steered molecular dynamics (SMD) and coarse-grained simulations to estimate the binding affinity of the monoclonal antibodies CR3022 and 4A8 to the SARS-CoV-2 receptor-binding domain (RBD) and SARS-CoV-2 N-terminal domain (NTD). Consistent with experiments, our SMD and coarse-grained simulations both indicate that CR3022 has a higher affinity for SARS-CoV-2 RBD than 4A8 for the NTD, and the coarse-grained simulations indicate the former binds three times stronger to its respective epitope. This finding shows that CR3022 is a candidate for Covid-19 therapy and is likely a better choice than 4A8. Energetic decomposition of the interaction energies between these two complexes reveals that electrostatic interactions explain the difference in the observed binding affinity between the two complexes. This result could lead to a new approach for developing anti-Covid-19 antibodies in which good candidates must contain charged amino acids in the area of contact with the virus.

Investigation of Binding Affinity between Potential Antiviral Agents and PB2 Protein of Influenza A: Non-equilibrium Molecular Dynamics Simulation Approach

The PB2 protein of the influenza virus RNA polymerase is a major virulence determinant of influenza viruses. It binds to the cap structure at the 5' end of host mRNA to generate short capped RNA fragments that are used as primers for viral transcription named cap-snatching. A large number of the compounds were shown to bind the minimal cap-binding domain of PB2 to inhibit the cap-snatching machinery. However, their binding in the context of an extended form of the PB2 protein has remained elusive. A previous study reported some promising compounds including azaindole and hydroxymethyl azaindole, which were analyzed here to predict binding affinity to PB2 protein using the steered molecular dynamics (SMD) and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) methods. The results show that the rupture force (F_max) value of three complexes is in agreement with the binding free energy value (ΔG_bind) estimated by the MM-PBSA method, whereas for the non-equilibrium pulling work (W_pull) value a small difference between A_PB2-4 and A_PB2-12 was observed. The binding affinity results indicate the A_PB2-12 complex is more favorable than the A_PB2-4 and A_PB2-16 complexes, which means the inhibitor (12) has the potential to be further developed as anti-influenza agents in the treatment of influenza A.

Conformational Response of Influenza A M2 Transmembrane Domain to Amantadine Drug Binding at Low pH (pH 5.5)

The M2 protein from influenza A plays important roles in its viral cycle. It contains a single transmembrane helix, which oligomerizes into a homotetrameric proton channel that conducts in the low-pH environment of the host-cell endosome and Golgi apparatus, leading to virion uncoating at an early stage of infection. We studied conformational rearrangements that occur in the M2 core transmembrane domain residing on the lipid bilayer, flanked by juxtamembrane residues (M2TMD21-49 fragment), upon its interaction with amantadine drug at pH 5.5 when M2 is conductive. We also tested the role of specific mutation and lipid chain length. Electron spin resonance (ESR) spectroscopy and electron microscopy were applied to M2TMD21-49, labeled at the residue L46C with either nitroxide spin-label or Nanogold® reagent, respectively. Electron microscopy confirmed that M2TMD21-49 reconstituted into DOPC/POPS at 1:10,000 peptide-to-lipid molar ratio (P/L) either with or without amantadine, is an admixture of monomers, dimers, and tetramers, confirming our model based on a dimer intermediate in the assembly of M2TMD21-49. As reported by double electron-electron resonance (DEER), in DOPC/POPS membranes amantadine shifts oligomer equilibrium to favor tetramers, as evidenced by an increase in DEER modulation depth for P/L's ranging from 1:18,000 to 1:160. Furthermore, amantadine binding shortens the inter-spin distances (for nitroxide labels) by 5-8 Å, indicating drug induced channel closure on the C-terminal side. No such effect was observed for the thinner membrane of DLPC/DLPS, emphasizing the role of bilayer thickness. The analysis of continuous wave (cw) ESR spectra of spin-labeled L46C residue provides additional support to a more compact helix bundle in amantadine-bound M2TMD 21-49 through increased motional ordering. In contrast to wild-type M2TMD21-49, the amantadine-bound form does not exhibit noticeable conformational changes in the case of G34A mutation found in certain drug-resistant influenza strains. Thus, the inhibited M2TMD21-49 channel is a stable tetramer with a closed C-terminal exit pore. This work is aimed at contributing to the development of structure-based anti-influenza pharmaceuticals.

Investigation of the free energy profiles of amantadine and rimantadine in the AM2 binding pocket

The purpose of this work was to study the mechanism of drug resistance of M2 channel proteins by analyzing the interactions between the drugs amantadine and rimantadine and M2 channel proteins (including the wild type and the three mutants V27A, S31N, and G34A) and the drug binding pathways, by use of a computational approach. Our results showed that multiple drug-binding sites were present in the M2 channel, and the trajectory of the drugs through the M2 channel was determined. A novel method was developed to investigate of free energy profiles of the ligand-protein complexes. Our work provides a new explanation of the large amount of experimental data on drug efficacy.

Resistance-Mutation (N31) Effects on Drug Orientation and Channel Hydration in Amantadine-Bound Influenza A M2

The mechanism of amantadine binding to the S31 variant of the M2 protein of Influenza A is well understood, but the reasons behind N31 M2 amantadine insensitivity remain under investigation. Many molecular dynamics studies have evaluated the influence of amantadine position within the channel pore on its ability to inhibit proton conductance in M2, but little is known about the influence of amantadine rotational orientation. Replica-exchange umbrella sampling, steered, and classic molecular dynamics simulations were performed on amantadine in the solid-state NMR structure of S31 M2 and an N31 M2 homologue, both in the homotetramer configuration, to explore the effects of the position and tilt angle of amantadine on inhibition of the M2 channel. Steered simulations show that amantadine rotates with the amine toward the bulk water as it passes into the hydrophobic entryway lined by Val27 side chains. Results from all simulation types performed indicate that amantadine has a strong, specific orientation with the amine turned inward toward the central cavity in the S31 M2 pore but has variable orientation and a strong propensity to remain outward pointing in N31 M2. Free energy profiles from umbrella sampling, measured relative to bulk water, show amantadine binds more strongly to the S31 M2 pore by 8 kcal/mol in comparison to amantadine in the N31 pore, suggesting that it can escape more readily from the N31 channel through the Val27 secondary gate, whereas it is captured by the S31 channel in the same region. Lower water density and distribution near amantadine in S31 M2 reveal that the drug inhibits proton conductance in S31 M2 because of its inward-pointing configuration, whereas in N31 M2, amantadine forms hydrogen bonds with an N31 side chain and does not widely occlude water occupancy in any configuration. Both amantadine's weaker binding to and weaker water occlusion in N31 M2 might contribute to its inefficacy as an inhibitor of the mutant protein.