Molecular interactions of type I and type II positive allosteric modulators with the human α7 nicotinic acetylcholine receptor: an in silico study

The binding site locations and structural components for type I and type II positive allosteric modulators (PAMs) of the α7 nicotinic acetylcholine receptor (nAChR) have not been fully characterized yet. In this regard, homology models of the human α7 nAChR and hα7/m5-HT_3A chimera, built using the crystal structure of the serotonin type 3A receptor (5-ΗΤ_3ΑR), were used for molecular docking and molecular dynamics simulations to study the molecular interactions of selected type I (5-hydroxyindol, NS-1738, and LY-2087101) and type II (PNU-120596, PAM-2, and TBS-516) PAMs. The docking results indicate: (1) a site located in the extracellular domain (ECD) for type I PAMs such as NS-1738 and LY-2087101, but not for 5-HI; (2) an overlapping site in the ECD–transmembrane domain (TMD) junction for all studied PAMs. Additional docking results on the hα7/m5-HT_3A chimera supported experimental results indicating that the ECD site might be relevant for type I PAM activity; and (3) two TMD sites, an intrasubunit site that recognizes type II PAMs, and an intersubunit pocket with high specificity for 5-HI (type I PAM). The in silico α7TSLMF mutant results support the view that M1–Ser223 and M3–Ile281 are key residues for the interaction of PAM-2 and PNU-120596 with the intrasubunit cavity. Our in silico results are in agreement with experimental data showing that the intrasubunit cavity is relevant for the activity of type II PAMs, and suggest that the ECD–TMD junction and intersubunit sites could be significant for the activity of type I PAMs.     

Molecular dynamics study of an α-cyclodextrin-phosphatidylinositol inclusion complex

Cyclodextrins are cyclic oligosaccharides known for their ability to include substrate molecules in their hydrophobic cavity. Moreover, cyclodextrins show a hemolytic activity when mm concentrations are added to blood. This hemolysis is commonly interpreted as a massive dissociation of phospholipids from the cell membrane due to the formation of complexes with the cyclodextrins. In the literature, a complexation between -cyclodextrin ( CD) and phosphatidylinositol (PI) specific to the inositol headgroup has been proposed. But the need for the detailed interaction mechanism between the two molecules motivated the present work based on molecular dynamics simulations. Investigation of long range electrostatic interactions shows that a mutual approach of the molecules is only possible when the primary hydroxyl side of CD faces the inositol headgroup of PI. This orientation is also the most favourable from adiabatic- and free-energy profiles calculated along a reaction coordinate that leads to an inclusion of PI into a CD. For free energy simulations, partial hydration of the model has been used. A study of glycosidic bond dihedral angles in CD shows an increase in dihedral fluctuations before complexation and a dihedral freezing once the complex is formed.     

Interaction of fast and slow dynamics in endocrine control systems with an application to β-cell dynamics

Endocrine dynamics spans a wide range of time scales, from rapid responses to physiological challenges to with slow responses that adapt the system to the demands placed on it. We outline a non-linear averaging procedure to extract the slower dynamics in a way that accounts properly for the non-linear dynamics of the faster time scale and is applicable to a hierarchy of more than two time scales, although we restrict our discussion to two scales for the sake of clarity. The procedure is exact if the slow time scale is infinitely slow (the dimensionless ε-quantity is the period of the fast time scale fluctuation times an upper bound to the slow time scale rate of change). However, even for an imperfect separation of time scales we find that this construction provides an excellent approximation for the slow-time dynamics at considerably reduced computational cost. Besides the computation advantage, the averaged equation provided a qualitative insight into the interaction of the time scales. We demonstrate the procedure and its advantages by applying the theory to the model described by Toli? et al. [I.M. Toli?, E. Mosekilde, J. Sturis, Modeling the insulin–glucose feedback system: the significance of pulsatile insulin secretion, J. Theor. Biol. 207 (2000) 361–375.] for ultradian dynamics of the glucose–insulin homeostasis feedback system, extended to include β-cell dynamics. We find that the dynamics of the β-cell mass are dependent not only on the glycemic load (amount of glucose administered to the system), but also on the way this load is applied (i.e. three meals daily versus constant infusion), effects that are lost in the inappropriate methods used by the earlier authors. Furthermore, we find that the loss of the protection against apoptosis conferred by insulin that occurs at elevated levels of insulin has a functional role in keeping the β-cell mass in check without compromising regulatory function. We also find that replenishment of β-cells from a rapidly proliferating pool of cells, as opposed to the slow turn-over which characterises fully differentiated β-cells, is essential to the prevention of type 1 diabetes.     

