Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.     

Modeling of the metallo-beta-lactamase from B. fragilis: structural and dynamic effects of inhibitor binding

The structure and dynamics of an inhibitor-bound complex of the metallo-beta-lactamase from Bacteroides fragilis are studied by using molecular dynamics. A search of the conformational space was performed to obtain three distinct models of the complex, which were then subjected to solvated molecular dynamics. A solvated molecular dynamics study of the apo protein was performed to serve as a baseline for comparison with the bound simulations. We find loop conformation changes due to binding as well as a decrease in flexibility of the protein as a whole and especially in the major loop of the beta-lactamase. We report the structural and dynamical features of the inhibitor-bound and apo models, as well as experimentally measurable quantities, which should be capable of distinguishing the two binding modes we have determined.     

Comparison of SARS and NL63 papain-like protease binding sites and binding site dynamics: inhibitor design implications

The human severe acute respiratory syndrome coronavirus (SARS-CoV) and the NL63 coronaviruses are human respiratory pathogens for which no effective antiviral treatment exists. The papain-like cysteine proteases encoded by the coronavirus (SARS-CoV: PLpro; NL63: PLP1 and PLP2) represent potential targets for antiviral drug development. Three recent inhibitor-bound PLpro structures highlight the role of an extremely flexible six-residue loop in inhibitor binding. The high binding site plasticity is a major challenge in computational drug discovery/design efforts. From conventional molecular dynamics and accelerated molecular dynamics (aMD) simulations, we find that with conventional molecular dynamics simulation, PLpro translationally samples the open and closed conformation of BL2 loop on a picosecond-nanosecond timescale but does not reproduce the peptide bond inversion between loop residues Tyr269 and Gln270 that is observed on inhibitor GRL0617 binding. Only aMD simulation, starting from the closed loop conformation, reproduced the 180° ?-ψ dihedral rotation back to the open loop state. The Tyr-Gln peptide bond inversion appears to involve a progressive conformational change of the full loop, starting at one side, and progressing to the other. We used the SARS-CoV apo X-ray structure to develop a model of the NL63-PLP2 catalytic site. Superimposition of the PLP2 model on the PLpro X-ray structure identifies binding site residues in PLP2 that contribute to the distinct substrate cleavage site specificities between the two proteases. The topological and electrostatic differences between the two protease binding sites also help explain the selectivity of non-covalent PLpro inhibitors.     

Effects of ligand binding on the dynamics of rice nonspecific lipid transfer protein 1: a model from molecular simulations

Plant nonspecific lipid transfer proteins (nsLTPs) are small, basic proteins constituted mainly of alpha-helices and stabilized by four conserved disulfide bridges. They are characterized by the presence of a tunnel-like hydrophobic cavity, capable of transferring various lipid molecules between lipid bilayers in vitro. In this study, molecular dynamics (MD) simulations were performed at room temperature to investigate the effects of lipid binding on the dynamic properties of rice nsLTP1. Rice nsLTP1, either in the free form or complexed with one or two lipids was subjected to MD simulations. The C-terminal loop was very flexible both before and after lipid binding, as revealed by calculating the root-mean-square fluctuation. After lipid binding, the flexibility of some residues that were not in direct contact with lipid molecules increased significantly, indicating an increase of entropy in the region distal from the binding site. Essential dynamics analysis revealed clear differences in motion between unliganded and liganded rice nsLTP1s. In the free form of rice nsLTP1, loop1 exhibited the largest directional motion. This specific essential motion mode diminished after binding one or two lipid molecules. To verify the origin of the essential motion observed in the free form of rice nsLTP1, we performed multiple sequence alignments to probe the intrinsic motion encoded in the primary sequence. We found that the amino acid sequence of loop1 is highly conserved among plant nsLTP1s, thus revealing its functional importance during evolution. Furthermore, the sequence of loop1 is composed mainly of amino acids with short side chains. In this study, we show that MD simulations, together with essential dynamics analysis, can be used to determine structural and dynamic differences of rice nsLTP1 upon lipid binding.     

