Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Epilepsy is a chronic brain disorder characterized by recurrent seizures due to abnormal, excessive and synchronous neuronal activities in the brain. It affects approximately 65 million people worldwide, one third of which are still estimated to suffer from refractory seizures. Glutamic acid decarboxylase (GAD) that converts glutamate into GABA is a key enzyme in the dynamic regulation of neural network excitability. Importantly, clinical evidence shows that lowered GAD activity is associated with several forms of epilepsy which are often treatment resistant. In the present study, we synthetized and explored the possibility of using ethyl ketopentenoate (EKP), a lipid-permeable GAD-inhibitor, to induce refractory seizures in zebrafish larvae. Our results demonstrate that EKP evoked robust convulsive locomotor activities, excessive epileptiform discharges and upregulated c-fos expression in zebrafish. Moreover, transgenic animals in which neuronal cells express apoaequorin, a Ca2+-sensitive bioluminescent photoprotein, displayed large luminescence signals indicating strong EKP-induced neuronal activation. Molecular docking data indicated that this proconvulsant activity resulted from the direct inhibition of both gad67 and gad65. Limited protective efficacy of tested anti-seizure drugs (ASDs) demonstrated a high level of treatment resistance of EKP-induced seizures. We conclude that the EKP zebrafish model can serve as a high-throughput platform for novel ASDs discovery.

Human papillomavirus (HPV) induced cervical cancer is the second most common cause of death, after breast cancer, in females. Three prophylactic vaccines by Merck Sharp & Dohme (MSD) and GlaxoSmithKline (GSK) have been confirmed to prevent high-risk HPV strains but these vaccines have been shown to be effective only in girls who have not been exposed to HPV previously. The constitutively expressed HPV oncoproteins E6 and E7 are usually used as target antigens for HPV therapeutic vaccines. These early (E) proteins are involved, for example, in maintaining the malignant phenotype of the cells. In this study, we predicted antigenic peptides of HPV types 16 and 18, encoded by E6 and E7 genes, using an immunoinformatics approach. To further evaluate the immunogenic potential of the predicted peptides, we studied their ability to bind to class I major histocompatibility complex (MHC-I) molecules in a computational docking study that was supported by molecular dynamics (MD) simulations and estimation of the free energies of binding of the peptides at the MHC-I binding cleft. Some of the predicted peptides exhibited comparable binding free energies and/or pattern of binding to experimentally verified MHC-I-binding epitopes that we used as references in MD simulations. Such peptides with good predicted affinity may serve as candidate epitopes for the development of therapeutic HPV peptide vaccines.

The Crimean-Congo Hemorrhagic Fever virus (CCHFV) is a segmented negative single-stranded RNA virus (−ssRNA) which causes severe hemorrhagic fever in humans with a mortality rate of ~50%. To date, no vaccine has been approved. Treatment is limited to supportive care with few investigational drugs in practice. Previous studies have identified viral RNA dependent RNA Polymerase (RdRp) as a potential drug target due to its significant role in viral replication and transcription. Since no crystal structure is available yet, we report the structural elucidation of CCHFV-RdRp by in-depth homology modeling. Even with low sequence identity, the generated model suggests a similar overall structure as previously reported RdRps. More specifically, the model suggests the presence of structural/functional conserved RdRp motifs for polymerase function, the configuration of uniform spatial arrangement of core RdRp sub-domains, and predicted positively charged entry/exit tunnels, as seen in sNSV polymerases. Extensive pharmacophore modeling based on per-residue energy contribution with investigational drugs allowed the concise mapping of pharmacophoric features and identified potential hits. The combination of pharmacophoric features with interaction energy analysis revealed functionally important residues in the conserved motifs together with in silico predicted common inhibitory binding modes with highly potent reference compounds.

Tumorigenesis in humans is a multistep progression that imitates genetic changes leading to cell transformation and malignancy. Oncogenic kinases play a central role in cancer progression, rendering them putative targets for the design of anti-cancer drugs. The presented work aims to identify the potential multi-target inhibitors of oncogenic receptor tyrosine kinases (RTKs) and serine/threonine kinases (STKs). For this, chemoinformatics and structure-based virtual screening approaches were combined with an in vitro validation of lead hits on both cancerous and non-cancerous cell lines. A total of 16 different kinase structures were screened against ~739,000 prefiltered compounds using diversity selection, after which the top hits were filtered for promising pharmacokinetic properties. This led to the identification of 12 and 9 compounds against RTKs and STKs, respectively. Molecular dynamics (MD) simulations were carried out to better comprehend the stability of the predicted hit kinase-compound complexes. Two top-ranked compounds against each kinase class were tested in vitro for cytotoxicity, with compound F34 showing the most promising inhibitory activity in HeLa, HepG2, and Vero cell lines with IC50 values of 145.46 μM, 175.48 μM, and 130.52 μM, respectively. Additional docking of F34 against various RTKs was carried out to support potential multi-target inhibition. Together with reliable MD simulations, these results suggest the promising potential of identified multi-target STK and RTK scaffolds for further kinase-specific anti-cancer drug development toward combinatorial therapies.

Ebola virus disease (EVD), caused by Ebola viruses, resulted in more than 11 500 deaths according to a recent 2018 WHO report. With mortality rates up to 90%, it is nowadays one of the most deadly infectious diseases. However, no Food and Drug Administration‐approved Ebola drugs or vaccines are available yet with the mainstay of therapy being supportive care. The high fatality rate and absence of effective treatment or vaccination make Ebola virus a category‐A biothreat pathogen. Fortunately, a series of investigational countermeasures have been developed to control and prevent this global threat. This review summarizes the recent therapeutic advances and ongoing research progress from research and development to clinical trials in the development of small‐molecule antiviral drugs, small‐interference RNA molecules, phosphorodiamidate morpholino oligomers, full‐length monoclonal antibodies, and vaccines. Moreover, difficulties are highlighted in the search for effective countermeasures against EVD with additional focus on the interplay between available in silico prediction methods and their evidenced potential in antiviral drug discovery.

The enoyl‐acyl carrier protein reductase enzyme FabI is essential for fatty acid biosynthesis in Staphylococcus aureus and represents a promising target for the development of novel, urgently needed anti‐staphylococcal agents. Here, we elucidate the mode of action of the kalimantacin antibiotics, a novel class of FabI inhibitors with clinically‐relevant activity against multidrug‐resistant S. aureus. By combining X‐ray crystallography with molecular dynamics simulations, in vitro kinetic studies and chemical derivatization experiments, we characterize the interaction between the antibiotics and their target, and we demonstrate that the kalimantacins bind in a unique conformation that differs significantly from the binding mode of other known FabI inhibitors. We also investigate mechanisms of acquired resistance in S. aureus and identify key residues in FabI that stabilize the binding of the antibiotics. Our findings provide intriguing insights into the mode of action of a novel class of FabI inhibitors that will inspire future anti‐staphylococcal drug development.

With no predicted indication of adverse effects on humans, this study unravels four active multi-target inhibitors against AChE with promising affinities and good ADMET profile for the potential use in AD treatment.