Network analysis of dynamically important residues in protein structures mediating ligand-binding conformational changes

Discovery of ellagic acid as a competitive inhibitor of Src homology phosphotyrosyl phosphatase 2 (SHP2) for cancer treatment: In vitro and in silico study

Targeting SHP2 has become a potential cancer treatment strategy. In this study, ellagic acid was first reported as a competitive inhibitor of SHP2, with an IC50 value of 0.69 ± 0.07 μM, and its inhibitory potency was 34.86 times higher that of the positive control NSC87877. Ellagic acid also had high inhibitory activity on the SHP2-E76K and SHP2-E76A mutants, with the IC50 values of 1.55 ± 0.17 μM and 0.39 ± 0.05 μM, respectively. Besides, the IC50 values of ellagic acid on homologous proteins SHP1, PTP1B, and TCPTP were 0.93 ± 0.08 μM, 2.04 ± 0.28 μM, and 11.79 ± 0.83 μM, with selectivity of 1.35, 2.96, and 17.09 times, respectively. The CCK8 proliferation experiment exhibited that ellagic acid would inhibit the proliferation of various cancer cells. It was worth noting that the combination of ellagic acid and KRASG12C inhibitor AMG510 would produce a strong synergistic effect in inhibiting NCI-H358 cells. Western blot experiment exhibited that ellagic acid would downregulate the phosphorylation levels of Erk and Akt in NCI-H358 and MDA-MB-468 cells. Molecular docking and molecular dynamics studies revealed the binding information between SHP2 and ellagic acid. In summary, this study provides new ideas for the development of SHP2 inhibitors.

Intermediate-aided allostery mechanism for α-glucosidase by Xanthene-11v as an inhibitor using residue interaction network analysis

Exploring allosteric inhibition and the discovery of new inhibitor binding sites are important studies in protein regulation mechanisms and drug discovery. Structural and network-based analyses of trajectories resulting from molecular dynamics (MD) simulations have been developed to discover protein dynamics, landscape, functions, and allosteric regions. Here, an experimentally suggested non-competitive inhibitor, xanthene-11v, was considered to explore its allosteric inhibition mechanism in α-glucosidase MAL12. Comparative structural and network analyses were applied to eight 250 ns independent MD simulations, four of which were performed in the free state and four of which were performed in ligand-bound forms. Projected two-dimensional free energy landscapes (FEL) were constructed from the probabilistic distribution of conformations along the first two principal components. The post-simulation analyses of the coordinates, side-chain torsion angles, non-covalent interaction networks, network communities, and their centralities were performed on α-glucosidase conformations and the intermediate sub-states. Important communities of residues have been found that connect the allosteric site to the active site. Some of these residues like Thr307, Arg312, TYR344, ILE345, Phe357, Asp406, Val407, Asp408, and Leu436 are the key messengers in the transition pathway between allosteric and active sites. Evaluating the probability distribution of distances between gate residues including Val407 in one community and Phe158, and Pro65 in another community depicted the closure of this gate due to the inhibitor binding. Six macro states of protein were deduced from the topology of FEL and analysis of conformational preference of free and ligand-bound systems to these macro states shows a combination of lock-and-key, conformational selection, and induced fit mechanisms are effective in ligand binding. All these results reveal structural states, allosteric mechanisms, and key players in the inhibition pathway of α-glucosidase by xanthene-11v.

Evolutionary progression of collective mutations in Omicron sub-lineages towards efficient RBD-hACE2: Allosteric communications between and within viral and human proteins

The interaction between the Spike (S) protein of SARS-CoV-2 and the human angiotensin converting enzyme 2 (hACE2) is essential for infection, and is a target for neutralizing antibodies. Consequently, selection of mutations in the S protein is expected to be driven by the impact on the interaction with hACE2 and antibody escape. Here, for the first time, we systematically characterized the collective effects of mutations in each of the Omicron sub-lineages (BA.1, BA.2, BA.3 and BA.4) on both the viral S protein receptor binding domain (RBD) and the hACE2 protein using post molecular dynamics studies and dynamic residue network (DRN) analysis. Our analysis suggested that Omicron sub-lineage mutations result in altered physicochemical properties that change conformational flexibility compared to the reference structure, and may contribute to antibody escape. We also observed changes in the hACE2 substrate binding groove in some sub-lineages. Notably, we identified unique allosteric communication paths in the reference protein complex formed by the DRN metrics betweenness centrality and eigencentrality hubs, originating from the RBD core traversing the receptor binding motif of the S protein and the N-terminal domain of the hACE2 to the active site. We showed allosteric changes in residue network paths in both the RBD and hACE2 proteins due to Omicron sub-lineage mutations. Taken together, these data suggest progressive evolution of the Omicron S protein RBD in sub-lineages towards a more efficient interaction with the hACE2 receptor which may account for the increased transmissibility of Omicron variants.

Allosteric pockets and dynamic residue network hubs of falcipain 2 in mutations including those linked to artemisinin resistance

