The folding and dimerization of HIV-1 protease: evidence for a stable monomer from simulations

Evolution of drug resistance drives destabilization of flap region dynamics in HIV-1 protease

The HIV-1 protease is one of several common key targets of combination drug therapies for human immunodeficiency virus infection and acquired immunodeficiency syndrome. During the progression of the disease, some individual patients acquire drug resistance due to mutational hotspots on the viral proteins targeted by combination drug therapies. It has recently been discovered that drug-resistant mutations accumulate on the "flap region" of the HIV-1 protease, which is a critical dynamic region involved in nonspecific polypeptide binding during invasion and infection of the host cell. In this study, we utilize machine learning-assisted comparative molecular dynamics, conducted at single amino acid site resolution, to investigate the dynamic changes that occur during functional dimerization and drug binding of wild-type and common drug-resistant versions of the main protease. We also use a multiagent machine learning model to identify conserved dynamics of the HIV-1 main protease that are preserved across simian and feline protease orthologs. We find that a key conserved functional site in the flap region, a solvent-exposed isoleucine (Ile50) that controls flap dynamics is functionally targeted by drug resistance mutations, leading to amplified molecular dynamics affecting the functional ability of the flap region to hold the drugs. We conclude that better long-term patient outcomes may be achieved by designing drugs that target protease regions that are less dependent upon single sites with large functional binding effects.

Mechanism of darunavir binding to monomeric HIV-1 protease: a step forward in the rational design of dimerization inhibitors

HIV protease (HIVPR) is a key target in AIDS therapeutics. All ten FDA-approved drugs that compete with substrates in binding to this dimeric enzyme's active site have become ineffective due to the emergence of drug resistant mutants. Blocking the dimerization interface of HIVPR is thus being explored as an alternate strategy. The latest drug, darunavir (DRV), which exhibited a high genetic barrier to viral resistance, is said to have a dual mode of action - (i) binding to the dimeric active site, and (ii) preventing the dimerization by binding to the HIVPR monomer. Despite several reports on DRV complexation with dimeric HIVPR, the mode and mechanism of the binding of DRV to the HIVPR monomer are poorly understood. In this study, we utilized all-atomic MD simulations and umbrella sampling techniques to identify the best possible binding mode of DRV to the monomeric HIVPR and its mechanism of association. The results suggest that DRV binds between the active site and the flap of the monomer, and the flap plays a crucial role in directing the drug to bind and driving the other protein domains to undergo induced fit changes for stronger complexation. The obtained binding mode of DRV was validated by comparing with various mutational data from clinical isolates to reported in vitro mutations. The identified binding pose was also able to successfully reproduce the experimental K_i value in the picomolar range. The residue-level information extracted from this study could accelerate the structure-based drug designing approaches targeting HIVPR dimerization.

The denatured state of HIV-1 protease under native conditions

The denatured state of several proteins has been shown to display transient structures that are relevant for folding, stability, and aggregation. To detect them by nuclear magnetic resonance (NMR) spectroscopy, the denatured state must be stabilized by chemical agents or changes in temperature. This makes the environment different from that experienced in biologically relevant processes. Using high‐resolution heteronuclear NMR spectroscopy, we have characterized several denatured states of a monomeric variant of HIV‐1 protease, which is natively structured in water, induced by different concentrations of urea, guanidinium chloride, and acetic acid. We have extrapolated the chemical shifts and the relaxation parameters to the denaturant‐free denatured state at native conditions, showing that they converge to the same values. Subsequently, we characterized the conformational properties of this biologically relevant denatured state under native conditions by advanced molecular dynamics simulations and validated the results by comparison to experimental data. We show that the denatured state of HIV‐1 protease under native conditions displays rich patterns of transient native and non‐native structures, which could be of relevance to its guidance through a complex folding process.

Experiments and simulation on ZIKV NS2B-NS3 protease reveal its complex folding

Zika virus has been identified in various body fluids such as semen, urine, saliva, cerebrospinal fluid, and vaginal secretion of an infected individual. The pH of these fluids varies from mildly acidic to mildly alkaline. So it is imperative to understand the impact of these conditions on viral protein functioning. We investigated the NS2B-NS3 protease stability and its activity in different denaturing environments. Finding indicates that NS2B-NS3 protease maintains stability at pH 4.8-8.7. Thus it suggests that the complex remains functionally active to hydrolyze the polyprotein within a diverse environmental condition such as variable pH. Despite a stable structure at a broad pH range, a change in environmental conditions dramatically influence its protease activity. Moreover, it is susceptible to structural transformation leading to increased β-strand or helix content in the presence of alcohol. This study may help further to understand the folding-function relationship of the general flaviviral protease complex.

Middle-way flexible docking: Pose prediction using mixed-resolution Monte Carlo in estrogen receptor α

