Metal Ions Can Modulate the Self-Assembly and Activity of Catalytic Peptide Amyloids.

Rational design of peptides has become a powerful tool to produce self-assembled nanostructures with the ability to catalyze different chemical reactions, paving the way to develop minimalistic enzyme-like nanomaterials. Catalytic amyloid-like assemblies have emerged among the most versatile and active, but they often require additional factors for activity. Elucidating how these factors influence the structure and activity is key for the design. Here, we showed that biologically relevant metal ions can guide and modulate the self-assembly of a small peptide into diverse amyloid architectures. The morphology and catalytic activity of the resulting fibrils were tuned by the specific metal ion decorating the surface, whereas X-ray structural analysis of the amyloids showed ion-dependent shape sizes. Molecular dynamics simulations showed that the metals can strongly affect the local conformational space, which can trigger major rearrangements of the fibrils. Our results demonstrate that the conformational landscape of catalytic amyloids is broad and tunable by external factors, which can be critical for future design strategies.

The prolactin receptor scaffolds Janus kinase 2 via co-structure formation with phosphoinositide-4,5-bisphosphate.

Class 1 cytokine receptors transmit signals through the membrane by a single transmembrane helix to an intrinsically disordered cytoplasmic domain that lacks kinase activity. While specific binding to phosphoinositides has been reported for the prolactin receptor (PRLR), the role of lipids in PRLR signaling is unclear. Using an integrative approach combining nuclear magnetic resonance spectroscopy, cellular signaling experiments, computational modeling, and simulation, we demonstrate co-structure formation of the disordered intracellular domain of the human PRLR, the membrane constituent phosphoinositide-4,5-bisphosphate (PI(4,5)P2) and the FERM-SH2 domain of the Janus kinase 2 (JAK2). We find that the complex leads to accumulation of PI(4,5)P2 at the transmembrane helix interface and that the mutation of residues identified to interact specifically with PI(4,5)P2 negatively affects PRLR-mediated activation of signal transducer and activator of transcription 5 (STAT5). Facilitated by co-structure formation, the membrane-proximal disordered region arranges into an extended structure. We suggest that the co-structure formed between PRLR, JAK2, and PI(4,5)P2 locks the juxtamembrane disordered domain of the PRLR in an extended structure, enabling signal relay from the extracellular to the intracellular domain upon ligand binding. We find that the co-structure exists in different states which we speculate could be relevant for turning signaling on and off. Similar co-structures may be relevant for other non-receptor tyrosine kinases and their receptors.

Structural model of tissue factor (TF) and TF-factor VIIa complex in a lipid membrane: A combined experimental and computational study.

Tissue factor (TF) is a membrane protein involved in blood coagulation. TF initiates a cascade of proteolytic reactions, ultimately leading to the formation of a blood clot. The first reaction consists of the binding of the coagulation factor VII and its conversion to the activated form, FVIIa. Here, we combined experimental, i.e. quartz crystal microbalance with dissipation monitoring and neutron reflectometry, and computational, i.e. molecular dynamics (MD) simulation, methods to derive a complete structural model of TF and TF/FVIIa complex in a lipid bilayer. This model shows that the TF transmembrane domain (TMD), and the flexible linker connecting the TMD to the extracellular domain (ECD), define the location of the ECD on the membrane surface. The average orientation of the ECD relative to the bilayer surface is slightly tilted towards the lipid headgroups, a conformation that we suggest is promoted by phosphatidylserine lipids, and favours the binding of FVIIa. On the other hand, the formation of the TF/FVIIa complex induces minor changes in the TF structure, and reduces the conformational freedom of both TF and FVIIA. Altogether we describe the protein-protein and protein-lipid interactions favouring blood coagulation, but also instrumental to the development of new drugs.

Collective mechanical responses of cadherin-based adhesive junctions as predicted by simulations.

