Statistical Characterization of Food-Derived α-Amylase Inhibitory Peptides: Computer Simulation and Partial Least Squares Regression Analysis

α-Amylase inhibitory peptides are used to treat diabetes, but few studies have statistically characterized their interaction with α-amylase. This study performed the molecular docking of α-amylase with inhibitory peptides from published papers. The key sites, side chain chargeability, and hydrogen bond distribution characteristics were analyzed. Molecular dynamics simulated the role of key sites in complex stability. Moreover, partial least squares regression (PLSR) was used to analyze the contribution of different amino acids in the peptides to inhibition. The results showed that, for the α-amylase molecule, His201 and Gln63, with the highest interaction numbers (INs, 15, 15) and hydrogen bond values (HBVs, 11.50, 10.33), are the key sites on α-amylase, and amino acids with positively charged side chains were important for inhibitory activity. For the inhibitory peptides, Asp and Arg had the highest HBVs, and amino acids with charged side chains were more likely to form hydrogen bonds and exert inhibitory activity. In molecular dynamics simulations, peptides involving key binding sites formed more stable complexes with α-amylase than α-amylase alone, suggesting enhanced inhibitory effects. Further, PLSR results showed that amino acids close to the N-terminus of the inhibitory peptide, located in the third and fifth positions, were significantly correlated with its inhibitory activity. In conclusion, this study provides a new approach to developing and screening α-amylase inhibitors.

Evaluation of action of steroid molecules on SARS-CoV-2 by inhibiting NSP-15, an endoribonuclease

Many countries in the world have recently experienced an outbreak of COVID-19, turned out to be a pandemic which significantly affected the world economy. Among many attempts to treat/control infection or to modulate host immunity, many small molecules including steroids were prescribed based on their use against other viral infection or inflammatory conditions. A recent report established the possibility of usage of a corticosteroid against the virus through inhibiting NSP-15; an mRNA endonuclease of SARS-CoV-2 and thereby viral replication. This study aimed to identify potential anti-viral agents for the virus through computational approaches and to validate binding properties with the protein target through molecular dynamics simulation. Unlike the conventional approaches, dedicated data base of steroid like compounds was used for initial screening along with dexamethasone and cortisone, which are used in the treatment of COVID-19 affected population in some countries. Molecular docking was performed for three compounds filtered from data base in addition to dexamethasone and Cortisone followed by molecular dynamics simulation analysis to validate the dynamics of binding at the active site. In addition, analysis of ADME properties established that these compounds have favorable drug-like properties. Based on docking, molecular dynamics simulation studies and various other trajectory analyses, compounds that are identified could be suggested as therapeutics or precursors towards designing new anti-viral agents against SARS-CoV-2, to combat COVID-19. Also, this is an attempt to study the impact of steroid compounds on NSP-15 of SARS-CoV-2, since many steroid like compounds are used during the treatment of COVID-19 patients.

Molecular mechanism of the association and dissociation of Deltarasin from the heterodimeric KRas4B-PDEδ complex

The formation of the KRas4B-PDEδ complex activates different signaling pathways required for the development and maintenance of cancer. Previous experimental and theoretical studies have allowed researchers to design an inhibitor of the KRas4B-PDEδ complex, "Deltarasin." This inhibitor binds to the prenyl-binding pocket of PDEδ and subsequently inhibits the proliferation of human pancreatic ductal adenocarcinoma cells that depend on oncogenic KRas4B. Nevertheless, structural and energetic information about the inhibitory effects of Deltarasin on the KRas4B-PDEδ complex are not available. In this study, we explore the properties of Deltarasin in inhibiting the formation of wild-type and mutant KRas4B-PDEδ complexes present in different cell lines expressing mutant RAS genes (G12D, G12C, G12V, G13D, Q61L, and Q61R) using 1.7 μs molecular dynamics (MD) simulations in combination with the MMGBSA approach. Our results revealed the energetic and structural mechanisms that suggest a higher affinity of Deltarasin for PDEδ than the farnesylated HVR. Moreover, Deltarasin exerts another dissociative effect by binding to the protein-protein dimeric interface of wild-type KRas4B-PDEδ, whereas associative and dissociative effects were observed for mutant KRas4B-PDEδ, providing a mechanistic explanation for the inhibitory effects of Deltarasin on different cancer cell lines.

Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine

