Fully Flexible Docking via Reaction-Coordinate-Independent Molecular Dynamics Simulations

Similarities and Differences in Ligand Binding to Protein and RNA Targets: The Case of Riboflavin

It is nowadays clear that RNA molecules can play active roles in several biological processes. As a result, an increasing number of RNAs are gradually being identified as potentially druggable targets. In particular, noncoding RNAs can adopt highly organized conformations that are suitable for drug binding. However, RNAs are still considered challenging targets due to their complex structural dynamics and high charge density. Thus, elucidating relevant features of drug-RNA binding is fundamental for advancing drug discovery. Here, by using Molecular Dynamics simulations, we compare key features of ligand binding to proteins with those observed in RNA. Specifically, we explore similarities and differences in terms of (i) conformational flexibility of the target, (ii) electrostatic contribution to binding free energy, and (iii) water and ligand dynamics. As a test case, we examine binding of the same ligand, namely riboflavin, to protein and RNA targets, specifically the riboflavin (RF) kinase and flavin mononucleotide (FMN) riboswitch. The FMN riboswitch exhibited enhanced fluctuations and explored a wider conformational space, compared to the protein target, underscoring the importance of RNA flexibility in ligand binding. Conversely, a similar electrostatic contribution to the binding free energy of riboflavin was found. Finally, greater stability of water molecules was observed in the FMN riboswitch compared to the RF kinase, possibly due to the different shape and polarity of the pockets.

Synthesis, and Insecticidal Activities of Propargyloxy-Diphenyl Oxide-Sulfonamide Derivatives

Agricultural pests are the primary contributing factor to crop yield reduction, particularly in underdeveloped regions. Despite the significant efficacy of pesticides in pest control, their extensive use has led to the drug-fast of insecticide resistance. Developing of new environmentally friendly plant-based pesticides is an urgent necessity. In this study, a series of diaryl ether compounds containing propargyloxy and sulfonamide groups were designed. The synthesis of these 36 compounds primarily relied on nuclear magnetic resonance for structure determination, while single-crystal X-ray diffraction was employed for certain compounds. Meanwhile, the insecticidal activities against Mythimna separata were also assessed. Some of the compounds exhibited significantly enhanced activity, the LC50 value of the highest activity compound TD8 (0.231 mg/mL) demonstrating respective increases by 100-fold compared to the plant pesticide celangulin V (23.9 mg/mL), and a 5-fold increase with the positive control L-1 (1.261 mg/mL). The interaction between the target compound and the target, as well as the consistency of the target, were verified through symptomological analysis and molecular docking. The structure-activity relationships were also conducted. This study offered a novel trajectory for the advancement and formulation of future pesticides.

Molecular Docking and Molecular Dynamics Simulations in Related to Leishmania donovani: An Update and Literature Review

Leishmaniasis, a disease caused by Leishmania parasites and transmitted via sandflies, presents in two main forms: cutaneous and visceral, the latter being more severe. With 0.7 to 1 million new cases each year, primarily in Brazil, diagnosing remains challenging due to diverse disease manifestations. Traditionally, the identification of Leishmania species is inferred from clinical and epidemiological data. Advances in disease management depend on technological progress and the improvement of parasite identification programs. Current treatments, despite the high incidence, show limited efficacy due to factors like cost, toxicity, and lengthy regimens causing poor adherence and resistance development. Diagnostic techniques have improved but a significant gap remains between scientific progress and application in endemic areas. Complete genomic sequence knowledge of Leishmania allows for the identification of therapeutic targets. With the aid of computational tools, testing, searching, and detecting affinity in molecular docking are optimized, and strategies that assess advantages among different options are developed. The review focuses on the use of molecular docking and molecular dynamics (MD) simulation for drug development. It also discusses the limitations and advancements of current treatments, emphasizing the importance of new techniques in improving disease management.

Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation

Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.

Combining Monte Carlo and Molecular Dynamics Simulations for Enhanced Binding Free Energy Estimation through Markov State Models

We present a multistep protocol, combining Monte Carlo and molecular dynamics simulations, for the estimation of absolute binding free energies, one of the most significant challenges in computer-aided drug design. The protocol is based on an initial short enhanced Monte Carlo simulation, followed by clustering of the ligand positions, which serve to identify the most relevant states of the unbinding process. From these states, extensive molecular dynamics simulations are run to estimate an equilibrium probability distribution obtained with Markov State Models, which is subsequently used to estimate the binding free energy. We tested the procedure on two different protein systems, the Plasminogen kringle domain 1 and Urokinase, each with multiple ligands, for an aggregated molecular dynamics length of 760 μs. Our results indicate that the initial sampling of the unbinding events largely facilitates the convergence of the subsequent molecular dynamics exploration. Moreover, the protocol is capable to properly rank the set of ligands examined, albeit with a significant computational cost for the, more realistic, Urokinase complexes. Overall, this work demonstrates the usefulness of combining enhanced sampling methods with regular simulation techniques as a way to obtain more reliable binding affinity estimates.

Exhaustive search of the configurational space of heat-shock protein 90 with its inhibitor by multicanonical molecular dynamics based dynamic docking

Multicanonical molecular dynamics based dynamic docking was used to exhaustively search the configurational space of an inhibitor binding to the N-terminal domain of heat-shock protein 90 (Hsp90). The obtained structures at 300 K cover a wide structural ensemble, with the top two clusters ranked by their free energy coinciding with the native binding site. The representative structure of the most stable cluster reproduced the experimental binding configuration, but an interesting conformational change in Hsp90 could be observed. The combined effects of solvation and ligand binding shift the equilibrium from a preferred loop-in conformation in the unbound state to an α-helical one in the bound state for the flexible lid region of Hsp90. Thus, our dynamic docking method is effective at predicting the native binding site while exhaustively sampling a wide configurational space, modulating the protein structure upon binding.

Selectivity Challenges in Docking Screens for GPCR Targets and Antitargets

To investigate large library docking's ability to find molecules with joint activity against on-targets and selectivity versus antitargets, the dopamine D2 and serotonin 5-HT2A receptors were targeted, seeking selectivity against the histamine H1 receptor. In a second campaign, κ-opioid receptor ligands were sought with selectivity versus the μ-opioid receptor. While hit rates ranged from 40% to 63% against the on-targets, they were just as good against the antitargets, even though the molecules were selected for their putative lack of binding to the off-targets. Affinities, too, were often as good or better for the off-targets. Even though it was occasionally possible to find selective molecules, such as a mid-nanomolar D2/5-HT2A ligand with 21-fold selectivity versus the H1 receptor, this was the exception. Whereas false-negatives are tolerable in docking screens against on-targets, they are intolerable against antitargets; addressing this problem may demand new strategies in the field.

Predicting how drug molecules bind to their protein targets

There have been substantial advances in the application of molecular modelling and simulation to drug discovery in recent years, as massive increases in computer power are coupled with continued development in the underlying methods and understanding of how to apply them. Here, we survey recent advances in one particular area-predicting how a known ligand binds to a particular protein. We focus on the four contributing classes of calculation: predicting where a binding site is on a protein; characterizing where chemical functional groups will bind to that site; molecular docking to generate a binding mode for a ligand and dynamics simulations to refine that pose and allow for protein conformation change. Examples of successful application are provided for each class.