Molecular modeling study on the dynamical structural features of human smoothened receptor and binding mechanism of antagonist LY2940680 by metadynamics simulation and free energy calculation

Studying noncovalent or covalent bond problem between smoothened and cholesterol by molecular dynamics simulation and Markov state model

Smoothened (SMO) is an attractive therapeutic target for the treatment and prevention of several malignant tumors of the nervous system. The crystal structure of SMO shows cholesterol interacts with residue Asp95 via the noncovalent bond. However, some studies indicate that cholesterol covalently binds to residue Asp95 of SMO. To study these contradictory results, we performed molecular dynamics (MD) simulations and Markov state model (MSM) on SMO in complex with noncovalent-bound and covalent-bound cholesterol. The MD simulated results showed that the noncovalent-bound cholesterol was extremely unstable around the position of residue Asp95 of SMO, while the covalent-bound cholesterol could keep the stable connection with residue Asp95 of SMO. The free energy landscape showed that noncovalent-bound cholesterol had more deep energy wells than covalent-bound cholesterol when it dynamically interacted with the extracellular domain of SMO crystal structure. The MSM results showed the noncovalent-bound cholesterol had more dynamic configuration transformation pathways than the covalent-bound cholesterol. These results theoretically revealed cholesterol should have a covalent bond with residue Asp95 if cholesterol could be stable in the near position of residue Asp95 of SMO. Our studies not only elucidate the covalent binding contradictory issue between cholesterol and residue Asp95 of SMO, but also supply helpful information for antagonists design of SMO.

The lac repressor hinge helix in context: The effect of the DNA binding domain and symmetry

The Lac system of genes has been an important model system in understanding gene regulation. When the dimer lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The hinge region is disordered when binding to nonoperator sequences. This region of the protein must pay a conformational entropic penalty to order when it is bound to operator DNA. Structural studies show that this region is flexible. Previous simulations showed that this region is disordered when free in solution without the DNA binding domain present. Our simulations corroborate that this region is extremely flexible in solution, but we find that the presence of the DNA binding domain proximal to the hinge helix and salt make the ordered conformation more favorable even without DNA present.

Uncovering the Resistance Mechanism of Mycobacterium tuberculosis to Rifampicin Due to RNA Polymerase H451D/Y/R Mutations From Computational Perspective

Tuberculosis is still one of the top 10 causes of deaths worldwide, especially with the emergence of multidrug-resistant tuberculosis. Rifampicin, as the most effective first-line antituberculosis drug, also develops resistance due to the mutation on Mycobacterium tuberculosis (Mtb) RNA polymerase. Among these mutations, three mutations at position 451 (H451D, H451Y, H451R) are associated with high-level resistance to rifampicin. However, the resistance mechanism of Mtb to rifampicin is still unclear. In this work, to explore the resistance mechanism of Mtb to rifampicin due to H451D/Y/R mutations, we combined the molecular dynamics simulation, molecular mechanics generalized-Born surface area calculation, dynamic network analysis, and residue interactions network analysis to compare the interaction change of rifampicin with wild-type RNA polymerase and three mutants. The results of molecular mechanics generalized-Born surface area calculations indicate that the binding free energy of rifampicin with three mutants decreases. In addition, the dynamic network analysis and residue interaction network analysis show that when H451 was mutated, the interactions of residue 451 with its adjacent residues such as Q438, F439, M440, D441, and S447 disappeared or weakened, increasing the flexibility of binding pocket. At the same time, the disappearance of hydrogen bonds between R613 and rifampicin caused by the flipping of R613 is another important reason for the reduction of binding ability of rifampicin in H451D/Y mutants. In H451R mutant, the mutation causes the binding pocket change too much so that the position of rifampicin has a large movement in the binding pocket. In this study, the resistance mechanism of rifampicin at the atomic level is proposed. The proposed drug-resistance mechanism will provide the valuable guidance for the design of antituberculosis drugs.

Discovery of novel natural compound inhibitors targeting estrogen receptor α by an integrated virtual screening strategy

