Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Synthesis of new non-natural L-glycosidic flavonoid derivatives and their evaluation as inhibitors of Trypanosoma cruzi ecto-nucleoside triphosphate diphosphohydrolase 1 (TcNTPDase1)

Trypanosoma cruzi is the pathogen of Chagas disease, a neglected tropical disease that affects more than 6 million people worldwide. There are no vaccines to prevent infection, and the therapeutic arsenal is very minimal and toxic. The unique E-NTPDase of T. cruzi (TcNTPDase1) plays essential roles in adhesion and infection and is a virulence factor. Quercetin is a flavonoid with antimicrobial, antiviral, and antitumor activities. Its potential as a partial inhibitor of NTPDases has also been demonstrated. In this work, we synthesized the non-natural L-glycoside derivatives of quercetin and evaluated them as inhibitors of recombinant TcNTPDase1 (rTcNTPDase1). These compounds, and quercetin and miquelianin, a natural quercetin derivative, were also tested. Compound 16 showed the most significant inhibitory effect (94%). Quercetin, miquelianin, and compound 14 showed inhibition close to 50%. We thoroughly investigated the inhibitory effect of 16. Our data suggested a competitive inhibition with a Ki of 8.39 μM (± 0.90). To better understand the interaction of compound 16 and rTcNTPDase1, we performed molecular dynamics simulations of the enzyme and docking analyses with the compounds. Our predictions show that compound 16 binds to the enzyme's catalytic site and interacts with important residues for NTPDase activity. As an inhibitor of a critical T. cruzi enzyme, (16) could be helpful as a starting point in the developing of a future treatment for Chagas disease. Furthermore, the discovery of (16) as an inhibitor of TcNTPDase1 may open new avenues in the study and development of new inhibitors of E-NTPDases.

Binding Mechanism of Riboswitch to Natural Ligand Elucidated by McMD-Based Dynamic Docking Simulations

Flavin mononucleotide riboswitches are common among many pathogenic bacteria and are therefore considered to be an attractive target for antibiotics development. The riboswitch binds riboflavin (RBF, also known as vitamin B2), and although an experimental structure of their complex has been solved with the ligand bound deep inside the RNA molecule in a seemingly unreachable state, the binding mechanism between these molecules is not yet known. We have therefore used our Multicanonical Molecular Dynamics (McMD)-based dynamic docking protocol to analyze their binding mechanism by simulating the binding process between the riboswitch aptamer domain and the RBF, starting from the apo state of the riboswitch. Here, the refinement stage was crucial to identify the native binding configuration, as several other binding configurations were also found by McMD-based docking simulations. RBF initially binds the interface between P4 and P6 including U61 and G62, which forms a gateway where the ligand lingers until this gateway opens sufficiently to allow the ligand to pass through and slip into the hidden binding site including A48, A49, and A85.

Improving protein-ligand docking results using the Semiempirical quantum mechanics: testing on the PDBbind 2016 core set

Molecular docking techniques are routinely employed for predicting ligand binding conformations and affinities in the in silico phase of the drug design and development process. In this study, a reliable semiempirical quantum mechanics (SQM) method, PM7, was employed for geometry optimization of top-ranked poses obtained from two widely used docking programs, AutoDock4 and AutoDock Vina. The PDBbind core set (version 2016), which contains high-quality crystal protein - ligand complexes with their corresponding experimental binding affinities, was used as an initial dataset in this research. It was shown that docking pose optimization improves the accuracy of pose predictions and is very useful for the refinement of docked complexes via removing clashes between ligands and proteins. It was also demonstrated that AutoDock Vina achieves a higher sampling power than AutoDock4 in generating accurate ligand poses (RMSD ≤ 2.0 Å), while AutoDock4 exhibits a better ranking power than AutoDock Vina. Finally, a new protocol based on a combination of the results obtained from the two docking programs was proposed for structure-based virtual screening studies, which benefits from the robust sampling abilities of AutoDock Vina and the reliable ranking performance of AutoDock4.Communicated by Ramaswamy H. Sarma.

Editors' Roundup: June 2023

This commentary article represents the latest edition of the Biophysical Reviews 'Editors' Roundup' Series - a platform made available to the editorial board members of any journal with a genuine interest in promoting biophysical content. Each journal associated editor is able to submit a short description of up to five articles recently appearing in their journals with an explanation of why these articles are of interest. This edition (Vol. 15 Issue 3 June 2023) carries contributions from editorial members associated with Biophysics and Physicobiology (Biophysical Society of Japan), Biophysics (Russian Academy of Sciences), Cell Biochemistry and Biophysics (Springer), and Biophysical Reviews (IUPAB-International Union for Pure and Applied Biophysics).

Interaction of Camptothecin Anticancer Drugs with Ribosomal Proteins L15 and L11: A Molecular Docking Study

The antitumor drug topotecan (TPT) is a potent inhibitor of topoisomerase I, triggering DNA breaks lethal for proliferating cancer cells. The mechanism is common to camptothecins SN38 (the active metabolite of irinotecan) and belotecan (BLT). Recently, TPT was shown to bind the ribosomal protein L15, inducing an antitumor immune activation independent of topoisomerase I. We have modeled the interaction of four camptothecins with RPL15 derived from the 80S human ribosome. Two potential drug-binding sites were identified at Ile135 and Phe129. SN38 can form robust RPL15 complexes at both sites, whereas BLT essentially gave stable complexes with site Ile135. The empirical energy of interaction (ΔE) for SN38 binding to RPL15 is similar to that determined for TPT binding to the topoisomerase I-DNA complex. Molecular models with the ribosomal protein L11 sensitive to topoisomerase inhibitors show that SN38 can form a robust complex at a single site (Cys25), much more stable than those with TPT and BLT. The main camptothecin structural elements implicated in the ribosomal protein interaction are the lactone moiety, the aromatic system and the 10-hydroxyl group. The study provides guidance to the design of modulators of ribosomal proteins L11 and L15, both considered anticancer targets.

Computational biophysics and structural biology of proteins-a Special Issue in honor of Prof. Haruki Nakamura's 70th birthday

Receiving his initial training jointly in theoretical and applied physics at the University of Tokyo, Professor Haruki Nakamura has had a long and eventful scientific career, along the way helping to shape the way that biophysics is carried out in Japan. Concentrating his research efforts on the simulation of protein structure and function, he has, over his career arc, acted as director of the Institute for Protein Research (Osaka, Japan), director of the Protein Data Bank of Japan (PDBj), president of the Biophysical Society of Japan (BSJ), president of the Protein Science Society of Japan (PSSJ), and group leader and professor of Bioinformatics and Computational Structural Biology at Osaka University. In 2022, Prof. Haruki Nakamura turned 70 years old, and to mark this occasion, his scientific colleagues from around the world have combined their efforts to produce this Festschrift Issue of the IUPAB Biophysical Reviews journal around the theme of the computational biophysics and structural biology of proteins.