Discovery and evaluation of potential sonic hedgehog signaling pathway inhibitors using pharmacophore modeling and molecular dynamics simulations

Sonic hedgehog (Shh) plays an important role in the activation of Shh signaling pathway that regulates preservation and rebirth of adult tissues. An abnormal activation of this pathway has been identified in hyperplasia and various tumorigenesis. Hence the inhibition of this pathway using a Shh inhibitor might be an efficient way to treat a wide range of malignancies. This study was done in order to develop a lead chemical candidate that has an inhibitory function in the Shh signaling pathway. We have generated common feature pharmacophore models using three-dimensional (3D) structural information of robotnikinin, an inhibitor of the Shh signaling pathway, and its analogs. These models have been validated with fit values of robotnikinin and its analogs, and the best model was used as a 3D structural query to screen chemical databases. The hit compounds resulted from the screening docked into a proposed binding site of the Shh named pseudo-active site. Molecular dynamics (MD) simulations were performed to investigate detailed binding modes and molecular interactions between the hit compounds and functional residues of the pseudo-active site. The results of the MD simulation analyses revealed that the hit compounds can bind the pseudo-active site with high affinity than robotnikinin. As a result of this study, a candidate inhibitor (GK 03795) was selected as a potential lead to be employed in future Shh inhibitor design.     

Hedgehog Pathway Antagonist 5E1 Binds Hedgehog at the Pseudo-active Site

Proper hedgehog (Hh) signaling is crucial for embryogenesis and tissue regeneration. Dysregulation of this pathway is associated with several types of cancer. The monoclonal antibody 5E1 is a Hh pathway inhibitor that has been extensively used to elucidate vertebrate Hh biology due to its ability to block binding of the three mammalian Hh homologs to the receptor, Patched1 (Ptc1). Here, we engineered a murine:human chimeric 5E1 (ch5E1) with similar Hh-binding properties to the original murine antibody. Using biochemical, biophysical, and x-ray crystallographic studies, we show that, like the regulatory receptors Cdon and Hedgehog-interacting protein (Hhip), ch5E1 binding to Sonic hedgehog (Shh) is enhanced by calcium ions. In the presence of calcium and zinc ions, the ch5E1 binding affinity increases 10–20-fold to tighter than 1 nm primarily because of a decrease in the dissociation rate. The co-crystal structure of Shh bound to the Fab fragment of ch5E1 reveals that 5E1 binds at the pseudo-active site groove of Shh with an epitope that largely overlaps with the binding site of its natural receptor antagonist Hhip. Unlike Hhip, the side chains of 5E1 do not directly coordinate the Zn²⁺ cation in the pseudo-active site, despite the modest zinc-dependent increase in 5E1 affinity for Shh. Furthermore, to our knowledge, the ch5E1 Fab-Shh complex represents the first structure of an inhibitor antibody bound to a metalloprotease fold.     

Macroautophagy supports Sonic Hedgehog signaling by promoting Patched1 degradation

Autophagy is a highly conservative self-digestion process to maintain intracellular homeostasis and to ensure the survival of cells under stress. Activation of Sonic Hedgehog (Shh) signaling depends on the normal endocytic degradation of pathway receptor Patched1 (Ptch1). It is unclear whether autophagy participates in the receptor endocytosis and modulates Shh signaling transduction. Here we found that blocking macroautophagy attenuates Shh signaling due to the failed transport of Smoothened (Smo) into primary cilia. At the upstream of Smo, Ptch1 was poly-ubiquitinated through K63-conjugated ubiquitin chains. Macroautophagy participates Shh-induced degradation of poly-ubiquitinated Ptch1, contributing to the activation of Shh signaling.     

Discovery of potential sclerostin inhibitors from plants with loop2 region of sclerostin inhibition by interacting with residues outside Pro-Asn-Ala-Ile-Gly motif

Abstract

