Combined molecular dynamics and continuum solvent approaches (MM-PBSA/GBSA) to predict noscapinoid binding to γ-tubulin dimer

An allosteric mechanism for potent inhibition of SARS-CoV-2 main proteinase

The main proteinase (Mpro) of SARS-CoV-2 plays a critical role in cleaving viral polyproteins into functional proteins required for viral replication and assembly, making it a prime drug target for COVID-19. It is well known that noncompetitive inhibition offers potential therapeutic options for treating COVID-19, which can effectively reduce the likelihood of cross-reactivity with other proteins and increase the selectivity of the drug. Therefore, the discovery of allosteric sites of Mpro has both scientific and practical significance. In this study, we explored the binding characteristics and inhibiting process of Mpro activity by two recently reported allosteric inhibitors, pelitinib and AT7519 which were obtained by the X-ray screening experiments, to probe the allosteric mechanism via molecular dynamic (MD) simulations. We found that pelitinib and AT7519 can stably bind to Mpro far from the active site. The binding affinity is estimated to be -24.37 ± 4.14 and - 26.96 ± 4.05 kcal/mol for pelitinib and AT7519, respectively, which is considerably stable compared with orthosteric drugs. Furthermore, the strong binding caused clear changes in the catalytic site of Mpro, thus decreasing the substrate accessibility. The community network analysis also validated that pelitinib and AT7519 strengthened intra- and inter-domain communication of Mpro dimer, resulting in a rigid Mpro, which could negatively impact substrate binding. In summary, our findings provide the detailed working mechanism for the two experimentally observed allosteric sites of Mpro. These allosteric sites greatly enhance the 'druggability' of Mpro and represent attractive targets for the development of new Mpro inhibitors.

Structural and energetic features of the dimerization of the main proteinase of SARS-CoV-2 using molecular dynamic simulations

The COVID-19 pandemic caused by SARS-CoV-2 has been declared a global health crisis. The development of anti-SARS-CoV-2 drugs heavily depends on the systematic study of the critical biological processes of key proteins of coronavirus among which the main proteinase (M^pro) dimerization is a key step for virus maturation. Because inhibiting the M^pro dimerization can efficiently suppress virus maturation, the key residues that mediate dimerization can be treated as targets of drug and antibody developments. In this work, the structure and energy features of the M^pro dimer of SARS-CoV-2 and SARS-CoV were studied using molecular dynamics (MD) simulations. The free energy calculations using the Generalized Born (GB) model showed that the dimerization free energy of the SARS-CoV-2 M^pro dimer (-107.5 ± 10.89 kcal mol^-1) is larger than that of the SARS-CoV M^pro dimer (-92.83 ± 9.81 kcal mol^-1), indicating a more stable and possibly a quicker formation of the M^pro dimer of SARS-CoV-2. In addition, the energy decomposition of each residue revealed 11 key attractive residues. Furthermore, Thr285Ala weakens the steric hindrance between the two protomers of SARS-CoV-2 that can form more intimate interactions. It is interesting to find 11 repulsive residues which effectively inhibit the dimerization process. At the interface of the M^pro dimer, we detected three regions that are rich in interfacial water which stabilize the SARS-CoV-2 M^pro dimer by forming hydrogen bonds with two protomers. The key residues and rich water regions provide important targets for the future design of anti-SARS-CoV-2 drugs through inhibiting M^pro dimerization.

Dysregulation of Microtubule Nucleating Proteins in Cancer Cells

Simple Summary

The dysfunction of microtubule nucleation in cancer cells changes the overall cytoskeleton organization and cellular physiology. This review focuses on the dysregulation of the γ-tubulin ring complex (γ-TuRC) proteins that are essential for microtubule nucleation. Recent research on the high-resolution structure of γ-TuRC has brought new insight into the microtubule nucleation mechanism. We discuss the effect of γ-TuRC protein overexpression on cancer cell behavior and new drugs directed to γ-tubulin that may offer a viable alternative to microtubule-targeting agents currently used in cancer chemotherapy.

Abstract