There is a vast gulf between the two primary strategies for simulating protein-ligand interactions. Docking methods significantly limit or eliminate protein flexibility to gain great speed at the price of uncontrolled inaccuracy, whereas fully flexible atomistic molecular dynamics simulations are expensive and often suffer from limited sampling. We have developed a flexible docking approach geared especially for highly flexible or poorly resolved targets based on mixed-resolution Monte Carlo (MRMC), which is intended to offer a balance among speed, protein flexibility, and sampling power. The binding region of the protein is treated with a standard atomistic force field, while the remainder of the protein is modeled at the residue level with a Gō model that permits protein flexibility while saving computational cost. Implicit solvation is used. Here we assess three facets of the MRMC approach with implications for other docking studies: (i) the role of receptor flexibility in cross-docking pose prediction; (ii) the use of non-equilibrium candidate Monte Carlo (NCMC) and (iii) the use of pose-clustering in scoring. We examine 61 co-crystallized ligands of estrogen receptor α, an important cancer target known for its flexibility. We also compare the performance of the MRMC approach with Autodock smina. Adding protein flexibility, not surprisingly, leads to significantly lower total energies and stronger interactions between protein and ligand, but notably we document the important role of backbone flexibility in the improvement. The improved backbone flexibility also leads to improved performance relative to smina. Somewhat unexpectedly, our implementation of NCMC leads to only modestly improved sampling of ligand poses. Overall, the addition of protein flexibility improves the performance of docking, as measured by energy-ranked poses, but we do not find significant improvements based on cluster information or the use of NCMC. We discuss possible improvements for the model including alternative coarse-grained force fields, improvements to the treatment of solvation, and adding additional types of NCMC moves.

Interactions between the Hepatitis C Virus Nonstructural 2 Protein and Host Adaptor Proteins 1 and 4 Orchestrate Virus Release

Hepatitis C virus (HCV) spreads via secreted cell-free particles or direct cell-to-cell transmission. Yet, virus-host determinants governing differential intracellular trafficking of cell-free- and cell-to-cell-transmitted virus remain unknown. The host adaptor proteins (APs) AP-1A, AP-1B, and AP-4 traffic in post-Golgi compartments, and the latter two are implicated in basolateral sorting. We reported that AP-1A mediates HCV trafficking during release, whereas the endocytic adaptor AP-2 mediates entry and assembly. We demonstrated that the host kinases AAK1 and GAK regulate HCV infection by controlling these clathrin-associated APs. Here, we sought to define the roles of AP-4, a clathrin-independent adaptor; AP-1A; and AP-1B in HCV infection. We screened for interactions between HCV proteins and the μ subunits of AP-1A, AP-1B, and AP-4 by mammalian cell-based protein fragment complementation assays. The nonstructural 2 (NS2) protein emerged as an interactor of these adaptors in this screening and by coimmunoprecipitations in HCV-infected cells. Two previously unrecognized dileucine-based motifs in the NS2 C terminus mediated AP binding and HCV release. Infectivity and coculture assays demonstrated that while all three adaptors mediate HCV release and cell-free spread, AP-1B and AP-4, but not AP-1A, mediate cell-to-cell spread. Live-cell imaging revealed HCV cotrafficking with AP-1A, AP-1B, and AP-4 and that AP-4 mediates HCV trafficking in a post-Golgi compartment. Lastly, HCV cell-to-cell spread was regulated by AAK1 and GAK and thus susceptible to treatment with AAK1 and GAK inhibitors. These data provide a mechanistic understanding of HCV trafficking in distinct release pathways and reveal a requirement for APs in cell-to-cell viral spread.

HCV spreads via cell-free infection or cell-to-cell contact that shields it from antibody neutralization, thereby facilitating viral persistence. Yet, factors governing this differential sorting remain unknown. By integrating proteomic, RNA interference, genetic, live-cell imaging, and pharmacological approaches, we uncover differential coopting of host adaptor proteins (APs) to mediate HCV traffic at distinct late steps of the viral life cycle. We reported that AP-1A and AP-2 mediate HCV trafficking during release and assembly, respectively. Here, we demonstrate that dileucine motifs in the NS2 protein mediate AP-1A, AP-1B, and AP-4 binding and cell-free virus release. Moreover, we reveal that AP-4, an adaptor not previously implicated in viral infections, mediates cell-to-cell spread and HCV trafficking. Lastly, we demonstrate cell-to-cell spread regulation by AAK1 and GAK, host kinases controlling APs, and susceptibility to their inhibitors. This study provides mechanistic insights into virus-host determinants that facilitate HCV trafficking, with potential implications for pathogenesis and antiviral agent design.

Mechanism of the formation of the RecA-ssDNA nucleoprotein filament structure: a coarse-grained approach

In prokaryotes, the RecA protein catalyzes the repair and strand exchange of double-stranded DNA. RecA binds to single-stranded DNA (ssDNA) and forms a presynaptic complex in which the protein polymerizes around the ssDNA to form a right-handed helical nucleoprotein filament structure. In the present work, the mechanism for the formation of the RecA-ssDNA filament structure is modeled using coarse-grained molecular dynamics simulations. Information from the X-ray structure was used to model the protein itself but not its interactions; the interactions between the protein and the ssDNA were modeled solely by electrostatic, aromatic, and repulsive energies. For the present study, the monomeric, dimeric, and trimeric units of RecA and 4, 8, and 11 NT-long ssDNA, respectively, were studied. Our results indicate that monomeric RecA is not sufficient for nucleoprotein filament formation; rather, dimeric RecA is the elementary binding unit, with higher multimeric units of RecA facilitating filament formation. Our results reveal that loop region flexibility at the primary binding site of RecA is essential for it to bind the incoming ssDNA, that the aromatic residues present in the loop region play an important role in ssDNA binding, and that ATP may play a role in guiding the ssDNA by changing the electrostatic potential of the RecA protein.

