Conformational analysis of nevirapine, a non-nucleoside HIV-1 reverse transcriptase inhibitor, based on quantum mechanical calculations

Theoretical investigation of the interactions in binding pocket of Reverse Transcriptase

Interactions in proteins have been studied using several chemical information techniques including quantum chemical methods that are applied to truncated systems composed of the ligand molecule and the surrounding amino acids of the receptor. In this work we adopt an approach to study these interactions accounting for as many as possible explicit solvent molecules and without the need of a fragmented calculation. Furthermore, we embed our quantum chemical calculations within a molecular dynamics framework that enables a fundamentally fast system for quantum molecular dynamic simulations (QCMD). Central to this new system for QCMD is the tight binding QC system, newly developed in our laboratories, and which combined with the MD paradigm results in an ultra accelerated QCMD method for protein–ligand interaction evaluations. We have applied our newly developed method to the Nevirapine (NVP)–Reverse Transcriptase (RT) system. We show how the proposed method leads us to new findings. The advanced QCMD was applied to a system of RT with NVP and it has led to the knowledge of specific groups and atoms that interact with surrounding amino acids of RT and help in drug binding. The information derived from this calculation may be used in designing drugs for NVP resistant virus strains that have binding capability like NVP.

Is nevirapine atropisomeric? Experimental and computational evidence for rapid conformational inversion

The non-nucleoside reverse transcriptase inhibitor nevirapine displays in its room temperature (1)H-NMR spectrum signals characteristic of a chiral compound. Following suggestions in the recent literature that nevirapine may display atropisomerism-and therefore be a chiral compound, due to slow interconversion between two enantiomeric conformers-we report the results of an NMR and computational study which reveal that while nevirapine does indeed possess two stable enantiomeric conformations, they interconvert with a barrier of about 76 kJ mol(-1) at room temperature. Nevirapine has a half life for enantiomerisation at room temperature of the order of seconds, is not atropisomeric, and cannot exist as separable enantiomers.

CONFORMATIONAL ANALYSIS OF NEVIRAPINE IN SOLUTIONS BASED ON NMR SPECTROSCOPY AND QUANTUM CHEMICAL CALCULATIONS

The conformational analysis of HIV-1 Reverse Transcriptase Inhibitor, nevirapine, 11-cyclopropyl-5,-11dihydro-4-methyl-6H-dipyrido[3,2-b2′,3′-e][1,4]diazepin-6-one, was investigated using ab initio and density functional theory calculations. The fully optimized structures and rotational potential energies of the nitrogen and carbon bonds in the cyclopropyl ring (C15-N11-C17-C19, α) were examined in detail. Geometries obtained from all applied calculations show similarities to the complex structure with HIV-1 reverse transcriptase. To obtain more information on the structure, conformational minima of nevirapine, optimized at the B3LYP/6-31G** level, were calculated for the 1H, 13C, and 15N-NMR chemical shifts at the B3LYP/6-311++G** level using the GIAO approach in DMSO and chloroform IEFPCM solvation models. The calculated 1H, 13C-NMR chemical shifts agree well with the experimental data, which indicates that the geometry of nevirapine in solution is similar to that of the molecule in the inhibition complex. Solvation free energies (ΔG sol ) of nevirapine in DMSO and chloroform were also obtained.

Bridge water mediates nevirapine binding to wild type and Y181C HIV-1 reverse transcriptase--evidence from molecular dynamics simulations and MM-PBSA calculations

The important role of the bridge water molecule in the binding of HIV-1 reverse transcriptase (RT) inhibitor complex was elucidated by molecular dynamics (MD) simulations using an MM-PBSA approach. Binding free energies and thermodynamic property differences for nevirapine bound to wild type and Y181C HIV-1 reverse transcriptase were investigated, and the results were compared with available experimental data. MD simulations over 3 ns revealed that the bridge water formed three characteristic hydrogen bonds to nevirapine and two residues, His235 and Leu234, in the binding pocket. The energetic derived model, which was determined from the consecutive addition of a water molecule, confirmed that only the contribution from the bridge water was essential in the binding configuration. Including this bridge water in the MM-PBSA calculations reoriented the binding energies from -32.20 to -37.65 kcal/mol and -28.07 to -29.82 kcal/mol in the wild type and Y181C HIV-1 RT, respectively. From the attractive interactions via the bridge water, His235 and Leu234 became major contributions. We found that the bridge water is the key in stabilizing the bound complex; however, in the Y181C RT complex this bridge water showed weaker hydrogen bond formation, lack of attractive force to nevirapine and lack of binding efficiency, leading to the failure of nevirapine against the Y181C HIV-1 RT. Moreover, the dynamics of Val179, Tyr181Cys, Gly190 and Leu234 in the binding pocket showed additional attractive energetic contributions in helping nevirapine binding. These findings that the presence of a water molecule in the hydrophobic binding site plays an important role are a step towards a quantitative understanding of the character of bridge water in enzyme-inhibitor binding. This can be helpful in developing designs for novel non-nucleoside HIV-1 RT inhibitors active against the mutant enzyme.