Dynamics of the WPD loop of the Yersinia protein tyrosine phosphatase

The bacterial protein tyrosine phosphatase YopH is an essential virulence determinant in Yersinia spp., causing gastrointestinal diseases and the plague. Like eukaryotic PTPases, YopH catalyzes the hydrolysis of the phosphate moiety of phosphotyrosine within a highly conserved binding pocket, which is also characterized by the closure of the so-called "WPD loop" upon ligand binding. In this study, we investigate the conformational changes and dynamics of the WPD loop by molecular dynamics simulations. Consistent with experimental observations, our simulations show that the WPD loop of YopH is intrinsically flexible and fluctuates between the open and closed conformation with a frequency of approximately 4 ns for the apo, native protein. The region of helix alpha4 spanning loop 384-392, which has been revealed experimentally as a second substrate-binding site in YopH, is found to be highly associated with the WPD loop, stabilizing it in the closed, active conformation, and providing a structural basis for the cooperation of the second-substrate binding site in substrate recognition. Loop L4 (residues 323-327) is shown to be involved in a parallel, correlated motion mode with the WPD loop that contributes the stabilization of a more extended open conformation. In addition, we have simulated the loop reopening in the ligand-bound protein complex by applying the locally enhanced sampling method. Finally, the dynamic behavior of the WPD loop for the C403S mutant differs from the wild-type YopH remarkably. These results shed light on the role of the WPD loop in PTPase-mediated catalysis, and are useful in structure-based design for novel, selective YopH inhibitors as antibacterial drugs.     

Molecular dynamics simulations of the bacterial periplasmic heme binding proteins ShuT and PhuT

ShuT and PhuT are two periplasmic heme binding proteins that shuttle heme between the outer and inner membranes of the Gram-negative bacteria. Periplasmic binding proteins (PBPs) generally exhibit considerable conformational changes during the ligand binding process, whereas ShuT and PhuT belong to a class of PBPs that do not show such behavior based on their apo and holo crystal structures. By employing a series of molecular dynamic simulations on the ShuT and the PhuT, the dynamics and functions of the two PBPs were investigated. Through monitoring the distance changes between the two conserved glutamates of ShuT and PhuT, it was found the two PBPs were more flexible than previously assumed, exhibiting obvious opening-closing motions which were more remarkable in the apo runs of ShuT. Based on the results of the domain motion analysis, large scale conformational transitions were found in all apo runs of ShuT and PhuT, hinting that the domain motions of the two PBPs may be intrinsic. On the basis of the results of the principle component analysis, distinct opening-closing and twisting motion tendencies were observed not only in the apo, but also in the holo simulations of the two PBPs. The Gaussian network model was applied in order to analyze the hinge bending regions. The most important bending regions of ShuT and PhuT are located around the midpoints of their respective connecting helixes. Finally, the flexibilities and the details of the simulations of ShuT and PhuT were discussed. Characterized by the remarkably large flexibilities, the loop constituted by Ala 169, Gly170 and Gly171 of ShuT and the beta-turn constituted by Ala176, Gly177 and Gly178 of PhuT may be important for the functions of the two PBPs. Furthermore, the Asn254 of ShuT and the Arg228 of PhuT may be indispensable for the binding or unbinding of heme, since it is involved in the important hydrogen bonding to the propionate side-chains of heme.     

Dynamics of cholesterol exchange in the oxysterol binding protein family