Evolution under Drug Pressure Remodels the Folding Free-Energy Landscape of Mature HIV-1 Protease

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.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Tabulation as a high-resolution alternative to coarse-graining protein interactions: Initial application to virus capsid subunits

Traditional coarse-graining based on a reduced number of interaction sites often entails a significant sacrifice of chemical accuracy. As an alternative, we present a method for simulating large systems composed of interacting macromolecules using an energy tabulation strategy previously devised for small rigid molecules or molecular fragments [S. Lettieri and D. M. Zuckerman, J. Comput. Chem. 33, 268-275 (2012); J. Spiriti and D. M. Zuckerman, J. Chem. Theory Comput. 10, 5161-5177 (2014)]. We treat proteins as rigid and construct distance and orientation-dependent tables of the interaction energy between them. Arbitrarily detailed interactions may be incorporated into the tables, but as a proof-of-principle, we tabulate a simple α-carbon Gō-like model for interactions between dimeric subunits of the hepatitis B viral capsid. This model is significantly more structurally realistic than previous models used in capsid assembly studies. We are able to increase the speed of Monte Carlo simulations by a factor of up to 6700 compared to simulations without tables, with only minimal further loss in accuracy. To obtain further enhancement of sampling, we combine tabulation with the weighted ensemble (WE) method, in which multiple parallel simulations are occasionally replicated or pruned in order to sample targeted regions of a reaction coordinate space. In the initial study reported here, WE is able to yield pathways of the final ∼25% of the assembly process.

Protein folding: understanding the role of water and the low Reynolds number environment as the peptide chain emerges from the ribosome and folds

The mechanism of protein folding during early stages of the process has three determinants. First, moving water molecules obey the rules of low Reynolds number physics without an inertial component. Molecular movement is instantaneous and size insensitive. Proteins emerging from the ribosome move and rotate without an external force if they change shape, forming and propagating helical structures that increases translocational efficiency. Forward motion ceases when the shape change or propelling force ceases. Second, application of quantum field theory to water structure predicts the spontaneous formation of low density coherent units of fixed size that expel dissolved atmospheric gases. Structured water layers with both coherent and non-coherent domains, form a sheath around the new protein. The surface of exposed hydrophobic amino acids is protected from water contact by small nanobubbles of dissolved atmospheric gases, 5 or 6 molecules on average, that vibrate, attracting even widely separated resonating nanobubbles. This force results from quantum effects, appearing only when the system is within and interacts with an oscillating electromagnetic field. The newly recognized quantum force sharply bends the peptide and is part of a dynamic field determining the pathway of protein folding. Third, the force initiating the tertiary folding of proteins arises from twists at the position of each hydrophobic amino acid, that minimizes surface exposure of the hydrophobic amino acids and propagates along the protein. When the total bend reaches 360°, the leading segment of water sheath intersects the trailing segment. This steric self-intersection expels water from overlapping segments of the sheath and by Newton׳s second law moves the polypeptide chain in an opposite direction. Consequently, with very few exceptions that we enumerate and discuss, tertiary structures are absent from proteins without hydrophobic amino acids, which control the early stages of protein folding and the overall shape of protein. Consequently, proteins only adopt a limited number of forms. The formation of quaternary structures is not necessarily prevented by the absence of hydrophobic amino acids.

Toward in silico biomolecular manipulation through static modes: atomic scale characterization of HIV-1 protease flexibility

Probing biomolecular flexibility with atomic-scale resolution is a challenging task in current computational biology for fundamental understanding and prediction of biomolecular interactions and associated functions. This paper makes use of the static mode method to study HIV-1 protease considered as a model system to investigate the full biomolecular flexibility at the atomic scale, the screening of active site biomechanical properties, the blind prediction of allosteric sites, and the design of multisite strategies to target deformations of interest. Relying on this single calculation run of static modes, we demonstrate that in silico predictive design of an infinite set of complex excitation fields is reachable, thanks to the storage of the static modes in a data bank that can be used to mimic single or multiatom contact and efficiently anticipate conformational changes arising from external stimuli. All along this article, we compare our results to data previously published and propose a guideline for efficient, predictive, and custom in silico experiments.

Natural polymorphisms and unusual mutations in HIV-1 protease with potential antiretroviral resistance: a bioinformatic analysis

Background

The correlations of genotypic and phenotypic tests with treatment, clinical history and the significance of mutations in viruses of HIV-infected patients are used to establish resistance mutations to protease inhibitors (PIs). Emerging mutations in human immunodeficiency virus type 1 (HIV-1) protease confer resistance to PIs by inducing structural changes at the ligand interaction site. The aim of this study was to establish an in silico structural relationship between natural HIV-1 polymorphisms and unusual HIV-1 mutations that confer resistance to PIs.