Estrogen receptor (ER) is a nuclear hormone receptor and plays an important role in mediating the cellular effects of estrogen. ER can be classified into two receptors: estrogen receptor alpha (ERα) and beta (ERβ), and the former is expressed in 50~80% of breast tumors and has been extensively investigated in breast cancer for decades. Excessive exposure to estrogen can obviously stimulate the growth of breast cancers primarily mediated by ERα, and thus anti-estrogen therapies by small molecules are of concern to clinicians and pharmaceutical industry in the treatment of ERα-positive breast cancers. Although a series of estrogen receptor modulators have been developed, these drugs can lead to resistance and side effects. Therefore, the development of small molecule inhibitors with high target specificity has been intensified. In this pursuit, an integrated computer-aided virtual screening technique, including molecular docking and pharmacophore model screening, was used to screen traditional Chinese medicine (TCM) databases. The compounds with high docking scores and fit values were subjected to ADME (adsorption, distribution, metabolism, excretion) and toxicity prediction, and ten hits were identified as potential inhibitors targeting ERα. Molecular docking was used to investigate the binding modes between ERα and three most potent hits, and molecular dynamic simulations were chosen to explore the stability of these complexes. The rank of the predicted binding free energies evaluated by MM/GBSA is consistent with the docking score. These novel scaffolds discovered in the present study can be used as critical starting point in the drug discovery process for treating ERα-positive breast cancer. Graphical abstract .

Safety and Tolerability of Sonic Hedgehog Pathway Inhibitors in Cancer

The hedgehog pathway, for which sonic hedgehog (Shh) is the most prominent ligand, is highly conserved and is tightly associated with embryonic development in a number of species. This pathway is also tightly associated with the development of several types of cancer, including basal cell carcinoma (BCC) and acute promyelocytic leukemia, among many others. Inactivating mutations in Patched-1 (PTCH1), leading to ligand-independent pathway activation, are frequent in several cancer types, but most prominent in BCC. This has led to the development of several compounds targeting this pathway as a cancer therapeutic. These compounds target the inducers of this pathway in Smoothened (SMO) and the GLI transcription factors, although targeting SMO has had the most success. Despite the many attempts at targeting this pathway, only three US FDA-approved drugs for cancers affect the Shh pathway. Two of these compounds, vismodegib and sonidegib, target SMO to suppress signaling from either PTCH1 or SMO mutations that lead to upregulation of the pathway. The other approved compound is arsenic trioxide, which can suppress this pathway at the level of the GLI proteins, although current evidence suggests it also has other targets. This review focuses on the safety and tolerability of these clinically approved drugs targeting the Shh pathway, along with a discussion on other Shh pathway inhibitors being developed.

Deciphering the Allosteric Effect of Antagonist Vismodegib on Smoothened Receptor Deactivation Using Metadynamics Simulation

The smoothened receptor (Smo) plays a key role in Hedgehog (Hh) signaling pathway and it has been regarded as an efficacious therapeutic target for basal cell carcinoma (BCC) and medulloblastoma (MB). Nevertheless, the resistance mutation and active mutants of Smo have put forward the requirement of finding more effective inhibitors. Herein, we performed metadynamics simulations on Smo bound with vismodegib (Smo-Vismod) and with cholesterol (Smo-CLR), respectively, to explore the inhibition mechanism of vismodegib. The simulation results indicated that vismodegib-induced shifts of TM5, TM6, and TM7, which permitted the extracellular extension of TM6 and extracellular loop3 (ECL3) to enter the extracellular cysteine-rich domain (CRD) groove. Therefore, an open CRD groove that has not been noticed previously was observed in Smo-Vismod complex. As a consequence, the occupied CRD groove prevents the binding of cholesterol. In addition, the HD and ECLs play crucial roles in the interaction of CRD and TMD. These results reveal that TM5, TM6, and TM7 play important roles in allosteric inhibition the activation of Smo and disrupting cholesterol binding by vismodegib binding. Our results are expected to contribute to understanding the allosteric inhibition mechanism of Smo by vismodegib. Moreover, the detailed conformational changes contribute to the development of novel Smo inhibitors against resistance mutation and active mutants of Smo.

Investigation of ECD conformational transition mechanism of GLP-1R by molecular dynamics simulations and Markov state model

