Compensatory Role of Double Mutation N348I/M184V on Nevirapine Binding Landscape: Insight from Molecular Dynamics Simulation

In Silico Repurposing of J147 for Neonatal Encephalopathy Treatment: Exploring Molecular Mechanisms of Mutant Mitochondrial ATP Synthase

BACKGROUND

Neonatal Encephalopathy (NE) is a mitochondrial ATP synthase (mATPase) disease, which results in the death of infants. The case presented here is reportedly caused by complex V deficiency as a result of mutation of Arginine to Cysteine at residue 329 in the mATPase. A recent breakthrough was the discovery of J147, which targets mATPase in the treatment of Alzheimer's disease. Based on the concepts of computational target-based drug design, this study investigated the possibility of employing J147 as a viable candidate in the treatment of NE.

OBJECTIVE/METHODS

The structural dynamic implications of this drug on the mutated enzyme are yet to be elucidated. Hence, integrative molecular dynamics simulations and thermodynamic calculations were employed to investigate the activity of J147 on the mutated enzyme in comparison to its already established inhibitory activity on the wild-type enzyme.

RESULTS

A correlated structural trend occurred between the wild-type and mutant systems whereby all the systems exhibited an overall conformational transition. Equal observations in favorable free binding energies further substantiated uniformity in the mobility, and residual fluctuation of the wild-type and mutant systems. The similarity in the binding landscape suggests that J147 could as well modulate mutant mATPase activity in addition to causing structural modifications in the wild-type enzyme.

CONCLUSION

Findings suggest that J147 can stabilize the mutant protein and restore it to a similar structural state as the wild-type which depicts functionality. These details could be employed in drug design for potential drug resistance cases due to mATPase mutations that may present in the future.

Exploring the Lapse in Druggability: Sequence Analysis, Structural Dynamics and Binding Site Characterization of K-RasG12C Variant, a Feasible Oncotherapeutics Target

BACKGROUND

The difficulty in druggability of K-Ras variant has presented a challenge in the treatment of cancer diseases associated with its dysfunctionality. Despite the identification of different binding sites, limited information exists in the literature about their characteristics. Therefore, identification, crossvalidation and characterization of its druggable sites would aid the design of chemical compounds that will arrest its dysfunctionality related oncogenesis.

OBJECTIVE

This study entails the identification, cross-validation and characterization of K-Ras G12C variant's binding sites for potential druggability, coupled with the elucidation of alterations in 3D conformations and dynamics.

METHOD

Molecular dynamics simulation was carried out on the inactive, the active and the hyperactive K-RasG12Cvariant using the amber software package. The SiteMap software was employed in identifying and characterizing the druggable binding sites while the validation of the binding sites was carried out with the SiteHound and MetaPocket servers.

RESULTS

Four druggable binding sites were identified, validated and characterized based on physicochemical attributes such as size, volume, degree of enclosure or exposure, degree of contact, hydrophobic/hydrophilic character, hydrophobic/hydrophilic balance and hydrogen-bonding features. Conformational studies also revealed that the K-Ras variant exhibited notable structural instability, increased flexibility and a strongly anticorrelated movement compared to the inactive and active wildtype forms.

CONCLUSION

The attributes of the characterized druggable sites will be useful in designing site-specific K-Ras inhibitors for the treatment of K-Ras variant associated cancer diseases.

Induced Mutation Proves a Potential Target for TB Therapy: A Molecular Dynamics Study on LprG

Molecular dynamics (MD) simulations of wild-type and V91W mutant Mycobacterium tuberculosis-LprG (Mtb-LprG) were performed with the goal to provide a comprehensive understanding of the Mtb-LprG as a potential antimycobacterial target. A long-range MD simulations and post-MD analyzes led us to various results that plainly explained the impact of V91W mutation on Mtb-LprG. Herein, the results revealed that the wild-type is less stable compared to V91W mutant. This was further supported by root mean square fluctuation, where the V91W mutant showed a higher degree of flexibility compared to the wild-type. Dynamic cross-correlation analysis revealed that induced mutation leads to higher residual flexibility in the mutant structure as compared to the wild-type structure thus resulting in the existence of negatively correlated motions. The difference in principal component analysis scatter plot across the first two normal modes suggests a greater mobility of the V91W mutant conformation compared to the wild-type. Thermodynamic calculations revealed that the van der Waals (E_vdw) forces contribute the most towards binding free energy in a case of the V91W mutant as compared to the wild-type LprG complex. In addition, the residue interaction networks revealed more of E_vdw interaction existence among residues in case of the V91W mutant. This study supports the Mtb-LprG as a potential antimycobacterial target and also serves as a cornerstone to identifying new potential targets that have no inhibitors.

Characterizing the ligand-binding landscape of Zika NS3 helicase-promising lead compounds as potential inhibitors

