Structural analysis and inhibition capacity of dioscorin protein yam: Theoretical study by molecular docking and dynamic simulations.

Proteases such as trypsins in the gut of Spodoptera frugiperda are responsible for breaking down dietary proteins into amino acids necessary for insect growth and development. In this study, we characterized the insecticidal potential of dioscorin, the storage protein of yam (Dioscorea alata), using molecular docking and molecular dynamics simulations to determine the interactions between trypsin enzymes and the protein inhibitor dioscorin. To achieve this, we used the three-dimensional structures of the trypsin-like digestive enzymes of S. frugiperda, a pest of corn and cotton, as receptors or target molecules. We performed protein-protein docking using Cluspro software, estimation of the binding free energy, and information on the dynamic and time-dependent behavior of dioscorin-trypsin complexes using the NAMD package. Our computational analysis showed that dioscorin can bind to the digestive trypsins of S. frugiperda, as confirmed by the affinity energy values (-1022.4 to -1236.9), stability of the complexes during the simulation trajectory, and binding free energy values between -57.3 and -66.9 kcal/mol. Additionally, dioscorin uses two reactive sites to bind trypsin, but the largest contribution to the interaction energy is made by amino acid residues between amino acid backbone positions 8-14 by hydrogen bonds, hydrophobic, and Van der Waals (VdW) interactions. VdW is the energy that makes the greatest contribution to the binding energy. Collectively, our findings demonstrate, for the first time, the binding capacity of the yam protein dioscorin to the digestive trypsin of S. frugiperda. These promising results suggest a possible bioinsecticide action of dioscorin.

Induced fit for cytochrome P450 3A4 based on molecular dynamics.

The present study aims at numerically describing to what extent substrate - enzyme complexes in solution may change over time as a natural process of conformational changes for a liganded enzyme in comparison to those movements which occur independently from substrate interaction, i.e. without a ligand. To this end, we selected structurally known pairs of liganded / unliganded CYP450 3A4 enzymes with different geometries hinting at induced fit events. We carried out molecular dynamics simulations (MD) comparing the trajectories in a "cross-over" protocol: (i) we added the ligand to the unliganded crystal form which should adopt geometries similar to the known geometry of the liganded crystal structure during MD, and - conversely - (ii) we removed the bound ligand form the known liganded complex to test if a geometry similar to the known unliganded (apo-) form can be adopted during MD. To compare continues changes we measured root means square deviations and frequencies. Results for case (i) hint at larger conformational changes required for accepting the substrate during its approach to final position - in contrast to case (ii) when mobility is fairly reduced by ligand binding (strain energy). In conclusion, a larger conformational sampling prior to ligand binding and the freezing-in (rigidity) of conformations for bound ligands can be interpreted as two conditions linked to induced-fit.

Cooperative hydrogen bonds and mobility of the non-aromatic ring as selectivity determinants for human acetylcholinesterase to similar anti-Alzheimer's galantaminics: a computational study.

Galantamine (Gnt) is a natural alkaloid inhibitor of acetylcholinesterase and is presently one of the most used drugs in the treatment against Alzheimer's disease during both the initial and intermediate stages. Among several natural Gnt derivatives, sanguinine (Sng) and lycoramine (Lyc) attract attention because of the way their subtle chemical differences from Gnt lead to drastic and opposite distinctions in inhibitory effects. However, to date, there is no solved structure for these natural derivatives. In the present study, we applied computational modeling and free energy calculation methods to better elucidate the molecular basis of the subtle distinctions between these derivatives and Gnt. The results showed that differences in the mobility of the non-aromatic ring carried by the Lyc-like sp2-sp3 modification display drastic conformational, vibrational, and entropic penalties at binding compared to Gnt. Additionally, the establishment of a stronger hydrogen bond network added enthalpic advantages for the linkage of the Sng-like methoxy-hydroxy substituted ligands. These results, which suggest an affinity ranking in agreement with that found in the literature, provided insights that are helpful for future planning and development of new anti-Alzheimer's disease drugs.

Molecular dynamics simulation and binding free energy studies of novel leads belonging to the benzofuran class inhibitors of Mycobacterium tuberculosis Polyketide Synthase 13.

In this work, the binding mechanism of new Polyketide Synthase 13 (Pks13) inhibitors has been studied through molecular dynamics simulation and free energy calculations. The drug Tam1 and its analogs, belonging to the benzofuran class, were submitted to 100 ns simulations, and according to the results obtained for root mean square deviation, all the simulations converged from approximately 30 ns. For the analysis of backbone flotation, the root mean square fluctuations were plotted for the Cα atoms; analysis revealed that the greatest fluctuation occurred in the residues that are part of the protein lid domain. The binding free energy value (ΔGbind) obtained for the Tam16 lead molecule was of -51.43 kcal/mol. When comparing this result with the ΔGbind values for the remaining analogs, the drug Tam16 was found to be the highest ranked: this result is in agreement with the experimental results obtained by Aggarwal and collaborators, where it was verified that the IC50 for Tam16 is the smallest necessary to inhibit the Pks13 (IC50 = 0.19 µM). The energy decomposition analysis suggested that the residues which most interact with inhibitors are: Ser1636, Tyr1637, Asn1640, Ala1667, Phe1670, and Tyr1674, from which the greatest energy contribution to Phe1670 was particularly notable. For the lead molecule Tam16, a hydrogen bond with the hydroxyl of the phenol not observed in the other analogs induced a more stable molecular structure. Aggarwal and colleagues reported this hydrogen bonding as being responsible for the stability of the molecule, optimizing its physic-chemical, toxicological, and pharmacokinetic properties.

