RESUMEN
We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, ß-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.
RESUMEN
Structure graphs, in which interacting amino acids/nucleotides correspond to linked nodes, represent cutting-edge tools to investigate macromolecular function. The graph-based approach defined as Protein Structure Network (PSN) was initially implemented in the Wordom software and subsequently in the webPSN server. PSNs are computed either on a molecular dynamics (MD) trajectory (PSN-MD) or on a single structure. In the latter case, information on atomic fluctuations is inferred from the Elastic Network Model-Normal Mode Analysis (ENM-NMA) (PSN-ENM). While Wordom performs both PSN-ENM and PSN-MD analyses but without output post-processing, the webPSN server performs only single-structure PSN-EMN but assisting the user in input setup and output analysis. Here we release for the first time the standalone software PSNtools, which allows calculation and post-processing of PSN analyses carried out either on single structures or on conformational ensembles. Relevant unique and novel features of PSNtools are either comparisons of two networks or computations of consensus networks on sets of homologous/analogous macromolecular structures or conformational ensembles. Network comparisons and consensus serve to infer differences in functionally different states of the same system or network-based signatures in groups of bio-macromolecules sharing either the same functionality or the same fold. In addition to the new software, here we release also an updated version of the webPSN server, which allows performing an interactive graphical analysis of PSN-MD, following the upload of the PSNtools output. PSNtools, the auxiliary binary version of Wordom software, and the WebPSN server are freely available at http://webpsn.hpc.unimo.it/wpsn3.php.
RESUMEN
A mixed Protein Structure Network (PSN) and Elastic Network Model-Normal Mode Analysis (ENM-NMA)-based strategy (i.e. PSN-ENM) was developed to investigate structural communication in bio-macromolecules. Protein Structure Graphs (PSGs) are computed on a single structure, whereas information on system dynamics is supplied by ENM-NMA. The approach was implemented in a webserver (webPSN), which was significantly updated herein. The webserver now handles both proteins and nucleic acids and relies on an internal upgradable database of network parameters for ions and small molecules in all PDB structures. Apart from the radical restyle of the server and some changes in the calculation setup, other major novelties concern the possibility to: a) compute the differences in nodes, links, and communication pathways between two structures (i.e. network difference) and b) infer links, hubs, communities, and metapaths from consensus networks computed on a number of structures. These new features are useful to identify commonalties and differences between two different functional states of the same system or structural-communication signatures in homologous or analogous systems. The output analysis relies on 3D-representations, interactive tables and graphs, also available for download. Speed and accuracy make this server suitable to comparatively investigate structural communication in large sets of bio-macromolecular systems. URL: http://webpsn.hpc.unimore.it.
Asunto(s)
Conformación Proteica , Programas Informáticos , Internet , Modelos Moleculares , Conformación de Ácido Nucleico , Ácidos Nucleicos/química , Proteínas/químicaRESUMEN
Uncontrolled activation of Rho signaling by RhoGEFs, in particular AKAP13 (Lbc) and its close homologs, is implicated in a number of human tumors with poor prognosis and resistance to therapy. Structure predictions and alanine scanning mutagenesis of Lbc identified a circumscribed hot region for RhoA recognition and activation. Virtual screening targeting that region led to the discovery of an inhibitor of Lbc-RhoA interaction inside cells. By interacting with the DH domain, the compound inhibits the catalytic activity of Lbc, halts cellular responses to activation of oncogenic Lbc pathways, and reverses a number of prostate cancer cell phenotypes such as proliferation, migration, and invasiveness. This study provides insights into the structural determinants of Lbc-RhoA recognition. This is a successful example of structure-based discovery of a small protein-protein interaction inhibitor able to halt oncogenic Rho signaling in cancer cells with therapeutic implications.
Asunto(s)
Proteínas de Anclaje a la Quinasa A/antagonistas & inhibidores , Neoplasias/tratamiento farmacológico , Proteínas Proto-Oncogénicas/antagonistas & inhibidores , Transducción de Señal/efectos de los fármacos , Bibliotecas de Moléculas Pequeñas/farmacología , Proteínas de Unión al GTP rho/antagonistas & inhibidores , Proteínas de Anclaje a la Quinasa A/metabolismo , Humanos , Antígenos de Histocompatibilidad Menor/metabolismo , Modelos Moleculares , Estructura Molecular , Neoplasias/metabolismo , Unión Proteica/efectos de los fármacos , Proteínas Proto-Oncogénicas/metabolismo , Bibliotecas de Moléculas Pequeñas/química , Proteínas de Unión al GTP rho/metabolismoRESUMEN
G protein coupled receptors (GPCRs) are allosteric proteins whose functioning fundamentals are the communication between the two poles of the helix bundle. Protein structure network (PSN) analysis is one of the graph theory-based approaches currently used to investigate the structural communication in biomolecular systems. Information on system's dynamics can be provided by atomistic molecular dynamics (MD) simulations or coarse grained elastic network models paired with normal mode analysis (ENM-NMA). The present review article describes the application of PSN analysis to uncover the structural communication in G protein coupled receptors (GPCRs). Strategies to highlight changes in structural communication upon misfolding, dimerization and activation are described. Focus is put on the ENM-NMA-based strategy applied to the crystallographic structures of rhodopsin in its inactive (dark) and signalling active (meta II (MII)) states, highlighting changes in structure network and centrality of the retinal chromophore in differentiating the inactive and active states of the receptor.
Asunto(s)
Receptores Acoplados a Proteínas G/fisiología , Simulación de Dinámica Molecular , Receptores Acoplados a Proteínas G/químicaRESUMEN
UNLABELLED: We developed a mixed Protein Structure Network (PSN) and Elastic Network Model-Normal Mode Analysis (ENM-NMA)-based strategy (i.e. PSN-ENM) to investigate structural communication in biomacromolecules. The approach starts from a Protein Structure Graph and searches for all shortest communication pathways between user-specified residues. The graph is computed on a single preferably high-resolution structure. Information on system's dynamics is supplied by ENM-NMA. The PSN-ENM methodology is made of multiple steps both in the setup and analysis stages, which may discourage inexperienced users. To facilitate its usage, we implemented WebPSN, a freely available web server that allows the user to easily setup the calculation, perform post-processing analyses and both visualize and download numerical and 3D representations of the output. Speed and accuracy make this server suitable to investigate structural communication, including allosterism, in large sets of bio-macromolecular systems. AVAILABILITY AND IMPLEMENTATION: The WebPSN server is freely available at http://webpsn.hpc.unimore.it.
Asunto(s)
Ensayos Analíticos de Alto Rendimiento , Internet , Proteínas/química , Análisis de Secuencia de Proteína/métodos , Programas Informáticos , Homología Estructural de Proteína , Humanos , Simulación de Dinámica Molecular , Dominios y Motivos de Interacción de Proteínas , Alineación de SecuenciaRESUMEN
In spite of the ever-increasing evidence that G protein-coupled receptors (GPCRs) form dimers/oligomers, the biological role(s) and structural architecture of homologous and heterologous receptor aggregation are, however, far from being clarified. This chapter reviews the insights gained so far, at multiscale levels of resolution, on GPCR dimerization/oligomerization from in vitro experiments, structure predictions, and structure determinations. Focus is put on the achievement by the FiPD-based approach, which proved effective in predicting the supramolecular organization of membrane proteins including GPCRs. The combination of FiPD-based quaternary structure predictions with molecular simulations and analyses can be a valuable tool to infer the effects of dimerization on the structural communication features of a receptor dimer/oligomer bound to functionally different ligands. Ultimately, the integration between atomistic and mesoscopic simulations is expected to be a promising tool to unveil functioning mechanisms that involve intricate protein networks.
Asunto(s)
Biología Computacional/métodos , Multimerización de Proteína , Receptores Acoplados a Proteínas G/química , Animales , Humanos , Modelos Moleculares , Estructura Cuaternaria de Proteína , Análisis de Secuencia de ProteínaRESUMEN
Graph theory is being increasingly used to study the structural communication in biomolecular systems. This requires incorporating information on the system's dynamics, which is time-consuming and not suitable for high-throughput investigations. We propose a mixed Protein Structure Network (PSN) and Elastic Network Model (ENM)-based strategy, i.e., PSN-ENM, for fast investigation of allosterism in biological systems. PSN analysis and ENM-Normal Mode Analysis (ENM-NMA) are implemented in the structural analysis software Wordom, freely available at http://wordom.sourceforge.net/ . The method performs a systematic search of the shortest communication pathways that traverse a protein structure. A number of strategies to compare the structure networks of a protein in different functional states and to get a global picture of communication pathways are presented as well. The approach was validated on the PDZ2 domain from tyrosine phosphatase 1E (PTP1E) in its free (APO) and peptide-bound states. PDZ domains are, indeed, the systems whose structural communication and allosteric features are best characterized both in vitro and in silico. The agreement between predictions by the PSN-ENM method and in vitro evidence is remarkable and comparable to or higher than that reached by more time-consuming computational approaches tested on the same biological system. Finally, the PSN-ENM method was able to reproduce the salient communication features of unbound and bound PTP1E inferred from molecular dynamics simulations. High speed makes this method suitable for high throughput investigation of the communication pathways in large sets of biomolecular systems in different functional states.
RESUMEN
Wordom is a versatile, user-friendly, and efficient program for manipulation and analysis of molecular structures and dynamics. The following new analysis modules have been added since the publication of the original Wordom paper in 2007: assignment of secondary structure, calculation of solvent accessible surfaces, elastic network model, motion cross correlations, protein structure network, shortest intra-molecular and inter-molecular communication paths, kinetic grouping analysis, and calculation of mincut-based free energy profiles. In addition, an interface with the Python scripting language has been built and the overall performance and user accessibility enhanced. The source code of Wordom (in the C programming language) as well as documentation for usage and further development are available as an open source package under the GNU General Purpose License from http://wordom.sf.net.
Asunto(s)
Modelos Moleculares , Péptidos/química , Proteínas/química , Programas Informáticos , Algoritmos , Animales , Cinética , Estructura Molecular , Propiedades de Superficie , Termodinámica , VibraciónRESUMEN
Disease-causing missense mutations in membrane proteins, such as rhodopsin mutations associated with the autosomal dominant form of retinitis pigmentosa (ADRP), are often linked to defects in folding and/or trafficking. The mechanical unfolding of wild-type rhodopsin was compared with that of 20 selected ADRP-linked mutants more or less defective in folding and retinal binding. Rhodopsin fold is characterized by networks of amino acids in the retinal and G-protein binding sites likely to play a role in the stability and function of the protein. The distribution of highly connected nodes in the network reflects the existence of a diffuse intramolecular communication inside and between the 2 poles of the helix bundle, which makes pathogenic mutations share similar phenotypes irrespective of topological and physicochemical differences between them. Because of this communication, the ADRP-linked rhodopsin mutations share a more or less marked ability to impair selected hubs in the protein structure network. The extent of this structural effect relates to the severity of the biochemical defect caused by mutation. The investigative strategy employed in this study is likely to apply to all structurally known membrane proteins particularly susceptible to misassembly-causing mutations.
Asunto(s)
Simulación de Dinámica Molecular , Mutación/genética , Retinitis Pigmentosa/metabolismo , Rodopsina/química , Rodopsina/genética , Secuencia de Aminoácidos , Animales , Bovinos , Modelos Moleculares , Datos de Secuencia Molecular , Pliegue de Proteína , Estructura Secundaria de Proteína , Retinitis Pigmentosa/genéticaRESUMEN
Intramolecular and intermolecular communication is a privileged issue in G protein-Coupled Receptor (GPCR) function as the prominent role of these receptors is to respond to extracellular signals by catalyzing nucleotide exchange in intracellular G proteins. In the last decade or so we have applied much effort in elaborating computational strategies to infer the mechanisms of intramolecular and intermolecular communication in a number of GPCRs of the rhodopsin family. In this article, we review the most relevant achievements on the matter. In summary, the receptor sites of activating mutations or ligand-binding communicate with a common allosteric site in the cytosolic domains. This was inferred from the observation that local perturbations by activating mutations or ligands correlate with increases in solvent accessibility of the neighborhoods of the highly conserved E/DRY receptor motif. The latter turned out to be the primary recognition point for the C-terminus of the G protein alpha-subunit, independent of the receptor or the G protein type. In spite of the highly composite nature of the receptor-G protein interface, receptor contacts with the C-terminus of the alpha5-helix seem to be the major players in the receptor-catalyzed formation of a nucleotide exit route. The latter would lie in between the alphaF-helix and the beta6/alpha5 loop, which detach from each other upon receptor binding, giving solvent accessibility to the nucleotide. A worthy inference of the studies is that GPCRs employ common pathways for the transfer of functionally relevant information.
Asunto(s)
Simulación por Computador , Modelos Moleculares , Receptores Acoplados a Proteínas G/química , Secuencias de Aminoácidos , Sitios de Unión , Ligandos , Mutación , Unión Proteica , Subunidades de Proteína/química , Subunidades de Proteína/genética , Receptores Acoplados a Proteínas G/genéticaRESUMEN
Multiple 4-micros molecular dynamics (MD) simulations are used to study the folding process of the cross-linked alpha-helical peptide Ac-EACAR(5)EAAAR(10)EAACR(15)Q-NH(2) (EAAAR peptide). The folding kinetics are single exponential at 330 K, while they are complex at 281 K with a clear deviation from single-exponential behavior, in agreement with time-resolved infrared (IR) spectroscopy measurements. Network analysis of the conformation space sampled by the MD simulations reveals four main folding channels which start from conformations with partially formed helical structure and non-native salt-bridges in a kinetically partitioned unfolded state. The independent folding pathways explain the comparable quality of models based on stretched exponential and multiexponential fitting of the kinetic traces at low temperature. The rearrangement of bulky side chains, and in particular their reorientation with respect to the cross-linker, makes the EAAAR peptide a slower folder at 281 K than a similar peptide devoid of the three glutamate side chains. On the basis of this simulation result, extracted from a total MD sampling of 1.0 ms, a mutant with additional bulky side chains (three methionines replacing alanines at positions 2, 7, and 12) is suggested to fold slower than the EAAAR peptide. This prediction is confirmed by time-resolved IR spectroscopy.
Asunto(s)
Péptidos/química , Péptidos/metabolismo , Pliegue de Proteína , Sales (Química)/química , Simulación por Computador , Cinética , Modelos Moleculares , Mutación/genética , Péptidos/genética , Estructura Secundaria de Proteína , Estructura Terciaria de Proteína , Espectrofotometría Infrarroja , Factores de TiempoRESUMEN
Retinitis pigmentosa (RP) refers to a group of debilitating, hereditary disorders that cause severe visual impairment in as many as 1.5 million patients worldwide. Rhodopsin mutations account for >25% of the autosomal dominant form of the disease (ADRP). Forty artificial and ADRP-associated mutations located in the second extracellular loop (EL2) that folds into a twisted ß-hairpin were screened through replica exchange molecular dynamics (REMD) simulations using the FACTS implicit solvent model. According to in vitro experiments, ADRP-linked mutants fail to express at the plasma membrane and/or to reconstitute with 11-cis-retinal, indicative of variable defects in protein folding and/or stability. The computational protocol was first probed on the protein G C-terminal ß-hairpin, proving the effectiveness of the implicit solvent model in reproducing the free energy landscape of ß-hairpin formation. Eight out of the 40 EL2 mutants resulted in misfolding effects on the native ß-hairpin structure, consistent with in vitro evidence that they all share severe impairments in folding/expression. Five mutants displayed moderate misfolding attitudes, whereas the remaining 27 mutants, overall characterized by milder effects on rhodopsin expression, did not perturb significantly the conformational behavior of the native ß-hairpin but are expected to exert variably disturbing effects on the native interactions of the loop with the chromophore and/or the surrounding receptor domains. Collectively, the results of this study add structural insight to the poorly resolved biochemical behavior of selected class II ADRP mutations, a fundamental step toward an understanding of the atomistic causes of the disease.
RESUMEN
This study represents the first attempt to couple, by computational experiments, the mechanisms of intramolecular and intermolecular communication concerning a guanidine nucleotide exchange factor (GEF), the thromboxane A2 receptor (TXA2R), and the cognate G protein (Gq) in its heterotrimeric GDP-bound state. Two-way pathways mediate the communication between the receptor-G protein interface and both the agonist binding site of the receptor and the nucleotide binding site of the G protein. The increase in solvent accessibility in the neighborhoods of the highly conserved E/DRY receptor motif, in response to agonist binding, is instrumental in favoring the penetration of the C-terminus of Gqalpha in between the cytosolic ends of H3, H5, and H6. The arginine of the E/DRY motif is predicted to be an important mediator of the intramolecular and intermolecular communication involving the TXA2R. The receptor-G protein interface is predicted to involve multiple regions from the receptor and the G protein alpha-subunit. However, receptor contacts with the C-terminus of the alpha5-helix seem to be the major players in the receptor-catalyzed motion of the alpha-helical domain with respect to the Ras-like domain and in the formation of a nucleotide exit route in between the alphaF-helix and beta6/alpha5 loop of Gqalpha. The inferences from this study are of wide interest, as they are expected to apply to the whole rhodopsin family, given also the considerable G protein promiscuity.
Asunto(s)
Proteínas de Unión al GTP/química , Receptores Acoplados a Proteínas G/química , Secuencia de Aminoácidos , Animales , Sitios de Unión , Simulación por Computador , Ratones , Modelos Químicos , Datos de Secuencia Molecular , Conformación Proteica , Estructura Terciaria de Proteína , Transducción de SeñalRESUMEN
Wordom is a versatile program for manipulation of molecular dynamics trajectories and efficient analysis of simulations. Original tools in Wordom include a procedure to evaluate significance of sampling for principal component analysis as well as modules for clustering multiple conformations and evaluation of order parameters for folding and aggregation. The program was developed with special emphasis on user-friendliness, effortless addition of new modules and efficient handling of large sets of trajectories.
Asunto(s)
Algoritmos , Bases de Datos de Proteínas , Almacenamiento y Recuperación de la Información/métodos , Modelos Químicos , Modelos Moleculares , Proteínas/química , Proteínas/ultraestructura , Programas Informáticos , Simulación por Computador , Sistemas de Administración de Bases de Datos , Conformación Proteica , Interfaz Usuario-ComputadorRESUMEN
In this computational study, we have investigated the implications of rhodopsin (Rho) oligomerization in transducin (Gt) recognition. The results of docking simulations between heterotrimeric Gt and monomeric, dimeric and tetrameric inactive Rho corroborate the hypothesis that Rho and Gt can be found coupled already in the dark. Moreover, our extensive computational analysis suggests that the most likely Rho:Gt stoichiometry is the 1:1 one. This means that the essential molecular determinants for Gt recognition and activation are contained in one Rho monomer. In this respect, the complex between one Rho molecule and one heterotrimeric Gt should be considered as the functional unit.
Asunto(s)
Rodopsina/química , Transducina/química , Animales , Sitios de Unión , Bovinos , Biología Computacional , Simulación por Computador , Oscuridad , Técnicas In Vitro , Modelos Moleculares , Complejos Multiproteicos , Estructura Cuaternaria de Proteína , Rodopsina/metabolismo , Transducina/metabolismoRESUMEN
Multiple molecular dynamics simulations of bacterioopsin pulling from its C-terminus show that its alpha-helices unfold individually. In the first metastable state observed in the simulations, helix G is unfolded at its C-terminal segment while the rest of helix G (residues 200-216) is folded and opposes resistance because of a salt-bridge network consisting of Asp-212 and Lys-216 on helix G and Arg-82 and Asp-85 on helix C. Helix G unfolds inside the bundle because the external force is applied to its C-terminal end in a direction perpendicular to the surface of the membrane. Inversely, helix F has to flip by 180 degrees to exit from the membrane because the applied force and the helical N-C axis point in opposite directions. At the highest peak of the force, which cannot be interpreted in single-molecule force spectroscopy experiments, helix F has a pronounced kink at Pro-186. Mutation of Pro-186 and/or the charged side chains mentioned above, which are involved in very favorable electrostatic interactions in the low-dielectric region of the membrane, are expected to reduce the highest peak of the force. Helices E and D unfold in a similar way to helices G and F, respectively. Hence, the force-distance profile and sequence of events during forced unfolding of bacterioopsin are influenced by the up-and-down topology of the seven-helix bundle. The sequential extraction of individual helices from the membrane suggests that the spontaneous (un)folding of bacterioopsin proceeds through metastable bundles of fewer than seven helices. The metastable states observed in the simulations provide atomic level evidence that corroborates the interpretation of very recent force spectroscopy experiments of bacteriorhodopsin refolding.
Asunto(s)
Bacteriorodopsinas/química , Bacteriorodopsinas/ultraestructura , Modelos Químicos , Modelos Moleculares , Membrana Púrpura/química , Simulación por Computador , Elasticidad , Cinética , Movimiento (Física) , Conformación Proteica , Desnaturalización Proteica , Pliegue de Proteína , Estrés MecánicoRESUMEN
Fragment complementation is gaining an increasing impact as a nonperturbing method to probe noncovalent interactions within protein supersecondary structures. In this study, the fast Fourier transform rigid-body docking algorithm ZDOCK has been employed for in silico reconstitution of the calcium binding protein calbindin D9k, from its two EF-hands subdomains, namely, EF1 (residues 1-43) and EF2 (residues 44-75). The EF1 fragment has been used both in its wild type and in nine mutant forms, in line with in vitro experiments. Consistent with in vitro data, ZDOCK reconstituted the proper fold of wild-type and mutated calbindin, locating the nativelike structures (i.e., holding a root-mean-square deviation < 1 A with respect to the X-ray structure) among the first 10 top-scored solutions out of 4000. Moreover, the three independent in silico reconstitutions of wild-type calbindin ranked a nativelike structure at the top of the output list, that is, the best scored one. The algorithm has been also successfully challenged in reconstituting the EF2 homodimer from two identical copies of the monomer. Furthermore, quantitative models consisting of linear correlations between thermodynamic data and ZDOCK scores were built, providing a tested tool for very fast in silico predictions of the free energy of association of protein-protein complexes solved at the atomic level and known to not undergo significant conformational changes upon binding.
Asunto(s)
Biología Computacional , Simulación por Computador , Proteína G de Unión al Calcio S100/química , Proteína G de Unión al Calcio S100/metabolismo , Calbindinas , Mutación , Conformación Proteica , Pliegue de Proteína , Subunidades de Proteína , Proteína G de Unión al Calcio S100/genética , Relación Estructura-ActividadRESUMEN
In the present study, more than 75 compounds structurally related to BMY 7378 have been designed and synthesized. Structural variations of each part of the reference molecule have been introduced, obtaining highly selective ligands for the alpha(1d) adrenergic receptor. The molecular determinants for selectivity at this receptor are essentially held by the phenyl substituent in the phenylpiperazine moiety. The integration of an extensive SAR analysis with docking simulations using the rhodopsin-based models of the three alpha(1)-AR subtypes and of the 5-HT(1A) receptor provides significant insights into the characterization of the receptor binding sites as well as into the molecular determinants of ligand selectivity at the alpha(1d)-AR and the 5-HT(1A) receptors. The results of multiple copies simultaneous search (MCSS) on the substituted phenylpiperazines together with those of manual docking of compounds BMY 7378 and 69 into the putative binding sites of the alpha(1a)-AR, alpha(1b)-AR, alpha(1d)-AR, and the 5-HT(1A) receptors suggest that the phenylpiperazine moiety would dock into a site formed by amino acids in helices 3, 4, 5, 6 and extracellular loop 2 (E2), whereas the spirocyclic ring of the ligand docks into a site formed by amino acids of helices 1, 2, 3, and 7. This docking mode is consistent with the SAR data produced in this work. Furthermore, the binding site of the imide moiety does not allow for the simultaneous involvement of the two carbonyl oxygen atoms in H-bonding interactions, consistent with the SAR data, in particular with the results obtained with the lactam derivative 128. The results of docking simulations also suggest that the second and third extracellular loops may act as selectivity filters for the substituted phenylpiperazines. The most potent and selective compounds for alpha(1d) adrenergic receptor, i.e., 69 (Rec 26D/038) and 128 (Rec 26D/073), are characterized by the presence of the 2,5-dichlorophenylpiperazine moiety.
Asunto(s)
Antagonistas de Receptores Adrenérgicos alfa 1 , Antagonistas Adrenérgicos alfa/síntesis química , Imidas/síntesis química , Compuestos de Espiro/síntesis química , Antagonistas Adrenérgicos alfa/química , Antagonistas Adrenérgicos alfa/farmacología , Secuencia de Aminoácidos , Animales , Sitios de Unión , Células CHO , Cricetinae , Cobayas , Células HeLa , Humanos , Imidas/química , Imidas/farmacología , Técnicas In Vitro , Ligandos , Masculino , Modelos Moleculares , Datos de Secuencia Molecular , Contracción Muscular/efectos de los fármacos , Músculo Liso Vascular/efectos de los fármacos , Músculo Liso Vascular/fisiología , Estructura Terciaria de Proteína , Ensayo de Unión Radioligante , Ratas , Ratas Sprague-Dawley , Receptor de Serotonina 5-HT1A/metabolismo , Receptores Adrenérgicos alfa 1/química , Receptores Adrenérgicos alfa 1/metabolismo , Alineación de Secuencia , Antagonistas del Receptor de Serotonina 5-HT1 , Compuestos de Espiro/química , Compuestos de Espiro/farmacología , Relación Estructura-ActividadRESUMEN
Comparative molecular dynamics simulations of the 5-HT(1A) receptor in its empty as well as agonist- (i.e. active) and antagonist-bound (i.e. nonactive) forms have been carried out. The agonists 5-HT and (R)-8-OH-DPAT as well as the antagonist WAY100635 have been employed. The results of this study strengthen the hypothesis that the receptor portions close to the E/DRY/W motif, with prominence to the cytosolic extensions of helices 3 and 6, are particularly susceptible to undergo structural modification in response to agonist binding. Despite the differences in the structural/dynamics behavior of the two agonists when docked into the 5-HT(1A) receptor, they both exert a destabilization of the intrahelical and interhelical interactions found in the empty and antagonist-bound receptor forms between the arginine of the E/DRY sequence and both D133(3.49) and E340(6.30). For both agonists, the chemical information transfer from the extracellular to the cytosolic domains is mediated by a cluster of aromatic amino acids in helix 6, following the ligand interaction with selected amino acids in the extracellular half of the receptor, such as D116(3.32), S199(5.42), Y195(5.38), and F361(6.51). A significant reduction in the bend at P360(6.50), as compared to the empty and the antagonist-bound receptor forms, is one of the features of the agonist-bound forms that is related to the breakage of the interhelical salt bridge between the E/DRY arginine and E340(6.30). Another structural feature, shared by the agonist-bound receptor forms and not by the empty and antagonist-bound forms, is the detachment of helices 2 and 4, as marked by the movement of W161(4.50) away from helix 2, toward the membrane space.