Design and development of high affinity dual anticancer peptide-inhibitors against p53-MDM2/X interaction

Insights into the Interaction Mechanisms of Peptide and Non-Peptide Inhibitors with MDM2 Using Gaussian-Accelerated Molecular Dynamics Simulations and Deep Learning

Inhibiting MDM2-p53 interaction is considered an efficient mode of cancer treatment. In our current study, Gaussian-accelerated molecular dynamics (GaMD), deep learning (DL), and binding free energy calculations were combined together to probe the binding mechanism of non-peptide inhibitors K23 and 0Y7 and peptide ones PDI6W and PDI to MDM2. The GaMD trajectory-based DL approach successfully identified significant functional domains, predominantly located at the helixes α2 and α2', as well as the β-strands and loops between α2 and α2'. The post-processing analysis of the GaMD simulations indicated that inhibitor binding highly influences the structural flexibility and collective motions of MDM2. Calculations of molecular mechanics-generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) not only suggest that the ranking of the calculated binding free energies is in agreement with that of the experimental results, but also verify that van der Walls interactions are the primary forces responsible for inhibitor-MDM2 binding. Our findings also indicate that peptide inhibitors yield more interaction contacts with MDM2 compared to non-peptide inhibitors. Principal component analysis (PCA) and free energy landscape (FEL) analysis indicated that the piperidinone inhibitor 0Y7 shows the most pronounced impact on the free energy profiles of MDM2, with the piperidinone inhibitor demonstrating higher fluctuation amplitudes along primary eigenvectors. The hot spots of MDM2 revealed by residue-based free energy estimation provide target sites for drug design toward MDM2. This study is expected to provide useful theoretical aid for the development of selective inhibitors of MDM2 family members.

The Synthetic Antimicrobial Peptide Derived From Melittin Displays Low Toxicity and Anti-infectious Properties

The low stability and nonspecific toxicity are the main limiting factors for the clinical applications of melittin (MLT). This study aimed to design and synthesize new analogs of MLT to increase stability, reduce toxicity, and retain their antimicrobial properties against bacterial pathogens. At first, peptide analogs were designed computationally by inducing single mutations in MLT peptides and evaluating their physicochemical properties. The stability of the analogs with the highest scores was determined by Gromacs software. In vitro assays were performed to examine the antimicrobial activity and toxicity of the selected analogs. Two peptide analogs, M1 and M2, were selected based on the SVM score in cell PPD. The M1 analog was created by replacing alanine with leucine on the 15th. The M2 analog was designed by substituting alanine with leucine and isoleucine with arginine at the 15th and 17th positions. According to the Gromacs results, the M2 peptide indicated more stability. RMSD and RMSF results showed no undesirable fluctuations during the 200 ns MD simulation. The MIC and MBC values for the M1 peptide were calculated in a range of 8-128 μg/ml, while the M2 peptide limited the bacterial growth to 32-128 μg/mL. Both peptides indicated less toxicity than natural MLT, based on MTT assay results. The hemolytic activity of the M1 analog was more than M2 at 16 μg/mL concentration. M1 peptide displayed the highest selectivity index against S. aureus and A. baumannii, which were approximately 5.27-fold improvements compared to MLT. In conclusion, we introduced two analogs of MLT with low toxicity, low hemolytic activity, and higher stability, along with retaining antimicrobial properties against gram-negative and positive bacteria.

Effect of flavonoids on the destabilization of α-synuclein fibrils and their conversion to amorphous aggregate: A molecular dynamics simulation and experimental study

The second most prevalent neurodegenerative disease, Parkinson's disease (PD), is caused by the accumulation and deposition of fibrillar aggregates of the α-Syn into the Lewy bodies. To create a potent pharmacological candidate to destabilize the preformed α-Syn fibril, it is important to understand the precise molecular mechanism underlying the destabilization of the α-Syn fibril. Through molecular dynamics simulations and experiments, we have examined the molecular mechanisms causing the destabilization and suppression of a newly discovered α-Syn fibril with a Greek-key-like shape and an aggregation prone state (APS) of α-Syn in the presence and absence of various Flvs. According to MD simulation and experimental evidence, morin, quercetin, and myricetin are the Flvs, most capable of destabilizing the fibrils and converting them into amorphous aggregates. Compared to galangin and kaempferol, they have more hydroxyl groups and form more hydrogen bonds with fibrils.The processes by which morin and myricetin prevent new fibril production from APS and destabilize the fibrils are different. According to linear interaction energy analysis, van der Waals interaction predominates with morin, and electrostatic interaction dominates with myricetin. Our MD simulation and experimental findings provide mechanistic insights into how Flvs destabilize α-Syn fibrils and change their morphology, opening the door to developing structure-based drugs for treating Parkinson's disease.

