Modeling the interaction between anti-cancer drug penicillamine and pristine and functionalized carbon nanotubes for medical applications: density functional theory investigation and a molecular dynamics simulation

Abstract

The present study focuses on the prediction and investigation of binding properties of penicillamine with pure (5,5) single-walled carbon nanotube (SWCNT) and functionalized SWCNT (f-SWCNT) through the B3LYP and M06-2X functionals using the 6-31G** basis set. The electronic and structural properties, adsorption energy and frontier molecular orbitals of various configurations are examined. Our theoretical results indicated that the interaction of the nanotubes with penicillamine molecule is weak so that the drug adsorption process is typically physisorption. Also, results of theoretical calculations show that the adsorption of the drug molecule on f-SWCNT is stronger with shorter intermolecular distances in comparison to pure SWCNT. The natural bond orbital (NBO) analysis of studied systems demonstrates that the charge is transferred from penicillamine molecule to the nanotubes. Furthermore, molecular dynamics (MD) simulation is employed to evaluate the dynamic and diffusion behavior of drug molecule on SWCNT and f-SWCNT. Energy results show that drug molecule spontaneously moves toward the carriers, and the van der Waals energy contributions in drug adsorption are more than electrostatic interaction. The obtained results from MD simulation confirm that the functionalization of SWCNT leads to increase in the solubility of the carrier in aqueous solution.

Communicated by Ramaswamy H. Sarma     

In silico studies on the interaction between bioactive ligands and ALK5, a biological target related to the cancer treatment

Studies have showed that there are many biological targets related to the cancer treatment, for example, TGF type I receptor (TGF-βRI or ALK5). The ALK5 inhibition is a strategy to treat some types of cancer, such as breast, lung, and pancreas. Here, we performed CoMFA and CoMSIA studies for 70 ligands with ALK5 inhibition. The internal validation for both models (q²_LOO = 0.887 and 0.822, respectively) showed their robustness, while the external validations showed their predictive power (CoMFA: r²_test = 0.998; CoMSIA: r²_test = 0.975). After all validations, CoMFA and CoMSIA maps indicated physicochemical evidences on the main factors involved in the interaction between bioactive ligands and ALK5. Therefore, these results suggest molecular modifications to design new ALK5 inhibitors.     

Temperature-induced structural changes of apo-lactoferrin and their functional implications: a molecular dynamics simulation study

Lactoferrin (Lf) is an iron-binding glycoprotein present in secretory fluids from human and bovine sources. Sequence alignment was employed to identify a region on the C-lobe of Lf capable of binding to bacterial cell surfaces, followed by all-atom explicit solvent molecular dynamics simulations to study the conformational changes of Lf after exposure to three processing temperatures: pasteurisation (72°C), spray drying (90°C) and ultra-high temperature (UHT) (127°C). Below 90°C, the simulations indicate relatively minor changes in overall protein structure. At UHT conditions (127°C), however, marked disruptions to protein structure were found as demonstrated by a substantial decrease in protein dimensions due to collapse in the inter-lobe region. There was also a marked increase in residue fluctuations in several regions of known functional importance, including antibacterial, iron-binding, and putative membrane binding regions, the latter of which is stabilised by a triplet of hydrophobic residues comprised of Leu446, Trp448 and Leu451 at low temperature, but which are disrupted under UHT conditions. A unique network analysis confirmed these results as demonstrated by large clusters of residues with increased dynamical correlation in the N-terminal lobe.     

Toward the understanding of the structure and dynamics of protein-carbohydrate interactions: molecular dynamics studies of the complexes between hevein and oligosaccharidic ligands