The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure-function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.

Understanding the microscopic binding mechanism of hydroxylated and sulfated polybrominated diphenyl ethers with transthyretin by molecular docking, molecular dynamics simulations and binding free energy calculations

Polybrominated diphenyl ethers (PBDEs), one typical type of persistent environmental contaminant, have toxicological effects such as disrupting thyroid homeostasis in the human body. The high binding affinities of hydroxylated metabolites of PBDEs (OH-PBDEs) with transthyretin (TTR) were considered to be one major reason for their extraordinary capacity of passing through the blood-brain barrier via competitive thyroid hormone (T4) transport protein binding. Recent findings showed that sulfated PBDEs can be formed in human liver cytosol as phase-II metabolites. However, experimentally determined data for the TTR binding potential of the sulfated PBDEs are still not available. Therefore, molecular docking and molecular dynamics (MD) simulations were employed in the present study to probe the molecular basis of TTR interacting with hydroxylated and sulfated PBDEs at the atomic level. The docking scores of LeDock were used to construct the structure-based predictive model. The calculated results showed that the sulfated PBDEs have stronger affinity for TTR than the corresponding OH-PBDEs. Further analysis of structural characteristics based on MD simulations indicated that upon the binding of PBDE metabolites, the stability of TTR was enhanced and the dissociation rate of the tetrameric protein structure was potentially decreased. Subsequent binding free energy calculations implied that van der Waals interactions are the dominant forces for the binding of these metabolites of PBDEs at the T4 site of TTR. The residues Ser117/Ser117' and Lys15/Lys15' were identified, by both residue energy decomposition and computational alanine-scanning mutagenesis methods, as key residues which play an important role in determining the binding orientations of the -OSO₃^- group of sulfated PBDEs by formation of either hydrogen bonds or electrostatic interactions, respectively. In general, the combination of docking calculations with MD simulations provided a theoretically toxicological assessment for the metabolites of PBDEs, deep insight into the recognition mechanism of TTR for these compounds, and thus more comprehensive understanding of the thyroid-related toxic effects of PBDEs as well.

Molecular modeling study of CP-690550 derivatives as JAK3 kinase inhibitors through combined 3D-QSAR, molecular docking, and dynamics simulation techniques

To develop more potent JAK3 kinase inhibitors, a series of CP-690550 derivatives were investigated using combined molecular modeling techniques, such as 3D-QSAR, molecular docking and molecular dynamics (MD). The leave-one-out correlation (q²) and non-cross-validated correlation coefficient (r²) of the best CoMFA model are 0.715 and 0.992, respectively. The q² and r² values of the best CoMSIA model are 0.739 and 0.995, respectively. The steric, electrostatic, and hydrophobic fields played important roles in determining the inhibitory activity of CP-690550 derivatives. Some new JAK3 kinase inhibitors were designed. Some of them have better inhibitory activity than the most potent Tofacitinib (CP-690550). Molecular docking was used to identify some key amino acid residues at the active site of JAK3 protein. 10ns MD simulations were successfully performed to confirm the detailed binding mode and validate the rationality of docking results. The calculation of the binding free energies by MMPBSA method gives a good correlation with the predicted biological activity. To our knowledge, this is the first report on MD simulations and free energy calculations for this series of compounds. The combination results of this study will be valuable for the development of potent and novel JAK3 kinase inhibitors.

Energetic and conformational features linked to the monomeric and dimeric states of bovine BLG

