stringtranslate.com

Ligando (bioquímica)

Mioglobina (azul) con su ligando hemo (naranja) unido. Basado en PDB : 1MBO

En bioquímica y farmacología , un ligando es una sustancia que forma un complejo con una biomolécula para cumplir un propósito biológico. La etimología proviene del latín ligare , que significa "unir". En la unión proteína-ligando, el ligando suele ser una molécula que produce una señal al unirse a un sitio de una proteína diana . La unión normalmente da como resultado un cambio de isomería conformacional (conformación) de la proteína diana. En los estudios de unión de ADN-ligando, el ligando puede ser una pequeña molécula, ion , [1] o proteína [2] que se une a la doble hélice del ADN . La relación entre el ligando y el compañero de unión es función de la carga, la hidrofobicidad y la estructura molecular.

La unión se produce mediante fuerzas intermoleculares , como enlaces iónicos , enlaces de hidrógeno y fuerzas de Van der Waals . La asociación o acoplamiento es en realidad reversible mediante disociación . El enlace covalente mensurablemente irreversible entre un ligando y una molécula objetivo es atípico en los sistemas biológicos. A diferencia de la definición de ligando en química metalorgánica e inorgánica , en bioquímica es ambiguo si el ligando generalmente se une a un sitio metálico , como es el caso de la hemoglobina . En general, la interpretación del ligando es contextual con respecto al tipo de unión que se ha observado.

La unión del ligando a una proteína receptora altera la conformación al afectar la orientación de la forma tridimensional. La conformación de una proteína receptora compone el estado funcional. Los ligandos incluyen sustratos , inhibidores , activadores , lípidos de señalización y neurotransmisores . La tasa de unión se llama afinidad y esta medida tipifica una tendencia o fuerza del efecto. La afinidad de unión se actualiza no solo por las interacciones huésped-huésped , sino también por los efectos del solvente que pueden desempeñar un papel estérico dominante que impulsa la unión no covalente en la solución. [3] El disolvente proporciona un entorno químico para que el ligando y el receptor se adapten y, por tanto, se acepten o rechacen mutuamente como socios.

Los radioligandos son compuestos marcados con radioisótopos que se utilizan in vivo como trazadores en estudios de PET y para estudios de unión in vitro .

Afinidad de unión al receptor/ligando

La interacción de ligandos con sus sitios de unión se puede caracterizar en términos de afinidad de unión. En general, la unión de ligandos de alta afinidad resulta de mayores fuerzas de atracción entre el ligando y su receptor , mientras que la unión de ligandos de baja afinidad implica menos fuerzas de atracción. En general, la unión de alta afinidad da como resultado una mayor ocupación del receptor por su ligando que en el caso de la unión de baja afinidad; el tiempo de residencia (vida útil del complejo receptor-ligando) no se correlaciona. La unión de alta afinidad de ligandos a receptores suele ser fisiológicamente importante cuando parte de la energía de unión se puede utilizar para provocar un cambio conformacional en el receptor, lo que da como resultado un comportamiento alterado, por ejemplo, de un canal iónico o enzima asociado .

Un ligando que puede unirse y alterar la función del receptor que desencadena una respuesta fisiológica se llama agonista del receptor . Los ligandos que se unen a un receptor pero no logran activar la respuesta fisiológica son antagonistas del receptor .

Dos agonistas con afinidad de unión similar

La unión del agonista a un receptor se puede caracterizar tanto en términos de cuánta respuesta fisiológica se puede desencadenar (es decir, la eficacia ) como en términos de la concentración del agonista que se requiere para producir la respuesta fisiológica (a menudo medida como CE 50 , la concentración requerida para producir la mitad de la respuesta máxima). La unión de ligando de alta afinidad implica que una concentración relativamente baja de un ligando es adecuada para ocupar al máximo un sitio de unión de ligando y desencadenar una respuesta fisiológica. La afinidad del receptor se mide mediante una constante de inhibición o valor Ki , la concentración necesaria para ocupar el 50% del receptor. Las afinidades de los ligandos se miden con mayor frecuencia indirectamente como un valor IC50 de un experimento de unión competitiva donde se determina la concentración de un ligando necesaria para desplazar el 50% de una concentración fija de ligando de referencia . El valor de Ki se puede estimar a partir de IC 50 mediante la ecuación de Cheng Prusoff . Las afinidades de los ligandos también se pueden medir directamente como una constante de disociación (Kd ) utilizando métodos como la extinción de la fluorescencia , la calorimetría de titulación isotérmica o la resonancia de plasmones superficiales . [4]