Insights into structural features of HDAC1 and its selectivity inhibition elucidated by Molecular dynamic simulation and Molecular Docking.

Histone deacetylases (HDACs) are a family of proteins whose main function is the removal of acetyl groups from lysine residues located on histone and non-histone substrates, which regulates gene transcription and other activities in cells. HDAC1 dysfunction has been implicated in cancer development and progression; thus, its inhibition has emerged as a new therapeutic strategy. Two additional metal binding sites (Site 1 and Site 2) in HDACs have been described that are primarily occupied by potassium ions, suggesting a possible structural role that affects HDAC activity. In this work, we explored the structural role of potassium ions in Site 1 and Site 2 and how they affect the interactions of compounds with high affinities for HDAC1 (AC1OCG0B, Chlamydocin, Dacinostat and Quisinostat) and SAHA (a pan-inhibitor) using molecular docking and molecular dynamics (MD) simulations in concert with a Molecular-Mechanics-Generalized-Born-Surface-Area (MMGBSA) approach. Four models were generated: one with a potassium ion (K+) in both sites (HDAC1k), a second with K+ only at site 1 (HDAC1ks1), a third with K+ only at site 2 (HDAC1ks2) and a fourth with no K+ (HDAC1wk). We found that the presence or absence of K+ not only impacted the structural flexibility of HDAC1, but also its molecular recognition, consistent with experimental findings. These results could therefore be useful for further structure-based drug design studies addressing new HDAC1 inhibitors.

Is It Reliable to Take the Molecular Docking Top Scoring Position as the Best Solution without Considering Available Structural Data?

Molecular docking is the most frequently used computational method for studying the interactions between organic molecules and biological macromolecules. In this context, docking allows predicting the preferred pose of a ligand inside a receptor binding site. However, the selection of the "best" solution is not a trivial task, despite the widely accepted selection criterion that the best pose corresponds to the best energy score. Here, several rigid-target docking methods were evaluated on the same dataset with respect to their ability to reproduce crystallographic binding orientations, to test if the best energy score is a reliable criterion for selecting the best solution. For this, two experiments were performed: (A) to reconstruct the ligand-receptor complex by performing docking of the ligand in its own crystal structure receptor (defined as self-docking), and (B) to reconstruct the ligand-receptor complex by performing docking of the ligand in a crystal structure receptor that contains other ligand (defined as cross-docking). Root-mean square deviation (RMSD) was used to evaluate how different the obtained docking orientation is from the corresponding co-crystallized pose of the same ligand molecule. We found that docking score function is capable of predicting crystallographic binding orientations, but the best ranked solution according to the docking energy is not always the pose that reproduces the experimental binding orientation. This happened when self-docking was achieved, but it was critical in cross-docking. Taking into account that docking is typically used with predictive purposes, during cross-docking experiments, our results indicate that the best energy score is not a reliable criterion to select the best solution in common docking applications. It is strongly recommended to choose the best docking solution according to the scoring function along with additional structural criteria described for analogue ligands to assure the selection of a correct docking solution.

The C-terminal extension of bacterial flavodoxin-reductases: involvement in the hydride transfer mechanism from the coenzyme.

To study the role of the mobile C-terminal extension present in bacterial class of plant type NADP(H):ferredoxin reductases during catalysis, we generated a series of mutants of the Rhodobacter capsulatus enzyme (RcFPR). Deletion of the six C-terminal amino acids beyond alanine 266 was combined with the replacement A266Y, emulating the structure present in plastidic versions of this flavoenzyme. Analysis of absorbance and fluorescence spectra suggests that deletion does not modify the general geometry of FAD itself, but increases exposure of the flavin to the solvent, prevents a productive geometry of FAD:NADP(H) complex and decreases the protein thermal stability. Although the replacement A266Y partially coats the isoalloxazine from solvent and slightly restores protein stability, this single change does not allow formation of active charge-transfer complexes commonly present in the wild-type FPR, probably due to restraints of C-terminus pliability. A proton exchange process is deduced from ITC measurements during coenzyme binding. All studied RcFPR variants display higher affinity for NADP(+) than wild-type, evidencing the contribution of the C-terminus in tempering a non-productive strong (rigid) interaction with the coenzyme. The decreased catalytic rate parameters confirm that the hydride transfer from NADPH to the flavin ring is considerably hampered in the mutants. Although the involvement of the C-terminal extension from bacterial FPRs in stabilizing overall folding and bent-FAD geometry has been stated, the most relevant contributions to catalysis are modulation of coenzyme entrance and affinity, promotion of the optimal geometry of an active complex and supply of a proton acceptor acting during coenzyme binding.