Sclerostin, an antagonist of the Wnt/β-catenin signaling pathway, was discovered as a potential therapeutic target for stimulating bone formation in osteoporosis. In this study, molecular docking was employed to predict the binding of 29 herbal compounds, which were reported as bone formation stimulators, to the loop2 region of sclerostin. Then, the 50 ns molecular dynamics (MD) simulation of the complexes between sclerostin and the top 10 hits obtained from molecular docking were carried out. Root mean square deviations (RMSDs) analysis of MD trajectories pointed out that all ligands-complexes remain stable throughout the duration of MD simulations. In addition, the molecular mechanics/generalized born surface area (MM/GBSA) binding free energy and energy decomposition analyses were determined. The results here suggested that baicalin is the most promising inhibitor of sclerostin. Interestingly, baicalin binds to sclerostin via the hydrophobic interaction with the amino acid residues on loop2 region but outside the Pro-Asn-Ala-Ile-Gly (PNAIG) motif, particularly the Arg-Gly-Lys-Trp-Trp-Arg (RGKWWR) motif. This finding could be a novel strategy for developing new sclerostin inhibitors in the future.

Communicated by Ramaswamy H. Sarma     

Exploration of Virtual Candidates for Human HMG-CoA Reductase Inhibitors Using Pharmacophore Modeling and Molecular Dynamics Simulations

3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is a rate-controlling enzyme in the mevalonate pathway which involved in biosynthesis of cholesterol and other isoprenoids. This enzyme catalyzes the conversion of HMG-CoA to mevalonate and is regarded as a drug target to treat hypercholesterolemia. In this study, ten qualitative pharmacophore models were generated based on chemical features in active inhibitors of HMGR. The generated models were validated using a test set. In a validation process, the best hypothesis was selected based on the statistical parameters and used for virtual screening of chemical databases to find novel lead candidates. The screened compounds were sorted by applying drug-like properties. The compounds that satisfied all drug-like properties were used for molecular docking study to identify their binding conformations at active site of HMGR. The final hit compounds were selected based on docking score and binding orientation. The HMGR structures in complex with the hit compounds were subjected to 10 ns molecular dynamics simulations to refine the binding orientation as well as to check the stability of the hits. After simulation, binding modes including hydrogen bonding patterns and molecular interactions with the active site residues were analyzed. In conclusion, four hit compounds with new structural scaffold were suggested as novel and potent HMGR inhibitors.     

Molecular insight into pseudolysin inhibition using the MM-PBSA and LIE methods

Pseudolysin, the extracellullar elastase of Pseudomonas aeruginosa (EC: 3.4.24.26) plays an important role in the pathogenesis of P. aeruginosa infections. In the present study, molecular dynamics simulations and theoretical affinity predictions were used to gain molecular insight into pseudolysin inhibition. Four low molecular weight inhibitors were docked at their putative binding sites and molecular dynamics (MD) simulations were performed for 5.0 ns, and the free energy of binding was calculated by the linear interaction energy method. The number and the contact surface area of stabilizing hydrophobic, aromatic, and hydrogen bonding interactions appears to reflect the affinity differences between the inhibitors. The proteinaceous inhibitor, Streptomyces metalloproteinase inhibitor (SMPI) was docked in three different binding positions and MD simulations were performed for 3.0 ns. The MD trajectories were used for molecular mechanics-Poisson-Boltzmann surface area analysis of the three binding positions. Computational alanine scanning of the average pseudolysin-SMPI complexes after MD revealed residues at the pseudolysin-SMPI interface giving the main contribution to the free energy of binding. The calculations indicated that SMPI interacts with pseudolysin via the rigid active site loop, but that also contact sites outside this loop contribute significantly to the free energy of association.     

Cardiac and CNS defects in a mouse with targeted disruption of suppressor of fused