As a member of the class B G protein-coupled receptors (GPCRs), the glucagon-like peptide-1 (GLP-1) can regulate the blood glucose level by binding to the glucagon-like peptide-1 receptor (GLP-1R). Since the extracellular domain (ECD) of GLP-1R is considered as one of the binding sites of GLP-1, the open and closed states of ECD play an important role in the binding process of GLP-1. To investigate the transition path of GLP-1R ECD, the crystal structures of GLP-1R in its bound and unbound states (apo-state) are chosen to perform a total of 1.6 μs of molecular dynamics simulations. The simulated results show that the ECD of GLP-1R closes in the GLP-1 bound state and opens in the GLP-1 unbound state. To determine the critical role that GLP-1 played in regulating the open and closed states of the ECD, we applied the independent gradient model (IGM) to the simulation trajectories. We found that the "hand-like" N-terminal of the GLP-1R ECD plays an important role in the GLP-1 binding. In contrast, the apo-state GLP-1R ECD opens and exposes the two ligand binding domains of GLP-1 after 200 ns of simulations. To elucidate the open and closed mechanisms of GLP-1R ECD in the apo-state and GLP-1 bound state, the Markov state model (MSM) is performed on the MD simulation trajectories. Our results provide possible transition pathways from the closed state to open state of the apo-state GLP-1R ECD. Each pathway contains several intermediate states that correspond to different local minima in deep wells. The dynamical relationships and the most possible conversion pathway between two states are detailed through the MSM analysis. Our results profile the conformation transition mechanism of the GLP-1R ECD and will help in hypoglycemic peptide design of GLP-1R.

Multiscale simulations on human Frizzled and Taste2 GPCRs

Recently, molecular dynamics simulations, from all atom and coarse grained to hybrid methods bridging the two scales, have provided exciting functional insights into class F (Frizzled and Taste2) GPCRs (about 40 members in humans). Findings include: (i) The activation of one member of the Frizzled receptors (FZD4) involves a bending of transmembrane helix TM7 far larger than that in class A GPCRs. (ii) The affinity of an anticancer drug targeting another member (Smoothened receptor) decreases in a specific drug-resistant variant, because the mutation ultimately disrupts the binding cavity and affects TM6. (iii) A novel two-state recognition mechanism explains the very large agonist diversity for at least one member of the Taste2 GPCRs, hTAS2R46.

State of the art of Smo antagonists for cancer therapy: advances in the target receptor and new ligand structures

Since the Hedgehog signaling pathway has been associated with cancer, it has emerged as a therapeutic target for cancer therapy. The main target among the key Hedgehog proteins is the GPCR-like Smo receptor. Therefore, some Smo antagonists that have entered clinical trials, including the US FDA-approved drugs vismodegib and sonidegib, to treat basal cell carcinoma and medulloblastoma. However, early resistance of these drugs has spawned the need to understand the molecular bases of this phenomena. We therefore reviewed details about Smo receptor structures and the best Smo antagonist chemical structures. In addition, we discussed strategies that should be considered to develop new, safer generations of Smo antagonists that avoid current clinical limitations.

Q817G mutation in phosphodiesterase type 5: Conformational analysis and dissociation profile of the inhibitor Tadalafil

Phosphodiesterase type 5 (PDE-5) is an important enzyme involved in the hydrolysis of cyclic guanosine monophosphate (cGMP) to guanosine monophosphate (GMP). The inhibition of this protein leads to the accumulation of cGMP in cells with various biological and therapeutic effects. Several PDE-5 inhibitors exist, with Tadalafil being one of the most commonly studied and used in clinical therapy. In this study, we applied Molecular Dynamics simulations coupled to the ABF (Adaptive Biasing Force) method to study the effect of the mutation on the Gln817 residue (Q817G). The results of the free energy profiles made clear that the affinity of the inhibitor for PDE-5 is dependent on the amino acid residue Gln817. The hydrogen bond made between the side chain of glutamine and the indole ring of Tadalafil results in the stabilization of the ligand in the catalytic site. Despite the prominent role of this interaction, it is important to highlight the contribution of other residues of the catalytic domain for the stabilization of the compound, due to the set of polar, hydrophobic and electrostatic interactions performed by specific amino acid residues.