In cells, microtubules typically nucleate from microtubule organizing centers, such as centrosomes. γ-Tubulin, which forms multiprotein complexes, is essential for nucleation. The γ-tubulin ring complex (γ-TuRC) is an efficient microtubule nucleator that requires additional centrosomal proteins for its activation and targeting. Evidence suggests that there is a dysfunction of centrosomal microtubule nucleation in cancer cells. Despite decades of molecular analysis of γ-TuRC and its interacting factors, the mechanisms of microtubule nucleation in normal and cancer cells remains obscure. Here, we review recent work on the high-resolution structure of γ-TuRC, which brings new insight into the mechanism of microtubule nucleation. We discuss the effects of γ-TuRC protein dysregulation on cancer cell behavior and new compounds targeting γ-tubulin. Drugs inhibiting γ-TuRC functions could represent an alternative to microtubule targeting agents in cancer chemotherapy.

Insights into interaction mechanism of inhibitors E3T, E3H and E3B with CREB binding protein by using molecular dynamics simulations and MM-GBSA calculations

CREB binding protein (CBP) and its paralog E1A binding protein (p300) are related to the development of inflammatory diseases, cancers and other diseases, and have been potential targets for the treatment of human diseases. In this work, interaction mechanism of three small molecules E3T, E3H, and E3B with CBP was investigated by employing molecular dynamics (MD) simulations, principal component analysis (PCA), and molecular mechanics/generalized born surface area (MM-GBSA) method. The results indicate that inhibitor bindings cause the changes of movement modes and structural flexibility of CBP, and van der Waals interactions mostly drive associations of inhibitors with CBP. In the meantime, the results based on inhibitor-residue interactions not only show that eight residues of CBP can strongly interact with E3T, E3H and E3B but also verify that the CH-CH, CH-π, and π-π interactions are responsible for vital contributions in associations of E3T, E3H and E3B with CBP. In addition, the H-O radial distribution functions (RDFs) were computed to assess the stability of hydrogen bonding interactions between inhibitors and CBP, and the obtained information identifies several key hydrogen bonds playing key roles in bindings of E3T, E3H and E3B to CBP. Potential new inhibitors have been proposed.

Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis

It is well known that bromodomain-containing protein 4 (BRD4) has been thought as a promising target utilized for treating various human diseases, such as inflammatory disorders, malignant tumours, acute myelogenous leukaemia (AML), bone diseases, etc. For this study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were integrated together to uncover binding modes of inhibitors 8P9, 8PU, and 8PX to BRD4(1). The results obtained from binding free energy calculations show that van der Waals interactions act as the main regulator in bindings of inhibitors to BRD4(1). The information stemming from PCA reveals that inhibitor associations extremely affect conformational changes, internal dynamics, and movement patterns of BRD4(1). Residue-based free energy decomposition method was wielded to unveil contributions of independent residues to inhibitor bindings and the data signify that hydrogen bonding interactions and hydrophobic interactions are decisive factors affecting bindings of inhibitors to BRD4(1). Meanwhile, eight residues Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asn140, and Ile146 are recognized as the common hot interaction spots of three inhibitors with BRD4(1). The results from this work are expected to provide a meaningfully theoretical guidance for design and development of effective inhibitors inhibiting of the activity of BRD4.

Molecular mechanism of inhibitor bindings to bromodomain-containing protein 9 explored based on molecular dynamics simulations and calculations of binding free energies

Recently, bromodomain-containing protein 9 (BRD9) has been a prospective therapeutic target for anticancer drug design. Molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were adopted to explore binding modes of three inhibitors (5SW, 5U2, and 5U6) to BRD9 and identify the hot spot of the inhibitor-BRD9 binding. The results indicate that the inhibitor 5SW has the strongest binding ability to BRD9 among the current three inhibitors. Furthermore, the rank of the binding free energies predicted by MM-GBSA approach agrees with that determined by the experimental values. In addition, inhibitor-residue interactions were computed by using residue-based free-energy decomposition method and the results suggest that residue His42 produces the CH-H interactions, residues Asn100, Ile53 and Val49 produce the CH-[Formula: see text] interactions with three inhibitors and Tyr106, Phe45 and Phe44 generate the π-π interactions with inhibitors. Notably, the residue Asn140 forms hydrogen bonding interactions with three inhibitors. This research is expected to provide useful molecular basis and dynamics information at atomic levels for the design of potent inhibitors inhibiting the activity of BRD9.