Cadherin-based adherens junctions and desmosomes help stabilize cell-cell contacts with additional function in mechano-signaling, while clustered protocadherin junctions are responsible for directing neuronal circuits assembly. Structural models for adherens junctions formed by epithelial cadherin (CDH1) proteins indicate that their long, curved ectodomains arrange to form a periodic, two-dimensional lattice stabilized by tip-to-tip trans interactions (across junction) and lateral cis contacts. Less is known about the exact architecture of desmosomes, but desmoglein (DSG) and desmocollin (DSC) cadherin proteins are also thought to form ordered junctions. In contrast, clustered protocadherin (PCDH)-based cell-cell contacts in neuronal tissues are thought to be responsible for self-recognition and avoidance, and structural models for clustered PCDH junctions show a linear arrangement in which their long and straight ectodomains form antiparallel overlapped trans complexes. Here, we report all-atom molecular dynamics simulations testing the mechanics of minimalistic adhesive junctions formed by CDH1, DSG2 coupled to DSC1, and PCDHγB4, with systems encompassing up to 3.7 million atoms. Simulations generally predict a favored shearing pathway for the adherens junction model and a two-phased elastic response to tensile forces for the adhesive adherens junction and the desmosome models. Complexes within these junctions first unbend at low tensile force and then become stiff to unbind without unfolding. However, cis interactions in both the CDH1 and DSG2-DSC1 systems dictate varied mechanical responses of individual dimers within the junctions. Conversely, the clustered protocadherin PCDHγB4 junction lacks a distinct two-phased elastic response. Instead, applied tensile force strains trans interactions directly, as there is little unbending of monomers within the junction. Transient intermediates, influenced by new cis interactions, are observed after the main rupture event. We suggest that these collective, complex mechanical responses mediated by cis contacts facilitate distinct functions in robust cell-cell adhesion for classical cadherins and in self-avoidance signaling for clustered PCDHs.

Elastic versus brittle mechanical responses predicted for dimeric cadherin complexes.

Cadherins are a superfamily of adhesion proteins involved in a variety of biological processes that include the formation of intercellular contacts, the maintenance of tissue integrity, and the development of neuronal circuits. These transmembrane proteins are characterized by ectodomains composed of a variable number of extracellular cadherin (EC) repeats that are similar but not identical in sequence and fold. E-cadherin, along with desmoglein and desmocollin proteins, are three classical-type cadherins that have slightly curved ectodomains and engage in homophilic and heterophilic interactions through an exchange of conserved tryptophan residues in their N-terminal EC1 repeat. In contrast, clustered protocadherins are straighter than classical cadherins and interact through an antiparallel homophilic binding interface that involves overlapped EC1 to EC4 repeats. Here we present molecular dynamics simulations that model the adhesive domains of these cadherins using available crystal structures, with systems encompassing up to 2.8 million atoms. Simulations of complete classical cadherin ectodomain dimers predict a two-phased elastic response to force in which these complexes first softly unbend and then stiffen to unbind without unfolding. Simulated α, ß, and γ clustered protocadherin homodimers lack a two-phased elastic response, are brittle and stiffer than classical cadherins and exhibit complex unbinding pathways that in some cases involve transient intermediates. We propose that these distinct mechanical responses are important for function, with classical cadherin ectodomains acting as molecular shock absorbers and with stiffer clustered protocadherin ectodomains facilitating overlap that favors binding specificity over mechanical resilience. Overall, our simulations provide insights into the molecular mechanics of single cadherin dimers relevant in the formation of cellular junctions essential for tissue function.

Order and disorder-An integrative structure of the full-length human growth hormone receptor.

Because of its small size (70 kilodalton) and large content of structural disorder (>50%), the human growth hormone receptor (hGHR) falls between the cracks of conventional high-resolution structural biology methods. Here, we study the structure of the full-length hGHR in nanodiscs with small-angle x-ray scattering (SAXS) as the foundation. We develop an approach that combines SAXS, x-ray diffraction, and NMR spectroscopy data obtained on individual domains and integrate these through molecular dynamics simulations to interpret SAXS data on the full-length hGHR in nanodiscs. The hGHR domains reorient freely, resulting in a broad structural ensemble, emphasizing the need to take an ensemble view on signaling of relevance to disease states. The structure provides the first experimental model of any full-length cytokine receptor in a lipid membrane and exemplifies how integrating experimental data from several techniques computationally may access structures of membrane proteins with long, disordered regions, a widespread phenomenon in biology.

Functional characterization of the ATPase-like activity displayed by a catalytic amyloid.

