Single Active Site Mutation Causes Serious Resistance of HIV Reverse Transcriptase to Lamivudine: Insight from Multiple Molecular Dynamics Simulations

'Polymorphism-aided' Selective Targeting and Inhibition of Caspase-6 by a Novel Allosteric Inhibitor Towards Efficient Alzheimer's Disease Treatment

The predominance of Alzheimer's disease (AD) among the aged remains a global challenge. As such, the search for alternative and effective therapeutic options continuous unabated. Among the therapeutic targets explored over the years toward impeding the progression of AD is caspase-6 (Casp6), although selectively targeting Casp6 remains a challenge due to high homology with other members of the caspase family. Methyl 3-[(2,3-dihydro-1-benzofuran-2-yl formamido) methyl]-5-(furan-2-amido) benzoate (C13), a novel allosteric inhibitor, is reportedly shown to exhibit selective inhibition against mutant human Casp6 variants (E35K). However, structural and atomistic insights accounting for the reported inhibitory prowess of C13 remains unresolved. In this study, we seek to unravel the mechanistic selectivity of C13 coupled with the complementary effects of E35K single-nucleotide polymorphism (SNP) relative to Casp6 inhibition. Analyses of binding dynamics revealed that the variant Lysine-35 mediated consistent high-affinity interactions with C13 at the allosteric site, possibly forming the molecular basis of the selectivity of C13 as well as its high binding free energy as estimated. Analysis of residue interaction network around Glu35 and Lys35 revealed prominent residue network distortions in the mutant Casp6 conformation evidenced by a decrease in node degree, reduced number of edges and an increase short in path length relative to a more compact conformation in the wild system. The relatively higher binding free energy of C13 coupled with the stronger intermolecular interactions elicited in the mutant conformation further suggests that the mutation E35K probably favours the inhibitory activity of C13. Further analysis of atomistic changes showed increased C-α atom deviations consistent with structural disorientations in the mutant Casp6. Structural Insights provided could open up a novel paradigm of structure-based design of selective allosteric inhibition of Casp6 towards the treatment of neurodegenerative 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.

Across the blood-brain barrier: Neurotherapeutic screening and characterization of naringenin as a novel CRMP-2 inhibitor in the treatment of Alzheimer's disease using bioinformatics and computational tools

The discovery and developmental processes of CNS drugs have been limited by the inability of potential drug molecules to pass through the blood-brain barrier (BBB). This presents a significant setback in the treatment of neurodegenerative disorders such as Alzheimer's disease (AD), hence the need for compounds that can adhere strictly to the selective criteria of suitable CNS drugs. Collapsin response mediator protein-2 (CRMP-2) has been recently identified as a viable target in neurotherapeutics due to its involvement in the etiology of AD. As shown in previous studies, Naringenin (NAR), a small molecule derivative of Drynaria rhizome (DR) extract, specifically binds CRMP-2 and reduces its phosphorylation. This was shown to facilitate axonal regrowth, with improvement in cognition and learning. Herein, we report the first account of the use of cheminformatics techniques to define the CNS drug-suitability of NAR using selective criteria, coupled with the prediction of possible biological activities and toxicities. Also, we evaluated the mechanistic activity of NAR by modeling its molecular interaction with human CRMP-2 (hCRMP-2). Physicochemical analyses revealed the suitability of NAR as a CNS drug and its ability to transverse the BBB. Possible neurogenic, anti-carcinogenic and cardioprotective activities were also predicted. NAR exhibited favorable binding to CRMP-2 and formed strong bonds with active site residues, which accounts for its stabilization and affinity. Moreover, NAR induced notable conformational changes in CRMP-2, an occurrence that could possibly disrupt kinase-mediated phosphorylation. These findings will aid in the optimization of NAR and improve its neurotherapeutic activities in the treatment of AD.

Non-active site mutations disturb the loop dynamics, dimerization, viral budding and egress of VP40 of the Ebola virus

