1
|
The structural, dynamic, and thermodynamic basis of darunavir resistance of a heavily mutated HIV-1 protease using molecular dynamics simulation. Front Mol Biosci 2022; 9:927373. [PMID: 36046605 PMCID: PMC9420863 DOI: 10.3389/fmolb.2022.927373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2022] [Accepted: 06/28/2022] [Indexed: 11/13/2022] Open
Abstract
The human immunodeficiency virus type 1 protease (HIV-1 PR) is an important enzyme in the life cycle of the HIV virus. It cleaves inactive pre-proteins of the virus and changes them into active proteins. Darunavir (DRV) suppresses the wild-type HIV-1 PR (WT-Pr) activity but cannot inhibit some mutant resistant forms (MUT-Pr). Increasing knowledge about the resistance mechanism can be helpful for designing more effective inhibitors. In this study, the mechanism of resistance of a highly MUT-Pr strain against DRV was investigated. For this purpose, complexes of DRV with WT-Pr (WT-Pr-D) and MUT-Pr (MUT-Pr-D) were studied by all-atom molecular dynamics simulation in order to extract the dynamic and energetic properties. Our data revealed that mutations increased the flap-tip flexibility due to the reduction of the flap-flap hydrophobic interactions. So, the protease’s conformation changed from a closed state to a semi-open state that can facilitate the disjunction of DRV from the active site. On the other hand, energy analysis limited to the final basins of the energy landscape indicated that the entropy of binding of DRV to MUT-Pr was more favorable than that of WT-Pr. However, the enthalpy penalty overcomes it and makes binding more unfavorable relative to the WT-Pr. The unfavorable interaction of DRV with R8, I50, I84, D25′, and A28′ residues in MUT-Pr-D relative to WT-Pr-D is the reason for this enthalpy penalty. Thus, mutations drive resistance to DRV. The hydrogen bond analysis showed that compared with WT-Pr, the hydrogen bonds between DRV and the active-site residues of MUT-Pr were disrupted.
Collapse
|
2
|
Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor. Phys Chem Chem Phys 2022; 24:3586-3597. [PMID: 35089990 PMCID: PMC8940534 DOI: 10.1039/d1cp05487b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ∼5-50 cm-1. This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effects of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal mol-1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.
Collapse
|
3
|
Multiscale Approach for Computing Gated Ligand Binding from Molecular Dynamics and Brownian Dynamics Simulations. J Chem Theory Comput 2021; 17:7912-7929. [PMID: 34739248 DOI: 10.1021/acs.jctc.1c00673] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We develop an approach to characterize the effects of gating by a multiconformation protein consisting of macrostate conformations that are either accessible or inaccessible to ligand binding. We first construct a Markov state model of the apo-protein from atomistic molecular dynamics simulations from which we identify macrostates and their conformations, compute their relative macrostate populations and interchange kinetics, and structurally characterize them in terms of ligand accessibility. We insert the calculated first-order rate constants for conformational transitions into a multistate gating theory from which we derive a gating factor γ that quantifies the degree of conformational gating. Applied to HIV-1 protease, our approach yields a kinetic network of three accessible (semi-open, open, and wide-open) and two inaccessible (closed and a newly identified, "parted") macrostate conformations. The parted conformation sterically partitions the active site, suggesting a possible role in product release. We find that the binding kinetics of drugs and drug-like inhibitors to HIV-1 protease falls in the slow gating regime. However, because γ = 0.75, conformational gating only modestly slows ligand binding. Brownian dynamics simulations of the diffusional association of eight inhibitors to the protease─having a wide range of experimental association constants (∼104-1010 M-1 s-1)─yields gated rate constants in the range of ∼0.5-5.7 × 108 M-1 s-1. This indicates that, whereas the association rate of some inhibitors could be described by the model, for many inhibitors either subsequent conformational transitions or alternate binding mechanisms may be rate-limiting. For systems known to be modulated by conformational gating, the approach could be scaled computationally efficiently to screen association kinetics for a large number of ligands.
Collapse
|
4
|
DEER-PREdict: Software for efficient calculation of spin-labeling EPR and NMR data from conformational ensembles. PLoS Comput Biol 2021; 17:e1008551. [PMID: 33481784 PMCID: PMC7857587 DOI: 10.1371/journal.pcbi.1008551] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2020] [Revised: 02/03/2021] [Accepted: 11/19/2020] [Indexed: 11/25/2022] Open
Abstract
Owing to their plasticity, intrinsically disordered and multidomain proteins require descriptions based on multiple conformations, thus calling for techniques and analysis tools that are capable of dealing with conformational ensembles rather than a single protein structure. Here, we introduce DEER-PREdict, a software program to predict Double Electron-Electron Resonance distance distributions as well as Paramagnetic Relaxation Enhancement rates from ensembles of protein conformations. DEER-PREdict uses an established rotamer library approach to describe the paramagnetic probes which are bound covalently to the protein.DEER-PREdict has been designed to operate efficiently on large conformational ensembles, such as those generated by molecular dynamics simulation, to facilitate the validation or refinement of molecular models as well as the interpretation of experimental data. The performance and accuracy of the software is demonstrated with experimentally characterized protein systems: HIV-1 protease, T4 Lysozyme and Acyl-CoA-binding protein. DEER-PREdict is open source (GPLv3) and available at github.com/KULL-Centre/DEERpredict and as a Python PyPI package pypi.org/project/DEERPREdict.
Collapse
|
5
|
Increasing the Sampling Efficiency of Protein Conformational Change by Combining a Modified Replica Exchange Molecular Dynamics and Normal Mode Analysis. J Chem Theory Comput 2020; 17:13-28. [PMID: 33351613 DOI: 10.1021/acs.jctc.0c00592] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Understanding conformational change at an atomic level is significant when determining a protein functional mechanism. Replica exchange molecular dynamics (REMD) is a widely used enhanced sampling method to explore protein conformational space. However, REMD with an explicit solvent model requires huge computational resources, immensely limiting its application. In this study, a variation of parallel tempering metadynamics (PTMetaD) with the omission of solvent-solvent interactions in exchange attempts and the use of low-frequency modes calculated by normal-mode analysis (NMA) as collective variables (CVs), namely ossPTMetaD, is proposed with the aim to accelerate MD simulations simultaneously in temperature and geometrical spaces. For testing the performance of ossPTMetaD, five protein systems with diverse biological functions and motion patterns were selected, including large-scale domain motion (AdK), flap movement (HIV-1 protease and BACE1), and DFG-motif flip in kinases (p38α and c-Abl). The simulation results showed that ossPTMetaD requires much fewer numbers of replicas than temperature REMD (T-REMD) with a reduction of ∼70% to achieve a similar exchange ratio. Although it does not obey the detailed balance condition, ossPTMetaD provides consistent results with T-REMD and experimental data. The high accessibility of the large conformational change of protein systems by ossPTMetaD, especially in simulating the very challenging DFG-motif flip of protein kinases, demonstrated its high efficiency and robustness in the characterization of the large-scale protein conformational change pathway and associated free energy profile.
Collapse
|
6
|
Highly drug-resistant HIV-1 protease reveals decreased intra-subunit interactions due to clusters of mutations. FEBS J 2020; 287:3235-3254. [PMID: 31920003 PMCID: PMC7343616 DOI: 10.1111/febs.15207] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2019] [Revised: 11/16/2019] [Accepted: 01/08/2020] [Indexed: 01/07/2023]
Abstract
Drug-resistance is a serious problem for treatment of the HIV/AIDS pandemic. Potent clinical inhibitors of HIV-1 protease show several orders of magnitude worse inhibition of highly drug-resistant variants. Hence, the structure and enzyme activities were analyzed for HIV protease mutant HIV-1 protease (EC 3.4.23.16) (PR) with 22 mutations (PRS5B) from a clinical isolate that was selected by machine learning to represent high-level drug-resistance. PRS5B has 22 mutations including only one (I84V) in the inhibitor binding site; however, clinical inhibitors had poor inhibition of PRS5B activity with kinetic inhibition value (Ki ) values of 4-1000 nm or 18- to 8000-fold worse than for wild-type PR. High-resolution crystal structures of PRS5B complexes with the best inhibitors, amprenavir (APV) and darunavir (DRV) (Ki ~ 4 nm), revealed only minor changes in protease-inhibitor interactions. Instead, two distinct clusters of mutations in distal regions induce coordinated conformational changes that decrease favorable internal interactions across the entire protein subunit. The largest structural rearrangements are described and compared to other characterized resistant mutants. In the protease hinge region, the N83D mutation eliminates a hydrogen bond connecting the hinge and core of the protease and increases disorder compared to highly resistant mutants PR with 17 mutations and PR with 20 mutations with similar hinge mutations. In a distal β-sheet, mutations G73T and A71V coordinate with accessory mutations to bring about shifts that propagate throughout the subunit. Molecular dynamics simulations of ligand-free dimers show differences consistent with loss of interactions in mutant compared to wild-type PR. Clusters of mutations exhibit both coordinated and antagonistic effects, suggesting PRS5B may represent an intermediate stage in the evolution of more highly resistant variants. DATABASES: Structural data are available in Protein Data Bank under the accession codes 6P9A and 6P9B for PRS5B/DRV and PRS5B/APV, respectively.
Collapse
|
7
|
Inhibition of the activity of HIV-1 protease through antibody binding and mutations probed by molecular dynamics simulations. Sci Rep 2020; 10:5501. [PMID: 32218488 PMCID: PMC7098958 DOI: 10.1038/s41598-020-62423-y] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2019] [Accepted: 03/12/2020] [Indexed: 01/31/2023] Open
Abstract
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.
Collapse
|
8
|
Protein NMR: Boundless opportunities. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2019; 306:187-191. [PMID: 31311710 PMCID: PMC6703950 DOI: 10.1016/j.jmr.2019.07.037] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/03/2019] [Revised: 04/22/2019] [Accepted: 07/08/2019] [Indexed: 05/04/2023]
Abstract
Over the past approximately three decades, isotope-directed NMR spectroscopy has become a powerful method for determining 3D structures of biological macromolecules and their complexes in solution. From a structural perspective NMR provides an invaluable tool for studying systems that are not amenable to crystallization, including intrinsically disordered proteins and weak complexes. In contrast to both X-ray crystallography and cryo-electron microscopy which afford a largely static view of the systems under consideration, the great power of NMR lies in its ability to quantitatively probe exchange dynamics between interconverting states, and to reveal and characterize at atomic resolution the existence of transient states that may be populated at levels as low as 1%. Such "excited" states play a key role in macromolecular recognition, allostery, signal transduction and macromolecular assembly, including the initial events involved in aggregation and amyloid formation. Optimal application of NMR to such systems of fundamental biological interest requires a sound footing of the physical underpinnings of today's and tomorrow's sophisticated NMR experiments.
Collapse
|
9
|
NMR and MD studies combined to elucidate inhibitor and water interactions of HIV-1 protease and their modulations with resistance mutations. JOURNAL OF BIOMOLECULAR NMR 2019; 73:365-374. [PMID: 31243634 PMCID: PMC6941145 DOI: 10.1007/s10858-019-00260-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Accepted: 06/19/2019] [Indexed: 06/09/2023]
Abstract
Over the last two decades, both the sensitivity of NMR and the time scale of molecular dynamics (MD) simulation have increased tremendously and have advanced the field of protein dynamics. HIV-1 protease has been extensively studied using these two methods, and has presented a framework for cross-evaluation of structural ensembles and internal dynamics by integrating the two methods. Here, we review studies from our laboratories over the last several years, to understand the mechanistic basis of protease drug-resistance mutations and inhibitor responses, using NMR and crystal structure-based parallel MD simulations. Our studies demonstrate that NMR relaxation experiments, together with crystal structures and MD simulations, significantly contributed to the current understanding of structural/dynamic changes due to HIV-1 protease drug resistance mutations.
Collapse
|
10
|
Binding kinetics and substrate selectivity in HIV-1 protease-Gag interactions probed at atomic resolution by chemical exchange NMR. Proc Natl Acad Sci U S A 2017; 114:E9855-E9862. [PMID: 29087351 DOI: 10.1073/pnas.1716098114] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
The conversion of immature noninfectious HIV-1 particles to infectious virions is dependent upon the sequential cleavage of the precursor group-specific antigen (Gag) polyprotein by HIV-1 protease. The precise mechanism whereby protease recognizes distinct Gag cleavage sites, located in the intrinsically disordered linkers connecting the globular domains of Gag, remains unclear. Here, we probe the dynamics of the interaction of large fragments of Gag and various variants of protease (including a drug resistant construct) using Carr-Purcell-Meiboom-Gill relaxation dispersion and chemical exchange saturation transfer NMR experiments. We show that the conformational dynamics within the flaps of HIV-1 protease that form the lid over the catalytic cleft play a significant role in substrate specificity and ordered Gag processing. Rapid interconversion between closed and open protease flap conformations facilitates the formation of a transient, sparsely populated productive complex between protease and Gag substrates. Flap closure traps the Gag cleavage sites within the catalytic cleft of protease. Modulation of flap opening through protease-Gag interactions fine-tunes the lifetime of the productive complex and hence the likelihood of Gag proteolysis. A productive complex can also be formed in the presence of a noncognate substrate but is short-lived owing to lack of optimal complementarity between the active site cleft of protease and the substrate, resulting in rapid flap opening and substrate release, thereby allowing protease to differentiate between cognate and noncognate substrates.
Collapse
|
11
|
Block-restraining of residual dipolar couplings to allow fluctuating relative alignments of molecular subdomains. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2017; 128:133-141. [PMID: 28223155 DOI: 10.1016/j.pbiomolbio.2017.02.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/12/2016] [Accepted: 02/07/2017] [Indexed: 11/19/2022]
Abstract
Residual dipolar couplings (RDCs), unlike most other types of NMR observables, provide orientational information, reporting on the alignment of inter-spin vectors (ISVs) relative to the magnetic field. A great challenge in using experimental RDCs to restrain molecular dynamics (MD) simulations is how to represent this alignment. An alignment tensor is often used to parameterise the contribution of molecular alignment to the angular dependence of RDCs. All ISVs that share the same tensor have fixed relative alignment, i.e. if just one tensor is used, the molecule is internally rigid. Here we propose and illustrate a method for subdividing molecules into individually aligned blocks during MD simulations restrained to fit RDCs. This allows the relative orientation of each block to vary during the simulation, which in turn ensures that the internal structure of each block is more realistically reproduced.
Collapse
|
12
|
Abstract
The study provided an integrated view of the transition pathway of the flap opening of HIV-1 protease using MD simulation.
Collapse
|
13
|
Structural Studies of a Rationally Selected Multi-Drug Resistant HIV-1 Protease Reveal Synergistic Effect of Distal Mutations on Flap Dynamics. PLoS One 2016; 11:e0168616. [PMID: 27992544 PMCID: PMC5161481 DOI: 10.1371/journal.pone.0168616] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2016] [Accepted: 12/02/2016] [Indexed: 12/20/2022] Open
Abstract
We report structural analysis of HIV protease variant PRS17 which was rationally selected by machine learning to represent wide classes of highly drug-resistant variants. Crystal structures were solved of PRS17 in the inhibitor-free form and in complex with antiviral inhibitor, darunavir. Despite its 17 mutations, PRS17 has only one mutation (V82S) in the inhibitor/substrate binding cavity, yet exhibits high resistance to all clinical inhibitors. PRS17 has none of the major mutations (I47V, I50V, I54ML, L76V and I84V) associated with darunavir resistance, but has 10,000-fold weaker binding affinity relative to the wild type PR. Comparable binding affinity of 8000-fold weaker than PR is seen for drug resistant mutant PR20, which bears 3 mutations associated with major resistance to darunavir (I47V, I54L and I84V). Inhibitor-free PRS17 shows an open flap conformation with a curled tip correlating with G48V flap mutation. NMR studies on inactive PRS17D25N unambiguously confirm that the flaps adopt mainly an open conformation in solution very similar to that in the inhibitor-free crystal structure. In PRS17, the hinge loop cluster of mutations, E35D, M36I and S37D, contributes to the altered flap dynamics by a mechanism similar to that of PR20. An additional K20R mutation anchors an altered conformation of the hinge loop. Flap mutations M46L and G48V in PRS17/DRV complex alter the Phe53 conformation by steric hindrance between the side chains. Unlike the L10F mutation in PR20, L10I in PRS17 does not break the inter-subunit ion pair or diminish the dimer stability, consistent with a very low dimer dissociation constant comparable to that of wild type PR. Distal mutations A71V, L90M and I93L propagate alterations to the catalytic site of PRS17. PRS17 exhibits a molecular mechanism whereby mutations act synergistically to alter the flap dynamics resulting in significantly weaker binding yet maintaining active site contacts with darunavir.
Collapse
|
14
|
Abstract
It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.
Collapse
|
15
|
Transient HIV-1 Gag-protease interactions revealed by paramagnetic NMR suggest origins of compensatory drug resistance mutations. Proc Natl Acad Sci U S A 2016; 113:12456-12461. [PMID: 27791180 DOI: 10.1073/pnas.1615342113] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Cleavage of the group-specific antigen (Gag) polyprotein by HIV-1 protease represents the critical first step in the conversion of immature noninfectious viral particles to mature infectious virions. Selective pressure exerted by HIV-1 protease inhibitors, a mainstay of current anti-HIV-1 therapies, results in the accumulation of drug resistance mutations in both protease and Gag. Surprisingly, a large number of these mutations (known as secondary or compensatory mutations) occur outside the active site of protease or the cleavage sites of Gag (located within intrinsically disordered linkers connecting the globular domains of Gag to one another), suggesting that transient encounter complexes involving the globular domains of Gag may play a role in guiding and facilitating access of the protease to the Gag cleavage sites. Here, using large fragments of Gag, as well as catalytically inactive and active variants of protease, we probe the nature of such rare encounter complexes using intermolecular paramagnetic relaxation enhancement, a highly sensitive technique for detecting sparsely populated states. We show that Gag-protease encounter complexes are primarily mediated by interactions between protease and the globular domains of Gag and that the sites of transient interactions are correlated with surface exposed regions that exhibit a high propensity to mutate in the presence of HIV-1 protease inhibitors.
Collapse
|
16
|
Evolution under Drug Pressure Remodels the Folding Free-Energy Landscape of Mature HIV-1 Protease. J Mol Biol 2016; 428:2780-92. [PMID: 27170547 PMCID: PMC4905781 DOI: 10.1016/j.jmb.2016.05.005] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2016] [Revised: 04/27/2016] [Accepted: 05/02/2016] [Indexed: 01/08/2023]
Abstract
Using high-pressure NMR spectroscopy and differential scanning calorimetry, we investigate the folding landscape of the mature HIV-1 protease homodimer. The cooperativity of unfolding was measured in the absence or presence of a symmetric active site inhibitor for the optimized wild type protease (PR), its inactive variant PRD25N, and an extremely multidrug-resistant mutant, PR20. The individual fit of the pressure denaturation profiles gives rise to first order, ∆GNMR, and second order, ∆VNMR (the derivative of ∆GNMR with pressure); apparent thermodynamic parameters for each amide proton considered. Heterogeneity in the apparent ∆VNMR values reflects departure from an ideal cooperative unfolding transition. The narrow to broad distribution of ∆VNMR spanning the extremes from inhibitor-free PR20D25N to PR-DMP323 complex, and distinctively for PRD25N-DMP323 complex, indicated large variations in folding cooperativity. Consistent with this data, the shape of thermal unfolding transitions varies from asymmetric for PR to nearly symmetric for PR20, as dimer-inhibitor ternary complexes. Lack of structural cooperativity was observed between regions located close to the active site, including the hinge and tip of the glycine-rich flaps, and the rest of the protein. These results strongly suggest that inhibitor binding drastically decreases the cooperativity of unfolding by trapping the closed flap conformation in a deep energy minimum. To evade this conformational trap, PR20 evolves exhibiting a smoother folding landscape with nearly an ideal two-state (cooperative) unfolding transition. This study highlights the malleability of retroviral protease folding pathways by illustrating how the selection of mutations under drug pressure remodels the free-energy landscape as a primary mechanism.
Collapse
|
17
|
ARTSY-J: Convenient and precise measurement of (3)JHNHα couplings in medium-size proteins from TROSY-HSQC spectra. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2016; 268:73-81. [PMID: 27179455 PMCID: PMC4915346 DOI: 10.1016/j.jmr.2016.05.001] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2016] [Revised: 04/22/2016] [Accepted: 05/02/2016] [Indexed: 06/05/2023]
Abstract
A new and convenient method, named ARTSY-J, is introduced that permits extraction of the (3)JHNHα couplings in proteins from the relative intensities in a pair of (15)N-(1)H TROSY-HSQC spectra. The pulse scheme includes (3)JHNHα dephasing of the narrower TROSY (1)H(N)-{(15)N} doublet component during a delay, integrated into the regular two-dimensional TROSY-HSQC pulse scheme, and compares the obtained intensity with a reference spectrum where (3)JHNHα dephasing is suppressed. The effect of passive (1)H(α) spin flips downscales the apparent (3)JHNHα coupling by a uniform factor that depends approximately linearly on both the duration of the (3)JHNHα dephasing delay and the (1)H-(1)H cross relaxation rate. Using such a correction factor, which accounts for the effects of both inhomogeneity of the radiofrequency field and (1)H(α) spin flips, agreement between prior and newly measured values for the small model protein GB3 is better than 0.3Hz. Measurement for the HIV-1 protease homodimer (22kDa) yields (3)JHNHα values that agree to better than 0.7Hz with predictions made on the basis of a previously parameterized Karplus equation. Although for Gly residues the two individual (3)JHNHα couplings cannot be extracted from a single set of ARTSY-J spectra, the measurement provides valuable ϕ angle information.
Collapse
|
18
|
Unraveling HIV protease flaps dynamics by Constant pH Molecular Dynamics simulations. J Struct Biol 2016; 195:216-226. [PMID: 27291071 DOI: 10.1016/j.jsb.2016.06.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Revised: 06/08/2016] [Accepted: 06/09/2016] [Indexed: 11/15/2022]
Abstract
The active site of HIV protease (HIV-PR) is covered by two flaps. These flaps are known to be essential for the catalytic activity of the HIV-PR, but their exact conformations at the different stages of the enzymatic pathway remain subject to debate. Understanding the correct functional dynamics of the flaps might aid the development of new HIV-PR inhibitors. It is known that, the HIV-PR catalytic efficiency is pH-dependent, likely due to the influence of processes such as charge transfer and protonation/deprotonation of ionizable residues. Several Molecular Dynamics (MD) simulations have reported information about the HIV-PR flaps. However, in MD simulations the protonation of a residue is fixed and thus it is not possible to study the correlation between conformation and protonation state. To address this shortcoming, this work attempts to capture, through Constant pH Molecular Dynamics (CpHMD), the conformations of the apo, substrate-bound and inhibitor-bound HIV-PR, which differ drastically in their flap arrangements. The results show that the HIV-PR flaps conformations are defined by the protonation of the catalytic residues Asp25/Asp25' and that these residues are sensitive to pH changes. This study suggests that the catalytic aspartates can modulate the opening of the active site and substrate binding.
Collapse
|
19
|
Abstract
It is now common knowledge that enzymes are mobile entities relying on complex atomic-scale dynamics and coordinated conformational events for proper ligand recognition and catalysis. However, the exact role of protein dynamics in enzyme function remains either poorly understood or difficult to interpret. This mini-review intends to reconcile biophysical observations and biological significance by first describing a number of common experimental and computational methodologies employed to characterize atomic-scale residue motions on various timescales in enzymes, and second by illustrating how the knowledge of these motions can be used to describe the functional behavior of enzymes and even act upon it. Two biologically relevant examples will be highlighted, namely the HIV-1 protease and DNA polymerase β enzyme systems.
Collapse
|
20
|
Binding of Clinical Inhibitors to a Model Precursor of a Rationally Selected Multidrug Resistant HIV-1 Protease Is Significantly Weaker Than That to the Released Mature Enzyme. Biochemistry 2016; 55:2390-400. [PMID: 27039930 DOI: 10.1021/acs.biochem.6b00012] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have systematically validated the activity and inhibition of a HIV-1 protease (PR) variant bearing 17 mutations (PR(S17)), selected to represent high resistance by machine learning on genotype-phenotype data. Three of five mutations in PR(S17) correlating with major drug resistance, M46L, G48V, and V82S, and five of 11 natural variations differ from the mutations in two clinically derived extreme mutants, PR20 and PR22 bearing 19 and 22 mutations, respectively. PR(S17), which forms a stable dimer (<10 nM), is ∼10- and 2-fold less efficient in processing the Gag polyprotein than the wild type and PR20, respectively, but maintains the same cleavage order. Isolation of a model precursor of PR(S17) flanked by the 56-amino acid transframe region (TFP-p6pol) at its N-terminus, which is impossible upon expression of an analogous PR20 precursor, allowed systematic comparison of inhibition of TFP-p6pol-PR(S17) and mature PR(S17). Resistance of PR(S17) to eight protease inhibitors (PIs) relative to PR (Ki) increases by 1.5-5 orders of magnitude from 0.01 to 8.4 μM. Amprenavir, darunavir, atazanavir, and lopinavir, the most effective of the eight PIs, inhibit precursor autoprocessing at the p6pol/PR site with IC50 values ranging from ∼7.5 to 60 μM. Thus, this process, crucial for stable dimer formation, shows inhibition ∼200-800-fold weaker than that of the mature PR(S17). TFP/p6pol cleavage, which occurs faster, is inhibited even more weakly by all PIs except darunavir (IC50 = 15 μM); amprenavir shows a 2-fold increase in IC50 (∼15 μM), and atazanavir and lopinavir show increased IC50 values of >42 and >70 μM, respectively.
Collapse
|
21
|
Abstract
The virally encoded protease is an important drug target for AIDS therapy. Despite the potency of the current drugs, infections with resistant viral strains limit the long-term effectiveness of therapy. Highly resistant variants of HIV protease from clinical isolates have different combinations of about 20 mutations and several orders of magnitude worse binding affinity for clinical inhibitors. Strategies are being explored to inhibit these highly resistant mutants. The existing inhibitors can be modified by introducing groups with the potential to form new interactions with conserved protease residues, and the flexible flaps. Alternative strategies are discussed, including designing inhibitors to bind to the open conformation of the protease dimer, and inhibition of the protease-catalyzed processing of the Gag-Pol precursor.
Collapse
|
22
|
Enhanced conformational sampling technique provides an energy landscape view of large-scale protein conformational transitions. Phys Chem Chem Phys 2016; 18:29170-29182. [DOI: 10.1039/c6cp05634b] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
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.
Collapse
|
23
|
Pressure-induced structural transition of mature HIV-1 protease from a combined NMR/MD simulation approach. Proteins 2015; 83:2117-23. [PMID: 26385843 DOI: 10.1002/prot.24931] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2015] [Revised: 09/01/2015] [Accepted: 09/01/2015] [Indexed: 11/10/2022]
Abstract
We investigate the pressure-induced structural changes in the mature human immunodeficiency virus type 1 protease dimer, using residual dipolar coupling (RDC) measurements in a weakly oriented solution. (1)DNH RDCs were measured under high-pressure conditions for an inhibitor-free PR and an inhibitor-bound complex, as well as for an inhibitor-free multidrug resistant protease bearing 20 mutations (PR20). While PR20 and the inhibitor-bound PR were little affected by pressure, inhibitor-free PR showed significant differences in the RDCs measured at 600 bar compared with 1 bar. The structural basis of such changes was investigated by MD simulations using the experimental RDC restraints, revealing substantial conformational perturbations, specifically a partial opening of the flaps and the penetration of water molecules into the hydrophobic core of the subunits at high pressure. This study highlights the exquisite sensitivity of RDCs to pressure-induced conformational changes and illustrates how RDCs combined with MD simulations can be used to determine the structural properties of metastable intermediate states on the folding energy landscape.
Collapse
|