BACKGROUND: Amyloids are highly ordered polypeptide aggregates stabilized by a beta-sheet structural core. Though classically associated to pathology, reports on novel functional roles of these proteins have increasingly emerged in the past decade. Moreover, the recent discovery that amyloids formed with rationally designed small peptides can exhibit catalytic reactivity has opened up new opportunities in both biology and biotechnology. The observed activities typically require the binding of divalent metals, giving rise to active metal-amyloid complexes. METHODS: Peptide (SDIDVFI) was aggregated in vitro. The structure of the self-assembled species was analyzed using fluorescence, transmission electron microscopy, circular dichroism and computational modeling. A kinetic characterization of the emerging catalytic activity was performed. RESULTS: The peptide self-assembled into canonical amyloids that exhibited catalytic activity towards hydrolysis of the phosphoanhydride bonds of adenosine triphosphate (ATP), partially mimicking an ATPase-like enzyme. Both amyloid formation and activity are shown to depend on manganese (Mn2+) binding. The activity was not restricted to ATP but also affected all other ribonucleotides (GTP, CTP and UTP). Peptides carrying a single aspartate exhibited a similar activity. CONCLUSIONS: The phosphoanhydride bonds appear as the main specificity target of the Mn2+-amyloid complex. A single aspartate per peptide is sufficient to enable the hydrolytic activity. GENERAL SIGNIFICANCE: Catalytic amyloids are shown for the first time to catalyze the hydrolysis of all four ribonucleotides. Our results should contribute towards understanding the biological implications of amyloid-mediated reactivity as well as in the design of future catalytic amyloids for biotechnological applications.

The intracellular lipid-binding domain of human Na+/H+ exchanger 1 forms a lipid-protein co-structure essential for activity.

Dynamic interactions of proteins with lipid membranes are essential regulatory events in biology, but remain rudimentarily understood and particularly overlooked in membrane proteins. The ubiquitously expressed membrane protein Na+/H+-exchanger 1 (NHE1) regulates intracellular pH (pHi) with dysregulation linked to e.g. cancer and cardiovascular diseases. NHE1 has a long, regulatory cytosolic domain carrying a membrane-proximal region described as a lipid-interacting domain (LID), yet, the LID structure and underlying molecular mechanisms are unknown. Here we decompose these, combining structural and biophysical methods, molecular dynamics simulations, cellular biotinylation- and immunofluorescence analysis and exchanger activity assays. We find that the NHE1-LID is intrinsically disordered and, in presence of membrane mimetics, forms a helical αα-hairpin co-structure with the membrane, anchoring the regulatory domain vis-a-vis the transport domain. This co-structure is fundamental for NHE1 activity, as its disintegration reduced steady-state pHi and the rate of pHi recovery after acid loading. We propose that regulatory lipid-protein co-structures may play equally important roles in other membrane proteins.

Structural determinants of protocadherin-15 mechanics and function in hearing and balance perception.

The vertebrate inner ear, responsible for hearing and balance, is able to sense minute mechanical stimuli originating from an extraordinarily broad range of sound frequencies and intensities or from head movements. Integral to these processes is the tip-link protein complex, which conveys force to open the inner-ear transduction channels that mediate sensory perception. Protocadherin-15 and cadherin-23, two atypically large cadherins with 11 and 27 extracellular cadherin (EC) repeats, are involved in deafness and balance disorders and assemble as parallel homodimers that interact to form the tip link. Here we report the X-ray crystal structure of a protocadherin-15 + cadherin-23 heterotetrameric complex at 2.9-Å resolution, depicting a parallel homodimer of protocadherin-15 EC1-3 molecules forming an antiparallel complex with two cadherin-23 EC1-2 molecules. In addition, we report structures for 10 protocadherin-15 fragments used to build complete high-resolution models of the monomeric protocadherin-15 ectodomain. Molecular dynamics simulations and validated crystal contacts are used to propose models for the complete extracellular protocadherin-15 parallel homodimer and the tip-link bond. Steered molecular dynamics simulations of these models suggest conditions in which a structurally diverse and multimodal protocadherin-15 ectodomain can act as a stiff or soft gating spring. These results reveal the structural determinants of tip-link-mediated inner-ear sensory perception and elucidate protocadherin-15's structural and adhesive properties relevant in disease.

Multidomain Convergence of Argonaute during RISC Assembly Correlates with the Formation of Internal Water Clusters.