Combined approach of homology modeling, molecular dynamics, and docking: computer-aided drug discovery

With the continuous development in software, algorithms, and increase in computer speed, the field of computer-aided drug design has been witnessing reduction in the time and cost of the drug designing process. Structure based drug design (SBDD), which is based on the 3D structure of the enzyme, is helping in proposing novel inhibitors. Although a number of crystal structures are available in various repositories, there are various proteins whose experimental crystallization is difficult. In such cases, homology modeling, along with the combined application of MD and docking, helps in establishing a reliable 3D structure that can be used for SBDD. In this review, we have reported recent works, which have employed these three techniques for generating structures and further proposing novel inhibitors, for cytoplasmic proteins, membrane proteins, and metal containing proteins. Also, we have discussed these techniques in brief in terms of the theory involved and the various software employed. Hence, this review can give a brief idea about using these tools specifically for a particular problem.

Noscapine, a Non-addictive Opioid and Microtubule-Inhibitor in Potential Treatment of Glioblastoma

Noscapine is a phthalide isoquinoline alkaloid that easily traverses the blood brain barrier and has been used for years as an antitussive agent with high safety. Despite binding opioid receptors, noscapine lacks significant hypnotic and euphoric effects rendering it safe in terms of addictive potential. In 1954, Hans Lettré first described noscapine as a mitotic poison. The drug was later tested for cancer treatment in the early 1960's, yet no effect was observed likely as a result of its short biological half-life and limited water solubility. Since 1998, it has regained interest thanks to studies from Emory University, which showed its anticancer activity in animal models with negligible toxicity. In contrast to other microtubule-inhibitors, noscapine does not affect the total intracellular tubulin polymer mass. Instead, it forces the microtubules to spend an increased amount of time in a paused state leading to arrest in mitosis and subsequently inducing mitotic slippage/mitotic catastrophe/apoptosis. In experimental models, noscapine does not induce peripheral neuropathy, which is common with other microtubule inhibitors. Noscapine also inhibits tumor growth and enhances cancer chemosensitivity via selective blockage of NF-κB, an important transcription factor in glioblastoma pathogenesis. Due to their anticancer activities and high penetration through the blood-brain barrier, noscapine analogues strongly deserve further study in various animal models of glioblastoma as potential candidates for future patient therapy.

Investigation of the probable homo-dimer model of the Xeroderma pigmentosum complementation group A (XPA) protein to represent the DNA-binding core

The Xeroderma pigmentosum complementation group A (XPA) protein functions as a primary damage verifier and as a scaffold protein in nucleotide excision repair (NER) in all higher organisms. New evidence of XPA's existence as a dimer and the redefinition of its DNA-binding domain (DBD) raises new questions regarding the stability and functional position of XPA in NER. Here, we have investigated XPA's dimeric status with respect to its previously defined DBD (XPA_98-219) as well as with its redefined DBD (XPA_98-239). We studied the stability of XPA_98-210 and XPA_98-239 homo-dimer systems using all-atom molecular dynamics simulation, and we have also characterized the protein-protein interactions (PPI) of these two homo-dimeric forms of XPA. After conducting the root mean square deviation (RMSD) analyses, it was observed that the XPA_98-239 homo-dimer has better stability than XPA_98-210. It was also found that XPA_98-239 has a larger number of hydrogen bonds, salt bridges, and hydrophobic interactions than the XPA_98-210 homo-dimer. We further found that Lys, Glu, Gln, Asn, and Arg residues shared the major contribution toward the intermolecular interactions in XPA homo-dimers. The binding free energy (BFE) analysis, which used the molecular mechanics Poisson-Boltzmann method (MM-PBSA) and the generalized Born and surface area continuum solvation model (GBSA) for both XPA homo-dimers, also substantiated the positive result in favor of the stability of the XPA_98-239 homo-dimer. Communicated by Ramaswamy H. Sarma.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Synthesis and biological evaluation of structurally simplified noscapine analogues as microtubule binding agents