La unión de baja afinidad (alto nivel de Ki ) implica que se requiere una concentración relativamente alta de un ligando antes de que el sitio de unión esté ocupado al máximo y se logre la máxima respuesta fisiológica al ligando. En el ejemplo que se muestra a la derecha, dos ligandos diferentes se unen al mismo sitio de unión al receptor. Sólo uno de los agonistas mostrados puede estimular al máximo el receptor y, por tanto, puede definirse como un agonista completo . Un agonista que sólo puede activar parcialmente la respuesta fisiológica se llama agonista parcial . En este ejemplo, la concentración a la que el agonista completo (curva roja) puede activar el receptor a la mitad del máximo es aproximadamente 5 x 10 −9 molar (nM = nanomolar ).

Dos ligandos con diferente afinidad de unión al receptor.

La afinidad de unión se determina más comúnmente utilizando un ligando radiomarcado , conocido como ligando marcado. Los experimentos de unión competitiva homóloga implican una competencia de unión entre un ligando marcado y un ligando no marcado. [5] Los métodos basados ​​en tiempo real, que a menudo no tienen etiquetas, como la resonancia de plasmón de superficie , la interferometría de doble polarización y la resonancia de plasmón de superficie multiparamétrica (MP-SPR), no solo pueden cuantificar la afinidad a partir de ensayos basados ​​en la concentración; sino también de la cinética de asociación y disociación y, en los últimos casos, del cambio conformacional inducido tras la unión. MP-SPR también permite mediciones en tampones de disociación con alto contenido salino gracias a una configuración óptica única. Se desarrolló la termoforesis a microescala (MST), un método sin inmovilización [6] . Este método permite la determinación de la afinidad de unión sin ninguna limitación al peso molecular del ligando. [7]

Para conocer el uso de la mecánica estadística en un estudio cuantitativo de la afinidad de unión ligando-receptor, consulte el artículo completo [8] sobre la función de partición configuracional .

Potencia de unión a fármacos u hormonas

Los datos de afinidad de unión por sí solos no determinan la potencia general de un fármaco o de una hormona producida naturalmente (biosintetizada). [9]

La potencia es el resultado de la compleja interacción tanto de la afinidad de unión como de la eficacia del ligando. [9]

Eficacia de unión a fármacos u hormonas

La eficacia del ligando se refiere a la capacidad del ligando para producir una respuesta biológica al unirse al receptor diana y la magnitud cuantitativa de esta respuesta. Esta respuesta puede ser como agonista , antagonista o agonista inverso , dependiendo de la respuesta fisiológica producida. [10]

Selectivo y no selectivo

Los ligandos selectivos tienden a unirse a tipos muy limitados de receptores, mientras que los ligandos no selectivos se unen a varios tipos de receptores. Esto juega un papel importante en farmacología , donde los fármacos que no son selectivos tienden a tener más efectos adversos , porque se unen a varios otros receptores además del que genera el efecto deseado.

Ligandos hidrófobos

Para ligandos hidrófobos (p. ej., PIP2) en complejo con una proteína hidrófoba (p. ej., canales iónicos activados por lípidos ), la determinación de la afinidad se complica por interacciones hidrófobas no específicas. Las interacciones hidrofóbicas no específicas pueden superarse cuando la afinidad del ligando es alta. [11] Por ejemplo, PIP2 se une con alta afinidad a los canales iónicos activados por PIP2.

Ligando bivalente