The hedgehog (Hh) pathway is conserved from Drosophila to humans and plays a key role in embryonic development. In addition, activation of the pathway in somatic cells contributes to cancer development in several tissues. Suppressor of fused is a negative regulator of Hh signaling. Targeted disruption of the murine suppressor of fused gene (Sufu) led to a phenotype that included neural tube defects and lethality at mid-gestation (9.0-10.5 dpc). This phenotype resembled that caused by loss of patched (Ptch1), another negative regulator of the Hh pathway. Consistent with this finding, Ptch1 and Sufu mutants displayed excess Hh signaling and resultant altered dorsoventral patterning of the neural tube. Sufu mutants also had abnormal cardiac looping, indicating a defect in the determination of left-right asymmetry. Marked expansion of nodal expression in 7.5 dpc embryos and variable degrees of node dysmorphology in 7.75 dpc embryos suggested that the pathogenesis of the cardiac developmental abnormalities was related to node development. Other mutants of the Hh pathway, such as Shh, Smo and Shh/Ihh compound mutants, also have laterality defects. In contrast to Ptch1 heterozygous mice, Sufu heterozygotes had no developmental defects and no apparent tumor predisposition. The resemblance of Sufu homozygotes to Ptch1 homozygotes is consistent with mouse Sufu being a conserved negative modulator of Hh signaling.     

FGF signaling enhances a sonic hedgehog negative feedback loop at the initiation of spinal cord ventral patterning

A prevalent developmental mechanism for the assignment of cell identities is the production of spatiotemporal concentration gradients of extracellular signaling molecules that are interpreted by the responding cells. One of such signaling systems is the Shh gradient that controls neuronal subtype identity in the ventral spinal cord. Using loss and gain of function approaches in chick and mouse embryos, we show here that the fibroblast growth factor (FGF) signaling pathway is required to restrict the domains of ventral gene expression as neuroepithelial cells become exposed to Shh during caudal extension of the embryo. FGF signaling activates the expression of the Shh receptor and negative pathway regulator Patched 2 (Ptch2) and therefore can enhance a negative feedback loop that restrains the activity of the pathway. Thus, we identify one of the mechanisms by which FGF signaling acts as a modulator of the onset of Shh signaling activity in the context of coordination of ventral patterning and caudal axis extension. © 2015 Wiley Periodicals, Inc. Develop Neurobiol 76: 956–971, 2016     

A molecular dynamics approach to identify an oxysterol-based hedgehog pathway inhibitor

BackgroundMultiple oxysterols (OHCs) have demonstrated the ability to act as agonists or antagonists of the hedgehog (Hh) signaling pathway, a developmental signaling pathway that has been implicated as a potential therapeutic target in a variety of human diseases. These OHCs are known to modulate Hh signaling through direct binding interactions with the N-terminal cysteine rich domain (CRD) of Smoothened, a key regulator of Hh signal transduction.MethodsHomology modeling, molecular dynamics simulations, and MM/GBSA energy calculations were utilized to explore binding interactions between the OHC scaffold and the human Smoothened CRD. Follow-up cellular assays explored the in vitro activity of potential Hh pathway modulators.ResultsStructural features that govern key molecular interactions between the Smoothened CRD and the OHC scaffold were identified. Orientation of the iso-octyl side chain as well as the overall entropy of the OHC-CRD complex are important for determining activity against the Hh pathway. OHC 9, which was previously thought to be inactive because it was not an Hh agonist, was identified as an inhibitor of Hh signal transmission.ConclusionsCalculated MM/GBSA binding energies for OHCs in complex with the CRD of Smoothened correlate well with in vitro Hh modulatory activity. Compounds with high affinity stabilize Smoothened and are antagonists, whereas compounds with reduced affinity allow a conformational change in Smoothened that results in pathway activation.General significanceComputational modeling and molecular dynamics simulations can be used to predict whether a small molecule that binds the Smoothened CRD will be an agonist or antagonist of the pathway.     

Shh signaling guides spatial pathfinding of raphespinal tract axons by multidirectional repulsion