Price of disorder in the lac repressor hinge helix

The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.

Insight into tetrodotoxin blockade and resistance mechanisms of Nav 1.2 sodium channel by theoretical approaches

Na_v 1.2, a member of voltage-gated sodium channels (Na_v s) that are responsible for the generation and propagation of action potentials along the cell membrane, and play a vital role in the process of information transmission within the nervous system and muscle contraction, is preferentially expressed in the central nervous system. As a potent and selective blocker of Na_v s, tetrodotoxin (TTX) has been extensively studied in biological and chemical sciences, whereas the detailed mechanism by which it blocks nine Na_v 1 channel subtypes remain elusive. Despite the high structural similarity, the TTX metabolite 4,9-anhydro-TTX is 161 times less effective toward the mammalian Na_v 1.2, which puzzled us to ask a question why such a subtle structural variation results in the largely binding affinity difference. In the current work, an integrated computational strategy, including homology modeling, induced fit docking, explicit-solvent MD simulations, and free energy calculations, was employed to investigate the binding mechanism and conformational determinants of TTX analogs. Based on the computational results, the H-bond interactions between C4-OH and C9-OH of TTX and the outer ring carboxylates of the selectivity-filter residues, and the cation-π interaction between the primary amine of guanidinium of TTX and Phe385 determine the difference of their binding affinities. Moreover, the computationally simulations were carried out for the D384N and E945K mutants of hNa_v 1.2-TTX, and the rank of the predicted binding free energies is in accordance with the experimental data. These observations provide a valuable model to design potent and selective neurotoxins of Na_v 1.2 and shed light on the blocking mechanism of TTX to sodium channels.

Importance of protein flexibility on molecular recognition: modeling binding mechanisms of aminopyrazine inhibitors to Nek2

NIMA-related kinase 2 (Nek2) plays a significant role in cell cycle regulation, and overexpression of Nek2 has been observed in several types of carcinoma, suggesting it is a potential target for cancer therapy. In this study, we attempted to gain more insight into the binding mechanisms of a series of aminopyrazine inhibitors of Nek2 through multiple molecular modeling techniques, including molecular docking, molecular dynamics (MD) simulations and free energy calculations. The simulation results showed that the induced fit docking and ensemble docking based on multiple protein structures yield better predictions than conventional rigid receptor docking, highlighting the importance of incorporating receptor flexibility into the accurate predictions of the binding poses and binding affinities of Nek2 inhibitors. Additionally, we observed that the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculations did not show better performance than the docking scoring to rank the binding affinities of the studied inhibitors, suggesting that MM/GBSA is system-dependent and may not be the best choice for the Nek2 systems. Moreover, the detailed information on protein-ligand binding was characterized by the MM/GBSA free energy decomposition, and a number of derivatives with improved docking scores were designed. It is expected that our study can provide valuable information for the future rational design of novel and potent inhibitors of Nek2.

Exploring the binding mechanisms of PDE5 with chromeno[2,3-c]pyrrol-9(2H)-one by theoretical approaches

Cyclic nucleotide phosphodiesterase type 5 (PDE5), exclusively specific for the cyclic guanosine monophosphate (cGMP), is an important drug target for the treatment of erectile dysfunction and pulmonary arterial hypertension (PAH). Although many PDE5 inhibitors have been approved, such as sildenafil, vardenafil, tadalafil and so on, extensive studies have reported some side effects, such as vision disturbance and hearing loss as a result of the amino acid sequence and the secondary structural similarity of other PDEs to the catalytic domain of PDE5. In this study, multiple docking strategies, molecular dynamics (MD) simulations, free energy calculations and decomposition were employed to explore the structural determinants of PDE5 with a series of chromeno[2,3-c]pyrrol-9(2H)-one derivatives. First, reliable docking results were obtained using quantum mechanics/molecular mechanics (QM/MM) docking. Then, MD simulations and MM/GBSA free energy calculations were used to explore the dynamic binding process and characterize the binding modes of the inhibitors with different activities. The predicted binding free energies are in good agreement with the experimental data, and the MM/GBSA free energy decomposition analysis sheds light on the importance of hydrogen bonds with Gln817, π–π stacks against Phe820 and hydrophobic residues for the PDE5 binding of the studied inhibitors. The structural and energetic insights obtained here are useful for understanding the molecular mechanism of ligand binding and designing novel potent and selective PDE5 inhibitors with new scaffolds.