Structure-function studies on jaburetox, a recombinant insecticidal peptide derived from jack bean (Canavalia ensiformis) urease.

BACKGROUND: Ureases are metalloenzymes involved in defense mechanisms in plants. The insecticidal activity of Canavalia ensiformis (jack bean) ureases relies partially on an internal 10kDa peptide generated by enzymatic hydrolysis of the protein within susceptible insects. A recombinant version of this peptide, jaburetox, exhibits insecticidal, antifungal and membrane-disruptive properties. Molecular modeling of jaburetox revealed a prominent ß-hairpin motif consistent with either neurotoxicity or pore formation. METHODS: Aiming to identify structural motifs involved in its effects, mutated versions of jaburetox were built: 1) a peptide lacking the ß-hairpin motif (residues 61-74), JbtxΔ-ß; 2) a peptide corresponding the N-terminal half (residues 1-44), Jbtx N-ter, and 3) a peptide corresponding the C-terminal half (residues 45-93), Jbtx C-ter. RESULTS: 1) JbtxΔ-ß disrupts liposomes, and exhibited entomotoxic effects similar to the whole peptide, suggesting that the ß-hairpin motif is not a determinant of these biological activities; 2) both Jbtx C-ter and Jbtx N-ter disrupted liposomes, the C-terminal peptide being the most active; and 3) while Jbtx N-ter persisted to be biologically active, Jbtx C-ter was less active when tested on different insect preparations. Molecular modeling and dynamics were applied to the urease-derived peptides to complement the structure-function analysis. MAJOR CONCLUSIONS: The N-terminal portion of the Jbtx carries the most important entomotoxic domain which is fully active in the absence of the ß-hairpin motif. Although the ß-hairpin contributes to some extent, probably by interaction with insect membranes, it is not essential for the entomotoxic properties of Jbtx. GENERAL SIGNIFICANCE: Jbtx represents a new type of insecticidal and membrane-active peptide.

Identification of a vaccine against schistosomiasis using bioinformatics and molecular modeling tools.

Schistosomiasis is a serious public health problem in Brazil and worldwide. Although the drugs used to treatment schistosomiasis are effective, the disease continues to expand in all endemic countries due to constant reinfection, poor sanitation, and the lack of effective programs for disease control. However, advances generated through genome projects have provided important information that has improved the understanding of the biology of this parasite. These advances, associated with the advent of bioinformatic analysis, are becoming an important tool in reverse vaccinology. Through database access to the DNA and protein sequences of Schistosoma mansoni and the use of bioinformatics programs, fourteen epitopes were identified. Five epitopes were obtained from proteins whose immunogenic potential had already been assessed in other studies (KP), and nine whose immunogenic potential is unknown (UP). To improve stimulation of the host immune system, the selected epitopes were modeled with a sugar moiety. After this addition, all of the epitopes showed structures similar to those observed in the native proteins, but only eleven of the peptides presented thermodynamically stable structures. Prediction analysis and molecular modeling showed that the glycopeptides presented here are important targets in the search for a vaccine against schistosomiasis. Additionally, they suggest that these molecules may be used in immunological assays to evaluate the level of protection, the effect on pathology reduction and the profile of cytokines and antibodies induced by them.

The structure of the CD3 ζζ transmembrane dimer in POPC and raft-like lipid bilayer: a molecular dynamics study.

Plasma membrane lipids significantly affect assembly and activity of many signaling networks. The present work is aimed at analyzing, by molecular dynamics simulations, the structure and dynamics of the CD3 ζζ dimer in palmitoyl-oleoyl-phosphatidylcholine bilayer (POPC) and in POPC/cholesterol/sphingomyelin bilayer, which resembles the raft membrane microdomain supposed to be the site of the signal transducing machinery. Both POPC and raft-like environment produce significant alterations in structure and flexibility of the CD3 ζζ with respect to nuclear magnetic resonance (NMR) model: the dimer is more compact, its secondary structure is slightly less ordered, the arrangement of the Asp6 pair, which is important for binding to the Arg residue in the alpha chain of the T cell receptor (TCR), is stabilized by water molecules. Different interactions of charged residues with lipids at the lipid-cytoplasm boundary occur when the two environments are compared. Furthermore, in contrast to what is observed in POPC, in the raft-like environment correlated motions between transmembrane and cytoplasmic regions are observed. Altogether the data suggest that when the TCR complex resides in the raft domains, the CD3 ζζ dimer assumes a specific conformation probably necessary to the correct signal transduction.