Herein we study, through all atom molecular dynamics simulations, the complex between hevein and two N-acetylated chitin oligomers, namely N,N(')-diacetylchitobiose and N,N('),N(")-triacetylchitotriose. The results of the simulations for two disaccharide complexes and one trisaccharide complex show that a carbohydrate oligomer is able to move on the surface of the relatively flat binding pocket of hevein, therefore occupying different binding subpockets. Statistical analysis methods were also applied in order to define the principal overall motions in the complexes, showing how the different ligands in the simulations modulate the protein motions. The oligosaccharide binding can be considered as defined by a subtle balance between enthalpic (formation of intermolecular interactions between the ligand and the receptor) and entropic (due mainly to the possibility for the sugar to move on the surface of the protein domain) effects, determining multiple binding conformations. This structural and dynamical view could parallel the results obtained by regularly used restrained MD simulations based on NOE NMR data that provide a well defined structure for both the disaccharide and trisaccharide complexes, and agrees with the observations for longer oligosaccharide chains.     

Molecular modeling study for the design of novel acetyl-CoA carboxylase inhibitors using 3D QSAR,molecular docking and dynamic simulations

Acetyl-CoA carboxylase (ACC) enzyme plays an important role in the regulation of biosynthesis and oxidation of fatty acids. ACC is a recognized drug target for the treatment of obesity and diabetes. Combination of ligand and structure-based in silico methods along with activity and toxicity prediction provides best lead compounds in the drug discovery process. In this study, a data-set of 100 ACC inhibitors were used for the development of comparative molecular field analysis (CoMFA) and comparative molecular similarity index matrix analysis (CoMSIA) models. The generated contour maps were used for the design of novel ACC inhibitors. CoMFA and CoMSIA models were used for the predication of activity of designed compounds. In silico toxicity risk prediction study was carried out for the designed compounds. Molecular docking and dynamic simulations studies were performed to know the binding mode of designed compounds with the ACC enzyme. The designed compounds showed interactions with key amino acid residues important for catalysis, and good correlation was observed between binding free energy and inhibition of ACC.     

3D-QSAR-aided design of potent c-Met inhibitors using molecular dynamics simulation and binding free energy calculation

Abstract

Mesenchymal-epithelial transition factor (c-Met) is a member of receptor tyrosine kinase. It involves in various cellular signaling pathways which includes proliferation, motility, migration, and invasion. Over-expression of c-Met has been reported in various cancers. Hence, it is an ideal therapeutic target for cancer. The main objective of the study is to identify crucial residues involved in the inhibition of c-Met kinase and to design a series of potent imidazo [4,5-b] pyrazine derivatives as c-Met inhibitors. Docking was used to identify important active site residues involved in the inhibition of c-Met kinase which was further validated by 100 ns of molecular dynamics simulation and free energy calculation using molecular mechanics generalized born surface area. Furthermore, binding energy decomposition identified that residues Tyr1230, Met1211, Asp1222, Tyr1159, Met1160, Val1092, Ala1108, and Leu1157 contributed favorably to the binding stability of compound 32. Receptor-guided Comparative Molecular Field Analysis (CoMFA) (q² = 0.751, NOC = 6, r² = 0.933) and Comparative Molecular Similarity Indices Analysis (COMSIA) (q² = 0.744, NOC = 6, r² = 0.950) models were generated based on the docked conformation of the most active compound 32. The robustness of these models was tested using various validation techniques and found to be predictive. The results of CoMFA and CoMSIA contour maps exposed the regions favorable to enhance the activity. Based on this information, 27 novel c-Met inhibitors were designed. These designed compounds exhibited potent activity than the most active compound of the existing dataset.

Communicated by Ramaswamy H. Sarma     

pH-induced structural changes of apo-lactoferrin and implications for its function: a molecular dynamics simulation study

The minor protein in milk, lactoferrin (Lf), is known for a variety of biological functions, and has been investigated as a protective encapsulant for probiotic bacteria in health-promoting food products. Lf is likely to be exposed to extreme pH conditions which are known to have disruptive influences on its functionality. The molecular mechanisms underlying these pH-dependent changes are not well-understood. To explore the potential of Lf as an encapsulant, molecular dynamics (MD) simulations were applied to study its conformational changes under extreme acidic (pH 1.0) or basic (pH 14.0) conditions, relative to neutral pH. Simulations indicate that the structure of apo-Lf is relatively stable at neutral pH, while acidic and basic pH result in substantially greater flexibility, partly induced by the loss of contacts between the N- and C-terminal lobes, causing them to undergo extensive relative bending and twisting motions. Basic pH causes greater structural disruption compared to acidic exposure. The latter has greater influence on the N-terminus, with increased fluctuations and disruptions of inter-residue contacts compared to those at neutral pH; while basic pH was found to more prominently disrupt contacts at the C-terminus. These results help elucidate possible functional consequences on Lf of exposure to extreme pH conditions.     

Understanding the molecular interaction of human argonaute-2 and miR-20a complex: A molecular dynamics approach

Argonaute-2 (AGO2), a member of the Argonaute family, is the only member possessing catalytic and RNA silencing activity. In here, a molecular dynamics (MDs) simulation was performed using the crystal structure of human AGO2 protein complex with miR-20a. miR-20a is involved with various kind of biological process like heart and lung development, oncogenic process, etc. In precise, MD simulation was carried out with AGO2 protein complex with wild type, two mutant sites and four mutant sites in guided microRNA (miRNA). It has been noted that root-mean-square deviation (RMSD) of atomic positions of nucleic acid for wild type and two mutant sites guided miRNA has the same pattern of fluctuations, which stabilizes around 0.27 nm after 2 ns. Cα atom of AGO2 protein in the complex shows that this complex with wild type and two mutant site mutation duplex has a stable RMSD value after 20 ns, ranging between 0.14 and 0.21 nm. From the root-mean-square fluctuation (RMSF), we observed an increased pattern of fluctuations for the atoms of four mutant complex of AGO2-miR-20a complex. This increased RMSF of non-mutated nucleic acids is contributed by U-A bond breaking at the site of the nucleotide of U2 of guided miRNA, as observed from the duplex structure taken at different time steps of the simulation. Superimposed structure of the miRNA-mRNA duplex for the three complexes depicts that the three miRNA-mRNA duplexes are stable during the simulation. Current work demonstrates the possible correlations between the conformational changes of this AGO2-miR-20a duplex structure and the interactions of different atoms.     

A study of the structural properties and thermal stability of human Bcl-2 by molecular dynamics simulations

The anti-apoptotic B-cell lymphoma 2 (Bcl-2) protein interacts with several proteins that regulate the apoptotic properties of cells. In this research, we conduct several all-atom molecular dynamics (MD) simulations under high-temperature unfolding conditions, from 400 to 800?K, for 25?ns. These simulations were performed using a model of an engineered Bcl-2 human protein (Bcl-2-Δ22Σ3), which lacks 22 C-terminal residues of the transmembrane domain. The aim of this study is to gain insight into the structural behavior of Bcl-2-Δ22Σ3 by mapping the conformational movements involved in Bcl-2 stability and its biological function. To build a Bcl-2-Δ22Σ3 three-dimensional model, the protein core was built by homology modeling and the flexible loop domain (FLD, residues 33-91) by ab initio methods. Further, the entire protein model was refined by MD simulations. Afterwards, the production MD simulations showed that the FLD at 400 and 500?K has several conformations reaching into the protein core, whereas at 600?K some of the alpha-helices were lost. At 800?K, the Bcl-2 core is destabilized suggesting a possible mechanism for protein unfolding, where the alpha helices 1 and 6 were the most stable, and a reduction in the number of hydrogen bonds initially occurs. In conclusion, the structural changes and the internal protein interactions suggest that the core and the FLD are crucial components of Bcl-2 in its function of regulate ng access to the recognition sites of kinases and caspases.     

2D-QSAR study,molecular docking,and molecular dynamics simulation studies of interaction mechanism between inhibitors and transforming growth factor-beta receptor I (ALK5)

Transforming growth factor type 1 receptor (ALK5) is kinase associated with a wide variety of pathological processes, and inhibition of ALK5 is a good strategy to treat many kinds of cancer and fibrotic diseases. Recently, a series of compounds have been synthesized as ALK5 inhibitors. However, the study of their selectivity against other potential targets remains elusive. In this research, a data-set of ALK5 inhibitors were collected and studied based on the combination of 2D-QSAR, molecular docking and molecular dynamics simulation. The quality of QSAR models were assessed statistically by F, R², and R²_ADJ, proved to be credible. The cross-validations for the models (q²_LOO = 0.571 and 0.629, respectively) showed their robustness, while the external validations (r²_test = 0.703 and 0.764, respectively) showed their predictive power. Besides, the predicted binding free energy results calculated by MM/GBSA method were in accordance with the experimental data, and the van der Waals energy term was the factor that had the most significant impact on ligand binding. What is more, several important residues were found to significantly affect the binding affinity. Finally, based on our analyses above, a proposed series of molecules were designed.     

Signaling pathways of PDZ2 domain: a molecular dynamics interaction correlation analysis

PDZ domains are found in many signaling proteins. One of their functions is to provide scaffolds for forming membrane-associated protein complexes by binding to the carboxyl termini of their partners. PDZ domains are thought also to play a signal transduction role by propagating the information that binding has occurred to remote sites. In this study, a molecular dynamics (MD) simulation-based approach, referred to as an interaction correlation analysis, is applied to the PDZ2 domain to identify the possible signal transduction pathways. A residue correlation matrix is constructed from the interaction energy correlations between all residue pairs obtained from the MD simulations. Two continuous interaction pathways, starting at the ligand binding pocket, are identified by a hierarchical clustering analysis of the residue correlation matrix. One pathway is mainly localized at the N-terminal side of helix alpha1 and the adjacent C-terminus of loop beta1-beta2. The other pathway is perpendicular to the central beta-sheet and extends toward the side of PDZ2 domain opposite to the ligand binding pocket. The results complement previous studies based on multiple sequence analysis, NMR, and MD simulations. Importantly, they reveal the energetic origin of the long-range coupling. The PDZ2 results, as well as the earlier rhodopsin analysis, show that the interaction correlation analysis is a robust approach for determining pathways of intramolecular signal transduction.     

Molecular dynamics investigations of structural and functional changes in Bcl-2 induced by the novel antagonist BDA-366

Apoptosis is a fundamental biological phenomenon, in which anti- or proapoptotic proteins of the Bcl-2 family regulate a committed step. Overexpression of Bcl-2, the prototypical antiapoptotic protein in this family, is associated with therapy resistance in various human cancers. Accordingly, Bcl-2 inhibitors intended for cancer therapy have been developed, typically against the BH3 domain. Recent experimental evidences have shown that the antiapoptotic function of Bcl-2 is not immutable, and that BDA-366, a novel antagonist of the BH4 domain, converts Bcl-2 from a survival molecule to an inducer of cell death. In this study, the underlying mechanisms of this functional conversion were investigated by accelerated molecular dynamics simulation. Results revealed that Pro127 and Trp30 in the BH4 domain rotate to stabilize BDA-366 via π-π interactions, and trigger a series of significant conformational changes of the α3 helix. This rearrangement blocks the hydrophobic binding site (HBS) in the BH3 domain and further prevents binding of BH3-only proteins, which consequently allows the BH3-only proteins to activate the proapoptotic proteins. Analysis of binding free energy confirmed that BDA-366 cross-inhibits BH3-only proteins, implying negative cooperative effects across separate binding sites. The newly identified blocked conformation of the HBS along with the open to closed transition pathway revealed by this study advances the understanding of the Bcl-2 transition from antiapoptotic to proapoptotic function, and yielded new structural insights for novel drug design against the BH4 domain.

Communicated by Ramaswamy H. Sarma     

Structure prediction and molecular dynamics simulations of a G-protein coupled receptor: human CCR2 receptor

CC chemokine receptor type-2 (CCR2) is a member of G-protein coupled receptors superfamily, expressed on the cell surface of monocytes and macrophages. It binds to the monocyte chemoattractant protein-1, a CC chemokine, produced at the sites of inflammation and infection. A homology model of human CCR2 receptor based on the recently available C-X-C chemokine recepor-4 crystal structure has been reported. Ligand information was used as an essential element in the homology modeling process. Six known CCR2 antagonists were docked into the model using simple and induced fit docking procedure. Docked complexes were then subjected to visual inspection to check their suitability to explain the experimental data obtained from site directed mutagenesis and structure-activity relationship studies. The homology model was refined, validated, and assessed for its performance in docking-based virtual screening on a set of CCR2 antagonists and decoys. The docked complexes of CCR2 with the known antagonists, TAK779, a dual CCR2/CCR5 antagonist, and Teijin-comp1, a CCR2 specific antagonist were subjected to molecular dynamics (MD) simulations, which further validated the binding modes of these antagonists. B-factor analysis of 20?ns MD simulations demonstrated that Cys190 is helpful in providing structural rigidity to the extracellular loop (EL2). Residues important for CCR2 antagonism were recognized using free energy decomposition studies. The acidic residue Glu291 from TM7, a conserved residue in chemokine receptors, is favorable for the binding of Teijin-comp1 with CCR2 by ΔG of ?11.4?kcal/mol. Its contribution arises more from the side chains than the backbone atoms. In addition, Tyr193 from EL2 contributes ?0.9?kcal/mol towards the binding of the CCR2 specific antagonist with the receptor. Here, the homology modeling and subsequent molecular modeling studies proved successful in probing the structure of human CCR2 chemokine receptor for the structure-based virtual screening and predicting the binding modes of CCR2 antagonists.     

Loop-driven conformational transition between the alternative and collapsed form of prethrombin-2: targeted molecular dynamics study

Two distinct crystal structures of prethrombin-2, the alternative and collapsed forms, are elucidated by X-ray crystallogrphy. We analyzed the conformational transition from the alternative to the collapsed form employing targeted molecular dynamics (TMD) simulation. Despite small RMSD difference in the two X-ray crystal structures, some hydrophobic residues (W60d, W148, W215, and F227) show a significant difference between the two conformations. TMD simulation shows that the four hydrophobic residues undergo concerted movement from dimer to trimer transition via tetramer state in the conformational change from the alternative to the collapsed form. We reveal that the concerted movement of the four hydrophobic residues is controlled by movement of specific loop regions behind. In this paper, we propose a sequential scenario for the conformational transition from the alternative form to the collapsed form, which is partially supported by the mutant W148A simulation.     

The investigation of the G-quadruplex aptamer selectivity to Pb2+ ion: a joint molecular dynamics simulation and density functional theory study

Abstract

The aptamers with the ability to form a G-quadruplex structure can be stable in the presence of some ions. Hence, study of the interactions between such aptamers and ions can be beneficial to determine the highest selective aptamer toward an ion. In this article, molecular dynamics (MD) simulations and quantum mechanics (QM) calculations have been applied to investigate the selectivity of the T30695 aptamer toward Pb²⁺ in comparison with some ions. The Free Energy Landscape (FEL) analysis indicates that Pb²⁺ has remained inside the aptamer during the MD simulation, while the other ions have left it. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding energies prove that the conformational stability of the aptamer is the highest in the presence of Pb²⁺. According to the compaction parameters, the greatest compressed ion-aptamer complex, and hence, the highest ion-aptamer interaction have been induced in the presence of Pb²⁺. The contact maps clarify the closer contacts between the nucleotides of the aptamer in the presence of Pb²⁺. The density functional theory (DFT) results show that Pb²⁺ forms the most stable complex with the aptamer, which is consistent with the MD results. The QM calculations reveal that the N-H bonds and the O…H distances are the longest and the shortest, respectively, in the presence of Pb²⁺. The obtained results verify that the strongest hydrogen bonds (HBs), and hence, the most compressed aptamer structure are induced by Pb²⁺. Besides, atoms in molecules (AIM) and natural bond orbital (NBO) analyses confirm the results.

Communicated by Ramaswamy H. Sarma     

Molecular docking, 3D-QSAR,molecular dynamics,synthesis and anticancer activity of tyrosine kinase 2 (TYK 2) inhibitors

Abstract

A therapeutic rationale is proposed by selectively targeting tyrosine kinase 2 (TYK 2) to obtain potent TYK 2 inhibitors by molecular modeling studies. In the present study, we have taken tyrosine kinase (TYK 2) inhibitors and carried out molecular docking, 3?D quantitative structure–activity relationship (3D-QSAR) analysis and molecular dynamics (MD). Based on the 3D-QSAR results thirteen new compounds (R-1 to R-13) were designed and synthesized in good yields. The synthesized molecules were evaluated for their in vitro anticancer activity against LnCap and A549 cell lines. The molecules R-1, R-3, R-5, R-7, and R-10 exhibited considerable anti cancer activity.     

Docking and molecular dynamics studies of new potential inhibitors of the human epidermal receptor 2

Compounds similar to lapatinib and gefitinib have been investigated as potential inhibitors of the intracellular receptor tyrosine kinase (RTK) domain of the human epidermal receptor 2 (HER2), which is a promising molecular target to the drug design of new chemotherapies for breast, lung, ovarian and colorectal cancers. In this study, we have searched potential HER2 inhibitors used for treatment of other illnesses such as hepatitis, bacterial infections and sexual impotence screened in the DrugBank. The compounds selected were subjected to virtual screening docking in order to evaluate the main interactions between them and the RTK domain of HER2. The selected compounds were investigated by flexible docking, molecular dynamics studies and ΔG _bind calculations. The results suggest that antrafenine, saprisartan, reserpine, irinotecan and udenafil are potential candidates to inhibit the RTK domain of HER2.     

Phosphorylation of KLHL3 at serine 433 impairs its interaction with the acidic motif of WNK4: a molecular dynamics study

Interaction between the acidic motif (AM) of protein kinase WNK4 and the Kelch domain of KLHL3 are involved in the pathogenesis of pseudohypoaldosteronism type II, a hereditary form of hypertension. This interaction is disrupted by some disease‐causing mutations in either WNK4 or KLHL3, or by angiotensin II‐ and insulin‐induced phosphorylation of KLHL3 at serine 433, which is also a site frequently mutated in patients. However, the mechanism by which this phosphorylation disrupts the interaction is unclear. In this study, we approached this problem using molecular dynamics simulation with structural, dynamical and energetic analyses. Results from independent simulations indicate that when S433 was phosphorylated, the electrostatic potential became more negative in the AM binding site of KLHL3 and therefore was unfavorable for binding with the negatively charged AM. In addition, the intermolecular hydrogen bond network that kept the AM stable in the binding site of KLHL3 was disrupted, and the forces for the hydrophobic interactions between the AM of WNK4 and KLHL3 were also reduced. As a result, the weakened interactions were no longer capable of holding the AM of WNK4 at its binding site in KLHL3. In conclusion, phosphorylation of KLHL3 at S433 disrupts the hydrogen bonds, hydrophobic and electrostatic interactions between the Kelch domain of KLHL3 and the AM of WNK4. This study provides a key molecular understanding of the KLHL3‐mediated regulation of WNK4, which is an integrative regulator of electrolyte homeostasis and blood pressure regulation in the kidney. Significances Statement: WNK4 is an integrative regulator of electrolyte homeostasis, which is important in the blood pressure regulation by the kidney. Interaction between WNK4 and KLHL3 is a key physiological process that is impaired in a hereditary form of hypertension. This study provides substantial new insights into the role of phosphorylation of KLHL3 in regulating the interaction with WNK4, and therefore advances our understanding of molecular pathogenesis of hypertension and the mechanism of blood pressure regulation.     

A molecular dynamics study of liquid aliphatic alcohols: simulation of density and self-diffusion coefficient using a modified OPLS force field

A set of 13 aliphatic alcohols was modelled by molecular dynamics simulations at temperatures from 288 to 338 K using the optimised potential for liquid simulations (OPLS) united-atom force field, the OPLS all-atom force field and the OPLS all-atom force field with modified partial charges of the hydroxyl group. The set includes primary and secondary alcohols, and mono-, di- and trialcohols, and covers a broad range of polarities from log P = ? 0.74 (methanol) to log P = 2.9 (octanol). The density, the radial distribution function, the self-diffusion coefficient and the dielectric constant were evaluated. A long equilibration time of at least 50 ns and a large size of the molecular system of more than 75,000 atoms were used. Except for glycerol, the OPLS all-atom force field reliably reproduced the experimentally determined density with deviations of less than 4% over the whole temperature range. In contrast, the modelled self-diffusion coefficient deviated from its experimental value by up to 55%. To modify the force field, the partial charges of the hydroxyl group were varied by up to 3%. Using the modified OPLS force field, the deviation of the self-diffusion coefficients from their experimental values decreased to less than 19%, while the densities changed by less than 1%.