Relatively little is known about the molecular mechanisms underlying spatial pathfinding in the descending serotonergic raphespinal tract (RST) in the developing spinal cord, one of the most important nerve pathways for pain, sensory and motor functions. We provide evidence that ventral floor plate-secreted Sonic hedgehog (Shh) is responsible for the establishment of decreasing gradients in both the anterior-to-posterior (A-P) and the medial-to-lateral (M-L) directions in the ventral spinal cord during serotonergic RST axon projection. Downstream components of the Shh pathway, Patched 1 (Ptch1) and Smoothened (Smo), were expressed in the serotonergic caudal raphe nuclei and enriched in the descending serotonergic RST axons. Diffusible Shh repulsion of serotonergic RST axons was shown to be mediated by Shh-Ptch1 interactions and derepression of Smo. Using a co-culture assay, we showed that A-P graded repulsion mediated by Shh signaling pushed the serotonergic axons caudally through the ventral spinal cord and M-L graded repulsion mediated by Shh signaling simultaneously restricted the serotonergic axons to the ventral and ventral-lateral funiculus. Prominent pathfinding errors of serotonergic RST axons were observed in various Shh, Ptch1 and Smo mutants. We conclude that Shh signaling-mediated multidirectional repulsion is required to push descending serotonergic RST axons in the A-P direction, and to restrict these axons to the ventral and ventral-lateral funiculus in the M-L direction. This is the first demonstration that Shh signaling-mediated multidirectional repulsion of serotonergic RST axons maintains spatial axon pathfinding in the developing spinal cord.     

In silico docking and molecular dynamics simulation of 3-dehydroquinate synthase (DHQS) from <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis>

The shikimate pathway is as an attractive target because it is present in bacteria, algae, fungi, and plants but does not occur in mammals. In Mycobacterium tuberculosis (MTB), the shikimate pathway is integral to the biosynthesis of naphthoquinones, menaquinones, and mycobactin. In these study, novel inhibitors of 3-dehydroquinate synthase (DHQS), an enzyme that catalyzes the second step of the shikimate pathway in MTB, were determined. 12,165 compounds were selected from two public databases through virtual screening and molecular docking analysis using PyRx 8.0 and Autodock 4.2, respectively. A total of 18 compounds with the best binding energies (?13.23 to ?8.22 kcal/mol) were then selected and screened for absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis, and nine of those compounds were found to satisfy all of the ADME and toxicity criteria. Among those nine, the three compounds—ZINC633887?(binding energy =??10.29 kcal/mol), ZINC08983432?(?9.34 kcal/mol), and PubChem73393?(?8.61 kcal/mol)—with the best binding energies were further selected for molecular dynamics (MD) simulation analysis. The results of the 50-ns MD simulations showed that the two compounds ZINC633887 and PubChem73393 formed stable complexes with DHQS and that the structures of those two ligands remained largely unchanged at the ligand-binding site during the simulations. These two compounds identified through docking and MD simulation are potential candidates for the treatment of TB, and should undergo validation in vivo and in vitro.     

Novel virtual lead identification in the discovery of hematopoietic cell kinase (HCK) inhibitors: application of 3D QSAR and molecular dynamics simulation

High level of hematopoietic cell kinase (Hck) is associated with drug resistance in chronic myeloid leukemia. Additionally, Hck activity has also been connected with the pathogenesis of HIV-1 and chronic obstructive pulmonary disease. In this study, three-dimensional (3D) QSAR pharmacophore models were generated for Hck based on experimentally known inhibitors. A best pharmacophore model, Hypo1, was developed with high correlation coefficient (0.975), Low RMS deviation (0.60) and large cost difference (49.31), containing three ring aromatic and one hydrophobic aliphatic feature. It was further validated by the test set (r?=?0.96) and Fisher’s randomization method (95%). Hypo 1 was used as a 3D query for screening the chemical databases, and the hits were further screened by applying Lipinski’s rule of five and ADMET properties. Selected hit compounds were subjected to molecular docking to identify binding conformations in the active site. Finally, the appropriate binding modes of final hit compounds were revealed by molecular dynamics (MD) simulations and free energy calculation studies. Hence, we propose the final three hit compounds as virtual candidates for Hck inhibitors.     

Targeting the NTPase site of Zika virus NS3 helicase for inhibitor discovery

Abstract