The oxysterol-binding protein-related protein (ORP) family is essential to sterol transfer and sterol-dependent signal transduction in eukaryotes. The crystal structure of one ORP family member, yeast Osh4, is known in apo and sterol-bound states. In the bound state, a 29 residue N-terminal lid region covers the opening of the cholesterol-binding tunnel, preventing cholesterol exchange. Equilibrium and steered molecular dynamics (MD) simulations of Osh4 were carried out to characterize the mechanism of cholesterol exchange. While most of the structural core was stable during the simulations, the lid was partly opened in the apo equilibrium MD simulation. Helix α7, which undergoes the largest conformational change in the crystallized bound and apo states, is conformationally coupled to the opening of the lid. The movement of α7 helps create a docking site for donor or acceptor membranes in the open state. In the steered MD simulations of cholesterol dissociation, we observed complete opening of the lid covering the cholesterol-binding tunnel. Cholesterol was found to exit the binding pocket in a step-wise process involving (i) the breaking of water-mediated hydrogen bonds and van der Waals contacts within the binding pocket, (ii) opening of the lid covering the binding pocket, and (iii) breakage of transient cholesterol contacts with the rim of the pocket and hydrophobic residues on the interior face of the lid.     

Conformational transitions in protein-protein association: binding of fasciculin-2 to acetylcholinesterase

The neurotoxin fasciculin-2 (FAS2) is a picomolar inhibitor of synaptic acetylcholinesterase (AChE). The dynamics of binding between FAS2 and AChE is influenced by conformational fluctuations both before and after protein encounter. Submicrosecond molecular dynamics trajectories of apo forms of fasciculin, corresponding to different conformational substates, are reported here with reference to the conformational changes of loop I of this three-fingered toxin. This highly flexible loop exhibits an ensemble of conformations within each substate corresponding to its functions. The high energy barrier found between the two major substates leads to transitions that are slow on the timescale of the diffusional encounter of noninteracting FAS2 and AChE. The more stable of the two apo substates may not be the one observed in the complex with AChE. It seems likely that the more stable apo form binds rapidly to AChE and conformational readjustments then occur in the resulting encounter complex.     

ATP hydrolysis at one of the two sites in ABC transporters initiates transport related conformational transitions

ABC transporters play important roles in all types of organisms by participating in physiological and pathological processes. In order to modulate the function of ABC transporters, detailed knowledge regarding their structure and dynamics is necessary. Available structures of ABC proteins indicate three major conformations, a nucleotide-bound "bottom-closed" state with the two nucleotide binding domains (NBDs) tightly closed, and two nucleotide-free conformations, the "bottom-closed" and the "bottom-open", which differ in the extent of separation of the NBDs. However, it remains a question how the widely open conformation should be interpreted, and whether hydrolysis at one of the sites can drive conformational transitions while the NBDs remain in contact. To extend our knowledge, we have investigated the dynamic properties of the Sav1866 transporter using molecular dynamics (MD) simulations. We demonstrate that the replacement of one ATP by ADP alters the correlated motion patterns of the NBDs and the transmembrane domains (TMD). The results suggest that the hydrolysis of a single nucleotide could lead to extracellular closure, driving the transport cycle. Essential dynamics analysis of simulations suggests that single nucleotide hydrolysis can drive the system toward a "bottom-closed" apo conformation similar to that observed in the structure of the MsbA transporter. We also found significant structural instability of the "bottom-open" form of the transporters in simulations. Our results suggest that ATP hydrolysis at one of the sites promotes transport related conformational changes leading to the "bottom-closed" apo conformation, which could thus be physiologically more relevant for describing the structure of the apo state.     

Loop dynamics and ligand binding kinetics in the reaction catalyzed by the Yersinia protein tyrosine phosphatase

The Yersinia protein tyrosine phosphatase (YopH) contains a loop of ten amino acids (the WPD loop) that covers the entrance of the active site of the enzyme during substrate binding. In this work the substrate mimicking competitive inhibitor p-nitrocatechol sulfate (PNC) is used as a probe of the active site. The dynamics of the WPD loop was determined by subjecting an equilibrated system containing YopH, PNC, and YopH bound to PNC to a laser induced temperature jump, and subsequently following the change in equilibrium due to the perturbation. Using this methodology the dynamics associated with substrate binding in YopH have been determined. These results indicate that substrate binding is coupled to the WPD loop motion, and WPD loop dynamics occur in the sub-millisecond time scale. The significance of these dynamic results is interpreted in terms of the catalytic cycle of the enzyme.     