Cyclic nucleotide phosphodiesterase type 5 (PDE5) is exclusively specific for the cyclic guanosine monophosphate (cGMP), and PDE5 is an important drug target for the treatment of erectile dysfunction and pulmonary arterial hypertension (PAH).

Affinity enhancement of nanobody binding to EGFR: in silico site-directed mutagenesis and molecular dynamics simulation approaches

Epidermal growth factor receptor (EGFR), a transmembrane glycoprotein, is overexpressed in many cancers such as head-neck, breast, prostate, and skin cancers for this reason it is a good target in cancer therapy and diagnosis. In nanobody-based cancer diagnosis and treatment, nanobodies with high affinity toward receptor (e.g. EGFR) results in effective treatment or diagnosis of cancer. In this regard, the main aim of this study is to develop a method based on molecular dynamic (MD) simulations for designing of 7D12 based nanobody with high affinity compared with wild-type nanobody. By surveying electrostatic and desolvation interactions between different residues of 7D12 and EGFR, the critical residues of 7D12 that play the main role in the binding of 7D12 to EGFR were elucidated and based on these residues, five logical variants were designed. Following the 50 ns MD simulations, pull and umbrella sampling simulation were performed for 7D12 and all its variants in complex with EGFR. Binding free energy of 7D12 (and all its variants) with EGFR was obtained by weighted histogram analysis method. According to binding free energy results, GLY101 to GLU mutation showed the highest binding affinity but this variant is unstable after 50 ns MD simulations. ALA100 to GLU mutation shows suitable binding enhancement with acceptable structural stability. Suitable position and orientation of GLU in residue 100 of 7D12 against related amino acids of EGFR formed some extra hydrogen and electrostatic interactions which resulted in binding enhancement.

Investigation of allosteric modulation mechanism of metabotropic glutamate receptor 1 by molecular dynamics simulations, free energy and weak interaction analysis

Metabotropic glutamate receptor 1 (mGlu₁), which belongs to class C G protein-coupled receptors (GPCRs), can be coupled with G protein to transfer extracellular signal by dimerization and allosteric regulation. Unraveling the dimer packing and allosteric mechanism can be of great help for understanding specific regulatory mechanism and designing more potential negative allosteric modulator (NAM). Here, we report molecular dynamics simulation studies of the modulation mechanism of FITM on the wild type, T815M and Y805A mutants of mGlu₁ through weak interaction analysis and free energy calculation. The weak interaction analysis demonstrates that van der Waals (vdW) and hydrogen bonding play an important role on the dimer packing between six cholesterol molecules and mGlu₁ as well as the interaction between allosteric sites T815, Y805 and FITM in wild type, T815M and Y805A mutants of mGlu₁. Besides, the results of free energy calculations indicate that secondary binding pocket is mainly formed by the residues Thr748, Cys746, Lys811 and Ser735 except for FITM-bound pocket in crystal structure. Our results can not only reveal the dimer packing and allosteric regulation mechanism, but also can supply useful information for the design of potential NAM of mGlu₁.

Unraveling the conformational determinants of LARP7 and 7SK small nuclear RNA by theoretical approaches

LARP7, a member of the La-related proteins (LARPs), shares a conserved La module comprising the La-motif (LAM) and the RNA-recognition motif (RRM1), binding exclusively to the non-coding RNA 7SK.

Computational study on the interaction between CCR5 and HIV-1 entry inhibitor maraviroc: insight from accelerated molecular dynamics simulation and free energy calculation