Molecular dynamics simulation of β₂-microglobulin in denaturing and stabilizing conditions

β₂‐Microglobulin has been a model system for the study of fibril formation for 20 years. The experimental study of β₂‐microglobulin structure, dynamics, and thermodynamics in solution, at atomic detail, along the pathway leading to fibril formation is difficult because the onset of disorder and aggregation prevents signal resolution in Nuclear Magnetic Resonance experiments. Moreover, it is difficult to characterize conformers in exchange equilibrium. To gain insight (at atomic level) on processes for which experimental information is available at molecular or supramolecular level, molecular dynamics simulations have been widely used in the last decade. Here, we use molecular dynamics to address three key aspects of β₂‐microglobulin, which are known to be relevant to amyloid formation: (1) 60 ns molecular dynamics simulations of β₂‐microglobulin in trifluoroethanol and in conditions mimicking low pH are used to study the behavior of the protein in environmental conditions that are able to trigger amyloid formation; (2) adaptive biasing force molecular dynamics simulation is used to force cis‐trans isomerization at Proline 32 and to calculate the relative free energy in the folded and unfolded state. The native‐like trans‐conformer (known as intermediate 2 and determining the slow phase of refolding), is simulated for 10 ns, detailing the possible link between cis‐trans isomerization and conformational disorder; (3) molecular dynamics simulation of highly concentrated doxycycline (a molecule able to suppress fibril formation) in the presence of β₂‐microglobulin provides details of the binding modes of the drug and a rationale for its effect. Proteins 2011. © 2010 Wiley‐Liss, Inc.     

Molecular cloning and characterization of the β-catenin gene from fine-wool sheep

β-Catenin is an evolutionarily conserved molecule that functions as a crucial effector in both cell-to-cell adhesion and Wnt signaling. To gain a better understanding of its role in the development of hair follicles, we cloned the cDNA sequence of the β-catenin gene from the skin of Aohan fine-wool sheep and performed a variety of bioinformatics analyses. We obtained the full-length sequence, which was 4573-bp long and contained a 2346-bp open reading frame encoding a protein of 781 amino acids. The protein had a predicted molecular weight of 85.4 kDa and a theoretical isoelectric point of 5.57. Domain architecture analysis of the β-catenin protein revealed an armadillo repeat region, which is a common feature of β-catenin in other species. The ovine β-catenin gene shares 97.91%, 94.25%, 94.59%, 83.89%, and 89.39% sequence identity with its homologs in Bos taurus, Homo sapiens, Sus scrofa, Gallus gallus, and Mus musculus, respectively, while the amino acid sequence is more than 99% identical with each of these species. The expression of β-catenin mRNA was detected in the heart, liver, spleen, lung, kidney, skin, muscle, and adipose tissue. Expression levels were maximal in the lung and minimal in the muscle, and the difference in expression in these tissues was significant (P < 0.01). Western blot analysis revealed the presence of the β-catenin protein in all tissues examined; expression was lowest in the skin and adipose tissues.     

Loop–loop interaction in an adenine-sensing riboswitch: A molecular dynamics study

Riboswitches are mRNA-based molecules capable of controlling the expression of genes. They undergo conformational changes upon ligand binding, and as a result, they inhibit or promote the expression of the associated gene. The close connection between structural rearrangement and function makes a detailed knowledge of the molecular interactions an important step to understand the riboswitch mechanism and efficiency. We have performed all-atom molecular dynamics simulations of the adenine-sensing add A-riboswitch to study the breaking of the kissing loop, one key tertiary element in the aptamer structure. We investigated the aptamer domain of the add A-riboswitch in complex with its cognate ligand and in the absence of the ligand. The opening of the hairpins was simulated using umbrella sampling using the distance between two loops as the reaction coordinate. A two-step process was observed in all the simulated systems. First, a general loss of stacking and hydrogen bond interactions is seen. The last interactions that break are the two base pairs G37-C61 and G38-C60, but the break does not affect the energy profile, indicating their pivotal role in the tertiary structure formation but not in the structure stabilization. The junction area is partially organized before the kissing loop formation and residue A24 anchors together the loop helices. Moreover, when the distance between the loops is increased, one of the hairpins showed more flexibility by changing its orientation in the structure, while the other conserved its coaxial arrangement with the rest of the structure.     