Bovine β-lactoglobulin (BLG) belong to the lipocalin family. This is a group of proteins involved in the binding and transporting of hydrophobic molecules. Experimental and theoretical reports have stated its complex structural behavior in solution, with coupled effects between homodimerization and ligand recognition. Nonetheless, structural evidence at the atomic level about the cause of this coupled effect has not been reported to date. To address this issue microsecond molecular dynamics (MD) simulations were combined with the molecular mechanics generalized Born surface area (MM/GBSA) approach, clustering analysis and principal component analysis (PCA), to explore the conformational complexity of BLG protein-protein self-association and palmitic acid (PLM) or dodecyl sulfate (SDS) ligand recognition in the monomeric and dimeric state. MD simulations, coupled to the MM/GBSA method, revealed that dimerization exerts contrasting effects on the ligand-binding capacity of BLG. Protein dimerization decreases PLM affinity, promoting dimer association. For SDS the dimeric state increases affinity, enhancing dimer dissociation. MD simulations based on PCA revealed that while few differences in the conformational subspace are observed between the free and bound monomer and dimer coupling for PLM, substantial changes are observed between the free and bound monomer and dimer coupling for SDS.

KECSA-Movable Type Implicit Solvation Model (KMTISM)

Computation of the solvation free energy for chemical and biological processes has long been of significant interest. The key challenges to effective solvation modeling center on the choice of potential function and configurational sampling. Herein, an energy sampling approach termed the “Movable Type” (MT) method, and a statistical energy function for solvation modeling, “Knowledge-based and Empirical Combined Scoring Algorithm” (KECSA) are developed and utilized to create an implicit solvation model: KECSA-Movable Type Implicit Solvation Model (KMTISM) suitable for the study of chemical and biological systems. KMTISM is an implicit solvation model, but the MT method performs energy sampling at the atom pairwise level. For a specific molecular system, the MT method collects energies from prebuilt databases for the requisite atom pairs at all relevant distance ranges, which by its very construction encodes all possible molecular configurations simultaneously. Unlike traditional statistical energy functions, KECSA converts structural statistical information into categorized atom pairwise interaction energies as a function of the radial distance instead of a mean force energy function. Within the implicit solvent model approximation, aqueous solvation free energies are then obtained from the NVT ensemble partition function generated by the MT method. Validation is performed against several subsets selected from the Minnesota Solvation Database v2012. Results are compared with several solvation free energy calculation methods, including a one-to-one comparison against two commonly used classical implicit solvation models: MM-GBSA and MM-PBSA. Comparison against a quantum mechanics based polarizable continuum model is also discussed (Cramer and Truhlar’s Solvation Model 12).

Influence of Chirality of Crizotinib on Its MTH1 Protein Inhibitory Activity: Insight from Molecular Dynamics Simulations and Binding Free Energy Calculations

As a promising target for the treatment of lung cancer, the MutT Homolog 1 (MTH1) protein can be inhibited by crizotinib. A recent work shows that the inhibitory potency of (S)-crizotinib against MTH1 is about 20 times over that of (R)-crizotinib. But the detailed molecular mechanism remains unclear. In this study, molecular dynamics (MD) simulations and free energy calculations were used to elucidate the mechanism about the effect of chirality of crizotinib on the inhibitory activity against MTH1. The binding free energy of (S)-crizotinib predicted by the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and Adaptive biasing force (ABF) methodologies is much lower than that of (R)-crizotinib, which is consistent with the experimental data. The analysis of the individual energy terms suggests that the van der Waals interactions are important for distinguishing the binding of (S)-crizotinib and (R)-crizotinib. The binding free energy decomposition analysis illustrated that residues Tyr7, Phe27, Phe72 and Trp117 were important for the selective binding of (S)-crizotinib to MTH1. The adaptive biasing force (ABF) method was further employed to elucidate the unbinding process of (S)-crizotinib and (R)-crizotinib from the binding pocket of MTH1. ABF simulation results suggest that the reaction coordinates of the (S)-crizotinib from the binding pocket is different from (R)-crizotinib. The results from our study can reveal the details about the effect of chirality on the inhibition activity of crizotinib to MTH1 and provide valuable information for the design of more potent inhibitors.

Molecular Dynamics Assisted Mechanistic Study of Isoniazid-Resistance against Mycobacterium tuberculosis InhA