C-C chemokine receptor type 5 (CCR5) is the co-receptor of human immunodeficiency virus type 1 (HIV-1) and plays an important role in HIV-1 virus infection. Maraviroc has been proved to be effective for anti-HIV-1 by targeting CCR5. Understanding the detailed interaction mechanism between CCR5 and Maraviroc will be of great help to the rational design of a more potential inverse agonist to block HIV-1 infection. Here, we performed molecular dynamics (MD) simulation and accelerated MD simulation (aMD) to study the interaction mechanism between CCR5 and Maraviroc based on a recently reported crystal structure. The results of MD simulation demonstrate that Maraviroc can form stable hydrogen bonds with residues Tyr37(1.39), Tyr251(6.51) and Glu283(7.39). The results of aMD simulation indicate that the carboxamide moiety is more flexible than the tropane group of Maraviroc in the pocket of CCR5. The electrostatic potential analysis proves that Maraviroc can escape from the pocket of CCR5 along the negative electrostatic potential pathway during the dissociation process. The free energy calculation illustrates that there exist three binding pockets during the dissociation process of Maraviroc. Our results will be useful for understanding the interaction mechanism between CCR5 and Maraviroc as well as for the rational design of a more potent inverse agonist.

Biochemical characterization of smoothened receptor antagonists by binding kinetics against drug-resistant mutant

Hedgehog (Hh) signaling critical for development, differentiation, and cell growth is involved in several cancers, including medulloblastoma and basal cell carcinoma. Although antagonism of the smoothened receptor (SMO), which mediates Hh signaling, is an attractive therapeutic target, a drug-resistant mutation in SMO (SMO-D473H) was identified in a clinical trial of the approved drug vismodegib. TAK-441 potently inhibits SMO-D473H, unlike vismodegib and another SMO antagonist, cyclopamine, whereas the differences in binding modes between these antagonists remain unknown. Here we report the biochemical characterization of TAK-441, vismodegib, and cyclopamine by binding kinetics. The association (kon) and dissociation (koff) rates were determined by kinetic binding studies using [(3)H]TAK-441, and dissociation was confirmed by label-free affinity selection-mass spectrometry (AS-MS). In the [(3)H]TAK-441 competition assay, TAK-441 but not vismodegib and cyclopamine showed time-dependent inhibition. Quantitative kinetic binding analysis revealed that koff of TAK-441 was >10-fold smaller than those of vismodegib and cyclopamine. To further assess the binding mode of antagonists, kinetic binding analysis was performed against SMO-D473H. The D473H mutation affected koff of TAK-441 but not kon. In contrast, only kon was changed by the D473H mutation in the case of vismodegib and cyclopamine. These results suggest that the difference in antagonist efficacy against D473H is associated with the binding mode of antagonists. These findings provide a new insight into the drug action of SMO antagonists and help develop potential therapeutics for drug-resistant mutants.

Ligand induced change of β2 adrenergic receptor from active to inactive conformation and its implication for the closed/open state of the water channel: insight from molecular dynamics simulation, free energy calculation and Markov state model analysis

The reported crystal structures of β2 adrenergic receptor (β2AR) reveal that the open and closed states of the water channel are correlated with the inactive and active conformations of β2AR. However, more details about the process by which the water channel states are affected by the active to inactive conformational change of β2AR remain illusive. In this work, molecular dynamics simulations are performed to study the dynamical inactive and active conformational change of β2AR induced by inverse agonist ICI 118,551. Markov state model analysis and free energy calculation are employed to explore the open and close states of the water channel. The simulation results show that inverse agonist ICI 118,551 can induce water channel opening during the conformational transition of β2AR. Markov state model (MSM) analysis proves that the energy contour can be divided into seven states. States S1, S2 and S5, which represent the active conformation of β2AR, show that the water channel is in the closed state, while states S4 and S6, which correspond to the intermediate state conformation of β2AR, indicate the water channel opens gradually. State S7, which represents the inactive structure of β2AR, corresponds to the full open state of the water channel. The opening mechanism of the water channel is involved in the ligand-induced conformational change of β2AR. These results can provide useful information for understanding the opening mechanism of the water channel and will be useful for the rational design of potent inverse agonists of β2AR.

The Hedgehog pathway: role in cell differentiation, polarity and proliferation