This paper reports on the results of chemical synthesis and biological assays performed on several new analogues of noscapine. We have successfully synthesized four noscapine analogues called 1a–4a, as well as their four corresponding enantiomers called 1b–4b. The chemical pathway consisted of three steps with yields in excess of 60% in each step. Subsequently, we have performed biological activity assays intended to reveal the mode of action of these compounds on microtubules in buffer and in cancer cell lines. We have assayed fluorescence quenching effects in microtubule polymerization experiments, cytotoxicity evaluation in breast cancer cell lines, as well as microtubule dynamicity assessments, for each of the synthesized compounds. Finally, we performed computational docking simulations to two binding sites on β-tubulin: (a) the colchicine binding site and (b) the noscapine binding site. Our results indicate that these compounds have relatively low cytotoxicity profile and less pronounced effects on microtubule dynamics compared with noscapine. Our computational results indicate that these compounds bind to both putative binding sites but have higher affinity for the colchicine site.

Characterizing the Hot Spots Involved in RON-MSPβ Complex Formation Using In Silico Alanine Scanning Mutagenesis and Molecular Dynamics Simulation

Purpose: Implication of protein-protein interactions (PPIs) in development of many diseases such as cancer makes them attractive for therapeutic intervention and rational drug design. RON (Recepteur d'Origine Nantais) tyrosine kinase receptor has gained considerable attention as promising target in cancer therapy. The activation of RON via its ligand, macrophage stimulation protein (MSP) is the most common mechanism of activation for this receptor. The aim of the current study was to perform in silico alanine scanning mutagenesis and to calculate binding energy for prediction of hot spots in protein-protein interface between RON and MSPβ chain (MSPβ). Methods: In this work the residues at the interface of RON-MSPβ complex were mutated to alanine and then molecular dynamics simulation was used to calculate binding free energy. Results: The results revealed that Gln¹⁹³, Arg²²⁰, Glu²⁸⁷, Pro²⁸⁸, Glu²⁸⁹, and His⁴²⁴ residues from RON and Arg⁵²¹, His⁵²⁸, Ser⁵⁶⁵, Glu⁶⁵⁸, and Arg⁶⁸³ from MSPβ may play important roles in protein-protein interaction between RON and MSP. Conclusion: Identification of these RON hot spots is important in designing anti-RON drugs when the aim is to disrupt RON-MSP interaction. In the same way, the acquired information regarding the critical amino acids of MSPβ can be used in the process of rational drug design for developing MSP antagonizing agents, the development of novel MSP mimicking peptides where inhibition of RON activation is required, and the design of experimental site directed mutagenesis studies.

Characterization of Promiscuous Binding of Phosphor Ligands to Breast-Cancer-Gene 1 (BRCA1) C-Terminal (BRCT): Molecular Dynamics, Free Energy, Entropy and Inhibitor Design

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.

Promiscuous proteins are commonly observed in biological systems, such as modular domains that recognize phosphopeptides during signal transduction. The use of phosphopeptides and compounds with phosphate groups as inhibitors to protein–protein interactions have attracted increasing interest for years. By using atomistic molecular dynamics simulations, we are able to perform detailed analyses of the dihedral space to explore protein fluctuation upon ligand binding to better understand promiscuous molecular recognition. Free energy calculation can further provide insights into the mechanism of binding, including both enthalpic and entropic contributions for molecular recognition, which assist in inhibitor design. Our calculation results show that pre-rigidifying a ligand is not always advantageous, suggesting the challenge in retaining optimized intermolecular interactions in pre-rigidified ligand. Instead, certain flexible ligands with multiple binding conformations can reduce entropic penalty, and therefore improves binding affinity. According to our computations, we can introduce new intermolecular interactions to flexible ligand to strengthen attractions while maintaining small entropic penalty by retaining its plasticity in the bound conformation. The study might cast light on a new general strategy for designing inhibitors targeting promiscuous modular domains and protein–protein interactions.