Los ligandos bivalentes constan de dos moléculas similares a fármacos (farmacóforos o ligandos) conectadas por un conector inerte. Hay varios tipos de ligandos bivalentes y, a menudo, se clasifican según el objetivo de los farmacóforos. Los ligandos homobivalentes se dirigen a dos de los mismos tipos de receptores. Los ligandos heterobivalentes se dirigen a dos tipos de receptores diferentes. [12] Los ligandos bitópicos se dirigen a sitios de unión ortostéricos y a sitios de unión alostéricos en el mismo receptor. [13] En la investigación científica, los ligandos bivalentes se han utilizado para estudiar los dímeros de los receptores e investigar sus propiedades. Philip S. Portoghese y sus compañeros de trabajo fueron pioneros en esta clase de ligandos mientras estudiaban el sistema de receptores de opioides. [14] [15] [16] Micheal Conn y sus compañeros de trabajo también informaron desde el principio sobre ligandos bivalentes para el receptor de la hormona liberadora de gonadotropina. [17] [18] Desde estos primeros informes, se han informado muchos ligandos bivalentes para varios sistemas de receptores acoplados a proteínas G (GPCR), incluidos cannabinoides, [19] serotonina, [20] [21] oxitocina, [22] y melanocortina. sistemas receptores, [23] [24] [25] y para sistemas GPCR - LIC (receptores D2 y nACh ). [12]

Los ligandos bivalentes generalmente tienden a ser más grandes que sus homólogos monovalentes y, por lo tanto, no "parecidos a fármacos" como en la regla de cinco de Lipinski . Muchos creen que esto limita su aplicabilidad en entornos clínicos. [26] [27] A pesar de estas creencias, ha habido muchos ligandos que han reportado estudios preclínicos exitosos en animales. [24] [25] [22] [28] [29] [30] Dado que algunos ligandos bivalentes pueden tener muchas ventajas en comparación con sus homólogos monovalentes (como selectividad tisular, mayor afinidad de unión y mayor potencia o eficacia), los bivalentes También puede ofrecer algunas ventajas clínicas.

Ligandos mono y polidésmicos

Los ligandos de proteínas también se pueden caracterizar por el número de cadenas de proteínas a las que se unen. Los ligandos "monodésmicos" (μόνος: único, δεσμός: unión) son ligandos que se unen a una sola cadena proteica, mientras que los ligandos "polidésmicos" (πολοί: muchos) [ 31] son ​​frecuentes en complejos proteicos y son ligandos que se unen a más de una proteína. cadena, típicamente en o cerca de las interfaces de proteínas. Investigaciones recientes muestran que el tipo de ligandos y la estructura del sitio de unión tienen profundas consecuencias para la evolución, función, alosterio y plegamiento de los complejos proteicos. [32] [33]

Andamio privilegiado

Un andamio privilegiado [34] es una estructura molecular o fracción química que es estadísticamente recurrente entre fármacos conocidos o entre una serie específica de compuestos biológicamente activos. Estos elementos privilegiados [35] pueden utilizarse como base para diseñar nuevos compuestos biológicos activos o bibliotecas de compuestos.

Métodos utilizados para estudiar la unión.

Los principales métodos para estudiar las interacciones proteína-ligando son las principales técnicas hidrodinámicas y calorimétricas, y los principales métodos espectroscópicos y estructurales, como

Otras técnicas incluyen: intensidad de fluorescencia, complementación de fluorescencia bimolecular, FRET (transferencia de energía por resonancia fluorescente) / resonancia de plasmón superficial de extinción FRET, interferometría de biocapa , ELISA indirecto de coinmunopreciptación, diálisis de equilibrio, electroforesis en gel, transferencia Far Western, anisotropía de polarización de fluorescencia, paramagnética electrónica. resonancia, termoforesis a microescala , switchSENSE .

El espectacular aumento de la potencia informática de los superordenadores y los ordenadores personales ha hecho posible estudiar las interacciones proteína-ligando también mediante química computacional . Por ejemplo, en el proyecto grid.org , que finalizó en abril de 2007 , se aprovechó una red mundial de más de un millón de PC comunes y corrientes para la investigación del cáncer. A Grid.org le han sucedido proyectos similares como World Community Grid , Human Proteome Folding Project. , Compute Against Cancer y Folding@Home .

Ver también

