Exploring ligand binding pathways on proteins using hypersound-accelerated molecular dynamics

生体分子の動きを効率的に捉えるシミュレーション技術を開発 --超高周波超音波照射によってタンパク質と医薬品の結合計算を加速--. 京都大学プレスリリース. 2021-05-28.Capturing the dynamic processes of biomolecular systems in atomistic detail remains difficult despite recent experimental advances. Although molecular dynamics (MD) techniques enable atomic-level observations, simulations of “slow” biomolecular processes (with timescales longer than submilliseconds) are challenging because of current computer speed limitations. Therefore, we developed a method to accelerate MD simulations by high-frequency ultrasound perturbation. The binding events between the protein CDK2 and its small-molecule inhibitors were nearly undetectable in 100-ns conventional MD, but the method successfully accelerated their slow binding rates by up to 10–20 times. Hypersound-accelerated MD simulations revealed a variety of microscopic kinetic features of the inhibitors on the protein surface, such as the existence of different binding pathways to the active site. Moreover, the simulations allowed the estimation of the corresponding kinetic parameters and exploring other druggable pockets. This method can thus provide deeper insight into the microscopic interactions controlling biomolecular processes

MOESM1 of Molmil: a molecular viewer for the PDB and beyond

Additional file 1. Appendices; Example renderings using Molmil & Overview of alternative WebGL based molecular viewers

Binding free-energy landscapes of small molecule binder and non-binder to FMN riboswitch: All-atom molecular dynamics

A small and flexible molecule, ribocil A (non-binder) or B (binder), binds to the deep pocket of the aptamer domain of the FMN riboswitch, which is an RNA molecule. This binding was studied by mD-VcMD, which is a generalized-ensemble simulation method. Ribocil A and B are structurally similar because they are optical isomers to each other. In the initial conformation of simulation, the ligands and the aptamer were completely dissociated in explicit solvent. The aptamer–ribocil B binding was stronger than the aptamer–ribocil A binding, which agrees with experiments. The computed free-energy landscape for the aptamer–ribocil B binding was funnel-like, whereas that for the aptamer–ribocil A binding was rugged. When passing through the gate (named “front gate”) of the binding pocket, each ligand interacted with bases of the riboswitch by non-native π-π stackings, and the stackings restrained the ligand’s orientation to be advantageous to reach the binding site smoothly. When the ligands reached the binding site in the pocket, the non-native stackings were replaced by the native stackings. The ligand’s orientation restriction is discussed referring to a selection mechanism reported in an earlier work on a drug–GPCR interaction. The present simulation showed another pathway leading the ligands to the binding site. The gate (“rear gate”) for this pathway was located completely opposite to the front gate on the aptamer’s surface. However, the approach from the rear gate required overcoming a free-energy barrier regarding ligand’s rotation before reaching the binding site

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

Accurate Prediction of Complex Structure and Affinity for a Flexible Protein Receptor and Its Inhibitor

In order to predict the accurate binding configuration as well as the binding affinity for a flexible protein receptor and its inhibitor drug, enhanced sampling with multicanonical molecular dynamics (McMD) simulation and thermodynamic integration (TI) were combined as a general drug docking method. CDK2, cyclin-dependent kinase 2, is involved in the cell cycle regulation. Malfunctions in CDK2 can cause tumorigenesis, and thus it is a potential drug target. Here, we performed a long McMD simulation for docking the inhibitor CS3 to CDK2 starting from the unbound structure. Subsequently, a potential binding/unbinding pathway was given from the multicanonical ensemble, and the binding free energy was readily computed by TI along the pathway. Using this combination, the correct binding configuration of CS3 to CDK2 was obtained, and its affinity coincided well with the experimental value