Despite the relevance of Argonaute proteins in RNA silencing, little is known about the structural steps of small RNA loading to form RNA-induced silencing complexes (RISCs). We report the 1.9 Å crystal structure of human Argonaute4 with guide RNA. Comparison with the previously determined apo structure of Neurospora crassa QDE2 revealed that the PIWI domain has two subdomains. Binding of guide RNA fastens the subdomains, thereby rearranging the active-site residues and increasing the affinity for TNRC6 proteins. We also identified two water pockets beneath the nucleic acid-binding channel that appeared to stabilize the mature RISC. Indeed, mutating the water-pocket residues of Argonaute2 and Argonaute4 compromised RISC assembly. Simulations predict that internal water molecules are exchangeable with the bulk solvent but always occupy specific positions at the domain interfaces. These results suggest that after guide RNA-driven conformational changes, water-mediated hydrogen-bonding networks tie together the converged domains to complete the functional RISC structure.

A Mechanically Weak Extracellular Membrane-Adjacent Domain Induces Dimerization of Protocadherin-15.

The cadherin superfamily of proteins is defined by the presence of extracellular cadherin (EC) "repeats" that engage in protein-protein interactions to mediate cell-cell adhesion, cell signaling, and mechanotransduction. The extracellular domains of nonclassical cadherins often have a large number of EC repeats along with other subdomains of various folds. Protocadherin-15 (PCDH15), a protein component of the inner-ear tip link filament essential for mechanotransduction, has 11 EC repeats and a membrane adjacent domain (MAD12) of atypical fold. Here we report the crystal structure of a pig PCDH15 fragment including EC10, EC11, and MAD12 in a parallel dimeric arrangement. MAD12 has a unique molecular architecture and folds as a ferredoxin-like domain similar to that found in the nucleoporin protein Nup54. Analytical ultracentrifugation experiments along with size-exclusion chromatography coupled to multiangle laser light scattering and small-angle x-ray scattering corroborate the crystallographic dimer and show that MAD12 induces parallel dimerization of PCDH15 near its membrane insertion point. In addition, steered molecular dynamics simulations suggest that MAD12 is mechanically weak and may unfold before tip-link rupture. Sequence analyses and structural modeling predict the existence of similar domains in cadherin-23, protocadherin-24, and the "giant" FAT and CELSR cadherins, indicating that some of them may also exhibit MAD-induced parallel dimerization.

Beyond Cell-Cell Adhesion: Sensational Cadherins for Hearing and Balance.

Cadherins form a large family of proteins often involved in calcium-dependent cellular adhesion. Although classical members of the family can provide a physical bond between cells, a subset of special cadherins use their extracellular domains to interlink apical specializations of single epithelial sensory cells. Two of these cadherins, cadherin-23 (CDH23) and protocadherin-15 (PCDH15), form extracellular "tip link" filaments that connect apical bundles of stereocilia on hair cells essential for inner-ear mechanotransduction. As these bundles deflect in response to mechanical stimuli from sound or head movements, tip links gate hair-cell mechanosensitive channels to initiate sensory perception. Here, we review the unusual and diverse structural properties of these tip-link cadherins and the functional significance of their deafness-related missense mutations. Based on the structural features of CDH23 and PCDH15, we discuss the elasticity of tip links and models that bridge the gap between the nanomechanics of cadherins and the micromechanics of hair-cell bundles during inner-ear mechanotransduction.

An elastic element in the protocadherin-15 tip link of the inner ear.

Tip link filaments convey force and gate inner-ear hair-cell transduction channels to mediate perception of sound and head movements. Cadherin-23 and protocadherin-15 form tip links through a calcium-dependent interaction of their extracellular domains made of multiple extracellular cadherin (EC) repeats. These repeats are structurally similar, but not identical in sequence, often featuring linkers with conserved calcium-binding sites that confer mechanical strength to them. Here we present the X-ray crystal structures of human protocadherin-15 EC8-EC10 and mouse EC9-EC10, which show an EC8-9 canonical-like calcium-binding linker, and an EC9-10 calcium-free linker that alters the linear arrangement of EC repeats. Molecular dynamics simulations and small-angle X-ray scattering experiments support this non-linear conformation. Simulations also suggest that unbending of EC9-10 confers some elasticity to otherwise rigid tip links. The new structure provides a first view of protocadherin-15's non-canonical EC linkers and suggests how they may function in inner-ear mechanotransduction, with implications for other cadherins.

Exploring the Membrane Potential of Simple Dual-Membrane Systems as Models for Gap-Junction Channels.