Results

Protease sequences isolated from 151 Mexican HIV-1 patients that were naïve to, or subjected to antiretroviral therapy, were examined. We identified 41 unrelated resistance mutations with a prevalence greater than 1%. Among these mutations, nine exhibited positive selection, three were natural polymorphisms (L63S/V/H) in a codon associated with drug resistance, and six were unusual mutations (L5F, D29V, L63R/G, P79L and T91V). The D29V mutation, with a prevalence of 1.32% in the studied population, was only found in patients treated with antiretroviral drugs. Using in silico modelling, we observed that D29V formed unstable protease complexes when were docked with lopinavir, saquinavir, darunavir, tipranavir, indinavir and atazanavir.

Conclusions

The structural correlation of natural polymorphisms and unusual mutations with drug resistance is useful for the identification of HIV-1 variants with potential resistance to PIs. The D29V mutation likely confers a selection advantage in viruses; however, in silico, presence of this mutation results in unstable enzyme/PI complexes, that possibly induce resistance to PIs.

Mapping conformational dynamics of proteins using torsional dynamics simulations

All-atom molecular dynamics simulations are widely used to study the flexibility of protein conformations. However, enhanced sampling techniques are required for simulating protein dynamics that occur on the millisecond timescale. In this work, we show that torsional molecular dynamics simulations enhance protein conformational sampling by performing conformational search in the low-frequency torsional degrees of freedom. In this article, we use our recently developed torsional-dynamics method called Generalized Newton-Euler Inverse Mass Operator (GNEIMO) to study the conformational dynamics of four proteins. We investigate the use of the GNEIMO method in simulations of the conformationally flexible proteins fasciculin and calmodulin, as well as the less flexible crambin and bovine pancreatic trypsin inhibitor. For the latter two proteins, the GNEIMO simulations with an implicit-solvent model reproduced the average protein structural fluctuations and sample conformations similar to those from Cartesian simulations with explicit solvent. The application of GNEIMO with replica exchange to the study of fasciculin conformational dynamics produced sampling of two of this protein's experimentally established conformational substates. Conformational transition of calmodulin from the Ca(2+)-bound to the Ca(2+)-free conformation occurred readily with GNEIMO simulations. Moreover, the GNEIMO method generated an ensemble of conformations that satisfy about half of both short- and long-range interresidue distances obtained from NMR structures of holo to apo transitions in calmodulin. Although unconstrained all-atom Cartesian simulations have failed to sample transitions between the substates of fasciculin and calmodulin, GNEIMO simulations show the transitions in both systems. The relatively short simulation times required to capture these long-timescale conformational dynamics indicate that GNEIMO is a promising molecular-dynamics technique for studying domain motion in proteins.

The maturation of HIV-1 protease precursor studied by discrete molecular dynamics

The equilibrium properties of a HIV-1-protease precursor are studied by means of an efficient molecular dynamics scheme, which allows for the simulation of the folding of the protein monomers and their dimerization into an active form and compare them with those of the mature protein. The results of the model provide, with atomic detail, an overall account of several experimental findings, including the NMR conformation of the mature dimer, the calorimetric properties of the system, the effects of the precursor tail on the dimerization constant, the secondary chemical shifts of the monomer, and the paramagnetic relaxation enhancement data associated with the conformations of the precursor. It is found that although the mature protein can dimerize in a unique, single way, the precursor populates several dimeric conformations in which monomers are always native-like, but their binding can be non-native.

Folding and association of a homotetrameric protein complex in an all-atom Go model

The 84-residue homotetrameric BBAT1 is one of the smallest stable protein complexes and therefore is a good test system to study the self-assembly of multimeric proteins. We have researched for this protein the interplay between the folding of monomers and their assembly into tetramers. Replica exchange molecular dynamics simulations relying on a Go model are compared with earlier simulations that use the physics-based coarse-grained UNRES model.

Experiments and comprehensive simulations of the formation of a helical turn

We investigate the kinetics and thermodynamics of a helical turn formation in the peptide Ac-WAAAH-NH(2). NMR measurements indicate that this peptide has significant tendency to form a structure of a helical turn, while temperature dependent CD establishes the helix fraction at different temperatures. Molecular dynamics and milestoning simulations agree with experimental observables and suggest an atomically detailed picture for the turn formation. Using a network representation, two alternative mechanisms of folding are identified: (i) a direct co-operative mechanism from the unfolded to the folded state without intermediate formation of hydrogen bonds and (ii) an indirect mechanism with structural intermediates with two residues in a helical conformation. This picture is consistent with kinetic measurements that reveal two experimental time scales of sub-nanosecond and several nanoseconds.

How Does Darunavir Prevent HIV-1 Protease Dimerization?