Referencias

  1. ^ Teif VB (octubre de 2005). "Condensación de ADN inducida por ligando: elección del modelo". Revista Biofísica . 89 (4): 2574–2587. Código Bib : 2005BpJ....89.2574T. doi : 10.1529/biophysj.105.063909. PMC  1366757 . PMID  16085765.
  2. ^ Teif VB, Rippe K (octubre de 2010). "Modelos de red estadístico-mecánico para la unión proteína-ADN en cromatina". Revista de Física: Materia Condensada . 22 (41): 414105. arXiv : 1004.5514 . Código Bib : 2010JPCM...22O4105T. doi :10.1088/0953-8984/22/41/414105. PMID  21386588. S2CID  103345.
  3. ^ Baron R, Setny P, McCammon JA (septiembre de 2010). "Reconocimiento de agua en cavidad-ligando". Revista de la Sociedad Química Estadounidense . 132 (34): 12091–12097. doi :10.1021/ja1050082. PMC 2933114 . PMID  20695475. 
  4. ^ "La diferencia entre los valores Ki, Kd, ​​IC50 y EC50". El caracol de la ciencia . 31 de diciembre de 2019.
  5. ^ Ver Curvas de vinculación competitiva homólogas Archivado el 19 de diciembre de 2007 en Wayback Machine , Una guía completa para la regresión no lineal, curvefit.com.
  6. ^ Baaske P, Wienken CJ, Reineck P, Duhr S, Braun D (marzo de 2010). "Termoforesis óptica para cuantificar la dependencia del tampón de la unión de aptámeros". Angewandte Chemie . 49 (12): 2238–2241. doi :10.1002/anie.200903998. PMID  20186894.
    • "Un camino caliente hacia nuevos fármacos". Phys.org . 24 de febrero de 2010.
  7. ^ Wienken CJ, Baaske P, Rothbauer U, Braun D, ​​Duhr S (octubre de 2010). "Ensayos de unión a proteínas en líquidos biológicos mediante termoforesis a microescala". Comunicaciones de la naturaleza . 1 (7): 100. Código Bib : 2010NatCo...1..100W. doi : 10.1038/ncomms1093 . PMID  20981028.
  8. ^ Vu-Quoc, L., Configuración integral (mecánica estadística), 2008. este sitio wiki no funciona; consulte este artículo en el archivo web del 28 de abril de 2012.
  9. ^ ab https://doi.org/10.1093/oso/9780199744121.003.0006
  10. ^ Kenakin TP (2006). Una introducción a la farmacología: teoría, aplicaciones y métodos. Prensa académica . pag. 79.ISBN _ 978-0-12-370599-0.
  11. ^ Cabañas, C; Wang, M; Han, X; Hansen, SB (8 de agosto de 2017). "Un ensayo de unión fluorescente soluble revela el antagonismo PIP2 de los canales TREK-1". Informes celulares . 20 (6): 1287–1294. doi :10.1016/j.celrep.2017.07.034. PMC 5586213 . PMID  28793254. 
  12. ^ ab Matera, Carlo; Pucci, Luca; Fiorentini, Chiara; Fucil, Sergio; Missale, Cristina; Grazioso, Giovanni; Clementi, Francisco; Zoli, Michele; De Amici, Marco (28 de agosto de 2015). "Compuestos bifuncionales dirigidos a receptores de nACh tanto D2 como no α7: diseño, síntesis y caracterización farmacológica". Revista europea de química medicinal . 101 : 367–383. doi :10.1016/j.ejmech.2015.06.039. PMID  26164842.
  13. ^ Matera, Carlo; Flamini, Lisa; Quadri, Marta; Vivo, Valentina; Ballabeni, Vigilio; Holzgrabe, Ulrike; Mohr, Klaus; De Amici, Marco; Barocelli, Elisabetta (21 de marzo de 2014). "Agonistas de tipo bis (amonio) alcano de los receptores muscarínicos de acetilcolina: síntesis, caracterización funcional in vitro y evaluación in vivo de su actividad analgésica". Revista europea de química medicinal . 75 : 222–232. doi :10.1016/j.ejmech.2014.01.032. PMID  24534538.
  14. ^ Erez M, Takemori AE, Portoghese PS (julio de 1982). "Potencia antagonista narcótica de ligandos bivalentes que contienen beta-naltrexamina. Evidencia de formación de puentes entre sitios de reconocimiento proximales". Revista de Química Medicinal . 25 (7): 847–849. doi :10.1021/jm00349a016. PMID  7108900.
  15. ^ Portoghese PS, Ronsisvalle G, Larson DL, Yim CB, Sayre LM, Takemori AE (1982). "Ligandos bivalentes agonistas y antagonistas de opioides como sondas receptoras". Ciencias de la vida . 31 (12-13): 1283-1286. doi :10.1016/0024-3205(82)90362-9. PMID  6292615.
  16. ^ Portoghese PS, Akgün E, Lunzer MM (enero de 2017). "Inducción de heterómeros: ¿un enfoque de farmacología única?". ACS Neurociencia Química . 8 (3): 426–428. doi : 10.1021/acschemneuro.7b00002 . PMID  28139906.
  17. ^ Blum JJ, Conn PM (diciembre de 1982). "Estimulación de la liberación de hormona luteinizante por hormona liberadora de gonadotropina: un modelo ligando-receptor-efector". Actas de la Academia Nacional de Ciencias de los Estados Unidos de América . 79 (23): 7307–7311. Código bibliográfico : 1982PNAS...79.7307B. doi : 10.1073/pnas.79.23.7307 . JSTOR  13076. PMC 347328 . PMID  6296828. 
  18. ^ Conn PM, Rogers DC, Stewart JM, Niedel J, Sheffield T (abril de 1982). "Conversión de un antagonista de la hormona liberadora de gonadotropina en un agonista". Naturaleza . 296 (5858): 653–655. Código Bib :1982Natur.296..653C. doi :10.1038/296653a0. PMID  6280058. S2CID  4303982.
  19. ^ Nimczick M, Pemp D, Darras FH, Chen X, Heilmann J, Decker M (agosto de 2014). "La síntesis y evaluación biológica de ligandos de receptores de cannabinoides bivalentes basados ​​en benzimidazoles selectivos de hCB₂R revelan propiedades intrínsecas inesperadas". Química bioorgánica y medicinal . 22 (15): 3938–3946. doi :10.1016/j.bmc.2014.06.008. PMID  24984935.
  20. ^ Russo O, Berthouze M, Giner M, Soulier JL, Rivail L, Sicsic S, Lezoualc'h F, Jockers R, Berque-Bestel I (septiembre de 2007). "Síntesis de sondas bivalentes específicas que interactúan funcionalmente con los dímeros del receptor 5-HT (4)". Revista de Química Medicinal . 50 (18): 4482–4492. doi :10.1021/jm070552t. PMID  17676726.
  21. ^ Soulier JL, Russo O, Giner M, Rivail L, Berthouze M, Ongeri S, Maigret B, Fischmeister R, Lezoualc'h F, Sicsic S, Berque-Bestel I (octubre de 2005). "Diseño y síntesis de sondas específicas para estudios de dimerización del receptor 5-HT4 humano" (PDF) . Revista de Química Medicinal . 48 (20): 6220–6228. doi :10.1021/jm050234z. PMID  16190749.
  22. ^ ab Busnelli M, Kleinau G, Muttenthaler M, Stoev S, Manning M, Bibic L, Howell LA, McCormick PJ, Di Lascio S, Braida D, Sala M, Rovati GE, Bellini T, Chini B (agosto de 2016). "Diseño y caracterización de ligandos bivalentes superpotentes dirigidos a dímeros del receptor de oxitocina a través de una estructura similar a un canal". Revista de Química Medicinal . 59 (15): 7152–7166. doi : 10.1021/acs.jmedchem.6b00564 . hdl : 2434/430357 . PMID  27420737.
  23. ^ Lensing CJ, Adank DN, Wilber SL, Freeman KT, Schnell SM, Speth RC, Zarth AT, Haskell-Luevano C (febrero de 2017). "Una comparación directa in vivo del agonista monovalente de melanocortina Ac-His-DPhe-Arg-Trp-NH2 versus el agonista bivalente Ac-His-DPhe-Arg-Trp-PEDG20-His-DPhe-Arg-Trp-NH2: un bivalente Ventaja". ACS Neurociencia Química . 8 (6): 1262-1278. doi :10.1021/acschemneuro.6b00399. PMC 5679024 . PMID  28128928. 
  24. ^ ab Xu L, Josan JS, Vagner J, Caplan MR, Hruby VJ, Mash EA, Lynch RM, Morse DL, Gillies RJ (diciembre de 2012). "Los ligandos heterobivalentes se dirigen a combinaciones de receptores de la superficie celular in vivo". Actas de la Academia Nacional de Ciencias de los Estados Unidos de América . 109 (52): 21295–21300. Código Bib : 2012PNAS..10921295X. doi : 10.1073/pnas.1211762109 . JSTOR  42553664. PMC 3535626 . PMID  23236171. 
  25. ^ ab Lensing CJ, Freeman KT, Schnell SM, Adank DN, Speth RC, Haskell-Luevano C (abril de 2016). "Una investigación in vitro e in vivo de ligandos bivalentes que muestran actividad funcional y de unión preferencial para diferentes homodímeros del receptor de melanocortina". Revista de Química Medicinal . 59 (7): 3112–3128. doi : 10.1021/acs.jmedchem.5b01894. PMC 5679017 . PMID  26959173. 
  26. ^ Shonberg J, Scammells PJ, Capuano B (junio de 2011). "Diseño de estrategias para ligandos bivalentes dirigidos a GPCR". ChemMedChem . 6 (6): 963–974. doi :10.1002/cmdc.201100101. PMID  21520422. S2CID  10561038.
  27. ^ Berque-Bestel I, Lezoualc'h F, Jockers R (diciembre de 2008). "Ligandos bivalentes como herramientas farmacológicas específicas para dímeros de receptores acoplados a proteína G". Tecnologías actuales de descubrimiento de fármacos . 5 (4): 312–318. doi :10.2174/157016308786733591. PMID  19075611.
  28. ^ Akgün E, Javed MI, Lunzer MM, Powers MD, Sham YY, Watanabe Y, Portoghese PS (noviembre de 2015). "Inhibición del dolor inflamatorio y neuropático dirigiéndose a un heterómero del receptor opioide Mu / receptor de quimiocinas 5 (MOR-CCR5)". Revista de Química Medicinal . 58 (21): 8647–8657. doi : 10.1021/acs.jmedchem.5b01245. PMC 5055304 . PMID  26451468. 
  29. ^ Daniels DJ, Lenard NR, Etienne CL, Law PY, Roerig SC, Portoghese PS (diciembre de 2005). "La tolerancia y dependencia inducida por opioides en ratones está modulada por la distancia entre farmacóforos en una serie de ligandos bivalentes". Actas de la Academia Nacional de Ciencias de los Estados Unidos de América . 102 (52): 19208–19213. Código Bib : 2005PNAS..10219208D. doi : 10.1073/pnas.0506627102 . JSTOR  4152590. PMC 1323165 . PMID  16365317. 
  30. ^ Smeester BA, Lunzer MM, Akgün E, Beitz AJ, Portoghese PS (noviembre de 2014). "Apuntar a los supuestos heterómeros del receptor 5 de glutamato opioide / metabotrópico mu produce una potente antinocicepción en un modelo de cáncer óseo murino crónico". Revista europea de farmacología . 743 : 48–52. doi :10.1016/j.ejphar.2014.09.008. PMC 4259840 . PMID  25239072. 
  31. ^ Abrusan G, Marsh JA (2019). "La estructura del sitio de unión del ligando da forma al plegamiento, ensamblaje y degradación de complejos de proteínas homoméricas". Revista de biología molecular . 431 (19): 3871–3888. doi :10.1016/j.jmb.2019.07.014. PMC 6739599 . PMID  31306664. 
  32. ^ Abrusan G, Marsh JA (2018). "La estructura del sitio de unión del ligando influye en la evolución de la función y la topología del complejo proteico". Informes celulares . 22 (12): 3265–3276. doi :10.1016/j.celrep.2018.02.085. PMC 5873459 . PMID  29562182. 
  33. ^ Abrusan G, Marsh JA (2019). "La estructura del sitio de unión del ligando da forma a la transducción de señales alostéricas y la evolución de la alosteria en complejos proteicos". Biología Molecular y Evolución . 36 (8): 1711-1727. doi :10.1093/molbev/msz093. PMC 6657754 . PMID  31004156. 
  34. ^ Welsch ME, Snyder SA, Stockwell BR (junio de 2010). "Andamios privilegiados para el diseño de bibliotecas y el descubrimiento de fármacos". Opinión actual en biología química . 14 (3): 347–61. doi :10.1016/j.cbpa.2010.02.018. PMC 2908274 . PMID  20303320. 
  35. ^ Kombarov R, Altieri A, Genis D, Kirpichenok M, Kochubey V, Rakitina N, Titarenko Z (febrero de 2010). "BioCores: identificación de un motivo estructural privilegiado basado en fármaco/producto natural para el descubrimiento de plomo de molécula pequeña". Diversidad Molecular . 14 (1): 193–200. doi :10.1007/s11030-009-9157-5. PMID  19468851. S2CID  23331540.

enlaces externos