The conductance of ion channels can be modulated by a transmembrane potential difference, due to alterations on ion-mobility and also by changes in the pore structure. Despite the vast knowledge regarding the influence of voltage on transport properties of ion channels, little attention has been paid to describe, with atomic detail, the modulation of ionic transport in gap-junction channels (GJCs). Hence, molecular dynamics simulations were performed to explore the conductance of simple dual-membrane systems that account for the very basic features of GJCs. In doing so, we studied the influence of different charge distributions in the channel surface on these idealized systems under external electric fields, paying attention to the behavior of the electrostatic potential, ion density, ion currents, and equilibrium properties. Our results demonstrate that the incorporation of a charge distribution akin GJCs decreased anionic currents, favoring the transport of cationic species. Moreover, a thermodynamic characterization of ionic transport in these systems demonstrate the existence of a kinetic barrier that hinders anionic currents, reinforcing the role played by the internal arrangement of charges in GJCs. Overall, our results provide insights at the atomic scale on the effects of charge distributions over ionic transport, constituting a step forward into a better understanding of GJCs.

Characterization of a novel water pocket inside the human Cx26 hemichannel structure.

Connexins (Cxs) are a family of vertebrate proteins constituents of gap junction channels (GJCs) that connect the cytoplasm of adjacent cells by the end-to-end docking of two Cx hemichannels. The intercellular transfer through GJCs occurs by passive diffusion allowing the exchange of water, ions, and small molecules. Despite the broad interest to understand, at the molecular level, the functional state of Cx-based channels, there are still many unanswered questions regarding structure-function relationships, perm-selectivity, and gating mechanisms. In particular, the ordering, structure, and dynamics of water inside Cx GJCs and hemichannels remains largely unexplored. In this work, we describe the identification and characterization of a believed novel water pocket-termed the IC pocket-located in-between the four transmembrane helices of each human Cx26 (hCx26) monomer at the intracellular (IC) side. Using molecular dynamics (MD) simulations to characterize hCx26 internal water structure and dynamics, six IC pockets were identified per hemichannel. A detailed characterization of the dynamics and ordering of water including conformational variability of residues forming the IC pockets, together with multiple sequence alignments, allowed us to propose a functional role for this cavity. An in vitro assessment of tracer uptake suggests that the IC pocket residue Arg-143 plays an essential role on the modulation of the hCx26 hemichannel permeability.

Dual inhibitory pathways of metallofullerenol Gd@C82(OH)22 on matrix metalloproteinase-2: molecular insight into drug-like nanomedicine.

Cancer metastasis is an important criterion to evaluate tumor malignancy. Matrix metalloproteinases (MMPs) play a crucial role in cancer proliferation and migration by virtue of their proteolytic functions in angiogenesis and extracelluar matrix (ECM) degradation, making them potential targets of anti-metastaic therapeutics. Recently we showed with both in vivo and in vitro experiments that metallofullerenol Gd@C82(OH)22 can effectively inhibit MMP-2 and MMP-9 with high antitumoral efficacy. Furthermore, our in silico study revealed that Gd@C82(OH)22 could indirectly inhibit the proteolysis of MMP-9 via allosteric modulation exclusively at the ligand specificity S1' loop. Here, we expand our study toward another gelatinase, MMP-2, using molecular dynamics simulations. Despite the high structural similarity with 64.3% sequence identity, their responses to Gd@C82(OH)22 were quite different. Toward MMP-2, Gd@C82(OH)22 could block either the Zn(2+)-catalylitic site directly or the S1' loop indirectly. Surface electrostatics uniquely determines the initial adsorption of Gd@C82(OH)22 on MMP-2, and then its further location of the most favorable binding site(s). These findings not only illustrated how the inhibitory mechanism of Gd@C82(OH)22 is distinguished between the two gelatinase MMPs with atomic details, but also shed light on the de novo design of anti-metastatic nanotherapeutics with enhanced target specificity.

Critical role of the first transmembrane domain of Cx26 in regulating oligomerization and function.