The drug Darunavir (DRV) is a potent inhibitor of HIV-1 protease (PR), a homodimeric essential enzyme of the AIDS virus. Recent experimental data suggest that DRV is able to prevent dimerization of HIV-1 PR, which, together with its high affinity for the mature enzyme, has been linked to the high genetic barrier to the development of viral resistance. The mechanism of dimerization inhibition and the binding mode(s) of DRV to monomeric HIV-1 PR are unknown. Here, multiple molecular dynamics simulations with explicit solvent (for a total of 11 μs with the CHARMM force field and 1 μs with the Amber force field) show that the monomer of HIV-1 PR is structurally stable and reveal a major binding mode of DRV. This binding mode is stabilized by favorable interactions between the apolar groups of DRV and the hydrophobic residues Ile32, Ile47, Ile50, Ile54, Pro79, Val82, and Ile84. The binding mode to monomeric HIV-1 PR identified by molecular dynamics is different from the two binding modes observed in the crystal structure of the complex with dimeric HIV-1 PR. As an example, there are no interactions between DRV and the catalytic Asp25 in the binding mode to monomeric HIV-1 PR revelead by the simulations. In contrast, the simulations show extensive and stable interactions between DRV and the flap (residues 46-55), which are likely to sterically hinder the formation of the flap interface as observed in the dimeric structure. Which of the two mechanisms of inhibition (dimerization inhibition by association with the flap or binding to the active site of the mature enzyme) dominates might depend on the HIV-1 PR mutations, and it is likely that dimerization inhibition is predominant for multiple mutations at the active site in multidrug resistant strains.

Visualization of early events in acetic acid denaturation of HIV-1 protease: a molecular dynamics study

Protein denaturation plays a crucial role in cellular processes. In this study, denaturation of HIV-1 Protease (PR) was investigated by all-atom MD simulations in explicit solvent. The PR dimer and monomer were simulated separately in 9 M acetic acid (9 M AcOH) solution and water to study the denaturation process of PR in acetic acid environment. Direct visualization of the denaturation dynamics that is readily available from such simulations has been presented. Our simulations in 9 M AcOH reveal that the PR denaturation begins by separation of dimer into intact monomers and it is only after this separation that the monomer units start denaturing. The denaturation of the monomers is flagged off by the loss of crucial interactions between the α-helix at C-terminal and surrounding β-strands. This causes the structure to transit from the equilibrium dynamics to random non-equilibrating dynamics. Residence time calculations indicate that denaturation occurs via direct interaction of the acetic acid molecules with certain regions of the protein in 9 M AcOH. All these observations have helped to decipher a picture of the early events in acetic acid denaturation of PR and have illustrated that the α-helix and the β-sheet at the C-terminus of a native and functional PR dimer should maintain both the stability and the function of the enzyme and thus present newer targets for blocking PR function.

Detailed atomistic analysis of the HIV-1 protease interface

HIV-1 protease is a very attractive target for the development of new anti-HIV drugs and has been extensively studied over the past decades. In this study, we present a detailed atomic level characterization of the dimer interface in the enzyme HIV-1 protease through computational alanine scanning mutagenesis and molecular dynamics simulations. In addition to a full mapping of the amino acid residues present at the subunit interface, in terms of the corresponding energetic contribution for dimer formation and of their classification as hot spots, warm spots, and null spots, we trace a dynamic analysis of the subunit interacting and solvent accessible surface areas and of the most important hydrogen bonds between subunits. The results presented illustrate the high energetic importance for dimer formation of a small set of five amino acid residue pairs at the subunit interface-Leu5, Ile50, Arg87, Leu97, and Phe99-and provide important clues on the most important structural and energetic determinants for dimer formation. In addition, the results presented suggest several key targets at the subunit interface for the development of new molecules that aim to inhibit HIV-1 protease (PR) activity through blocking the formation of the fully active PR homodimeric form, providing important clues for drug design.

Interdomain communication revealed in the diabetes drug target mitoNEET

MitoNEET is a recently identified drug target for a commonly prescribed diabetes drug, Pioglitazone. It belongs to a previously uncharacterized ancient family of proteins for which the hallmark is the presence of a unique 39 amino acid CDGSH domain. In order to characterize the folding landscape of this novel fold, we performed thermodynamic simulations on MitoNEET using a structure-based model. Additionally, we implement a method of contact map clustering to partition out alternate pathways in folding. This cluster analysis reveals a detour late in folding and enables us to carefully examine the folding mechanism of each pathway rather than the macroscopic average. We observe that tightness in a region distal to the iron-sulfur cluster creates a constraint in folding and additionally appears to mediate communication in folding between the two domains of the protein. We demonstrate that by making changes at this site we are able to tweak the order of folding events in the cluster binding domain as well as decrease the barrier to folding.

Autocatalytic maturation, physical/chemical properties, and crystal structure of group N HIV-1 protease: relevance to drug resistance

