Clarifying allosteric control of flap conformations in the 1TW7 crystal structure of HIV-1 protease

Inhibition of the activity of HIV-1 protease through antibody binding and mutations probed by molecular dynamics simulations

HIV-1 protease is an essential enzyme in the life cycle of the HIV-1 virus. The conformational dynamics of the flap region of the protease is critical for the ligand binding mechanism, as well as for the catalytic activity. The monoclonal antibody F11.2.32 raised against HIV-1 protease inhibits its activity on binding. We have studied the conformational dynamics of protease in its free, inhibitor ritonavir and antibody bound forms using molecular dynamics simulations. We find that upon Ab binding to the epitope region (residues 36-46) of protease, the overall flexibility of the protease is decreased including the flap region and the active site, which is similar to the decrease in flexibility observed by inhibitor binding to the protease. This suggests an allosteric mechanism to inhibit protease activity. Further, the protease mutants G40E and G40R are known to have decreased activity and were also subjected to MD simulations. We find that the loss of flexibility in the mutants is similar to that observed in the protease bound to the Ab/inhibitor. These insights highlight the role played by dynamics in the function of the protease and how control of flexibility through Ab binding and site specific mutations can inhibit protease activity.

Enhanced conformational sampling technique provides an energy landscape view of large-scale protein conformational transitions

A novel in silico approach (NMA–ITS) is introduced to rapidly and effectively sample the configuration space and give quantitative data for exploring the conformational changes of proteins.

Potential selective inhibitors against Rv0183 of Mycobacterium tuberculosis targeting host lipid metabolism

Tuberculosis is the second leading infectious killer with 9 million new cases in 2009. Extensive use of pathogen's lipid metabolism especially in utilizing the host lipids and virulence highlights the importance of exported lipid-catabolizing enzymes. Current study aims to emphasize the importance of Rv0183, an exported monoacylglycerol lipase, involved in metabolizing the host cell membrane lipids. Sequence analysis and homology modeling shows Rv0183 is highly conserved throughout mycobacterial species even in Mycobacterium leprae and also significantly divergent from mammalian lipases. Additionally, employing virtual screening using NCI diversity set and ZINC database with criteria of molecules with higher predicted free energy of binding toward Rv0183 than human lipase, potential inhibitors have been identified for Rv0183. A tautomer of ZINC13451138, known inhibitor for HIV-1 integrase is the best hit with difference in free energy of binding of 8.72 kcal/mol. The sequence and structure analysis were helpful in identifying the ligand binding sites and molecular function of the mycobacterial specific monoacylglycerol lipase. Rv0183 represents a suitable and promising drug target and is also a step towards understanding dormancy development and reactivation, thereby addressing pathogen's drug resistance. Experimental studies on the discovered potential inhibitors in this virtual screen should further validate the therapeutic utility of Rv0183.

Contribution of the 80s loop of HIV-1 protease to the multidrug-resistance mechanism: crystallographic study of MDR769 HIV-1 protease variants

The flexible flaps and the 80s loops (Pro79-Ile84) of HIV-1 protease are crucial in inhibitor binding. Previously, it was reported that the crystal structure of multidrug-resistant 769 (MDR769) HIV-1 protease shows a wide-open conformation of the flaps owing to conformational rigidity acquired by the accumulation of mutations. In the current study, the effect of mutations on the conformation of the 80s loop of MDR769 HIV-1 protease variants is reported. Alternate conformations of Pro81 (proline switch) with a root-mean-square deviation of 3-4.8 Å in the C(α) atoms of the I10V mutant and a side chain with a `flipped-out' conformation in the A82F mutant cause distortion in the S1/S1' binding pockets that affects inhibitor binding. The A82S and A82T mutants show local changes in the electrostatics of inhibitor binding owing to the mutation from nonpolar to polar residues. In summary, the crystallographic studies of four variants of MDR769 HIV-1 protease presented in this article provide new insights towards understanding the drug-resistance mechanism as well as a basis for design of future protease inhibitors with enhanced potency.

Insights into the dynamics of HIV-1 protease: a kinetic network model constructed from atomistic simulations

The conformational dynamics in the flaps of HIV-1 protease plays a crucial role in the mechanism of substrate binding. We develop a kinetic network model, constructed from detailed atomistic simulations, to determine the kinetic mechanisms of the conformational transitions in HIV-1 PR. To overcome the time scale limitation of conventional molecular dynamics (MD) simulations, our method combines replica exchange MD with transition path theory (TPT) to study the diversity and temperature dependence of the pathways connecting functionally important states of the protease. At low temperatures the large-scale flap opening is dominated by a small number of paths; at elevated temperatures the transition occurs through many structurally heterogeneous routes. The expanded conformation in the crystal structure 1TW7 is found to closely mimic a key intermediate in the flap-opening pathways at low temperature. We investigated the different transition mechanisms between the semi-open and closed forms. The calculated relaxation times reveal fast semi-open ↔ closed transitions, and infrequently the flaps fully open. The ligand binding rate predicted from this kinetic model increases by 38-fold from 285 to 309 K, which is in general agreement with experiments. To our knowledge, this is the first application of a network model constructed from atomistic simulations together with TPT to analyze conformational changes between different functional states of a natively folded protein.

Structure-based design of carbon nanotubes as HIV-1 protease inhibitors: atomistic and coarse-grained simulations

Nanoparticles such as fullerenes and carbon nanotubes have been extensively studied for biomedical applications. In this paper, we report the design of carbon nanotubes as HIV-1 protease inhibitors. Docking and molecular dynamics calculations are performed using an atomistic model to explore the optimal interaction structure and free energy between the nanotube and HIV-1 protease. A coarse-grained model is then developed based on the atomistic model, allowing us to investigate the dynamic behaviors of the protease in the bound and unbound states. The dynamic process reveals that the carbon nanotube is able to bind to the active site of the protease and prevent the active flaps from opening up, thus blocking the function of the protease. This process is strongly influenced by the size of the nanotube. The binding of carbon nanotubes to an alternative binding site other than the active site is also explored. Therefore, carbon nanotube-based inhibitors have great potential for application as HIV-1 protease inhibitors.

Structure of the unbound form of HIV-1 subtype A protease: comparison with unbound forms of proteases from other HIV subtypes

The crystal structure of the unbound form of HIV-1 subtype A protease (PR) has been determined to 1.7 A resolution and refined as a homodimer in the hexagonal space group P6(1) to an R(cryst) of 20.5%. The structure is similar in overall shape and fold to the previously determined subtype B, C and F PRs. The major differences lie in the conformation of the flap region. The flaps in the crystal structures of the unbound subtype B and C PRs, which were crystallized in tetragonal space groups, are either semi-open or wide open. In the present structure of subtype A PR the flaps are found in the closed position, a conformation that would be more anticipated in the structure of HIV protease complexed with an inhibitor. The amino-acid differences between the subtypes and their respective crystal space groups are discussed in terms of the differences in the flap conformations.

Method to Predict Crowding Effects by Postprocessing Molecular Dynamics Trajectories: Application to the Flap Dynamics of HIV-1 Protease

The internal dynamics of proteins inside of cells may be affected by the crowded intracellular environments. Here, we test a novel approach to simulations of crowding, in which simulations in the absence of crowders are postprocessed to predict crowding effects, against the direct approach of simulations in the presence of crowders. The effects of crowding on the flap dynamics of HIV-1 protease predicted by the postprocessing approach are found to agree well with those calculated by the direct approach. The postprocessing approach presents distinct advantages over the direct approach in terms of accuracy and speed and is expected to have broad impact on atomistic simulations of macromolecular crowding.