To identify motifs involved in oligomerization of the gap junction protein Cx26, we studied individual transmembrane (TM) domains and the full-length protein. Using the TOXCAT assay for interactions of isolated TM α-helices, we found that TM1, a Cx26 pore domain, had a strong propensity to homodimerize. We identified amino acids Val-37-Ala-40 (VVAA) as the TM1 motif required for homodimerization. Two deafness-associated Cx26 mutations localized in this region, Cx26V37I and Cx26A40G, differentially affected dimerization. TM1-V37I dimerized only weakly, whereas TM1-A40G did not dimerize. When the full-length mutants were expressed in HeLa cells, both Cx26V37I and Cx26A40G formed oligomers less efficiently than wild-type Cx26. A Cx26 cysteine substitution mutant, Cx26V37C formed dithiothreitol-sensitive dimers. Substitution mutants of Val-37 formed intercellular channels with reduced function, while mutants of Ala-40 did not form functional gap junction channels. Unlike wild-type Cx26, neither Cx26V37I nor Cx26A40G formed functional hemichannels in low extracellular calcium. Thus the VVAA motif of Cx26 is critical for TM1 dimerization, hexamer formation, and channel function. The differential effects of VVAA mutants on hemichannels and gap junction channels imply that inter-TM interactions can differ in unapposed and docked hemichannels. Moreover, Cx26 oligomerization appears dependent on transient TM1 dimerization as an intermediate step.

Molecular modeling of Trypanosoma cruzi glutamate cysteine ligase and investigation of its interactions with glutathione.

Trypanosoma cruzi glutamate cysteine ligase (TcGCL) is considered a potential drug target to develop novel antichagasic drugs. We have used a variety of computational methods to investigate the interactions between TcGCL with Glutathione (GSH). The three-dimensional structure of TcGCL was constructed by comparative modeling methods using the Saccharomyces cerevisiae glutamate cysteine ligase as template. Molecular dynamics simulations were used to validate the TcGCL model and to analyze the molecular interactions with GSH. Using RMSD clustering, the most prevalent GSH binding modes were identified paying attention to the residues involved in the molecular interactions. The GSH binding modes were used to propose pharmacophore models that can be exploited in further studies to identify novel antichagasic compounds.

Molecular dynamics study of the archaeal aquaporin AqpM.

BACKGROUND: Aquaporins are a large family of transmembrane channel proteins that are present throughout all domains of life and are implicated in human disorders. These channels, allow the passive but selective movement of water and other small neutral solutes across cell membranes. Aquaporins have been classified into two sub-families: i) strict aquaporins that only allow the passage of water and ii) the less selective aquaglyceroporins that transport water and other neutral solutes, such as glycerol, CO2 or urea. Recently, the identification and characterization of a number of archaeal and bacterial aquaporins suggested the existence of a third sub-family; one that is neither a strict aquaporin nor an aquaglyceroporin. The function and phylogeny of this third family is still a matter of debate. RESULTS: Twenty nanosecond molecular dynamics (MD) simulation of a fully hydrated tetramer of AqpM embedded in a lipid bilayer permitted predictions to be made of key biophysical parameters including: single channel osmotic permeability constant (pf), single channel diffusive permeability constant (pd), channel radius, potential water occupancy of the channel and water orientation inside the pore. These properties were compared with those of well characterized representatives of the two main aquaporin sub-families. Results show that changes in the amino acid composition of the aromatic/arginine region affect the size and polarity of the selectivity filter (SF) and could help explain the difference in water permeability between aquaporins. In addition, MD simulation results suggest that AqpM combines characteristics of strict aquaporins, such as the narrow SF and channel radius, with those of aquaglyceroporins, such as a more hydrophobic and less polar SF. CONCLUSIONS: MD simulations of AqpM extend previous evidence that this archaeal aquaporin exhibits hybrid features intermediate between the two known aquaporin sub-families, supporting the idea that it may constitute a member of a novel class of aquaporins.

Computational modeling study of functional microdomains in cannabinoid receptor type 1.

The seven transmembrane helices (TMH) G-protein-coupled receptors (GPCRs) constitute one of the largest superfamily of signaling proteins found in mammals. Some of its members, in which the cannabinoid (CB) receptors are included, stand out because their functional states can be modulated by a broad spectrum of effector molecules. The relative ligand promiscuity exhibited by these receptors could be related with particular attributes conferred by their molecular architecture and represents a motivating issue to be explored. In this regard, this study represents an effort to investigate the cannabinoid receptor type 1 (CB1) ligand recognition plasticity, using comparative modeling, molecular dynamics (MD) simulations and docking. Our results suggest that a cooperative set of subtle structural rearrangements within the TMHs provide to the CB1 protein the plasticity to reach alternate configurations. These changes include the relaxation of intramolecular constraints, the rotations, translations and kinks of the majority of TMHs and the reorganization of the ligand binding cavities.