The mature protease from Group N human immunodeficiency virus Type 1 (HIV-1) (PR1(N)) differs in 20 amino acids from the extensively studied Group M protease (PR1(M)) at positions corresponding to minor drug-resistance mutations (DRMs). The first crystal structure (1.09 Å resolution) of PR1(N) with the clinical inhibitor darunavir (DRV) reveals the same overall structure as PR1(M), but with a slightly larger inhibitor-binding cavity. Changes in the 10s loop and the flap hinge propagate to shift one flap away from the inhibitor, whereas L89F and substitutions in the 60s loop perturb inhibitor-binding residues 29-32. However, kinetic parameters of PR1(N) closely resemble those of PR1(M), and calorimetric results are consistent with similar binding affinities for DRV and two other clinical PIs, suggesting that minor DRMs coevolve to compensate for the detrimental effects of drug-specific major DRMs. A miniprecursor (TFR 1-61-PR1(N)) comprising the transframe region (TFR) fused to the N-terminus of PR1(N) undergoes autocatalytic cleavage at the TFR/PR1(N) site concomitant with the appearance of catalytic activity characteristic of the dimeric, mature enzyme. This cleavage is inhibited at an equimolar ratio of precursor to DRV (∼6 μM), which partially stabilizes the precursor dimer from a monomer. However, cleavage at L34/W35 within the TFR, which precedes the TFR 1-61/PR1(N) cleavage at pH ≤ 5, is only partially inhibited. Favorable properties of PR1(N) relative to PR1(M) include its suitability for column fractionation by size under native conditions and >10-fold higher dimer dissociation constant (150 nM). Exploiting these properties may facilitate testing of potential dimerization inhibitors that perturb early precursor processing steps.

Multiple routes and milestones in the folding of HIV-1 protease monomer

Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV–1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV–1 protease.

Post-translational modification of glutamine and lysine residues of HIV-1 aspartyl protease by transglutaminase increases its catalytic activity

The human immunodeficiency virus type 1 aspartyl protease (HIV-1 PR) is a homodimeric aspartyl endopeptidase that is required for virus replication. HIV-1 PR was shown to act invitro as acyl-donor and -acceptor for both guinea pig liver transglutaminase (TG, EC 2.3.2.13) and human Factor XIIIa. These preliminary evidences suggested that the HIV-1 PR contains at least three TG-reactive glutaminyl and one lysyl residues. We report here that the incubation of HIV-1 PR with TG increases its catalytic activity. This increase is dependent upon the time of incubation, the concentration of TG and the presence of Ca2+. Identification of epsilon-(gamma-glutamyl)lysine in the proteolytic digest of the TG-modified HIV-1 PR suggested intramolecular covalent cross-linking of this protease which may promote a non-covalent dimerization and subsequent activation of this enzyme via a conformational change. This hypothesis is supported by the observation that the TG-catalyzed activation of HIV-1 PR was completely abolished by spermidine (SPD) which acts as a competitive inhibitor of epsilon-(gamma-glutamyl)lysine formation. Indeed, in the presence of 1mM SPD the formation of the isopeptide was decreased of about 80%. The main products of the TG-catalyzed modification of HIV-1 PR in the presence of SPD were N(1)-mono(gamma-glutamyl)SPD and N(8)-mono(gamma-glutamyl)SPD. Negligible amount of N(1),N(8)-bis(gamma-glutamyl)SPD were found. The significance of these results is discussed with respect to the activation of the protease by post-translational modification and design of potential inhibitors.

Capillarity theory for the fly-casting mechanism

Biomolecular folding and function are often coupled. During molecular recognition events, one of the binding partners may transiently or partially unfold, allowing more rapid access to a binding site. We describe a simple model for this fly-casting mechanism based on the capillarity approximation and polymer chain statistics. The model shows that fly casting is most effective when the protein unfolding barrier is small and the part of the chain which extends toward the target is relatively rigid. These features are often seen in known examples of fly casting in protein-DNA binding. Simulations of protein-DNA binding based on well-funneled native-topology models with electrostatic forces confirm the trends of the analytical theory.

The folding free-energy surface of HIV-1 protease: insights into the thermodynamic basis for resistance to inhibitors

Spontaneous mutations at numerous sites distant from the active site of human immunodeficiency virus type 1 protease enable resistance to inhibitors while retaining enzymatic activity. As a benchmark for probing the effects of these mutations on the conformational adaptability of this dimeric beta-barrel protein, the folding free-energy surface of a pseudo-wild-type variant, HIV-PR(*), was determined by a combination of equilibrium and kinetic experiments on the urea-induced unfolding/refolding reactions. The equilibrium unfolding reaction was well described by a two-state model involving only the native dimeric form and the unfolded monomer. The global analysis of the kinetic folding mechanism reveals the presence of a fully folded monomeric intermediate that associates to form the native dimeric structure. Independent analysis of a stable monomeric version of the protease demonstrated that a small-amplitude fluorescence phase in refolding and unfolding, not included in the global analysis of the dimeric protein, reflects the presence of a transient intermediate in the monomer folding reaction. The partially folded and fully folded monomers are only marginally stable with respect to the unfolded state, and the dimerization reaction provides a modest driving force at micromolar concentrations of protein. The thermodynamic properties of this system are such that mutations can readily shift the equilibrium from the dimeric native state towards weakly folded states that have a lower affinity for inhibitors but that could be induced to bind to their target proteolytic sites. Presumably, subsequent secondary mutations increase the stability of the native dimeric state in these variants and, thereby, optimize the catalytic properties of the resistant human immunodeficiency virus type 1 protease.

Toolbox for analyzing finite two-state trajectories

