Hermes Luís Neubauer de Amorim

Linear Interaction Energy (LIE) Method in Lead Discovery and Optimization

Currently, in order to accelerate the process of drug development and also reduce costs, many of the experimental assays related to lead discovery and lead optimization processes are being replaced by computational, in silico, methods. In this context, the LIE (linear interaction energy) method has been used to calculate binding free energies for widely different compounds by averaging interaction energies obtained from molecular dynamics (MD) or Monte Carlo (MC) simulations. In particular, the combination of docking and affinity predictions with the LIE method can thus save valuable resources in lead discovery and optimization projects. This review presents a description of LIE methodology and some recent studies that illustrate the importance and utility of the method in the field of pharmaceutical research.

Thrombin allosteric modulation revisited: a molecular dynamics study.

The regulatory properties of thrombin are derived predominantly from its capacity to produce different functional conformations. Functional studies have revealed that two antagonistic thrombin conformations exist in equilibrium: the fast (procoagulant) and slow (anticoagulant) forms. The mechanisms whereby thrombin activity is regulated by the binding of different effectors remain among the most enigmatic and controversial subjects in the field of protein function. In order to obtain more detailed information on the dynamic events originating from the interaction with the Na+ effector and ligand binding at the active site and anion binding exosite 1 (ABE1), we carried out molecular dynamics simulations of thrombin in different bound states. The results indicated that Na+ release results in a more closed conformation of thrombin, which can be compared to the slow form. The conformational changes induced by displacement of the sodium ion from the Na-binding site include: (1) distortion of the 220- and 186-loops that constitute the Na-binding site; (2) folding back of the Trp148 loop towards the body of the protein, (3) a 180° rotation of the Asp189 side-chain, and (4) projection of the Trp60D loop toward the solvent accompanied by the rearrangement of the Trp215 side chain toward the 95–100 loop. Our findings correlate well with the known structural and recognition properties of the slow and fast forms of thrombin, and are in accordance with the hypothesis that there is communication between the diverse functional domains of thrombin. The theoretical models generated from our MD simulations complement and advance the structural information currently available, leading to a more detailed understanding of thrombin structure and function.

Sequence and structural characterization of tbnat gene in isoniazid-resistant Mycobacterium tuberculosis: identification of new mutations.

The present study was carried out to investigate the presence of polymorphism in the N-acetyltransferase gene of 41 clinical isolates of Mycobacterium tuberculosis, that were resistant to isoniazid (INH) with no mutations in the hot spots of the genes previously described to be involved in INH resistance (katG, inhA and ahpC). We observed single nucleotide polymorphisms (SNPs) in ten of these, including the G619A SNP in five isolates and an additional four so far un-described mutations in another five isolates. Among the latter SNPs, two were synonymous (C276T, n=1 and C375G, n=3), while two more non-synonymous SNPs were composed of C373A (Leu→Met) and T503G (Met→Arg) were observed in respectively one and two isolates. Molecular modeling and structural analysis based in a constructed full length 3D models of wild type TBNAT (TBNAT_H37Rv) and the isoforms (TBNAT_L125M and TBNAT_M168R) were also performed. The refined models show that, just as observed in human NATs, the carboxyl terminus extends deep within the folded enzyme, into close proximity to the buried catalytic triad. Analysis of tbnat that present non-synonymous mutations indicates that both substitutions are plausible to affect enzyme specificity or acetyl-CoA binding capacity. The results contribute to a better understanding of structure-function relationships of NATs. However, further investigation including INH-sensitive strains as a control group is needed to get better understanding of the possible role of these new mutations on tuberculosis control.

Interaction of wild type, G68R and L125M isoforms of the arylamine- N -acetyltransferase from Mycobacterium tuberculosis with isoniazid: a computational study on a new possible mechanism of resistance

AbstractIsoniazid (INH) is a front-line drug used in the treatment of tuberculosis (TB), a disease that remains a major cause of death worldwide. Isoniazid is a prodrug, requiring activation in the mycobacterial cell by the catalase-peroxidase (CP) enzyme. Recent studies have suggested that acetylation of INH by the arylamine-N-acetyltransferase from Mycobacterium tuberculosis (TBNAT) may be a possible cause of inactivation of the drug thus resulting in resistant strains. In this study, computational techniques were applied to investigate the binding of isoniazid to three TBNAT isoforms: wild type, G68R and L125M. Since there is no experimental structure available, molecular dynamics (MD) simulations were initially used for the refinement of TBNAT homology models. Distinct conformations of the models were selected during the production stage of MD simulations for molecular docking experiments with the drug. Finally, each mode of binding was refined by new molecular MD simulations. Essential dynamics (ED) analysis and linear interaction energy calculations (LIE) were used to evaluate the impact of amino acid substitutions on the structural and binding properties of the enzymes. The results suggest that the wild type and the G68R TBNATs have a similar pattern of affinity to INH. On the other hand, the calculated enzyme-INH dissociation constant (KD) was estimated 33 times lower for L125M isoform in comparison with wild type enzyme. This last finding is consistent with the hypothesis that isolated mutations in the tbnat gene can produce M. tuberculosis strains resistant to isoniazid. FigureGlobal structure of the complexes and binding mode of INH calculated for the last frame of MD simulations