The major threats linked to Zika virus (ZIKV) are microcephaly, Guillain-Barre syndrome, and the ability to transfer through sexual transmission. Despite these threats, Zika specific FDA approved drugs or vaccines are not available as of yet. Additionally, the involvement of pregnant women makes the drug screening process lengthy and complicated in terms of safety and minimum toxicity of the molecules. Since NS3 helicase of ZIKV performs the critical function of unwinding double-stranded RNA during replication, it is considered as a promising drug target to block ZIKV replication. In the present study, we have exploited the NTPase site of ZIKV NS3 helicase for screening potential inhibitor compounds by molecular docking, and molecular dynamics (MD) simulation approaches. NS3 helicase hydrolyzes the ATP to use its energy for unwinding RNA. We have chosen twenty natural compounds from ZINC library with known antiviral properties and a helicase focused library (HFL) of small molecules from Life Chemicals compounds. After going through docking, the top hit molecules from ZINC and HFL library were further analysed by MD simulations to find out stable binding poses. Finally, we have reported the molecules with potential of binding at NTPase pocket of ZIKV NS3 helicase, which could be further tested on virus through in vitro experiments to check their efficacy.

Communicated by Ramaswamy H. Sarma     

Distinct Effects of the mesenchymal dysplasia Gene Variant of Murine Patched-1 Protein on Canonical and Non-canonical Hedgehog Signaling Pathways

Hedgehog (Hh) signaling requires regulation of the receptor Patched-1 (Ptch1), which, in turn, regulates Smoothened activity (canonical Hh signaling) as well as other non-canonical signaling pathways. The mutant Ptch1 allele mesenchymal dysplasia (mes), which truncates the Ptch1 C terminus, produces a limited spectrum of developmental defects in mice as well as deregulation of canonical Hh signaling in some, but not all, affected tissues. Paradoxically, mes suppresses canonical Hh signaling and binds to Hh ligands with an affinity similar to wild-type mouse Ptch1 (mPtch1). We characterized the distinct activities of the mes variant of mPtch1 mediating Hh signaling through both canonical and non-canonical pathways. We demonstrated that mPtch1 bound c-src in an Hh-regulated manner. Stimulation with Sonic Hedgehog (Shh) of primary mammary mesenchymal cells from wild-type and mes animals activated Erk1/2. Although Shh activated c-src in wild-type cells, c-src was constitutively activated in mes mesenchymal cells. Transient assays showed that wild-type mPtch1, mes, or mPtch1 lacking the C terminus repressed Hh signaling in Ptch1-deficient mouse embryo fibroblasts and that repression was reversed by Shh, revealing that the C terminus was dispensable for mPtch1-dependent regulation of canonical Hh signaling. In contrast to these transient assays, constitutively high levels of mGli1 but not mPtch1 were present in primary mammary mesenchymal cells from mes mice, whereas the expression of mPtch1 was similarly induced in both mes and wild-type cells. These data define a novel signal transduction pathway involving c-src that is activated by the Hh ligands and reveals the requirement for the C terminus of Ptch in regulation of canonical and non-canonical Hh signaling pathways.     

Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models?

The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D₃ dopamine receptor (D₃R) in complex with an antagonist. Close to 60 μs of simulation time was generated and the resulting MD refined models were compared to a D₃R crystal structure. In the MD simulations, the receptor models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.     

The p53 inhibitor MDM2 facilitates Sonic Hedgehog-mediated tumorigenesis and influences cerebellar foliation

Disruption of cerebellar granular neuronal precursor (GNP) maturation can result in defects in motor coordination and learning, or in medulloblastoma, the most common childhood brain tumor. The Sonic Hedgehog (Shh) pathway is important for GNP proliferation; however, the factors regulating the extent and timing of GNP proliferation, as well as GNP differentiation and migration are poorly understood. The p53 tumor suppressor has been shown to negatively regulate the activity of the Shh effector, Gli1, in neural stem cells; however, the contribution of p53 to the regulation of Shh signaling in GNPs during cerebellar development has not been determined. Here, we exploited a hypomorphic allele of Mdm2 (Mdm2(puro)), which encodes a critical negative regulator of p53, to alter the level of wild-type MDM2 and p53 in vivo. We report that mice with reduced levels of MDM2 and increased levels of p53 have small cerebella with shortened folia, reminiscent of deficient Shh signaling. Indeed, Shh signaling in Mdm2-deficient GNPs is attenuated, concomitant with decreased expression of the Shh transducers, Gli1 and Gli2. We also find that Shh stimulation of GNPs promotes MDM2 accumulation and enhances phosphorylation at serine 166, a modification known to increase MDM2-p53 binding. Significantly, loss of MDM2 in Ptch1(+/-) mice, a model for Shh-mediated human medulloblastoma, impedes cerebellar tumorigenesis. Together, these results place MDM2 at a major nexus between the p53 and Shh signaling pathways in GNPs, with key roles in cerebellar development, GNP survival, cerebellar foliation, and MB tumorigenesis.     