Aim: This study aims to provide insight into the binding features of the ATPase and ssRNA sites of the NS3 helicase. Methods: Clinically approved Flavivirus inhibitors were docked to the corresponding active sites of the protein, and the three best compounds were validated with molecular dynamic simulations. Result: Binding of Ivermectin to ssRNA site and Lapachol and HMC-HO1α to the ATPase site allowed for conformational rigidity of the Zika NS3 helicase, thus stabilizing residue fluctuations and allowing for protein stability. Favorable free binding energies were also noted between compounds and the helicase, thus supporting the intermolecular forces at the helicase active site. Conclusion: The pharmacophoric characteristics found in Lapachol, HMC-HO1α and Ivermectin may be utilized in the design of a potent hybrid drug that is able to show efficient inhibition of a multitude of diseases including the detrimental co-infection of Zika virus, dengue and chikungunya.

VP40 of the Ebola Virus as a Target for EboV Therapy: Comprehensive Conformational and Inhibitor Binding Landscape from Accelerated Molecular Dynamics

The first account of the dynamic features of the loop region of VP40 of the Ebola virus was studied using accelerated molecular dynamics simulations and reported herein. Among the proteins of the Ebola virus, the matrix protein (VP40) plays a significant role in the virus lifecycle thereby making it a promising therapeutic target. Of interest is the newly elucidated N-terminal domain loop region of VP40 comprising residues K127, T129, and N130 which when mutated to alanine have demonstrated an unrecognized role for N-terminal domain-plasma membrane interaction for efficient VP40-plasma membrane localization, oligomerization, matrix assembly, and egress. The molecular understanding of the conformational features of VP40 in complex with a known inhibitor still remains elusive. Using accelerated molecular dynamics approaches, we conducted a comparative study on VP40 apo and bound systems to understand the conformational features of VP40 at the molecular level and to determine the effect of inhibitor binding with the aid of a number of post-dynamic analytical tools. Significant features were seen in the presence of an inhibitor as per molecular mechanics/generalized born surface area binding free energy calculations. Results revealed that inhibitor binding to VP40 reduces the flexibility and mobility of the protein as supported by root mean square fluctuation and root mean square deviation calculations. The study revealed a characteristic "twisting" motion and coiling of the loop region of VP40 accompanied by conformational changes in the dimer interface upon inhibitor binding. We believe that results presented in this study will ultimately provide useful insight into the binding landscape of VP40 which could assist researchers in the discovery of potent Ebola virus inhibitors for anti-Ebola therapies.

Sliding Clamp of DNA Polymerase III as a Drug Target for TB Therapy: Comprehensive Conformational and Binding Analysis from Molecular Dynamic Simulations

Tuberculosis, caused by Mycobacterium tuberculosis, is one of the most common causes of death in the world. Mycobacterium tuberculosis -sliding clamp is a protein essential for many important DNA transactions including replication and DNA repair proteins, thus, a potential drug target for tuberculosis. Further investigation is needed in understanding DNA polymerase sliding clamp structure, especially from a computational perspective. In this study, we employ a wide-range of comparative molecular dynamic analyses on two systems: Mycobacterium tuberculosis - sliding clamp enzyme in its apo and bound form. The results reported in this study shows apo conformation to be less stable, as compared to bound conformation with an average radius of gyration of 25.812 and 25.459 Å, respectively. This was further supported by root mean square fluctuation, where an apo enzyme showed a higher degree of flexibility. However, the presence of the ligand lowers radius of gyration and root mean square fluctuation and also leads to an existence of negative correlated motions. Principal component analysis further justifies the same findings, whereby the apo enzyme exhibits a higher fluctuation compared to the bound complex. In addition, a stable 3₁₀ helix located at the binding site appears to be unstable in the presence of the ligand. Hence, it is possible that the binding of the ligand may have caused a rearrangement of the structure, leading to a change in the unwinding of 3₁₀ helix. Findings reported in this study further enhance the understanding of Mycobacterium tuberculosis -DnaN and also give a lead to the development of potent tuberculosis drugs.

Structural Characterizations of the Fas Receptor and the Fas-Associated Protein with Death Domain Interactions

The Fas receptor is a representative death receptor, and the Fas-associated protein with death domain (FADD) is a crucial adapter protein needed to support the Fas receptor's activity. The Fas-FADD interactions constitute an important signaling pathway that ultimately induces apoptosis or programmed cell death in biological systems. The interactions responsible for this cell-death process are governed by the binding process of the Fas ligand to the Fas, followed by the caspase cascade activation. Using a computational approach, the present communication explores certain essential structural aspects of the Fas-FADD death domains and their interfacial interactions.

Insights on the structural perturbations in human MTHFR Ala222Val mutant by protein modeling and molecular dynamics