Molecular docking and dynamics simulations on the interaction of cationic porphyrin–anthraquinone hybrids with DNA G-quadruplexes

A series of cationic porphyrin–anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.     

Distinct interactions between the human adrenergic β2 receptor and Gαs—an in silico study

The aim of this study was to perform an in silico analysis of the interaction of the human β₂ adrenergic receptor with Gα_s. In a first step, a systematic surface-interaction-scan between the inactive or active human β₂ adrenergic receptor and Gα_s was performed in order to gain knowledge about energetically preferred areas on the potential energy surface. Subsequently, two energetically favored regions for the active human β₂ adrenergic receptor–Gα_s complex were identified. Two representative complex structures were put into a POPC (1-palmitoyl-2-oleoyl-phosphatidylcholine) bilayer and solvated in order to perform molecular dynamic simulations. The simulations revealed that both conformations, which have comparable potential energy, are stable. A mean number of about 14 hydrogen bonds was observed between the active receptor and Gα_s for both conformations. Based on these results, two energetically favored β₂–Gα_scomplexes can be proposed.     

Prognostic values of osteopontin-c,E-cadherin and β-catenin in breast cancer

ObjectiveTo determine the correlation of cell adhesion molecules (osteopontin-c, E-cadherin and β-catenin) with clinicopathological characteristics in breast cancer.MethodsImmunostaining of osteopontin-c, E-cadherin and β-catenin were conducted in 170 samples of breast cancer and 30 samples of adjacent normal breast tissues. The correlation of osteopontin-c, E-cadherin and β-catenin expression level with clinicopathological characteristics was evaluated by Pearson's chi-square and Wilcoxon rank-sum test. Univariate and multivariate Cox hazard regression model was used to assess the prognostic values of osteopontin-c, E-cadherin and β-catenin in clinical outcome of breast cancer.ResultsA higher level of osteopontin-c whereas lower levels of E-cadherin and β-catenin were observed in breast cancer as compared with the normal breast tissues. The expression of osteopontin-c was negatively associated with the expression of E-cadherin and β-catenin. The expression of osteopontin-c correlated with lymph node metastasis, and advanced TNM stage and histologic grade. The expression of E-cadherin correlated with low histologic grade; and β-catenin with low TNM stage and histological grade. Moreover, high osteopontin-c level correlated with tumor recurrence or metastasis as well as triple negative subtype. The expression of osteopontin-c was an independent prognostic factor for both disease-free and overall survival of breast cancer patients.ConclusionThe data suggest that the expression of osteopontin-c could serve as a prognostic factor of breast cancer.     

Nonenzymatic hydrolysis reactions of adenosine 5′-triphosphate and related compounds: II. Kinetic studies of the hydrolysis of ATP and related compounds with various polyamines

Various structural features of polyamines which are responsible for the acceleration of the hydrolysis of ATP to ADP and P_i at pH 3–4 were surveyed by means of kinetic studies, leading to the following conclusions: 1) The ethyleneimine chain of the polyamines should be as long as possible; 2) the number of methylene carbon atoms between the two adjacent nitrogen atoms of the polyamine has to be two; 3) the terminal groups of the ethyleneimine chain should be primary amino groups.The rate of ATP hydrolysis in the presence of pentaethylenehexamine (pentaen), which possesses the above properties, was found to be 15 times as high as that of hydrolysis at pH 3.5 in the absence of amines. The kinetic data support the previous assumption that there is formation of an ATP-pentaen complex in the hydrolysis reaction. The formation constant of the complex has been calculated to be K = 1.9 × 10⁴M⁻¹ at pH 3.5 and 50°C from the kinetic data. From the temperature dependence of the rates for pentaen or tetraethylenepentamine, the thermodynamic data of these reactions have been obtained.On the other hand, it has been found that pentaen enhances the hydrolysis of GTP and UTP as well as ATP. No phosphate ester bonds of AMP, p-nitrophenylphosphate and α-d-glucose-1-phosphate were hydrolyzed. Therefore, it may be concluded that hydrolysis of the phosphate ester bond with the polyamine is characteristic of ATP, GTP and UTP.     

