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.

Recently emerged SARS-CoV-2 caused a major outbreak of coronavirus disease 2019 (COVID-19) and instigated a widespread fear, threatening global health safety. To date, no licensed antiviral drugs or vaccines are available against COVID-19 although several clinical trials are under way to test possible therapies. During this urgent situation, computational drug discovery methods provide an alternative to tiresome high-throughput screening, particularly in the hit-to-lead-optimization stage. Identification of small molecules that specifically target viral replication apparatus has indicated the highest potential towards antiviral drug discovery. In this work, we present potential compounds that specifically target SARS-CoV-2 vital proteins, including the main protease, Nsp12 RNA polymerase and Nsp13 helicase. An integrative virtual screening and molecular dynamics simulations approach has facilitated the identification of potential binding modes and favourable molecular interaction profile of corresponding compounds. Moreover, the identification of structurally important binding site residues in conserved motifs located inside the active site highlights relative importance of ligand binding based on residual energy decomposition analysis. Although the current study lacks experimental validation, the structural information obtained from this computational study has paved way for the design of targeted inhibitors to combat COVID-19 outbreak.

Study of monogenic forms of obesity has demonstrated the pivotal role of the central leptin-melanocortin pathway in controlling energy balance, appetite and body weight 1 . The majority of loss-of-function mutations (mostly recessive or co-dominant) have been identified in genes that are directly involved in leptin-melanocortin signaling. These genes, however, only explain obesity in <5% of cases, predominantly from outbred populations 2 . We previously showed that, in a consanguineous population in Pakistan, recessive mutations in known obesity-related genes explain ~30% of cases with severe obesity3-5. These data suggested that new monogenic forms of obesity could also be identified in this population. Here we identify and functionally characterize homozygous mutations in the ADCY3 gene encoding adenylate cyclase 3 in children with severe obesity from consanguineous Pakistani families, as well as compound heterozygous mutations in a severely obese child of European-American descent. These findings highlight ADCY3 as an important mediator of energy homeostasis and an attractive pharmacological target in the treatment of obesity.

The recent outbreak of Zika virus (ZIKV) infection in Brazil has developed to a global health concern due to its likely association with birth defects (primary microcephaly) and neurological complications. Consequently, there is an urgent need to develop a vaccine to prevent or a medicine to treat the infection. In this study, immunoinformatics approach was employed to predict antigenic epitopes of Zika viral proteins to aid in development of a peptide vaccine against ZIKV. Both linear and conformational B-cell epitopes as well as cytotoxic T-lymphocyte (CTL) epitopes were predicted for ZIKV Envelope (E), NS3 and NS5 proteins. We further investigated the binding interactions of altogether 15 antigenic CTL epitopes with three class I major histocompatibility complex (MHC I) proteins after docking the peptides to the binding groove of the MHC I proteins. The stability of the resulting peptide-MHC I complexes was further studied by molecular dynamics simulations. The simulation results highlight the limits of rigid-body docking methods. Some of the antigenic epitopes predicted and analyzed in this work might present a preliminary set of peptides for future vaccine development against ZIKV.

The SARS-CoV-2 was confirmed to cause the global pandemic of coronavirus disease 2019 (COVID-19). The 3-chymotrypsin-like protease (3CLpro), an essential enzyme for viral replication, is a valid target to combat SARS-CoV and MERS-CoV. In this work, we present a structure-based study to identify potential covalent inhibitors containing a variety of chemical warheads. The targeted Asinex Focused Covalent (AFCL) library was screened based on different reaction types and potential covalent inhibitors were identified. In addition, we screened FDA-approved protease inhibitors to find candidates to be repurposed against SARS-CoV-2 3CLpro. A number of compounds with significant covalent docking scores were identified. These compounds were able to establish a covalent bond (C–S) with the reactive thiol group of Cys145 and to form favorable interactions with residues lining the substrate-binding site. Moreover, paritaprevir and simeprevir from FDA-approved protease inhibitors were identified as potential inhibitors of SARS-CoV-2 3CLpro. The mechanism and dynamic stability of binding between the identified compounds and SARS-CoV-2 3CLpro were characterized by molecular dynamics (MD) simulations. The identified compounds are potential inhibitors worthy of further development as COVID-19 drugs. Importantly, the identified FDA-approved anti-hepatitis-C virus (HCV) drugs paritaprevir and simeprevir could be ready for clinical trials to treat infected patients and help curb COVID-19. Communicated by Ramaswamy H. Sarma