In many experiments, the aim is to deduce an underlying multisubstate on-off kinetic scheme (KS) from the statistical properties of a two-state trajectory. However, a two-state trajectory that is generated from an on-off KS contains only partial information about the KS, and so, in many cases, more than one KS can be associated with the data. We recently showed that the optimal way to solve this problem is to use canonical forms of reduced dimensions (RDs). RD forms are on-off networks with connections only between substates of different states, where the connections can have nonexponential waiting time probability density functions (WT-PDFs). In theory, only a single RD form can be associated with the data. To utilize RD forms in the analysis of the data, a RD form should be associated with the data. Here, we give a toolbox for building a RD form from a finite time, noiseless, two-state trajectory. The methods in the toolbox are based on known statistical methods in data analysis, combined with statistical methods and numerical algorithms designed specifically for the current problem. Our toolbox is self-contained-it builds a mechanism based only on the information it extracts from the data, and its implementation is fast (analyzing a 10;{6}cycle trajectory from a 30-parameter mechanism takes a couple of hours on a PC with a 2.66GHz processor). The toolbox is automated and is freely available for academic research upon electronic request.

Characterization of the unfolded state of repeat proteins

The unfolded state ensemble of proteins has been described as a structurally featureless state. While this approach is supported by the fact that many unfolded proteins follow the scaling law behavior of a random coil, there is evidence that the unfolded states of various proteins are stabilized by native or non-native interactions. Recently, the existence of extensive non-native structure was reported for a repeat protein, which resulted in a scaling law exponent that is significantly smaller than that of a random polymer [Cortajarena et al., J. Mol. Biol. 382(1), 203-212 (2008)]. It was concluded that the high compactness of this protein stems from a significant fraction of interacting PP(II) helical segments in the unfolded state. In this study, we aim at providing possible molecular understanding of this anomalous compactness of the unfolded state and to investigate its origin. Using a hierarchy of computational models, we ask whether in general the unfolded state of a repeat protein is likely to be intrinsically more compact than the unfolded state of globular proteins, or whether this phenomenon depends mostly on the occurrence of a specific sequence that promotes PP(II) conformations. Our results suggest that the formation of the PP(II) conformation is indeed essential, yet the recurring sequence of repeat proteins promotes the interactions between these PP(II) segments and the formation of non-native interactions in the unfolded state.

Folding of elongated proteins: conventional or anomalous?

The complexity of the mechanisms by which proteins fold has been shown by many studies to be governed by their native-state topologies. This was manifested in the ability of the native topology-based model to capture folding mechanisms and the success of folding rate predictions based on various topological measures, such as the contact order. However, while the finer details of topological complexity have been thoroughly examined and related to folding kinetics, simpler characteristics of the protein, such as its overall shape, have been largely disregarded. In this study, we investigated the folding of proteins with an unusual elongated geometry that differs substantially from the common globular structure. To study the effect of the elongation degree on the folding kinetics, we used repeat proteins, which become more elongated as they include more repeating units. Some of these have apparently anomalous experimental folding kinetics, with rates that are often less than expected on the basis of rates for globular proteins possessing similar topological complexity. Using experimental folding rates and a larger set of rates obtained from simulations, we have shown that as the protein becomes increasingly elongated, its folding kinetics becomes slower and deviates more from the rate expected on the basis of topology measures fitted for globular proteins. The observed slow kinetics is a result of a more complex pathway in which stable intermediates composed of several consecutive repeats can appear. We thus propose a novel measure, an elongation-sensitive contact order, that takes into account both the extent of elongation and the topological complexity of the protein. This new measure resolves the apparent discrimination between the folding of globular and elongated repeat proteins. Our study extends the current capabilities of folding-rate predictions by unifying the kinetics of repeat and globular proteins.

Atomistic simulations of the HIV-1 protease folding inhibition

Biochemical experiments have recently revealed that the p-S8 peptide, with an amino-acid sequence identical to the conserved fragment 83-93 (S8) of the HIV-1 protease, can inhibit catalytic activity of the enzyme by interfering with protease folding and dimerization. In this study, we introduce a hierarchical modeling approach for understanding the molecular basis of the HIV-1 protease folding inhibition. Coarse-grained molecular docking simulations of the flexible p-S8 peptide with the ensembles of HIV-1 protease monomers have revealed structurally different complexes of the p-S8 peptide, which can be formed by targeting the conserved segment 24-34 (S2) of the folding nucleus (folding inhibition) and by interacting with the antiparallel termini beta-sheet region (dimerization inhibition). All-atom molecular dynamics simulations of the inhibitor complexes with the HIV-1 PR monomer have been independently carried out for the predicted folding and dimerization binding modes of the p-S8 peptide, confirming the thermodynamic stability of these complexes. Binding free-energy calculations of the p-S8 peptide and its active analogs are then performed using molecular dynamics trajectories of the peptide complexes with the HIV-1 PR monomers. The results of this study have provided a plausible molecular model for the inhibitor intervention with the HIV-1 PR folding and dimerization and have accurately reproduced the experimental inhibition profiles of the active folding inhibitors.