Methylenetetrahydrofolate reductase (MTHFR) protein catalyzes the only biochemical reaction which produces methyltetrahydrofolate, the active form of folic acid essential for several molecular functions. The Ala222Val polymorphism of human MTHFR encodes a thermolabile protein associated with increased risk of neural tube defects and cardiovascular disease. Experimental studies have shown that the mutation does not affect the kinetic properties of MTHFR, but inactivates the protein by increasing flavin adenine dinucleotide (FAD) loss. The lack of completely solved crystal structure of MTHFR is an impediment in understanding the structural perturbations caused by the Ala222Val mutation; computational modeling provides a suitable alternative. The three-dimensional structure of human MTHFR protein was obtained through homology modeling, by taking the MTHFR structures from Escherichia coli and Thermus thermophilus as templates. Subsequently, the modeled structure was docked with FAD using Glide, which revealed a very good binding affinity, authenticated by a Glide XP score of -10.3983 (kcal mol(-1)). The MTHFR was mutated by changing Alanine 222 to Valine. The wild-type MTHFR-FAD complex and the Ala222Val mutant MTHFR-FAD complex were subjected to molecular dynamics simulation over 50 ns period. The average difference in backbone root mean square deviation (RMSD) between wild and mutant variant was found to be ~.11 Å. The greater degree of fluctuations in the mutant protein translates to increased conformational stability as a result of mutation. The FAD-binding ability of the mutant MTHFR was also found to be significantly lowered as a result of decreased protein grip caused by increased conformational flexibility. The study provides insights into the Ala222Val mutation of human MTHFR that induces major conformational changes in the tertiary structure, causing a significant reduction in the FAD-binding affinity.

Flap flexibility amongst plasmepsins I, II, III, IV, and V: Sequence, structural, and molecular dynamics analyses

Herein, for the first time, we comparatively report the opening and closing of apo plasmepsin I - V. Plasmepsins belong the aspartic protease family of enzymes, and are expressed during the various stages of the P. falciparum lifecycle, the species responsible for the most lethal and virulent malaria to infect humans. Plasmepsin I, II, IV and HAP degrade hemoglobin from infected red blood cells, whereas plasmepsin V transport proteins crucial to the survival of the malaria parasite across the endoplasmic reticulum. Flap-structures covering the active site of aspartic proteases (such as HIV protease) are crucial to the conformational flexibility and dynamics of the protein, and ultimately control the binding landscape. The flap-structure in plasmepsins is made up of a flip tip in the N-terminal lying perpendicular to the active site, adjacent to the flexible loop region in the C-terminal. Using molecular dynamics, we propose three parameters to better describe the opening and closing of the flap-structure in apo plasmepsins. Namely, the distance, d1, between the flap tip and the flexible region; the dihedral angle, ϕ, to account for the twisting motion; and the TriCα angle, θ1. Simulations have shown that as the flap-structure twists, the flap and flexible region move apart opening the active site, or move toward each other closing the active site. The data from our study indicate that of all the plasmepsins investigated in the present study, Plm IV and V display the highest conformational flexibility and are more dynamic structures versus Plm I, II, and HAP.

Heat-shock protein 90 (Hsp90) as anticancer target for drug discovery: an ample computational perspective

There are over 100 different types of cancer, and each is classified based on the type of cell that is initially affected. If left untreated, cancer can result in serious health problems and eventually death. Recently, the paradigm of cancer chemotherapy has evolved to use a combination approach, which involves the use of multiple drugs each of which targets an individual protein. Inhibition of heat-shock protein 90 (Hsp90) is one of the novel key cancer targets. Because of its ability to target several signaling pathways, Hsp90 inhibition emerged as a useful strategy to treat a wide variety of cancers. Molecular modeling approaches and methodologies have become 'close counterparts' to experiments in drug design and discovery workflows. A wide range of molecular modeling approaches have been developed, each of which has different objectives and outcomes. In this review, we provide an up-to-date systematic overview on the different computational models implemented toward the design of Hsp90 inhibitors as anticancer agents. Although this is the main emphasis of this review, different topics such as background and current statistics of cancer, different anticancer targets including Hsp90, and the structure and function of Hsp90 from an experimental perspective, for example, X-ray and NMR, are also addressed in this report. To the best of our knowledge, this review is the first account, which comprehensively outlines various molecular modeling efforts directed toward identification of anticancer drugs targeting Hsp90. We believe that the information, methods, and perspectives highlighted in this report would assist researchers in the discovery of potential anticancer agents.

Flap dynamics of plasmepsin proteases: insight into proposed parameters and molecular dynamics

Herein, for the first time, we report the flap opening and closing in Plasmepsin proteases – plasmepsin II (PlmII) was used as a prototype model.