Dynamics of cellular retinoic acid binding protein I on multiple time scales with implications for ligand binding

Cellular retinoic acid binding protein I (CRABPI) belongs to the family of intracellular lipid binding proteins (iLBPs), all of which bind a hydrophobic ligand within an internal cavity. The structures of several iLBPs reveal minimal structural differences between the apo (ligand-free) and holo (ligand-bound) forms, suggesting that dynamics must play an important role in the ligand recognition and binding processes. Here, a variety of nuclear magnetic resonance (NMR) spectroscopy methods were used to systematically study the dynamics of both apo and holo CRABPI at various time scales. Translational and rotational diffusion constant measurements were used to study the overall motions of the proteins. Both apo and holo forms of CRABPI tend to self-associate at high (1.2 mM) concentrations, while at low concentrations (0.2 mM), they are predominantly monomeric. Rapid amide exchange rate and laboratory frame relaxation rate measurements at two spectrometer field strengths (500 and 600 MHz) were used to probe the internal motions of the individual residues. Several residues in the apo form, notably within the ligand recognition region, exhibit millisecond time scale motions that are significantly arrested in the holo form. In contrast, no significant differences in the high-frequency motions were observed between the two forms. These results provide direct experimental evidence for dynamics-induced ligand recognition and binding at a specifically defined time scale. They also exemplify the importance of dynamics in providing a more comprehensive understanding of how a protein functions.     

Effect of metal ion binding on the secondary structure of bovine alpha-lactalbumin as examined by infrared spectroscopy

We have examined the influence of monovalent and divalent cations on the secondary structure of bovine alpha-lactalbumin at neutral pH using Fourier-transform infrared spectroscopy. Our present studies are based on previously reported amide I' component band assignments for this protein [Prestrelski, S. J., Byler, D. M., & Thompson, M. P. (1991) Int. J. Pept. Protein Res. 37, 508-512]. The results indicate that upon dissolution, alpha-lactalbumin undergoes a small, but significant, time-dependent conformational change, regardless of the ions present. Additionally, these studies provide the first quantitative measure of the well-known secondary structural change which accompanies calcium binding. Results indicate that removal of Ca2+ from holo alpha-lactalbumin results in local unfolding of the Ca(2+)-binding loop; the spectra indicate that approximately 16% of the backbone chain changes from a rigid coordination complex to an unordered loop. We have also examined the effects of binding of several other metal ions. Our studies have revealed that binding of Mn2+ to apo alpha-lactalbumin (Ca(2+)-free), while inducing a small, but significant, conformational change, does not cause the alpha-lactalbumin backbone conformation to change to that of the holo (Ca(2+)-bound) form as characterized by infrared spectroscopy. Similar changes to those induced by Mn2+ are observed upon binding of Na+ to apo alpha-lactalbumin, and furthermore, even at very high concentrations (0.2 M), Na+ does not stabilize a structure similar to the holo form. Binding of Zn2+ to the apo form of alpha-lactalbumin does not result in significant backbone conformational changes, suggesting a rigid Zn(2+)-binding site.(ABSTRACT TRUNCATED AT 250 WORDS)     

Conformational dynamics of the membrane enzyme LspA upon antibiotic and substrate binding

Lipoprotein signal peptidase (LspA) is an aspartyl protease that cleaves the transmembrane helix signal peptide of lipoproteins as part of the lipoprotein-processing pathway. Members of this pathway are excellent targets for the development of antibiotic therapeutics because they are essential in Gram-negative bacteria, are important for virulence in Gram-positive bacteria, and may not develop antibiotic resistance. Here, we report the conformational dynamics of LspA in the apo state and bound to the antibiotic globomycin determined using molecular dynamics simulations and electron paramagnetic resonance. The periplasmic helix fluctuates on the nanosecond timescale and samples unique conformations in the different states. In the apo state, the dominant conformation is the most closed and occludes the charged active site from the lipid bilayer. With antibiotic bound there are multiple binding modes with the dominant conformation of the periplasmic helix in a more open conformation. The different conformations observed in both bound and apo states indicate a flexible and adaptable active site, which explains how LspA accommodates and processes such a variety of substrates.     

