Imidazole[1,5-a]pyridine derivatives as EGFR tyrosine kinase inhibitors unraveled by umbrella sampling and steered molecular dynamics simulations.

Although the use of the tyrosine kinase inhibitors (TKIs) has been proved that it can save live in a cancer treatment, the currently used drugs bring in many undesirable side-effects. Therefore, the search for new drugs and an evaluation of their efficiency are intensively carried out. Recently, a series of eighteen imidazole[1,5-a]pyridine derivatives were synthetized by us, and preliminary analyses pointed out their potential to be an important platform for pharmaceutical development owing to their promising actions as anticancer agents and enzyme (kinase, HIV-protease,…) inhibitors. In the present theoretical study, we further analyzed their efficiency in using a realistic scenario of computational drug design. Our protocol has been developed to not only observe the atomistic interaction between the EGFR protein and our 18 novel compounds using both umbrella sampling and steered molecular dynamics simulations, but also determine their absolute binding free energies. Calculated properties of the 18 novel compounds were in detail compared with those of two known drugs, erlotinib and osimertinib, currently used in cancer treatment. Inspiringly the simulation results promote three imidazole[1,5-a]pyridine derivatives as promising inhibitors into a further step of clinical trials.

Treatment of flexibility of protein backbone in simulations of protein-ligand interactions using steered molecular dynamics.

To ensure that an external force can break the interaction between a protein and a ligand, the steered molecular dynamics simulation requires a harmonic restrained potential applied to the protein backbone. A usual practice is that all or a certain number of protein's heavy atoms or Cα atoms are fixed, being restrained by a small force. This present study reveals that while fixing both either all heavy atoms and or all Cα atoms is not a good approach, while fixing a too small number of few atoms sometimes cannot prevent the protein from rotating under the influence of the bulk water layer, and the pulled molecule may smack into the wall of the active site. We found that restraining the Cα atoms under certain conditions is more relevant. Thus, we would propose an alternative solution in which only the Cα atoms of the protein at a distance larger than 1.2 nm from the ligand are restrained. A more flexible, but not too flexible, protein will be expected to lead to a more natural release of the ligand.

Novel design of amine and metal hydroxide functional group modified onto sludge biochar for arsenic removal.

This study involved novel-designed sludge biochar (SB) adsorbed for arsenic removal with lower operating costs and higher adsorption efficiency properties. Generally, biochar only relies on micropores for pollutant adsorption, but physical adsorption is not highly efficient for arsenic removal. Therefore, in order to improve the removal efficiency of arsenic by SB, diethylenetriamine (DETA) and FeCl3 were used in this study to modify the surface of SB by an immersion method. The objectives of this research are to obtain optimum operation conditions by assessing the effect of different Fe content, pH and initial concentration on adsorbing arsenic. This study is the first to use Density Functional Theory (DFT) to simulate and verify the adsorption mechanism of arsenic by SB. Results showed the presence of amine/iron oxyhydroxides functional groups greatly promoted SB surface activity and its arsenic adsorption potential. The surface area, pore volume and pore size of the SB were estimated to be 525 m2 g-1, 0.35 cm3 g-1 and 8.71 nm, respectively. The DFT model result is the same as the result of arsenic adsorption performance with high adsorption energy (-246.3 kJmol-1) and shorter bond distances (1.42 Å), indicating strong chemical adsorption between arsenic and material. The reaction mechanism is divided into four pathways, including oxidation-reduction, complexation, electrostatic adsorption and pore adsorption.

Dependence of Work on the Pulling Speed in Mechanical Ligand Unbinding.

In single-molecule force spectroscopy, the rupture force Fmax required for mechanical unfolding of a biomolecule or for pulling a ligand out of a binding site depends on the pulling speed V and, in the linear Bell-Evans regime, Fmax ∼ ln(V). Recently, it has been found that non-equilibrium work W is better than Fmax in describing relative ligand binding affinity, but the dependence of W on V remains unknown. In this paper, we developed an analytical theory showing that in the linear regime, W ∼ c1 ln(V) + c2 ln2(V), where c1 and c2 are constants. This quadratic dependence was also confirmed by all-atom steered molecular dynamics simulations of protein-ligand complexes. Although our theory was developed for ligand unbinding, it is also applicable to other processes, such as mechanical unfolding of proteins and other biomolecules, due to its universality.

Remdesivir Strongly Binds to Both RNA-Dependent RNA Polymerase and Main Protease of SARS-CoV-2: Evidence from Molecular Simulations.

The outbreak of a new coronavirus SARS-CoV-2 (severe acute respiratory syndrome-coronavirus 2) has caused a global COVID-19 (coronavirus disease 2019) pandemic, resulting in millions of infections and thousands of deaths around the world. There is currently no drug or vaccine for COVID-19, but it has been revealed that some commercially available drugs are promising, at least for treating symptoms. Among them, remdesivir, which can block the activity of RNA-dependent RNA polymerase (RdRp) in old SARS-CoV and MERS-CoV viruses, has been prescribed to COVID-19 patients in many countries. A recent experiment showed that remdesivir binds to SARS-CoV-2 with an inhibition constant of µM, but the exact target has not been reported. In this work, combining molecular docking, steered molecular dynamics, and umbrella sampling, we examined its binding affinity to two targets including the main protease (Mpro), also known as 3C-like protease, and RdRp. We showed that remdesivir binds to Mpro slightly weaker than to RdRp, and the corresponding inhibition constants, consistent with the experiment, fall to the µM range. The binding mechanisms of remdesivir to two targets differ in that the electrostatic interaction is the main force in stabilizing the RdRp-remdesivir complex, while the van der Waals interaction dominates in the Mpro-remdesivir case. Our result indicates that remdesivir can target not only RdRp but also Mpro, which can be invoked to explain why this drug is effective in treating COVID-19. We have identified residues of the target protein that make the most important contribution to binding affinity, and this information is useful for drug development for this disease.