Rational design of an anti-cancer peptide inhibiting CD147 / Cyp A interaction

The CD147 / Cyp A interaction is a critical pathway in cancer types and an essential factor in entering the COVID-19 virus into the host cell. Melittin acts as an inhibitory peptide in cancer types by blocking the CD147/ Cyp A interaction. The clinical application of Melittin is limited due to weak penetration into cancer cells. TAT is an arginine-rich peptide with high penetration ability into cells widely used in drug delivery systems. This study aimed to design a hybrid peptide derived from Melittin and TAT to inhibit CD147 /Cyp A interaction. An amino acid region with high anti-cancer activity in Melittin was selected based on the physicochemical properties. Based on the results, a truncated Melittin peptide with 15 amino acids by the GGGS linker was fused to a TAT peptide (nine amino acids) to increase the penetration rate into the cell. A new hybrid peptide analog(TM) was selected by replacing the glycine with serine based on random point mutation. Docking results indicated that the TM peptide acts as an inhibitory peptide with high binding energy when interacting with CD147 and the CypA proteins. RMSD and RMSF results confirmed the high stability of the TM peptide in interaction with CD147. Also, the coarse-grained simulation showed the penetration potential of TM peptide into the DOPS-DOPC model membrane. Our findings indicated that the designed multifunctional peptide could be an attractive therapeutic candidate to halter tumor types and COVID-19 infection.

Rational Design of Anti-Angiogenic Peptides to Inhibit VEGF/VEGFR2 Interactions for Cancer Therapeutics

BACKGROUND

Angiogenesis is a critical physiological process that plays a key role in tumor progression, metastatic dissemination, and invasion. In the last two decades, the vascular endothelial growth factor (VEGF) signaling pathway has been the area of extensive researches. VEGF executes its special effects by binding to vascular endothelial growth factor receptors (VEGFRs), particularly VEGFR-2.

OBJECTIVE

The inhibition of VEGF/VEGFR2 interaction is known as an effective cancer therapy strategy. The current study pointed to design and model an anti-VEGF peptide based on VEGFR2 binding regions.

METHOD

The large-scale peptide mutation screening was used to achieve a potent peptide with high binding affinity to VEGF for possible application in inhibition of VEGF/VEGFR2 interaction. The AntiCP and Peptide Ranker servers were used to generate the possible peptides library with anticancer activities and prediction of peptides bioactivity. Then, the interaction of VEGF and all library peptides were analyzed using Hex 8.0.0 and ClusPro tools. A number of six peptides with favorable docking scores were achieved. All of the best docking scores of peptides in complexes with VEGF were evaluated to confirm their stability, using molecular dynamics simulation (MD) with the help of the GROMACS software package.

RESULTS

As a result, two antiangiogenic peptides with 13 residues of PepA (NGIDFNRDFFLGL) and PepC (NGIDFNRDKFLFL) were achieved and introduced to inhibit VEGF/VEGFR2 interactions.

CONCLUSIONS

In summary, this study provided new insights into peptide-based therapeutics development for targeting VEGF signaling pathway in tumor cells. PepA and PepC are recommended as potentially promising anticancer agents for further experimental evaluations.

The p53 Protein Family in the Response of Tumor Cells to Ionizing Radiation: Problem Development

Survival mechanisms are activated in tumor cells in response to therapeutic ionizing radiation. This reduces a treatment's effectiveness. The p53, p63, and p73 proteins belonging to the family of proteins that regulate the numerous pathways of intracellular signal transduction play a key role in the development of radioresistance. This review analyzes the p53-dependent and p53-independent mechanisms involved in overcoming the resistance of tumor cells to radiation exposure.

