Comparing Molecular Modeling and Dynamics with Crystallography: Brazilian HIV-1 Protease mutants
Batista, P.R.1, Loureiro, P.A.G.L.1, Sanches, M 2, Polikarpov, I 2 and Pascutti, PG1
1 IBCCF – UFRJ – RJ; 2 Physics Department, USP, Sao Carlos, SP
Introduction:
An important question relates to effectiveness of anti-HIV therapy
against different subtypes. All HIV-1 inhibitors were developed for B
subtype, the prevalent in developed countries but not in the world. In Brazil,
B subtype is the major but F presents an important role in drug
resistance. The differences between subtypes are present in many
targets of anti-HIV drugs, as the protease (PR). There's no structural
information of non-B PR in the Protein Data Bank and
the first crystals were done by Polikarpov's group however these
structures were not published yet.TL-3 inhibits efficiently PR
containing the known resistance mutations G48V and V82F against PR
inhibitors (PI). Objectives:
To better understand resistance and TL-3 efficiency against different
subtypes, we exploit Molecular Dynamics (MD) to perform a structural
study with PR of Bwt, Fwt, B and F mutants. In addition, with the
crystallographic unpublished structures we played MD simulations in
order to compare with the modeled ones, making possible analyze the
efficacy of Molecular Modeling as a structure predictor tool. Methods:
We applied Comparative Modeling using as template the B subtype PR
structure (3THL) complexed with TL‑3, obtaining Fwt, Bmut and Fmut
models. The sequences utilized were based on HIV‑1 drug-resistant
Brazilian patients' isolates. Then we carried out MD simulations (10ns)
in a periodic boundary water box, using GROMACS package. TL-3 topology
was built to GROMOS96 FF, with aid of PRODRG server and ab initio calculations. Results:
The MD estimated Cα's B-factor was in agreement with the
crystallographic one. Essential Dynamics analyses reveal
destabilization of active site and flaps in Fmut trajectories,
confirmed by RMS fluctuations around the catalytic ASP 25 and ILE50.
High mobility of P4 and P4' groups were observed. Only in Fmut, P1
group presents large deviation lowering its affinity, confirmed by ΔG
calculations. Now, are in progress MD analyzes with the crystallography
trajectories. Conclusions:
Our MD results of Fmut PR showed a stability and affinity decrease when
compared to B subtype. Analyzes of the crystallography structure's MD
could provide an important contribution for validation of this modeling
technique.
|