Hedgehog (Hh) is first described as a genetic mutation that has "spiked" phenotype in the cuticles of Drosophila in later 1970s. Since then, Hh signaling has been implicated in regulation of differentiation, proliferation, tissue polarity, stem cell population and carcinogenesis. The first link of Hh signaling to cancer was established through discovery of genetic mutations of Hh receptor gene PTCH1 being responsible for Gorlin syndrome in 1996. It was later shown that Hh signaling is associated with many types of cancer, including skin, leukemia, lung, brain and gastrointestinal cancers. Another important milestone for the Hh research field is the FDA approval for the clinical use of Hh inhibitor Erivedge/Vismodegib for treatment of locally advanced and metastatic basal cell carcinomas. However, recent clinical trials of Hh signaling inhibitors in pancreatic, colon and ovarian cancer all failed, indicating a real need for further understanding of Hh signaling in cancer. In this review, we will summarize recent progress in the Hh signaling mechanism and its role in human cancer.

Search for β2 adrenergic receptor ligands by virtual screening via grid computing and investigation of binding modes by docking and molecular dynamics simulations

We designed a program called MolGridCal that can be used to screen small molecule database in grid computing on basis of JPPF grid environment. Based on MolGridCal program, we proposed an integrated strategy for virtual screening and binding mode investigation by combining molecular docking, molecular dynamics (MD) simulations and free energy calculations. To test the effectiveness of MolGridCal, we screened potential ligands for β2 adrenergic receptor (β2AR) from a database containing 50,000 small molecules. MolGridCal can not only send tasks to the grid server automatically, but also can distribute tasks using the screensaver function. As for the results of virtual screening, the known agonist BI-167107 of β2AR is ranked among the top 2% of the screened candidates, indicating MolGridCal program can give reasonable results. To further study the binding mode and refine the results of MolGridCal, more accurate docking and scoring methods are used to estimate the binding affinity for the top three molecules (agonist BI-167107, neutral antagonist alprenolol and inverse agonist ICI 118,551). The results indicate agonist BI-167107 has the best binding affinity. MD simulation and free energy calculation are employed to investigate the dynamic interaction mechanism between the ligands and β2AR. The results show that the agonist BI-167107 also has the lowest binding free energy. This study can provide a new way to perform virtual screening effectively through integrating molecular docking based on grid computing, MD simulations and free energy calculations. The source codes of MolGridCal are freely available at http://molgridcal.codeplex.com.

Structure and dynamics analysis on plexin-B1 Rho GTPase binding domain as a monomer and dimer

Plexin-B1 is a single-pass transmembrane receptor. Its Rho GTPase binding domain (RBD) can associate with small Rho GTPases and can also self-bind to form a dimer. In total, more than 400 ns of NAMD molecular dynamics simulations were performed on RBD monomer and dimer. Different analysis methods, such as root mean squared fluctuation (RMSF), order parameters (S(2)), dihedral angle correlation, transfer entropy, principal component analysis, and dynamical network analysis, were carried out to characterize the motions seen in the trajectories. RMSF results show that after binding, the L4 loop becomes more rigid, but the L2 loop and a number of residues in other regions become slightly more flexible. Calculating order parameters (S(2)) for CH, NH, and CO bonds on both backbone and side chain shows that the L4 loop becomes essentially rigid after binding, but part of the L1 loop becomes slightly more flexible. Backbone dihedral angle cross-correlation results show that loop regions such as the L1 loop including residues Q25 and G26, the L2 loop including residue R61, and the L4 loop including residues L89-R91, are highly correlated compared to other regions in the monomer form. Analysis of the correlated motions at these residues, such as Q25 and R61, indicate two signal pathways. Transfer entropy calculations on the RBD monomer and dimer forms suggest that the binding process should be driven by the L4 loop and C-terminal. However, after binding, the L4 loop functions as the motion responder. The signal pathways in RBD were predicted based on a dynamical network analysis method using the pathways predicted from the dihedral angle cross-correlation calculations as input. It is found that the shortest pathways predicted from both inputs can overlap, but signal pathway 2 (from F90 to R61) is more dominant and overlaps all of the routes of pathway 1 (from F90 to P111). This project confirms the allosteric mechanism in signal transmission inside the RBD network, which was in part proposed in the previous experimental study.