Plausible blockers of Spike RBD in SARS-CoV2-molecular design and underlying interaction dynamics from high-level structural descriptors

Abstract

COVID-19 is characterized by an unprecedented abrupt increase in the viral transmission rate (SARS-CoV-2) relative to its pandemic evolutionary ancestor, SARS-CoV (2003). The complex molecular cascade of events related to the viral pathogenicity is triggered by the Spike protein upon interacting with the ACE2 receptor on human lung cells through its receptor binding domain (RBD_Spike). One potential therapeutic strategy to combat COVID-19 could thus be limiting the infection by blocking this key interaction. In this current study, we adopt a protein design approach to predict and propose non-virulent structural mimics of the RBD_Spike which can potentially serve as its competitive inhibitors in binding to ACE2. The RBD_Spike is an independently foldable protein domain, resilient to conformational changes upon mutations and therefore an attractive target for strategic re-design. Interestingly, in spite of displaying an optimal shape fit between their interacting surfaces (attributed to a consequently high mutual affinity), the RBD_Spike–ACE2 interaction appears to have a quasi-stable character due to a poor electrostatic match at their interface. Structural analyses of homologous protein complexes reveal that the ACE2 binding site of RBD_Spike has an unusually high degree of solvent-exposed hydrophobic residues, attributed to key evolutionary changes, making it inherently “reaction-prone.” The designed mimics aimed to block the viral entry by occupying the available binding sites on ACE2, are tested to have signatures of stable high-affinity binding with ACE2 (cross-validated by appropriate free energy estimates), overriding the native quasi-stable feature. The results show the apt of directly adapting natural examples in rational protein design, wherein, homology-based threading coupled with strategic “hydrophobic ↔ polar” mutations serve as a potential breakthrough.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s00894-021-04779-0.

The role of miRNAs in MDMX-p53 interplay

MicroRNAs (miRNAs) are endogenous noncoding RNAs of 19-24 nucleotides in length and are tightly related to tumorigenesis and progression. Recent studies have demonstrated that the tumor suppressor p53 and its negative controller MDMX are regulated by miRNAs in different ways. Some miRNAs directly target p53 and regulate its expression and function, whereas some miRNAs target MDMX and regulate p53's activity indirectly. The overexpression of several miRNAs can restore the activity of p53 by negatively regulating MDMX in cancer cells. Therefore, a better understanding of the miRNAs-MDMX-p53 network will put forward potential research directions for developing anticancer therapeutics. In the present review, we mainly focus on the regulatory effects of miRNAs on the MDMX-p53 interplay as well as the role of the miRNAs-MDMX-p53 network in human cancer.

Probing biased activation of mu-opioid receptor by the biased agonist PZM21 using all atom molecular dynamics simulation

Morphine is a commonly used opioid drug to treat acute pain by binding to the mu-opioid receptor (MOR), but its effective analgesic efficacy via triggering of the heterotrimeric G_i protein pathway is accompanied by a series of adverse side effects via triggering of the β-arrestin pathway. Recently, PZM21, a recently developed MOR biased agonist, shows preferentially activating the G protein pathway over β-arrestin pathway. However, there is no high-resolution receptor structure in complex with PZM21 and its action mechanism remains elusive. In this study, PZM21 and Morphine were docked to the active human MOR-1 homology structure and then subjected to the molecular dynamics (MD) simulations in two different situations (i.e., one situation includes the crystal waters but another does not). Detailed comparisons between the two systems were made to characterize the differences in protein-ligand interactions, protein secondary and tertiary structures and dynamics networks. PZM21 could strongly interact with Y328^7.43 of TM7, besides the residues (Asp149^3.32 and Tyr150^3.33) of TM3. The two systems' network paths to the intracellular end of TM6 were roughly similar but the paths to the end of TM7 were different. The PZM21-bound MOR's intracellular ends of TM5-7 bent outward more along with the distance changes of the three key molecular switches (ionic lock, transmission and Tyr toggle) and the distance increase of some conserved inter-helical residue pairs. The larger intracellular opening of the receptor could potentially facilitate G protein binding.

Simulations on the dual effects of flavonoids as suppressors of Aβ42 fibrillogenesis and destabilizers of mature fibrils

