Accelerated Molecular Dynamics Simulation for Helical Proteins Folding in Explicit Water

Preventing the Interaction between Coronaviruses Spike Protein and Angiotensin I Converting Enzyme 2: An In Silico Mechanistic Case Study on Emodin as a Potential Model Compound

Emodin, a widespread natural anthraquinone, has many biological activities including health-protective and adverse effects. Amongst beneficial effects, potential antiviral activity against coronavirus responsible for the severe acute respiratory syndrome outbreak in 2002–2003 has been described associated with the inhibition of the host cells target receptors recognition by the viral Spike protein. However, the inhibition mechanisms have not been fully characterized, hindering the rational use of emodin as a model compound to develop more effective analogues. This work investigates emodin interaction with the Spike protein to provide a mechanistic explanation of such inhibition. A 3D molecular modeling approach consisting of docking simulations, pharmacophoric analysis and molecular dynamics was used. The plausible mechanism is described as an interaction of emodin at the protein–protein interface which destabilizes the viral protein-target receptor complex. This analysis has been extended to the Spike protein of the coronavirus responsible for the current pandemic hypothesizing emodin’s functional conservation. This solid knowledge-based foothold provides a possible mechanistic rationale of the antiviral activity of emodin as a future basis for the potential development of efficient antiviral cognate compounds. Data gaps and future work on emodin-related adverse effects in parallel to its antiviral pharmacology are explored.

Ab-initio binding of barnase–barstar with DelPhiForce steered Molecular Dynamics (DFMD) approach

Receptor–ligand interactions are involved in various biological processes, therefore understanding the binding mechanism and ability to predict the binding mode are essential for many biological investigations. While many computational methods exist to predict the 3D structure of the corresponding complex provided the knowledge of the monomers, here we use the newly developed DelPhiForce steered Molecular Dynamics (DFMD) approach to model the binding of barstar to barnase to demonstrate that first-principles methods are also capable of modeling the binding. Essential component of DFMD approach is enhancing the role of long-range electrostatic interactions to provide guiding force of the monomers toward their correct binding orientation and position. Thus, it is demonstrated that the DFMD can successfully dock barstar to barnase even if the initial positions and orientations of both are completely different from the correct ones. Thus, the electrostatics provides orientational guidance along with pulling force to deliver the ligand in close proximity to the receptor.

Effects of Disulfide Bonds on Binding of Inhibitors to β-Amyloid Cleaving Enzyme 1 Decoded by Multiple Replica Accelerated Molecular Dynamics Simulations

The β-amyloid cleaving enzyme 1 (BACE1) has been thought to be an efficient target for treatment of Alzheimer's disease (AD). Deep insight into inhibitor-BACE1 binding mechanism is of significance for design of potent drugs toward BACE1. In this work, multiple replica accelerated molecular dynamics (MR-aMD) simulations, principal component (PC) analysis, and free energy landscapes were integrated to decode the effect of disulfide bonds (SSBs) in BACE1 on bindings of three inhibitors 3KO, 3KT, and 779 to BACE1. The results from cross-correlation analysis suggest that the breaking of SSBs exerts significant influence on structural flexibility and internal dynamics of inhibitor-bound BACE1. PC analysis and free energy landscapes reveal that the breaking of SSBs not only evidently induces the conformational rearrangement of BACE1 but also highly changes binding poses of three inhibitors in BACE1 and leads to more disordered binding of three inhibitors to BACE1, which is further supported by the increase in binding entropy of inhibitors to BACE1 due to the breaking of SSBs. Residue-based free energy decomposition method was utilized to evaluate contributions of separate residues to inhibitor-BACE1 binding. The results suggest that although the breaking of SSBs in BACE1 does not destroy the interaction network of inhibitors with BACE1 it changes interaction strength of some residues with inhibitors. Meanwhile, the information from residue-based free energy decomposition indicates that residues L91, S96, V130, Y132, Q134, W137, F169, I171, and I179 can be used as efficient targets of drug design toward BACE1.

Structural dynamic studies on identification of EGCG analogues for the inhibition of Human Papillomavirus E7

High risk human papillomaviruses are highly associated with the cervical carcinoma and the other genital tumors. Development of cervical cancer passes through the multistep process initiated from benign cyst to increasingly severe premalignant dysplastic lesions in an epithelium. Replication of this virus occurs in the fatal differentiating epithelium and involves in the activation of cellular DNA replication proteins. The oncoprotein E7 of human papillomavirus expressed in the lower epithelial layers constrains the cells into S-phase constructing an environment favorable for genome replication and cell proliferation. To date, no suitable drug molecules exist to treat HPV infection whereas anticipation of novel anti-HPV chemotherapies with distinctive mode of actions and identification of potential drugs are crucial to a greater extent. Hence, our present study focused on identification of compounds analogue to EGCG, a green tea molecule which is considered to be safe to use for mammalian systems towards treatment of cancer. A three dimensional similarity search on the small molecule library from natural product database using EGCG identified 11 potential small molecules based on their structural similarity. The docking strategies were implemented with acquired small molecules and identification of the key interactions between protein and compounds were carried out through binding free energy calculations. The conformational changes between the apoprotein and complexes were analyzed through simulation performed thrice demonstrating the dynamical and structural effects of the protein induced by the compounds signifying the domination. The analysis of the conformational stability provoked us to describe the features of the best identified small molecules through electronic structure calculations. Overall, our study provides the basis for structural insights of the identified potential identified small molecules and EGCG. Hence, the identified analogue of EGCG can be potent inhibitors against the HPV 16 E7 oncoprotein.

Rational Drug Design of Axl Tyrosine Kinase Type I Inhibitors as Promising Candidates Against Cancer

The high level of Axl tyrosine kinase expression in various cancer cell lines makes it an attractive target for the development of anti-cancer drugs. In this study, we carried out several sets of in silico screening for the ATP-competitive Axl kinase inhibitors based on different molecular docking protocols. The best drug-like candidates were identified, after parental structure modifications, by their highest affinity to the target protein. We found that our newly designed compound R5, a derivative of the R428 patented analog, is the most promising inhibitor of the Axl kinase according to the three molecular docking algorithms applied in the study. The molecular docking results are in agreement with the molecular dynamics simulations using the MM-PBSA/GBSA implicit solvation models, which confirm the high affinity of R5 toward the protein receptor. Additionally, the selectivity test against other kinases also reveals a high affinity of R5 toward ABL1 and Tyro3 kinases, emphasizing its promising potential for the treatment of malignant tumors.

Interactions of a multifunctional di-triazole derivative with Alzheimer's Aβ42monomer and Aβ42protofibril: a systematic molecular dynamics study

The molecular dynamics simulations results highlighted that the multi-target-directed ligand6nstabilizes the native α-helix conformation of the Aβ₄₂monomer and induces a sizable destabilization in the Aβ₄₂protofibril structure.