Molecular Docking Assessment of Environmental Pesticides as Potential Endocrine Disruptors of Human Estrogen Receptor Alpha

기술 요약: 인간 에스트로겐 수용체 알파(ERα)의 잠재적 내분비 교란 물질로서 환경 농약에 대한 분자 도킹 평가

문제 제기
환경 농약은 호르몬 작용을 모방하거나 차단할 수 있는 내분비 교란 물질(EDCs)로서 인체와 동물 건강에 상당한 위험을 초래한다는 점이 점점 더 널리 인식되고 있다. 신경독성, 간독성, 생식 독성을 포함하여 농약의 급성 및 만성 독성에 관한 광범위한 문헌이 존재하지만, 이러한 화합물이 특히 에스트로겐 수용체 알파(ERα)를 방해함으로써 생식 건강을 저해하는 구체적인 분자적 메커니즘에 대한 포괄적인 이해에는 여전히 공백이 존재한다. 본 연구는 다섯 가지 널리 사용되는 농약(DDT, 클로르피리포스, 엔도설판, 아트라진, 글리포세이트트)의 결합 친화도와 상호작용 프로파일을 인간 ERα에 대해 평가하여, 이들이 에스트로겐성 교란 물질로서 작용할 가능성을 평가하고자 한다.

연구 방법론
본 연구는 분자 도킹 및 상호작용 프로파일링을 활용한 in silico 계산적 접근 방식을 채택하였다:

• 단백질 준비: 인간 ERα의 리간드 결합 도메인(PDB ID: 3ERT)의 3D 구조를 Protein Data Bank에서 추출하였다. 도킹 정확도를 보장하기 위해 고해상도이며 천연 리간드와 공동 결정화된 구조를 선택하였다.
• 리간드 준비: 다섯 가지 농약(아트라진, DDT, 클로르피리포스, 엔도설판, 글리포세이트)과 천연 리간드인 에스트라디올(대조군으로 사용)을 선정하였다. 3D 컨포머(Conformer)는 PubChem에서 .sdf 형식으로 다운로드하여 검증하였다.
• 도킹 프로토콜: 시뮬레이션은 포즈 예측 및 결합 에너지 계산을 위해 Auto Dock Vina를 사용하는 CB-Dock 웹 서버를 사용하여 수행되었다. 해당 서버는 잠재적인 결합 공동(Cavity)을 자동으로 감지하며, 상위 5개의 예측된 포켓에 집중하였다. 에스트라디올을 먼저 도킹하여 천연 결합 부위를 검증하였다.
• 상호작용 분석: 단백질-리간드 상호작용 프로파일러(PLIP) 서버를 사용하여 수소 결합, 소수성 접촉, 염다리(Salt bridge), π–카티온(π–cation) 상호작용을 포함한 구체적인 분자 상호작용을 분석하였다.
• 한계점: 저자들은 서버 제한으로 인해 전체 규모의 분자 역학(MD) 시뮬레이션 및 종합적인 ADMET(흡수, 분포, 대사, 배설 및 독성) 예측을 수행하지 못했음을 명시하였다. 대신, 본 연구는 ADMET 파라미터에 대해 기존 문헌 데이터를 활용하였다.

주요 결과
테스트된 모든 리간드는 ERα의 천연 호르몬 결합 부위에 해당하는 "Pocket 1"에 우선적으로 도킹되는 것으로 나타났다. 농약들의 결합 친화도와 상호작용 프로파일은 다음과 같이 유의미하게 달랐다:

• DDT: 농약 중 가장 강력한 결합 친화도(-8.5 kcal/mol)를 보였으며, 이는 대조군인 에스트라디올(-9.3 kcal/mol)에 근접한 수치이다. DDT는 비극성 잔기(LEU346, ALA350, TRP383, LEU387, PHE404, LEU525)와 광범위한 소수성 접촉(6개 접촉)을 형성하였으나 수소 결합은 형성하지 않았다.
• 클로르피리포스(Chlorpyrifos): 주로 잔기 LEU346, LEU391, ILE424와의 소수성 접촉(3개 접촉)에 의해 구동되는 중간 정도의 결합(-6.3 kcal/mol)을 보였다.
• 엔도설판(Endosulfan): 3개의 수소 결합, 2개의 소수성 접촉, 1개의 π–카티온 상호작용을 포함한 혼합 모드를 통해 중간 정도의 결합(-6.0 kcal/mol)을 나타냈다.
• 아트라진(Atrazine): 4개의 수소 결합과 2개의 소수성 접촉을 특징으로 하는 혼합 결합 모드와 함께 중간 정도의 결합(-5.9 kcal/mol)을 보였다.
• 글리포세이트(Glyphosate): 소수성 상호작용 없이 가장 약한 결합(-4.7 kcal/mol)을 나타냈다. 안정화는 LYS449와의 3개 수소 결합 및 1개의 염다리를 포함한 극성 상호작용을 통해 주로 이루어졌다.

주요 기여

• 비교 결합 평가: 본 연구는 다섯 가지 서로 다른 환경 농약을 ERα 결합 포켓에 대해 직접적인 계산적 비교를 제공하여, 결합 에너지를 기반으로 한 잠재적 내분비 교란의 계층 구조를 확립하였다.
• 상호작용 프로파일링: PLIP을 활용하여 각 농약의 결합을 주도하는 구체적인 분자력(소수성 vs 극성)을 상세히 기술함으로써, DDT가 글리포세이트와 같은 약한 결합체에 비해 어떻게 천연 에스트로겐 신호 전달을 경쟁적으로 억제할 수 있는지에 대한 메커니즘적 통찰을 제공하였다.
• 스크리닝 유용성: 본 연구는 환경 화학물질 중 잠재적 EDC를 식별하기 위한 예측 스크리닝 도구로서 분자 도킹(특히 CB-Dock)의 유용성을 입증하였다.

의의 및 주장
본 논문은 계산 결과가 특정 농약, 특히 DDT, 클로리피리포스, 엔도설판이 ERα의 천연 호르몬 포켓 내에서 경쟁적으로 결합할 수 있다는 가설을 뒷받end하는 분자 수준의 증거를 제공한다고 주장한다. 저자들은 결합 에너지가 -6.0 kcal/mol보다 강한 농약은 에스트로겐 신호 전달을 모방하거나 방해함으로써 잠재적인 내분비 교란제로 작용할 수 있다고 제안한다.

본 연구는 이러한 결과가 DDT 및 엔도설판과 같은 화합물의 지속성 및 생물 농축에 관한 기존 독성학적 데이터를 뒷받침하지만, 결과는 예비적인 단계임을 강조한다. 저자들은 본 연구가 복잡한 안정성을 평가하기 위한 장기 MD 시뮬레이션과 in silico 예측을 확인하고 실제 호르몬 경로에 미치는 영향을 평가하기 위한 실험실 실험(예: 수용체 결합 분석 및 유전자 발현 연구)을 포함한 추가적인 실험적 검증의 필요성을 강조하고 있음을 명시하였다.