Syntheses of aminoalcohol-derived macrocycles leading to a small-molecule binder to and inhibitor of Sonic Hedgehog

We report the synthesis and biological activity of a library of aminoalcohol-derived macrocycles from which robotnikinin (17), a binder to and inhibitor of Sonic Hedgehog, was derived. Using an asymmetric alkylation to set a key stereocenter and an RCM reaction to close the macrocycle, we were able to synthesize compounds for testing. High-throughput screening via small-molecule microarray (SMM) technology led to the discovery of a compound capable of binding ShhN. Follow-up chemistry led to a library of macrocycles with enhanced biological activity relative to the original hit compounds. Differences in ring size and stereochemistry, leading to alterations in the mode of binding, may account for differences in the degree of biological activity. These compounds are the first ones reported that inhibit Shh signaling at the ShhN level.     

Combination of Cyclopamine and Tamoxifen Promotes Survival and Migration of MCF-7 Breast Cancer Cells – Interaction of Hedgehog-Gli and Estrogen Receptor Signaling Pathways

Hedgehog-Gli (Hh-Gli) signaling pathway is one of the new molecular targets found upregulated in breast tumors. Estrogen receptor alpha (ERα) signaling has a key role in the development of hormone-dependent breast cancer. We aimed to investigate the effects of inhibiting both pathways simultaneously on breast cancer cell survival and the potential interactions between these two signaling pathways. ER-positive MCF-7 cells show decreased viability after treatment with cyclopamine, a Hh-Gli pathway inhibitor, as well as after tamoxifen (an ERα inhibitor) treatment. Simultaneous treatment with cyclopamine and tamoxifen on the other hand, causes short-term survival of cells, and increased migration. We found upregulated Hh-Gli signaling under these conditions and protein profiling revealed increased expression of proteins involved in cell proliferation and migration. Therefore, even though Hh-Gli signaling seems to be a good potential target for breast cancer therapy, caution must be advised, especially when combining therapies. In addition, we also show a potential direct interaction between the Shh protein and ERα in MCF-7 cells. Our data suggest that the Shh protein is able to activate ERα independently of the canonical Hh-Gli signaling pathway. Therefore, this may present an additional boost for ER-positive cells that express Shh, even in the absence of estrogen.     

Biophysical and computational studies of membrane penetration by the GRP1 pleckstrin homology domain

The pleckstrin homology (PH) domain of the general receptor for phosphoinositides 1 (GRP1) exhibits specific, high-affinity, reversible binding to phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P(3)) at?the plasma membrane, but the nature and extent of the interaction between this bound complex and the surrounding membrane environment remains unclear. Combining equilibrium and nonequilibrium molecular dynamics (MD) simulations, NMR spectroscopy, and monolayer penetration experiments, we characterize the membrane-associated state of?GRP1-PH. MD simulations show loops flanking the binding site supplement the interaction with PI(3,4,5)P(3) through multiple contacts with the lipid bilayer. NMR data show large perturbations in chemical shift for these loop regions on binding to PI(3,4,5)P(3)-containing DPC micelles. Monolayer penetration experiments and further MD simulations demonstrate that mutating hydrophobic residues to polar residues in the flanking loops reduces membrane penetration. This supports a "dual-recognition" model of binding, with specific GRP1-PH-PI(3,4,5)P(3) interactions supplemented by interactions of loop regions with the lipid bilayer.