The role of encapsulation by β-cyclodextrin in the interaction of raloxifene with macromolecular targets: a study by spectroscopy and molecular modeling

We report the binding of the drug raloxifene with Calf thymus DNA (ctDNA) and bovine serum albumin (BSA) in the presence and absence of β-cyclodextrin (β-CD) and explain the influence of β-cyclodextrin on the binding of the drug to macromolecules. UV-Vis absorption, fluorescence, proton nuclear magnetic resonance and two-dimensional rotating-frame nuclear overhauser effect spectroscopic techniques are used to study the stoichiometry and the binding strength of the complexes. Molecular modeling is used in combination with other techniques to propose the structure of the inclusion complex and the interaction with ctDNA. The Stern–Volmer quenching constants of the interaction of raloxifene with ctDNA in aqueous and in β-CD solution are compared. The competition for binding of ctDNA with raloxifene and Methylene Blue is studied. The apparent binding constant and the number of binding sites for the binding of raloxifene with BSA in aqueous solution are significantly different from those in the presence of β-CD. The influence of β-CD on the binding of the small molecules with biological macromolecules is discussed. We infer that the binding strengths between raloxifene and macromolecules, viz., ctDNA and BSA are influenced by the β-CD encapsulation. These results may suggest new ways to tune the drug binding to biomacromolecules by encapsulating specific moieties of drugs.     

Differential interaction of β2e with phosphoinositides: A comparative study between β2e and MARCKS

Voltage-gated calcium (Ca_V) channels are responsible for Ca²⁺ influx in excitable cells. As one of the auxiliary subunits, the Ca_V β subunit plays a pivotal role in the membrane expression and receptor modulation of Ca_V channels. In particular, the subcellular localization of the β subunit is critical for determining the biophysical properties of Ca_V channels. Recently, we showed that the β2e isotype is tethered to the plasma membrane. Such a feature of β2e is due to the reversible electrostatic interaction with anionic membrane phospholipids. Here, we further explored the membrane interaction property of β2e by comparing it with that of myristoylated alanine-rich C kinase substrate (MARCKS). First, the charge neutralization of the inner leaf of the plasma membrane induced the translocation of both β2e and MARCKS to the cytosol, while the transient depletion of poly-phosphoinositides (poly-PIs) by translocatable pseudojanin (PJ) systems induced the cytosolic translocation of β2e but not MARCKS. Second, the activation of protein kinase C (PKC) induced the translocation of MARCKS but not β2e. We also found that after the cytosolic translocation of MARCKS by receptor activation, depletion of poly-PIs slowed the recovery of MARCKS to the plasma membrane. Together, our data demonstrate that both β2e and MARCKS bind to the membrane through electrostatic interaction but with different binding affinity, and thus, they are differentially regulated by enzymatic degradation of membrane PIs.     

Oncogenic function of DACT1 in colon cancer through the regulation of β-catenin

The Wnt/β-catenin signaling pathway plays important roles in the progression of colon cancer. DACT1 has been identified as a modulator of Wnt signaling through its interaction with Dishevelled (Dvl), a central mediator of both the canonical and noncanonical Wnt pathways. However, the functions of DACT1 in the WNT/β-catenin signaling pathway remain unclear. Here, we present evidence that DACT1 is an important positive regulator in colon cancer through regulating the stability and sublocation of β-catenin. We have shown that DACT1 promotes cancer cell proliferation in vitro and tumor growth in vivo and enhances the migratory and invasive potential of colon cancer cells. Furthermore, the higher expression of DACT1 not only increases the nuclear and cytoplasmic fractions of β-catenin, but also increases its membrane-associated fraction. The overexpression of DACT1 leads to the increased accumulation of nonphosphorylated β-catenin in the cytoplasm and particularly in the nuclei. We have demonstrated that DACT1 interacts with GSK-3β and β-catenin. DACT1 stabilizes β-catenin via DACT1-induced effects on GSK-3β and directly interacts with β-catenin proteins. The level of phosphorylated GSK-3β at Ser9 is significantly increased following the elevated expression of DACT1. DACT1 mediates the subcellular localization of β-catenin via increasing the level of phosphorylated GSK-3β at Ser9 to inhibit the activity of GSK-3β. Taken together, our study identifies DACT1 as an important positive regulator in colon cancer and suggests a potential strategy for the therapeutic control of the β-catenin-dependent pathway.     