Structural studies of the aggregation inhibition of the amyloid-β peptide (Aβ) by different natural compounds are of the utmost importance due to their great potential as neuroprotective and therapeutic agents for Alzheimer’s disease. We provided the simulation of molecular dynamics for two different states of Aβ42, including “monomeric aggregation-prone state (APS)” and “U-shaped pentamers of amyloidogenic protofilament intermediates” in the absence and presence of polyphenolic flavonoids (Flvs, myricetin and morin) in order to verify the possible mechanism of Flvs fibrillogenesis suppression. Data showed that Flvs directly bind into Aβ42 species in both states of “monomeric APS β-sheets” and “pentameric amyloidogenic intermediates”. Binding of Flvs with amyloidogenic protofilament intermediates caused the attenuation of some inter-chains H-bonds, salt bridges, van der Waals and interpeptide interaction energies without interfering with their secondary β-sheets. Therefore, Flvs redirect oligomeric amyloidogenic intermediates into unstructured aggregates by significant disruption of the "steric zipper" motif of fibrils—pairs of self-complementary β-sheets—without changing the amount of β-sheets. It is while Flvs completely destruct the disadvantageous secondary β-sheets of monomeric APS conformers by converting them into coil/helix structures. It means that Flvs suppress the fibrillogenesis process of the monomeric APS structures by converting their β-sheets into proper soluble coil/helices structures. The different actions of Flvs in contact with two different states of Aβ conformers are related to high interaction tendency of Flvs with additional H-bonds for monomeric APS β-sheet, rather than oligomeric protofilaments. Linear interaction energy (LIE) analysis confirmed the strong binding of monomeric Aβ-Flvs with more negative ∆G_binding, rather than oligomeric Aβ-Flvs system. Therefore, atomic scale computational evaluation of Flvs actions demonstrated different dual functions of Flvs, concluded from the application of two different monomeric and pentameric Aβ42 systems. The distinct dual functions of Flvs are proposed as suppressing the aggregation by converting β-sheets of monomeric APS to proper soluble structures and disrupting the "steric zipper" fibril motifs of oligomeric intermediate by converting on-pathway into off-pathway. Taken together, our data propose that Flvs exert dual and more effective functions against monomeric APS (fibrillogenesis suppression) and remodel the Aβ aggregation pathway (fibril destabilization).

Target molecular treatment markers in Intrahepatic Cholangiocarcinoma based on Chinese population

BACKGROUND

As more and more molecular markers have been identified in Intrahepatic Cholangiocarcinoma (ICC), target treatments are promising all around the world. However, geographical and ethnic variations in the ICC epidemiology suggest different genetic variance prevalence in western and eastern countries.

METHODS

Six genetic variations in Chinese ICC populations were analyzed by fluorescent in situ hybridization (FISH) or Sanger sequencing, listed as IDH1/2 mutation, FGFR2 translocation, NTRK1 amplification, MDM2 amplification, HER2 amplification and MET amplification, all of which have corresponding target drugs; meanwhile, they were compared with these gene prevalence in Spanish population based on the cBioPortal database.

RESULTS

The incidences of IDH1/2 mutation, FGFR2 translocation, NTRK1 amplification, MDM2 amplification, HER2 amplification and MET amplification were 29.5 %, 12.9 %, 1.51 %, 2.27 %, 3.03 % and 0.75 %, respectively, in the Spanish population and 7.14 %, 5.71 %, 7.86 %, 5.71 %,4.29 % and 2.14 %, respectively, in the Chinese population. For the gene NTRK1, 11 samples showed signal apart companied amplified using FISH break-apart probes but none of them demonstrated genetic fusion by next-generation sequencing. As to clinicopathological characteristics, patients carrying IDH1/2 mutation showed longer overall survival in the Chinese population, while those carrying FGFR2 translocation tended to be younger in the Spanish population. For HER2, MDM2 and MET, gene amplification predicted protein high-expression, whereas FGFR2 translocation and NTRK1 amplification did not predict protein high-expression.

CONCLUSIONS

Although many target drugs have been speeded up for approval such as BGJ398 for FGFR2 fusion positive ICC patients in western countries, the beneficiary populations are very small in China. The regular target drug such as trastuzumab for HER2 amplification and Crizotinib for MET amplification may be potential candidates in target treatment based on the Chinese population.