Todos os reagentes eram de grau analítico e usados sem purificação adicional. O progresso da reação foi monitorado por cromatografia de camada fina (TLC) em gel de sílica usando n-hexano/acetato de etila (7:3–6:4, v/v), com pontos visualizados sob luz ultravioleta. Salvo indicação em contrário, todas as reações foram realizadas sob atmosfera ambiente. Os compostos sintetizados foram caracterizados por espectroscopia de transformada de Fourier no infravermelho (FTIR), RMN de 1H e RMN de 13C. Os espectros de RMN de 1H e 13C foram registrados em CDCl₃ a 600 MHz e 150 MHz, respectivamente, e deslocamentos químicos foram relatados em ppm em relação ao sinal residual do solvente. A via sintética é representada na Figura 1. Os dados completos de caracterização fisicoquímica para todos os compostos sintetizados são fornecidos no Arquivo Suplementar 1.
Síntese do benzoato de etil (1)
Ácido benzóico (50 g, 0,41 mol) foi colocado em um frasco de fundo arredondado de 500 mL, e etanol absoluto (150 mL) foi adicionado. Ácido sulfúrico concentrado (6 mL) foi adicionado gota a gota como catalisador sob agitação magnética contínua. A mistura de reação era aquecida sob refluxo a 78 °C por 4–6 horas usando um condensador resfriado a água. O progresso da reação foi monitorado pela TLC. Após a conclusão, a mistura era resfriada até a temperatura ambiente e neutralizada com carbonato de sódio aquoso saturado até que o pH atingisse 7–8. A mistura de reação foi transferida para um funil separador, e a camada orgânica foi coletada, lavada com água destilada e seca sobre sulfato de sódio anidro. Após a filtração, o benzoato de etil (1) foi usado diretamente na etapa seguinte, sem purificação adicional (Figura 1).
Síntese de benzohidrazida (2)
O benzoato de etílico (1; 0,41 mol) foi dissolvido em etanol absoluto (100 mL), e hidrazina hidratada (80%, 25 mL, excesso) foi adicionada sob agitação contínua. A mistura de reação era aquecida sob refluxo a 78 °C por 4–5 horas em ar usando um condensador resfriado a água. O progresso da reação foi monitorado pela TLC. Após a conclusão, o precipitado cristalino branco da benzohidrazida (2) foi coletado por filtração, lavado com etanol frio, secado até peso constante em temperatura ambiente e usado em etapas sintéticas subsequentes (Figura 1).
Síntese de 5-fenil-1,3,4-oxadiazol-2-tiol(3)
Benzohidrazida (2; 10 g, 0,073 mol) e hidróxido de potássio (4,12 g, 0,073 mol) foram colocados em um frasco de fundo arredondado de 250 mL contendo etanol (30 mL). A mistura foi aquecida sob refluxo a 78 °C por 10 minutos usando um condensador refrigerado a água. O dissulfeto de carbono (CS₂; 5 mL, 0,083 mol) foi então adicionado, e o refluxo foi continuado por mais 3–4 horas. O progresso da reação foi monitorado pela TLC. Após a conclusão, a mistura de reação era transferida para água destilada e aquecida com agitação contínua. Ácido clorídrico diluído foi adicionado gota a gota até que o pH atingisse 5, resultando na precipitação do composto (3) (Figura 1). O precipitado era coletado por filtração, lavado cuidadosamente com água destilada fria, seco até peso constante e pesado para determinar o rendimento. O composto foi caracterizado por espectroscopia FTIR, RMN de 1H e RMN de 13C.
Síntese de 4,5-difenil-4H-1,2,4-triazola-3-tiol(4)
Benzidraida (2; 5 g, 0,037 mol) e fenilisotiocianato (4 mL, 0,037 mol) foram colocados em um frasco de fundo arredondado de 100 mL equipado com um condensador refrigerado a água. A mistura de reação foi aquecida sob refluxo a 78 °C até que a formação do intermediário tiosemicarbazida estivesse completa, conforme confirmado pelo TLC usando n-hexano/acetato de etila (7:3, v/v). Hidróxido de potássio (2 g, 0,036 mol), dissolvido em etanol (20 mL), foi então adicionado, e o refluxo continuou até a formação do composto (4), conforme confirmado pela TLC. A mistura de reação foi transferida para água destilada, aquecida com agitação contínua e acidificada com ácido clorídrico diluído até pH 5 para precipitar o composto (4) (Figura 1). O precipitado foi coletado por filtração, lavado cuidadosamente com água destilada fria, seco até peso constante e caracterizado por espectroscopia FTIR, RMN de 1H e RMN de 13°C.
Síntese de intermediários de bromoacetamida 7(a–f)
Cada derivado substituído da anilina 5(a–f) (1 equiv.) foi dissolvido em clorofórmio e resfriado em um banho de gelo para manter a temperatura da reação entre 0 e 5 °C. Bromocetila brometo (6, ≥98%, 1 equivalente) foi adicionado gota em gota com mexer contínuo. Durante a reação, foi adicionada uma solução aquosa de carbonato de sódio a 18% conforme necessário para manter o pH da reação entre 9 e 10, neutralizando o ácido bromídico gerado. Após a adição completa do brometo de bromoacetil, a mistura de reação era deixada aquecer até a temperatura ambiente e misturada até que a TLC confirmasse o consumo completo do material inicial. Os intermediários resultantes de bromoacetamida 7(a–f) foram coletados por filtração, lavados com água destilada fria, secos até peso constante e usados diretamente na etapa sintética subsequente, sem purificação adicional (Figura 1).
Síntese de derivados de acetamida à base de triazol 8(a–c)
Derivado de triazol (4; 0,5 g, 0,002 mol) foi dissolvido em N, N-dimetilformamida (10 mL) em um frasco de fundo arredondado de 50 mL. O intermediário apropriado de bromoacetamida (7a, 7b ou 7c; 0,002 mol, 1 equiv.) foi adicionado com agitação contínua, seguido por hidreto de sódio (0,05 g, 0,002 mol) para promover a alquilação S. A mistura de reação foi agitada em temperatura ambiente por 24–48 horas, com o progresso da reação monitorado pelo TLC. Ao concluir, a mistura de reação era despejada em água gelada para induzir precipitação. O precipitado foi coletado por filtração, lavado cuidadosamente com água fria, seco e caracterizado por espectroscopia FTIR, RMN de 1H e RMN de 13°C (Figura 1).
Síntese de derivados de acetamida à base de oxadiazol 9(d–f)
Derivado de oxadiazol (3; 0,5 g, 0,002 mol) foi dissolvido em N, N-dimetilformamida (10 mL) em um frasco de fundo arredondado de 50 mL. O intermediário correspondente de bromoacetamida, 7d, 7e ou 7f (0,002 mol, 1 equiv.), foi adicionado, seguido pelo hidreto de sódio (0,05 g, 0,002 mol) para facilitar a alquilação S. A mistura de reação foi agitada em temperatura ambiente por 24–48 horas, com progresso monitorado pelo TLC. Ao finalizar, a mistura era despejada em água gelada para precipitar o produto. O precipitado foi coletado por filtração, lavado cuidadosamente com água fria, seco e caracterizado por espectroscopia FTIR, RMN de 1H e RMN de 13°C (Figura 1).
Ensaio de inibição de α-glucosidase
A atividade inibitória da α-glicosidase foi avaliada misturando 50 μL da solução amostrada de teste com 100 μL de solução de α-glucosidase (0,35 U/mL) preparada em tampão fosfato (pH 6,8). A mistura de reação foi incubada a 37 °C por 10 minutos, depois 100 μL de 1,5 mM de p-nitrofenil-α-D-glipirospinósido foram adicionados como substrato. Após nova incubação a 37 °C por 20 minutos, a reação foi encerrada com a adição de 1 mL de carbonato de sódio (Na₂CO₃). A absorvância foi medida em 400nm 14. Acarbose foi usado como controle positivo. A mistura de reação sem o composto de teste servia como controle negativo, enquanto a mistura de reação sem a enzima era usada como blank.
Todos os experimentos foram realizados em triplicado, e os resultados são expressos como a média ± desvio padrão (DS). A porcentagem de inibição da atividade α-glucosidase foi calculada usando a seguinte equação(1)14:
(1)
onde As representa a absorvância da solução da amostra de teste, Ab é a absorvância do reagente em branco (sem α-glucosidase), e A0 é a absorvância do controle negativo (sem a amostra). Os valores do IC50 foram determinados testando concentrações compostas variando de 3,125 a 100 μM. A análise de regressão não linear foi realizada usando software padrão de gráficos.
Ensaio de inibição da urease
A atividade inibitória da urease foi determinada usando um método modificado de Berthelot. A mistura de reação consistia em 300 μL da solução amostrada de teste, 200 μL de solução de urease de feijão jack e 200 μL de 25 mM de ureia. A mistura foi incubada em banho-maria a 37 °C por 30 minutos. Após a incubação, foram adicionados 500 μL de uma solução reveladora de cor contendo 10 g/L de fenol e 50 g/L de nitroprussídeo de sódio, seguidos por 500 μL de uma segunda solução contendo 5% (v/v) de hipoclorito de sódio e 5 g/L de hidróxido de sódio. A mistura de reação foi então incubada por mais 15–30 minutos para permitir o desenvolvimento da cor.
A absorção da solução foi medida em 625 nm contra um reagenteem branco 15. A tioureia foi usada como controle positivo, a mistura de reação sem o composto de teste como controle negativo, e a mistura sem a enzima urease como branco. Todos os experimentos foram realizados em triplicado (n = 3), e os resultados são expressos como valores médios ± desvio padrão (S.D.). A porcentagem de inibição foi calculada usando a seguinte fórmula matemática (2)15:
(2)
onde A0 é a absorvância do controle (sem inibidor) e As é a absorvância da amostra de teste (com inibidor). Gráficos e análise de regressão não linear para determinar os valores IC50 foram realizados usando software padrão de plotagem de dados.
Ensaio de ligação BSA
As interações de ligação entre a albumina sérica bovina (BSA) e os compostos sintetizados foram investigadas usando espectroscopia de extinção por fluorescência. Uma solução de estoque de BSA de 10 μM foi preparada em tampão de fosfato de 20 mM (pH 7,4). Para cada experimento de titulação, 1,0 mL da solução de BSA foi misturado com 2,0 mL de tampão de fosfato em uma cuveta de quartzo. A solução foi ajustada por adições sucessivas de alíquotas de 5 μL do composto de teste (1,5 mM em metanol). Espectros de emissão de fluorescência foram registrados na faixa de comprimentos de onda de 300–400 nm após excitação em 295 nm, com a emissão máxima observada em aproximadamente 336 nm. As medições de fluorescência foram realizadas usando um espectrofotômetro de fluorescência equipado com software padrão de aquisição de fluorescência, com larguras de fendas de excitação e emissão definidas para 10 nm e 2,5 nm, respectivamente. Para determinar os parâmetros termodinâmicos da ligação, incluindo a energia livre de Gibbs (ΔG), entalpia (ΔH) e entropia (ΔS), os experimentos de titulação foram repetidos em 298, 308 e 313 K.
Para investigar o mecanismo de têmpera por fluorescência, os dados experimentais foram analisados usando a equação de Stern–Volmer (Equação 3)16:
(3)
onde F, 0 e F representam, respectivamente, as intensidades de fluorescência da BSA na ausência e presença do quencher; Kq é a constante da taxa de têmpera bimolecular; τ0 é a vida média de vida de BSA na ausência de um quencher; KSV é a constante de têmpera de Stern-Volmer; e [Q] é a concentração do quencher. A constante de ligação (Kb) e o número de sítios de ligação (n) foram posteriormente determinados ao plotar o logaritmo duplo de acordo com a Equação (4)16:
(4)
Além disso, a variação padrão de energia livre de Gibbs (ΔG°), a mudança padrão da entalpia (ΔH°) e a mudança padrão de entropia (ΔS°) que regem a interação de ligação foram calculadas usando as relações de Van't Hoff e Gibbs-Helmholtz expressas nas Equações (5) e (6)17:
ΔG∘ = -RT Ln Kb (5)
(6)
onde R é a constante universal do gás e T é a temperatura absoluta em Kelvin.
Acoplamento molecular
Todas as simulações de modelagem molecular e acoplamento foram realizadas usando a Schrödinger Suite (interface Maestro), com acoplamento feito usando o módulo Glide. As estruturas cristalinas da α-glucosidase (PDB ID: 5NN8) e da urease (PDB ID: 4H9M) foram recuperadas do Banco de Dados de Proteínas RCSB. As estruturas proteicas foram preparadas usando o Assistente de Preparação de Proteínas, removendo ligandos cocristalizados, moléculas de água e heteroátomos não essenciais, adicionando resíduos e laços ausentes, atribuindo ordens corretas de ligações e adicionando átomos de hidrogênio polar. Os estados de protonação e tautoméricos foram atribuídos no pH fisiológico (7,0 ± 2,0) usando o módulo Epik. As estruturas proteicas preparadas foram subsequentemente submetidas a minimização de energia contida com o campo de força OPLS para aliviar os confrontos estéricos.
As estruturas bidimensionais dos compostos 3, 4, 8(a–c) e 9(d–f), juntamente com os inibidores de referência acarbose e tioureia, foram preparadas usando o módulo LigPrep. Ordens de ligação, estereoquímica e estados de ionização no pH fisiológico foram otimizados, e as conformações tridimensionais resultantes dos ligantes foram minimizadas em energia usando o campo de força OPLS.
Grades receptoras foram geradas usando a ferramenta Geração de Grade de Receptores, com as caixas de grade centralizadas nos ligantes nativos cocristalizados para definir os sítios ativos. As dimensões da grade foram definidas para 25 × 25 × 25 Å para ambas as proteínas. A grade da urease estava centrada em x = −38,76, y = −44,14 e z = 66,88, enquanto a grade da α-glicosidase estava centrada em x = 21,14, y = −7,56 e z = 24,37.
Os cálculos de acoplamento eram realizados usando o modo de extra-precisão (XP) do Glide com amostragem flexível de ligandos. Múltiplas poses de ligação foram geradas para cada ligante, e a pose com a menor pontuação de acoplamento Glide foi selecionada para análise posterior. O protocolo de acoplamento foi validado reacoplando os ligantes cocristalizados em suas respectivas grades receptoras. Valores de desvio raiz-médio-quadrático (RMSD) de ≤2,0 Å confirmaram a confiabilidade e reprodutibilidade do protocolo de acoplamento.
As conformações de ligação dos complexos acoplados foram visualizadas usando o PyMOL (versão 3.1)18. Interações proteína–ligante, incluindo ligações de hidrogênio, contatos hidrofóbicos e interações eletrostáticas, foram analisadas usando o Visualizador19 do Discovery Studio.
Análise da teoria do funcional da densidade (DFT)
Todos os cálculos da teoria da função da densidade (DFT) foram realizados usando a Gaussiana 09, e as geometrias moleculares otimizadas e os orbitais moleculares de fronteira foram visualizados usando o GaussView 5. Os derivados sintetizados à base de acetamida à base de triazol e oxadiazol, 8(a–c) e 9(d–f), foram submetidos à otimização geométrica na fase aquosa usando o funcional híbrido de três parâmetros Lee–Yang–Parr (B3LYP) de Becke, em conjunto com o conjunto de base 6-311+G(d,p) 20. Os efeitos dos solventes foram incorporados usando o Modelo de Contínuo Polarizável Semelhante ao Condutor (CPCM)21. Cálculos de frequência vibracional foram posteriormente realizados no mesmo nível teórico para confirmar que todas as estruturas otimizadas correspondiam aos verdadeiros mínimos locais na superfície de energia potencial, como indicado pela ausência de frequências imaginárias.
Para entender as transições eletrônicas e a estabilidade cinética dos compostos sintetizados, foram calculadas as energias do orbital molecular mais ocupado (EHOMO) e do orbital molecular desocupado mais baixo (ELUMO). A diferença de energia orbital molecular fronteira (ΔE = ELUMO - EHOMO) FOI DETERMINADA PARA CADA DERIVADA PARA AVALIAR SUA REATIVIDADE QUÍMICA E ESTABILIDADE CINÉTICA, COM UMA DIFERENÇA DE ENERGIA MENOR INDICANDO MAIOR REATIVIDADE.
Além disso, descritores globais de reatividade química — incluindo dureza química (η), suavidade (S), potencial químico (μ), eletronegatividade (X) e índice de eletrofilicidade (w) — foram calculados a partir das energias orbitais de fronteira usando as seguintes equações operacionaispadrão 22:
(7)
(8)
(9)
(10)
As propriedades eletrônicas computadas e os descritores globais foram sistematicamente comparados entre as séries para estabelecer relações estrutura-atividade (SARs) que refletam seus perfis observados de α-glicosidase e inibidores da ureasa.
Além disso, mapas de potencial eletrostático molecular (MESP) foram gerados e visualizados usando VESTA para mapear a distribuição de carga entre as estruturas moleculares otimizadas. As superfícies MESP permitiram a identificação de regiões eletrofílicas (deficientes em elétrons, preferidas por ataque nucleofílico) e nucleofílicas (ricas em elétrons, favorecidas por ataques eletrofílicos). Essa distribuição de potencial eletrostático foi, em última análise, correlacionada com as interações específicas de ligação não covalente, como ligações de hidrogênio e contatos eletrostáticos, observadas nas simulações de acoplamento molecular.
Análise do ADME
Para avaliar a semelhança com medicamentos e as propriedades biofarmacêuticas dos derivados sintetizados à base de triazol e oxadiazol da acetamida 8(a–c) e 9(d–f), foi realizada uma análise in silico de absorção, distribuição, metabolismo e excreção (ADME). As representações simplificadas do sistema de entrada molecular de entrada de linha (SMILES) de todos os compostos sintetizados, juntamente com os inibidores de referência acarbose e tioureia, foram geradas usando um pacote de software de modelagem molecular e submetidas ao servidor webSwissADME 23.
Os parâmetros físico-químicos calculados incluíam peso molecular (MW), lipofilicidade (logP), área superficial polar topológica (TPSA), número de doadores de ligações de hidrogênio (HBD), número de aceitadores de ligações de hidrogênio (HBA) e número de ligações rotativas (RB).
As propriedades farmacocinéticas foram avaliadas prevendo a absorção gastrointestinal (GI) e a permeabilidade da barreira hematoencefálica (BBB). Além disso, os compostos foram avaliados quanto à conformidade com a regra dos cinco de Lipinski e outros critérios estabelecidos de semelhança com medicamentos para estimar sua biodisponibilidade oral. Por fim, a pontuação de acessibilidade sintética (SA) foi calculada para cada composto em uma escala que varia de 1 (fácil de sintetizar) a 10 (muito difícil de sintetizar). Os perfis previstos de ADME e semelhança de fármacos dos compostos sintetizados foram comparados com os dos inibidores de referência para avaliar seu potencial como candidatos principais ao desenvolvimento de fármacos.