Continually emerging resistant strains of malarial parasites to current drugs present challenges. Understanding the underlying resistance mechanisms, especially those linked to allostery is, thus, highly crucial for drug design. This forms the main concern of the paper through a case study of falcipain 2 (FP-2) and its mutations, some of which are linked to artemisinin (ART) drug resistance. Here, we applied a variety of in silico approaches and tools that we developed recently, together with existing computational tools. This included novel essential dynamics and dynamic residue network (DRN) analysis algorithms. We identified six pockets demonstrating dynamic differences in the presence of some mutations. We observed striking allosteric effects in two mutant proteins. In the presence of M245I, a cryptic pocket was detected via a unique mechanism in which Pocket 2 fused with Pocket 6. In the presence of the A353T mutation, which is located at Pocket 2, the pocket became the most rigid among all protein systems analyzed. Pocket 6 was also highly stable in all cases, except in the presence of M245I mutation. The effect of ART linked mutations was more subtle, and the changes were at residue level. Importantly, we identified an allosteric communication path formed by four unique averaged BC hubs going from the mutated residue to the catalytic site and passing through the interface of three identified pockets. Collectively, we established and demonstrated that we have robust tools and a pipeline that can be applicable to the analysis of mutations.

New Insight Into the Structure-Activity Relationship of Sweet-Tasting Proteins: Protein Sector and Its Role for Sweet Properties

Sweet-tasting protein is a kind of biomacromolecule that has remarkable sweetening power and is regarded as the promising sugar replacer in the future. Some sweet-tasting proteins has been used in foods and beverages. However, the structure and function relationship of these proteins is still elusive, and guidelines for their protein engineering is limited. It is well-known that the sweet-tasting proteins bind to and activate the sweet taste receptor T1R2/T1R3, thus eliciting their sweetness. The “wedge-model” for describing the interaction between sweet-tasting proteins and sweet taste receptor to elucidate their sweetness has been reported. In this perspective article, we revealed that the intramolecular interaction forces in sweet-tasting proteins is directly correlated to their properties (sweetness and stability). This intramolecular interaction pattern, named as “protein sector,” refers to a small subset of residues forming physically connections, which cooperatively affect the function of the proteins. Based on the analysis of previous experimental data, we suggest that “protein sector” of sweet-tasting proteins is pivotal for their sweet properties, which are meaningful guidelines for the future protein engineering.

Structural and conformational changes induced by missense variants in the zinc finger domains of GATA3 involved in breast cancer

Breast cancer (BC) is the main cancer in women having multiple receptor based tumour subtypes. Large scale genome sequencing studies of BC have identified several genes among which GATA3 is reported as a highly mutated gene followed by TP53 and PIK3CA. GATA3 is a crucial transcription factor, and was initially identified as a DNA-binding protein involved in the regulation of immune cell functions. Different missense mutations in the region of the DNA-binding domain of GATA3 are associated with BC and other neoplastic disorders. In this study, computational based approaches have been exploited to reveal associations of various mutations on structure, stability, conformation and function of GATA3. Our findings have suggested that, all analysed missense mutations were deleterious and highly pathogenic in nature. A molecular dynamics simulation study showed that all mutations led to structural destabilisation by reducing protein globularity and flexibility, by altering secondary structural configuration and decreasing protein ligand stability. Essential dynamics analysis indicated that mutations in GATA3 decreased protein mobility and increased its conformational instability. Furthermore, residue network analysis showed that the mutations affected the signal transduction of important residues that potentially influenced GATA3-DNA binding. The present study highlights the importance of different variants of GATA3 which have potential impact on neoplastic progression in breast cancer and may facilitate development of precise and personalized therapeutics.

Mutations in the N- and C-finger domains of GATA3 lead to breast cancer.

Identification of novel inhibitor of protein tyrosine phosphatases delta: structure-based pharmacophore modeling, virtual screening, flexible docking, molecular dynamics simulation, and post-molecular dynamics analysis

Owing to their unique functions in regulating the synapse activity of protein tyrosine phosphatases delta (PTPδ) that has drawn special attention for developing drugs to autism spectrum disorders (ASDs). In this study, the PTPδ pharmacophore was first established by the structure-based pharmacophore method. Subsequently, 10 compounds contented Lipinski's rule of five was acquired by the virtual screening of the PTPδ pharmacophore against ZINC and PubChem databases. Then, the 10 identified molecules were discovered that had better binding affinity than a known PTPδ inhibitors compound SCHEMBL16375396. Two compounds SCHEMBL16375408 and ZINC19796658 with high binding score, low toxicity were gained. They were observed by docking analysis and molecular dynamics simulations that the novel potential inhibitors not only possessed the same function as SCHEMBL16375396 did in inhibiting PTPδ, but also had more favorable conformation to bind with the catalytic active regions. This study provides a new method for identify PTPδ inhibitor for the treatment of ASDs disease.Communicated by Ramaswamy H. Sarma.