How Good is Jarzynski's Equality for Computer-Aided Drug Design?

Accurate determination of the binding affinity of the ligand to the receptor remains a difficult problem in computer-aided drug design. Here, we study and compare the efficiency of Jarzynski's equality (JE) combined with steered molecular dynamics and the linear interaction energy (LIE) method by assessing the binding affinity of 23 small compounds to six receptors, including ß-lactamase, thrombin, factor Xa, HIV-1 protease (HIV), myeloid cell leukemia-1, and cyclin-dependent kinase 2 proteins. It was shown that Jarzynski's nonequilibrium binding free energy ΔGneqJar correlates with the available experimental data with the correlation levels R = 0.89, 0.86, 0.83, 0.80, 0.83, and 0.81 for six data sets, while for the binding free energy ΔGLIE obtained by the LIE method, we have R = 0.73, 0.80, 0.42, 0.23, 0.85, and 0.01. Therefore, JE is recommended to be used for ranking binding affinities as it provides accurate and robust results. In contrast, LIE is not as reliable as JE, and it should be used with caution, especially when it comes to new systems.

Probing the Binding Affinity by Jarzynski's Nonequilibrium Binding Free Energy and Rupture Time.

Binding affinity of a small ligand to a receptor is the important quantity in drug design, and it might be characterized by different quantities. The most popular one is the binding free energy, which can be estimated by several methods in conventional molecular dynamics simulation. So far in steered molecular dynamics (SMD), one can use either the rupture force or nonequilibrium pulling work as a measure for binding affinity. In this paper, we have shown that the nonequilibrium binding free energy Δ GneqJar, obtained by Jarzynski's equality at a finite pulling speed, has good correlation with experimental data on inhibition constants, implying that this quantity can be used as a good scoring function for binding affinity. A similar correlation has also been disclosed for binding and unbinding free energy barriers. Applying the SMD method to unbinding of 23 small compounds from the binding site of ß-lactamase protein, a bacteria-produced enzyme, we have demonstrated that the rupture or unbinding time strongly correlates with experimental data with correlation level R ≈ 0.84. As follows from the Jarzynski's equality, the rupture time depends on the unbinding barrier exponentially. We show that Δ GneqJar, the rupture time, and binding and unbinding free energy barriers are good descriptors for binding affinity. Our observation may be useful for fast screening of potential leads as the SMD simulation is not time-consuming. On the basis of nonequilibrium simulation, we disclosed that, in agreement with the experiment, the binding time is much longer than the unbinding one.

EGCG inhibits the oligomerization of amyloid beta (16-22) hexamer: Theoretical studies.

An extensive replica exchange molecular dynamics (REMD) simulation was performed to investigate the progress patterns of the inhibition of (-)-epigallocatechin-3-gallate (EGCG) on the Aß16-22 hexamer. Structural variations of the oligomers without and with EGCG were monitored and analyzed in detail. It has been found that EGCG prevents the formation of Aß oligomer through two different ways by either accelerating the Aß oligomerization or reducing the ß-content of the hexamer. It also decreases the potential "highly toxic" conformations of Aß oligomer, which is related to the conformations having high order ß-sheet sizes. Both electrostatic and van der Waals interaction energies are found to be involved to the binding process. Computed results using quantum chemical methods show that the π-π stacking is a critical factor of the interaction between EGCG and the peptides. As a result, the binding free energy of the EGCG to the Aß peptides is slightly larger than that of the curcumin.

Replica exchange molecular dynamics study of the truncated amyloid beta (11-40) trimer in solution.

Amyloid beta (Aß) oligomers are neurotoxic compounds that destroy the brain of Alzheimer's disease patients. Recent studies indicated that the trimer is one of the most cytotoxic forms of low molecular weight Aß oligomers. As there was limited information about the structure of the Aß trimer, either by experiment or by computation, we determined in this work the structure of the 3Aß11-40 oligomer for the first time using the temperature replica exchange molecular dynamics simulations in the presence of an explicit solvent. More than 20.0 µs of MD simulations were performed. The probability of the ß-content and random coil structure of the solvated trimer amounts to 42 ± 6 and 49 ± 7% which is in good agreement with experiments. Intermolecular interactions in central hydrophobic cores play a key role in stabilizing the oligomer. Intermolecular polar contacts between D23 and residues 24-29 replace the salt bridge D23-K28 to secure the loop region. The hydrophilic region of the N-terminus is maintained by the intermolecular polar crossing contacts H13A-Q15B and H13B-Q15C. The difference in the free energy of binding between the constituting monomers and the others amounts to -36 ± 8 kcal mol-1. The collision cross section of the representative structures of the trimer was computed to be 1330 ± 47 Å2, which is in good agreement with previous experiments.