Examination of InhA mutants I16T, I21V, I47T, S94A, and I95P showed that direct and water mediated H-bond interactions between NADH and binding site residues reduced drastically. It allowed conformational flexibility to NADH, particularly at the pyrophosphate region, leading to weakening of its binding at dinucleotide binding site. The highly scattered distribution of pyrophosphate dihedral angles and chi1 side chain dihedral angles of corresponding active site residues therein confirmed weak bonding between InhA and NADH. The average direct and water mediated bridged H-bond interactions between NADH and mutants were observed weaker as compared to wild type. Further, estimated NADH binding free energy in mutants supported the observed weakening of InhA-NADH interactions. Similarly, per residue contribution to NADH binding was also found little less as compared to corresponding residues in wild type. This investigation clearly depicted and supported the effect of mutations on NADH binding and can be accounted for isoniazid resistance as suggested by previous biochemical and mutagenic studies. Further, structural analysis of InhA provided the crucial points to enhance the NADH binding affinity towards InhA mutants in the presence of direct InhA inhibitors to combat isoniazid drug resistance. This combination could be a potential alternative for treatment of drug resistant tuberculosis.

Interaction mechanism exploration of R-bicalutamide/S-1 with WT/W741L AR using molecular dynamics simulations

R-Bicalutamide is a first generation antiandrogen used to treat prostate cancer, which inhibits androgen action by competitively binding to the androgen receptor (AR). However, R-bicalutamide was discovered to exhibit some agonistic properties in clinical application. According to reports, the W741L AR mutation may lead to resistance towards R-bicalutamide. But the mechanism of the R-bicalutamide switch from an antagonist to an agonist due to the mutation of AR W741L is still not so clear. Another molecule, S-1, owing to a very similar structure to R-bicalutamide, is always agonistic to both the wild type and W741L AR. The main difference between these two chemicals is that S-1 has an ether linkage while R-bicalutamide has a sulfonyl group. To study the drug-resistant mechanism caused by W741L mutation and the opposite effects arising from subtle structure differences, molecular dynamics (MD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) calculations were employed to explore the interaction mechanisms between R-bicalutamide/S-1 and WT/W741L AR. The calculated binding free energies are in accordance with the reported experimental values. The obtained results indicate that M895 and W741 are vital amino acids in the antagonism of R-bicalutamide. The bulkier substitution of sulfonyl and tryptophan push aside M895, together with helix 12 (H12), to expose the ligand-binding domain resulting in the antagonistic conformation of the AR. If W741 is mutated to L741, the B-ring of these two chemicals would shift toward L741. At the same time, M895 dragging helix H12, would also move closer to L741. So H12 tends to cover the AR ligand-binding domain to a certain degree, changing the androgen receptor from an antagonistic to an agonistic conformation, which may explain the agonism of R-bicalutamide to the mutant W741L AR.

Factors affecting the interactions between beta-lactoglobulin and fatty acids as revealed in molecular dynamics simulations

Beta-lactoglobulin (BLG), a bovine dairy protein, is a promiscuously interacting protein that can bind multiple hydrophobic ligands. Fatty acids (FAs), common hydrophobic molecules bound to BLG, are important sources of fuel for life because they yield large quantities of ATP when metabolized. The binding affinity increases with the length of the ligands, indicating the importance of the van der Waals (vdW) interactions between the hydrocarbon tail and the hydrophobic calyx of BLG. An exception to this rule is caprylic acid (OCA) which is two-carbon shorter but has a stronger binding affinity than capric acid. Theoretical calculations in the current literature are not accurate enough to shed light on the underlying physics of this exception. The computed affinity values are greater for longer fatty acids without respect for the caprylic exception and those values are generally several orders of magnitude away from the experimental data. In this work, we used hybrid steered molecular dynamics to accurately compute the binding free energies between BLG and the five saturated FAs of 8 to 16 carbon atoms. The computed binding free energies agree well with experimental data not only in rank but also in absolute values. We gained insights into the exceptional behavior of caprylic acid in the computed values of entropy and electrostatic interactions. We found that the electrostatic interaction between the carboxyl group of caprylic acid and the two amino groups of K60/69 in BLG is much stronger than the vdW force between the OCA's hydrophobic tail and the BLG calyx. This pulls OCA to the top of the beta barrel where it is easier to fluctuate, giving rise to greater entropy of OCA at the binding site.