Conformation of a coarse-grained protein chain (an aspartic acid protease) model in effective solvent by a bond-fluctuating Monte Carlo simulation

In a coarse-grained description of a protein chain, all of the 20 amino acid residues can be broadly divided into three groups: Hydrophobic (H) , polar (P) , and electrostatic (E) . A protein can be described by nodes tethered in a chain with a node representing an amino acid group. Aspartic acid protease consists of 99 residues in a well-defined sequence of H , P , and E nodes tethered together by fluctuating bonds. The protein chain is placed on a cubic lattice where empty lattice sites constitute an effective solvent medium. The amino groups (nodes) interact with the solvent (S) sites with appropriate attractive (PS) and repulsive (HS) interactions with the solvent and execute their stochastic movement with the Metropolis algorithm. Variations of the root mean square displacements of the center of mass and that of its center node of the protease chain and its gyration radius with the time steps are examined for different solvent strength. The structure of the protease swells on increasing the solvent interaction strength which tends to enhance the relaxation time to reach the diffusive behavior of the chain. Equilibrium radius of gyration increases linearly on increasing the solvent strength: A slow rate of increase in weak solvent regime is followed by a faster swelling in stronger solvent. Variation of the gyration radius with the time steps suggests that the protein chain moves via contraction and expansion in a somewhat quasiperiodic pattern particularly in strong solvent.

Molecular dynamics simulations of HIV‐1 protease monomer: Assembly of N‐terminus and C‐terminus into β‐sheet in water solution

HIV-1 protease (HIV-PR) consists of two identical subunits that are united together through a four-stranded antiparallel beta-sheet formed of the peptide termini of each monomer. Since the active site exists only in the dimer, a strategy that is attracting more and more attention in inhibitor design and which may overcome the serious drug resistance caused by competitive inhibitors is to block the peptide termini of the monomer, thereby interfering with formation of the active dimer. In the present work, we performed several extensive molecular dynamics (MD) simulations of the HIV-PR monomer in water to illustrate its solvated conformation and dynamics behavior. We found that the peptide termini usually assembled into beta-sheet after several nanoseconds' simulation, and became much less flexible. This beta-sheet is stabilized by intramolecular interactions and is not easily disaggregated under the present MD simulation conditions. This transformation may be an important transition during the relaxing and equilibrating of the HIV-PR monomer in aqueous solution, and the terminal beta-sheet may be one of the major conformations of the solvated HIV-PR monomer termini in water. This work may provide new insights into the dynamics behavior and dimerization mechanism of HIV-PR, and more significantly, offer a more rational receptor model for the design and discovery of novel dimerization inhibitors than crystalline structures.

Mapping folding energy landscapes with theory and experiment

The detailed characterization of the overall free energy landscape associated with the folding process of a protein is the ultimate goal in protein folding studies. Modern experimental techniques and all-atom simulations provide a way to obtain accurate thermodynamic and kinetic measurements, but they are oftentimes restricted to probe limited regions of a protein landscape. Although simplified protein models can access larger regions of the landscape, they are built on assumptions and approximations that can affect the accuracy of the results. We review here recent promising approaches that allow to combine the complementary strengths of theory and experiment for a more complete characterization of a protein folding landscape at multiple resolutions. Recent results and possible applications are discussed.

Large multiprotein structures modeling and simulation: the need for mesoscopic models

Recent observational techniques based upon confocal microscopy make it possible to observe cells at a scale that has never been probed before: the mesoscopic scale. In the eukaryotic cell nucleus, many objects demonstrating phenomena occurring at this scale, such as nuclear bodies, are current subjects of investigations. But from a modeling perspective, this scale has not been widely explored, and hence there is a lack of suitable models for such studies. By reviewing higher and lower scale modeling techniques, we analyze their relevance in the context of mesoscale phenomena. We emphasize important characteristics that should be included in a mesoscopic model: an explicit continuous three-dimensional space with discrete simplified molecules that still have the characteristics of steric volume exclusion and realistic distant interaction forces. Then we present 3DSPI, a model dedicated to studies of nuclear bodies based on a simple formalism inspired from molecular dynamics and coarse-grained models: particles interacting through a potential energy function and driven by an overdamped Langevin equation. Finally, we present the features expected to be included in the model, pointing out the difficulties that might arise.

In silico protein fragmentation reveals the importance of critical nuclei on domain reassembly

Protein complementation assays (PCAs) based on split protein fragments have become powerful tools that facilitate the study and engineering of intracellular protein-protein interactions. These assays are based on the observation that a given protein can be split into two inactive fragments and these fragments can reassemble into the original properly folded and functional structure. However, one experimentally observed limitation of PCA systems is that the folding of a protein from its fragments is dramatically slower relative to that of the unsplit parent protein. This is due in part to a poor understanding of how PCA design parameters such as split site position in the primary sequence and size of the resulting fragments contribute to the efficiency of protein reassembly. We used a minimalist on-lattice model to analyze how the dynamics of the reassembly process for two model proteins was affected by the location of the split site. Our results demonstrate that the balanced distribution of the "folding nucleus," a subset of residues that are critical to the formation of the transition state leading to productive folding, between protein fragments is key to their reassembly.