Dynamic Fluctuations Provide the Basis of a Conformational Switch Mechanism in Apo Cyclic AMP Receptor Protein

Escherichia coli cyclic AMP Receptor Protein (CRP) undergoes conformational changes with cAMP binding and allosterically promotes CRP to bind specifically to the DNA. In that, the structural and dynamic properties of apo CRP prior to cAMP binding are of interest for the comprehension of the activation mechanism. Here, the dynamics of apo CRP monomer/dimer and holo CRP dimer were studied by Molecular Dynamics (MD) simulations and Gaussian Network Model (GNM). The interplay of the inter-domain hinge with the cAMP and DNA binding domains are pre-disposed in the apo state as a conformational switch in the CRP''s allosteric communication mechanism. The hinge at L134-D138 displaying intra- and inter-subunit coupled fluctuations with the cAMP and DNA binding domains leads to the emergence of stronger coupled fluctuations between the two domains and describes an on state. The flexible regions at K52-E58, P154/D155 and I175 maintain the dynamic coupling of the two domains. With a shift in the inter-domain hinge position towards the N terminus, nevertheless, the latter correlations between the domains loosen and become disordered; L134-D138 dynamically interacts only with the cAMP and DNA binding domains of its own subunit, and an off state is assumed. We present a mechanistic view on how the structural dynamic units are hierarchically built for the allosteric functional mechanism; from apo CRP monomer to apo-to-holo CRP dimers.     

A dynamic loop provides dual control over the catalytic and membrane binding activity of a bacterial serine hydrolase

The bacterial acyl protein thioesterase (APT) homologue FTT258 from the gram-negative pathogen Francisella tularensis exists in equilibrium between a closed and open state. Interconversion between these two states is dependent on structural rearrangement of a dynamic loop overlapping its active site. The dynamics and structural properties of this loop provide a simple model for how the catalytic activity of FTT258 could be spatiotemporally regulated within the cell. Herein, we characterized the dual roles of this dynamic loop in controlling its catalytic and membrane binding activity. Using a comprehensive library of loop variants, we determined the relative importance of each residue in the loop to these two biological functions. For the catalytic activity, a centrally located tryptophan residue (Trp66) was essential, with the resulting alanine variant showing complete ablation of enzyme activity. Detailed analysis of Trp66 showed that its hydrophobicity in combination with spatial arrangement defined its essential role in catalysis. Substitution of other loop residues congregated along the N-terminal side of the loop also significantly impacted catalytic activity, indicating a critical role for this loop in controlling catalytic activity. For membrane binding, the centrally located hydrophobic residues played a surprising minor role in membrane binding. Instead general electrostatic interactions regulated membrane binding with positively charged residues bracketing the dynamic loop controlling membrane binding. Overall for FTT258, this dynamic loop dually controlled its biological activities through distinct residues within the loop and this regulation provides a new model for the spatiotemporal control over FTT258 and potentially homologous APT function.     

Dynamic behavior of the post‐SET loop region of NSD1: Implications for histone binding and drug development

NSD1 is a SET‐domain histone methyltransferase that methylates lysine 36 of histone 3. In the crystal structure of NSD1, the post‐SET loop is in an autoinhibitory position that blocks binding of the histone peptide as well as the entrance to the lysine‐binding channel. The conformational dynamics preceding histone binding and the mechanism by which the post‐SET loop moves to accommodate the target lysine is currently unknown, although potential models have been proposed. Using molecular dynamics simulations, we have identified potential conformations of the post‐SET loop differing from those of previous studies, as well as proposed a model of peptide‐bound NSD1. Our simulations illustrate the dynamic behavior of the post‐SET loop and the presence of a few distinct conformations. In every case, the post‐SET loop remains in an autoinhibitory position blocking the peptide‐binding cleft, suggesting that another interaction is required to optimally position NSD1 in an active conformation. This finding provides initial evidence for a mechanism by which NSD1 preferentially binds nucleosomal substrates.     