The first account of the dynamic features of the loop region of VP40 of the Ebola virus (EboV) using accelerated molecular dynamics (aMD) simulations is reported herein. Due to its major role in the Ebola life cycle, VP40 is considered a promising therapeutic target. The available experimental data on the N-terminal domain (NTD) loop indicates that mutations K127A, T129A and N130A demonstrate an unrecognized role for NTD-plasma membrane (PM) interaction for efficient VP40-PM localization, oligomerization, matrix assembly and egress. Despite experimental results, the molecular description of VP40 and the information it can provide still remain vague. Therefore, to gain further molecular insight into the effect of mutations on the loop region of VP40 and its effects on the overall protein conformation and VP40 dimerization, aMD simulations and post-dynamic analyses were employed for wildtype (WT) and mutant systems. The results showed significant variations in the presence of mutations as per RMSF, RMSD, R_g, PCA and distance calculations in comparison to the WT. These results could provide researchers with insight with regards to the conformational aspects concerning VP40 and its close relation to the experimental data. We believe that the results presented in this study will ultimately provide a useful understanding of the structural landscape of the loop region of VP40, which would contribute towards the discovery of novel EboV inhibitors.

Variability of the Cyclin-Dependent Kinase 2 Flexibility Without Significant Change in the Initial Conformation of the Protein or Its Environment; a Computational Study

BACKGROUND

Protein flexibility, which has been referred as a dynamic behavior has various roles in proteins' functions. Furthermore, for some developed tools in bioinformatics, such as protein-protein docking software, considering the protein flexibility, causes a higher degree of accuracy. Through undertaking the present work, we have accomplished the quantification plus analysis of the variations in the human Cyclin Dependent Kinase 2 (hCDK2) protein flexibility without affecting a significant change in its initial environment or the protein per se.

OBJECTIVES

The main goal of the present research was to calculate variations in the flexibility for each residue of the hCDK2, analysis of their flexibility variations through clustering, and to investigate the functional aspects of the residues with high flexibility variations.

MATERIALS AND METHODS

Using Gromacs package (version 4.5.4), three independent molecular dynamics (MD) simulations of the hCDK2 protein (PDB ID: 1HCL) was accomplished with no significant changes in their initial environments, structures, or conformations, followed by Root Mean Square Fluctuations (RMSF) calculation of these MD trajectories. The amount of variations in these three curves of RMSF was calculated using two formulas.

RESULTS

More than 50% of the variation in the flexibility (the distance between the maximum and the minimum amount of the RMSF) was found at the region of Val-154. As well, there are other major flexibility fluctuations in other residues. These residues were mostly positioned in the vicinity of the functional residues. The subsequent works were done, as followed by clustering all hCDK2 residues into four groups considering the amount of their variability with respect to flexibility and their position in the RMSF curves.

CONCLUSIONS

This work has introduced a new class of flexibility aspect of the proteins' residues. It could also help designing and engineering proteins, with introducing a new dynamic aspect of hCDK2, and accordingly, for the other similar globular proteins. In addition, it could provide a better computational calculation of the protein flexibility, which is, especially important in the comparative studies of the proteins' flexibility.

Molecular Dynamics Simulations of Ligand-Induced Flap Conformational Changes in Cathepsin-D-A Comparative Study

The flap region in aspartic proteases is a unique structural feature to this class of enzymes, and found to have a profound impact on protein overall structure, function, and dynamics. Understanding the structure and dynamic behavior of the flap regions is crucial in the design of selective inhibitors against aspartic proteases. Cathepsin-D, an aspartic protease enzyme, has been implicated in a long list of degenerative diseases as well as breast cancer progression. Presented herein, for the first time, is a comprehensive description of the conformational flap dynamics of cathepsin-D using a comparative 50 ns "multiple" molecular dynamics simulations. Diverse collective metrics were proposed to accurately define flap dynamics. These are distance d1 between the flap tips residues (Gly79 and Met301); dihedral angle ϕ; in addition to TriCα angles Gly79-Asp33-Asp223, θ1 , and Gly79-Asp223-Met301, θ2 . The maximum distance attained throughout the simulation was 17.42 and 11.47 Å for apo and bound cathepsin-D, respectively, while the minimum distance observed was 8.75 and 6.32 Å for apo and bound cathepsin-D, respectively. The movement of the flap as well as the twist of the active pocket can properly be explained by measuring the angle, θ1 , between Gly79-Asp33-Met301 and correlating it with the distance Cα of the flap tip residues. The asymmetrical opening of the binding cavity was best described by the large shift of -6.26° to +20.94° in the dihedral angle, ϕ, corresponding to the full opening of the flap at a range of 31-33 ns. A wide-range of post-dynamic analyses was also applied in this report to supplement our findings. We believe that this report would augment current efforts in designing potent structure-based inhibitors against cathepsin-D in the treatment of breast cancer and other degenerative diseases. J. Cell. Biochem. 117: 2643-2657, 2016. © 2016 Wiley Periodicals, Inc.