Complexation of the mycotoxin zearalenone with β-cyclodextrin: Study of the interaction and first promising applications

This work reports the study of the interactions between native and substituted β-cyclodextrins and zearalenone and its derivatives α- and β-zearelonol. The data obtained by fluorescence and NMR experiments suggested that zearalenone, α- and β-zearalenol and cyclodextrins give rise to host-guest complexation, with the inclusion of the phenolic moiety inside the cyclodextrin cavity. The high stability of these complexes induces a high fluorescence enhancement upon complexation. These results have been successfully applied to the spectrofluorimetric determination of zearalenone in maize raw samples, without any chromatographic separation. Presented at the 29^th Mykotoxin-Workshop, Fellbach, Germany, May 14–16, 2007 Financial support: Emilia-Romagna region (project SIQUAL)     

Structural insights into the Aedes aegypti aquaporins and aquaglyceroporins – an in silico study

There has been a fair bit of understanding on the structure–function relationship of Aquaporins (AQPs) from plants and vertebrates obtained from available X-ray crystallography data. However, there is a lacuna in understanding the structure of AQPs from sanguinivorous insects like the mosquito where it plays a crucial role in survival. In this study, we have built homology models for the Aedes aegypti AQPs, identified key channel lining residues and compared the structure and sequence with orthodox AQPs. Although Ar/R filter residues of AaAQP1 were exactly similar to orthodox AQPs, AaAQP2 has a substitution at LE1position possibly making it less efficient in high capacity water transport. The huge difference in the selectivity filter region of AaAQP3 suggests a different transport property for this channel. The changes observed in the H5 position of the filter of AaAQP4 and AaAQP5 may explain the presence of a larger pore aperture to permit the passage of larger solute molecules. AaAQP6 possesses a completely hydrophobic filter like that in mammalian super aquaporins. The identified key residues are pivotal in understanding the mechanism of action and gating of these channels.     

Molecular dynamics simulations of the effect of the G-protein and diffusible ligands on the β2-adrenergic receptor

G-protein-coupled receptors have extraordinary therapeutic potential as targets for a broad spectrum of diseases. Understanding their function at the molecular level is therefore essential. A variety of crystal structures have made the investigation of the inactive receptor state possible. Recently released X-ray structures of opsin and the β₂-adrenergic receptor (β₂AR) have provided insight into the active receptor state. In addition, we have contributed to the crystal structure of an irreversible agonist-β₂ adrenoceptor complex. These extensive studies and biophysical investigations have revealed that agonist binding leads to a low-affinity conformation of the active state that is suggested to facilitate G-protein binding. The high-affinity receptor state, which promotes signal transduction, is only formed in the presence of both agonist and G-protein. Despite numerous crystal structures, it is not yet clear how ligands tune receptor dynamics and G-protein binding. We have now used molecular dynamics simulations to elucidate the distinct impact of agonist and inverse agonist on receptor conformation and G-protein binding by investigating the influence of the ligands on the structure and dynamics of a complex composed of β₂AR and the C-terminal end of the Gα_s subunit (GαCT). The simulations clearly showed that the agonist isoprenaline and the inverse agonist carazolol influence the ligand-binding site and the interaction between β₂AR and GαCT differently. Isoprenaline induced an inward motion of helix 5, whereas carazolol blocked the rearrangement of the extracellular part of the receptor. Moreover, in the presence of isoprenaline, β₂AR and GαCT form a stable interaction that is destabilized by carazolol.