Dynamic origins of differential RNA binding function in two dsRBDs from the miRNA "microprocessor" complex

MicroRNAs (miRNAs) affect gene regulation by base pairing with mRNA and contribute to the control of cellular homeostasis. The first step in miRNA maturation is conducted in the nucleus by the "microprocessor" complex made up of an RNase III enzyme, Drosha, that contains one dsRNA binding domain (dsRBD), and DGCR8, that contains two dsRBDs in tandem. The crystal structure of DGCR8-Core (493-720), containing both dsRBDs, and the NMR solution structure of Drosha-dsRBD (1259-1337) have been reported, but the solution dynamics have not been explored for any of these dsRBDs. To better define the mechanism of dsRNA binding and thus the nuclear maturation step of miRNA processing, we report NMR spin relaxation and MD simulations of Drosha-dsRBD (1259-1337) and DGCR8-dsRBD1 (505-583). The study was motivated by electrophoretic mobility shift assays (EMSAs) of the two dsRBDs, which showed that Drosha-dsRBD does not bind a representative miRNA but isolated DGCR8-dsRBD1 does (K(d) = 9.4 ± 0.4 μM). Our results show that loop 2 in both dsRBDs is highly dynamic but the pattern of the correlations observed in MD is different for the two proteins. Additionally, the extended loop 1 of Drosha-dsRBD is more flexible than the corresponding loop in DGCR8-dsRBD1 but shows no correlation with loop 2, which potentially explains the lack of dsRNA binding by Drosha-dsRBD in the absence of the RNase III domains. The results presented in this study provide key structural and dynamic features of dsRBDs that contribute to the binding mechanism of these domains to dsRNA.     

Structures of Aplysia AChBP complexes with nicotinic agonists and antagonists reveal distinctive binding interfaces and conformations

Upon ligand binding at the subunit interfaces, the extracellular domain of the nicotinic acetylcholine receptor undergoes conformational changes, and agonist binding allosterically triggers opening of the ion channel. The soluble acetylcholine-binding protein (AChBP) from snail has been shown to be a structural and functional surrogate of the ligand-binding domain (LBD) of the receptor. Yet, individual AChBP species display disparate affinities for nicotinic ligands. The crystal structure of AChBP from Aplysia californica in the apo form reveals a more open loop C and distinctive positions for other surface loops, compared with previous structures. Analysis of Aplysia AChBP complexes with nicotinic ligands shows that loop C, which does not significantly change conformation upon binding of the antagonist, methyllycaconitine, further opens to accommodate the peptidic antagonist, alpha-conotoxin ImI, but wraps around the agonists lobeline and epibatidine. The structures also reveal extended and nonoverlapping interaction surfaces for the two antagonists, outside the binding loci for agonists. This comprehensive set of structures reflects a dynamic template for delineating further conformational changes of the LBD of the nicotinic receptor.     

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand–solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand–DNA binding.     

Fully flexible docking models of the complex between alpha7 nicotinic receptor and a potent heptapeptide inhibitor of the beta-amyloid peptide binding

The heptapeptide IQTTWSR (IQ), recently reported as inhibitor of the beta-amyloid (Abeta) binding to nicotinic acetylcholine receptors (nAChrs), was docked to the homology model of the alpha7 nicotinic acetylcholine receptor. The most representative models were further subjected to molecular dynamics simulations. The data obtained here suggest that Abeta needs highly specific structural motifs to bind to the alpha7nAChR. These structural motifs are located principally in the upper and lower surroundings of loop C, including loop F and sheets beta1, beta2, beta6, beta9, and beta10 of the receptor. Overall, these results suggest that IQ can be mimicked by more bioavailable, stable compounds that would be helpful for the understanding of the Abeta binding site and its dynamics, and for the design of novel agents to be used as an effective alternative against Alzheimer's disease.