GenoTriplo: A SNP genotype calling method for triploids.

Triploidy is very useful in both aquaculture and some cultivated plants as the induced sterility helps to enhance growth and product quality, as well as acting as a barrier against the contamination of wild populations by escapees. To use genetic information from triploids for academic or breeding purposes, an efficient and robust method to genotype triploids is needed. We developed such a method for genotype calling from SNP arrays, and we implemented it in the R package named GenoTriplo. Our method requires no prior information on cluster positions and remains unaffected by shifted luminescence signals. The method relies on starting the clustering algorithm with an initial higher number of groups than expected from the ploidy level of the samples, followed by merging groups that are too close to each other to be considered as distinct genotypes. Accurate classification of SNPs is achieved through multiple thresholds of quality controls. We compared the performance of GenoTriplo with that of fitPoly, the only published method for triploid SNP genotyping with a free software access. This was assessed by comparing the genotypes generated by both methods for a dataset of 1232 triploid rainbow trout genotyped for 38,033 SNPs. The two methods were consistent for 89% of the genotypes, but for 26% of the SNPs, they exhibited a discrepancy in the number of different genotypes identified. For these SNPs, GenoTriplo had >95% concordance with fitPoly when fitPoly genotyped better. On the contrary, when GenoTriplo genotyped better, fitPoly had less than 50% concordance with GenoTriplo. GenoTriplo was more robust with less genotyping errors. It is also efficient at identifying low-frequency genotypes in the sample set. Finally, we assessed parentage assignment based on GenoTriplo genotyping and observed significant differences in mismatch rates between the best and second-best couples, indicating high confidence in the results. GenoTriplo could also be used to genotype diploids as well as individuals with higher ploidy level by adjusting a few input parameters.

Structure elucidation of the elusive Enzyme I monomer reveals the molecular mechanisms linking oligomerization and enzymatic activity.

Enzyme I (EI) is a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate. This reaction initiates a five-step phosphorylation cascade in the bacterial phosphotransferase (PTS) transduction pathway. Under physiological conditions, EI exists in an equilibrium between a functional dimer and an inactive monomer. The monomer-dimer equilibrium is a crucial factor regulating EI activity and the phosphorylation state of the overall PTS. Experimental studies of EI's monomeric state have yet been hampered by the dimer's high thermodynamic stability, which prevents its characterization by standard structural techniques. In this study, we modified the dimerization domain of EI (EIC) by mutating three amino acids involved in the formation of intersubunit salt bridges. The engineered variant forms an active dimer in solution that can bind and hydrolyze PEP. Using hydrostatic pressure as an additional perturbation, we were then able to study the complete dissociation of the variant from 1 bar to 2.5 kbar in the absence and the presence of EI natural ligands. Backbone residual dipolar couplings collected under high-pressure conditions allowed us to determine the conformational ensemble of the isolated EIC monomeric state in solution. Our calculations reveal that three catalytic loops near the dimerization interface become unstructured upon monomerization, preventing the monomeric enzyme from binding its natural substrate. This study provides an atomic-level characterization of EI's monomeric state and highlights the role of the catalytic loops as allosteric connectors controlling both the activity and oligomerization of the enzyme.

Intricate coupling between the transactivation and basic-leucine zipper domains governs phosphorylation of transcription factor ATF4 by casein kinase 2.

Most transcription factors possess at least one long intrinsically disordered transactivation domain that binds to a variety of coactivators and corepressors and plays a key role in modulating the transcriptional activity. Despite the crucial importance of these domains, the structural and functional basis of transactivation remains poorly understood. Here, we focused on activating transcription factor 4 (ATF4)/cAMP response element-binding protein-2, an essential transcription factor for cellular stress adaptation. Bioinformatic sequence analysis of the ATF4 transactivation domain sequence revealed that the first 125 amino acids have noticeably less propensity for structural disorder than the rest of the domain. Using solution nuclear magnetic resonance spectroscopy complemented by a range of biophysical methods, we found that the isolated transactivation domain is predominantly yet not fully disordered in solution. We also observed that a short motif at the N-terminus of the transactivation domain has a high helical propensity. Importantly, we found that the N-terminal region of the transactivation domain is involved in transient long-range interactions with the basic-leucine zipper domain involved in DNA binding. Finally, in vitro phosphorylation assays with the casein kinase 2 show that the presence of the basic-leucine zipper domain is required for phosphorylation of the transactivation domain. This study uncovers the intricate coupling existing between the transactivation and basic-leucine zipper domains of ATF4, highlighting its potential regulatory significance.

MbnC Is Not Required for the Formation of the N-Terminal Oxazolone in the Methanobactin from Methylosinus trichosporium OB3b.

Methanobactins (MBs) are ribosomally synthesized and posttranslationally modified peptides (RiPPs) produced by methanotrophs for copper uptake. The posttranslational modification that defines MBs is the formation of two heterocyclic groups with associated thioamines from X-Cys dipeptide sequences. Both heterocyclic groups in the MB from Methylosinus trichosporium OB3b (MB-OB3b) are oxazolone groups. The precursor gene for MB-OB3b is mbnA, which is part of a gene cluster that contains both annotated and unannotated genes. One of those unannotated genes, mbnC, is found in all MB operons and, in conjunction with mbnB, is reported to be involved in the formation of both heterocyclic groups in all MBs. To determine the function of mbnC, a deletion mutation was constructed in M. trichosporium OB3b, and the MB produced from the ΔmbnC mutant was purified and structurally characterized by UV-visible absorption spectroscopy, mass spectrometry, and solution nuclear magnetic resonance (NMR) spectroscopy. MB-OB3b from the ΔmbnC mutant was missing the C-terminal Met and was also found to contain a Pro and a Cys in place of the pyrrolidinyl-oxazolone-thioamide group. These results demonstrate MbnC is required for the formation of the C-terminal pyrrolidinyl-oxazolone-thioamide group from the Pro-Cys dipeptide, but not for the formation of the N-terminal 3-methylbutanol-oxazolone-thioamide group from the N-terminal dipeptide Leu-Cys. IMPORTANCE A number of environmental and medical applications have been proposed for MBs, including bioremediation of toxic metals and nanoparticle formation, as well as the treatment of copper- and iron-related diseases. However, before MBs can be modified and optimized for any specific application, the biosynthetic pathway for MB production must be defined. The discovery that mbnC is involved in the formation of the C-terminal oxazolone group with associated thioamide but not for the formation of the N-terminal oxazolone group with associated thioamide in M. trichosporium OB3b suggests the enzymes responsible for posttranslational modification(s) of the two oxazolone groups are not identical.

Thermodynamic stability of hnRNP A1 low complexity domain revealed by high-pressure NMR.

We have investigated the pressure- and temperature-induced conformational changes associated with the low complexity domain of hnRNP A1, an RNA-binding protein able to phase separate in response to cellular stress. Solution NMR spectra of the hnRNP A1 low-complexity domain fused with protein-G B1 domain were collected from 1 to 2500 bar and from 268 to 290 K. While the GB1 domain shows the typical pressure-induced and cold temperature-induced unfolding expected for small globular domains, the low-complexity domain of hnRNP A1 exhibits unusual pressure and temperature dependences. We observed that the low-complexity domain is pressure sensitive, undergoing a major conformational transition within the prescribed pressure range. Remarkably, this transition has the inverse temperature dependence of a typical folding-unfolding transition. Our results suggest the presence of a low-lying extended and fully solvated state(s) of the low-complexity domain that may play a role in phase separation. This study highlights the exquisite sensitivity of solution NMR spectroscopy to observe subtle conformational changes and illustrates how pressure perturbation can be used to determine the properties of metastable conformational ensembles.

Dissociation of the trimeric gp41 ectodomain at the lipid-water interface suggests an active role in HIV-1 Env-mediated membrane fusion.

The envelope glycoprotein gp41 mediates the process of membrane fusion that enables entry of the HIV-1 virus into the host cell. The actual fusion process involves a switch from a homotrimeric prehairpin intermediate conformation, consisting of parallel coiled-coil helices, to a postfusion state where the ectodomains are arranged as a trimer of helical hairpins, adopting a six-helix bundle (6HB) state. Here, we show by solution NMR spectroscopy that a water-soluble 6HB gp41 ectodomain binds to zwitterionic detergents that contain phosphocholine or phosphatidylcholine head groups and phospholipid vesicles that mimic T-cell membrane composition. Binding results in the dissociation of the 6HB and the formation of a monomeric state, where its two α-helices, N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR), become embedded in the lipid-water interface of the virus and host cell. The atomic structure of the gp41 ectodomain monomer, based on NOE distance restraints and residual dipolar couplings, shows that the NHR and CHR helices remain mostly intact, but they completely lose interhelical contacts. The high affinity of the ectodomain helices for phospholipid surfaces suggests that unzippering of the prehairpin intermediate leads to a state where the NHR and CHR helices become embedded in the host cell and viral membranes, respectively, thereby providing a physical force for bringing these membranes into close juxtaposition before actual fusion.

Monomeric Aß(1-40) and Aß(1-42) Peptides in Solution Adopt Very Similar Ramachandran Map Distributions That Closely Resemble Random Coil.

The pathogenesis of Alzheimer's disease is characterized by the aggregation and fibrillation of amyloid peptides Aß(1-40) and Aß(1-42) into amyloid plaques. Despite strong potential therapeutic interest, the structural pathways associated with the conversion of monomeric Aß peptides into oligomeric species remain largely unknown. In particular, the higher aggregation propensity and associated toxicity of Aß(1-42) compared to that of Aß(1-40) are poorly understood. To explore in detail the structural propensity of the monomeric Aß(1-40) and Aß(1-42) peptides in solution, we recorded a large set of nuclear magnetic resonance (NMR) parameters, including chemical shifts, nuclear Overhauser effects (NOEs), and J couplings. Systematic comparisons show that at neutral pH the Aß(1-40) and Aß(1-42) peptides populate almost indistinguishable coil-like conformations. Nuclear Overhauser effect spectra collected at very high resolution remove assignment ambiguities and show no long-range NOE contacts. Six sets of backbone J couplings ((3)JHNHα, (3)JC'C', (3)JC'Hα, (1)JHαCα, (2)JNCα, and (1)JNCα) recorded for Aß(1-40) were used as input for the recently developed MERA Ramachandran map analysis, yielding residue-specific backbone ϕ/ψ torsion angle distributions that closely resemble random coil distributions, the absence of a significantly elevated propensity for ß-conformations in the C-terminal region of the peptide, and a small but distinct propensity for αL at K28. Our results suggest that the self-association of Aß peptides into toxic oligomers is not driven by elevated propensities of the monomeric species to adopt ß-strand-like conformations. Instead, the accelerated disappearance of Aß NMR signals in D2O over H2O, particularly pronounced for Aß(1-42), suggests that intermolecular interactions between the hydrophobic regions of the peptide dominate the aggregation process.

Accurate measurement of (3)J(HNHα) couplings in small or disordered proteins from WATERGATE-optimized TROSY spectra.

Provided that care is taken in adjusting the WATERGATE element of a (1)H-(15)N TROSY-HSQC experiment, such that neither the water magnetization nor the (1)H(α) protons are inverted by its final 180° pulse, (3)JHNHα couplings can be measured directly from splittings in the (1)H dimension of the spectrum. With band-selective (1)H decoupling, very high (15)N resolution can be achieved. A complete set of (3)JHNHα values, ranging from 3.4 to 10.1 Hz was measured for the 56-residue third domain of IgG-binding protein G (GB3). Using the H-N-C(α)-H(α) dihedral angles extracted from a RDC-refined structure of GB3, (3)JHNHα values predicted by a previously parameterized Karplus equation agree to within a root-mean-square deviation (rmsd) of 0.37 Hz with the experimental data. Values measured for the Alzheimer's implicated Aß(1-40) peptide fit to within an rmsd of 0.45 Hz to random coil (3)JHNHα values.

Exploring the Protein Folding Pathway with High-Pressure NMR: Steady-State and Kinetics Studies.

Defining the physical-chemical determinants of protein folding and stability, under normal and pathological conditions has constituted a major subfield in biophysical chemistry for over 50 years. Although a great deal of progress has been made in recent years towards this goal, a number of important questions remain. These include characterizing the structural, thermodynamic and dynamic properties of the barriers between conformational states on the protein energy landscape, understanding the sequence dependence of folding cooperativity, defining more clearly the role of solvation in controlling protein stability and dynamics and probing the high energy thermodynamic states in the native state basin and their role in misfolding and aggregation. Fundamental to the elucidation of these questions is a complete thermodynamic parameterization of protein folding determinants. In this chapter, we describe the use of high-pressure coupled to Nuclear Magnetic Resonance (NMR) spectroscopy to reveal unprecedented details on the folding energy landscape of proteins.

Pressure-induced structural transition of mature HIV-1 protease from a combined NMR/MD simulation approach.

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.

Evolutionarily Conserved Pattern of Interactions in a Protein Revealed by Local Thermal Expansion Properties.

The way in which the network of intramolecular interactions determines the cooperative folding and conformational dynamics of a protein remains poorly understood. High-pressure NMR spectroscopy is uniquely suited to examine this problem because it combines the site-specific resolution of the NMR experiments with the local character of pressure perturbations. Here we report on the temperature dependence of the site-specific volumetric properties of various forms of staphylococcal nuclease (SNase), including three variants with engineered internal cavities, as measured with high-pressure NMR spectroscopy. The strong temperature dependence of pressure-induced unfolding arises from poorly understood differences in thermal expansion between the folded and unfolded states. A significant inverse correlation was observed between the global thermal expansion of the folded proteins and the number of strong intramolecular hydrogen bonds, as determined by the temperature coefficient of the backbone amide chemical shifts. Comparison of the identity of these strong H-bonds with the co-evolution of pairs of residues in the SNase protein family suggests that the architecture of the interactions detected in the NMR experiments could be linked to a functional aspect of the protein. Moreover, the temperature dependence of the residue-specific volume changes of unfolding yielded residue-specific differences in expansivity and revealed how mutations impact intramolecular interaction patterns. These results show that intramolecular interactions in the folded states of proteins impose constraints against thermal expansion and that, hence, knowledge of site-specific thermal expansivity offers insight into the patterns of strong intramolecular interactions and other local determinants of protein stability, cooperativity, and potentially also of function.

Complete dissociation of the HIV-1 gp41 ectodomain and membrane proximal regions upon phospholipid binding.

The envelope glycoprotein gp41 mediates the process of membrane fusion that enables entry of the HIV-1 virus into the host cell. Strong lipid affinity of the ectodomain suggests that its heptad repeat regions play an active role in destabilizing membranes by directly binding to the lipid bilayers and thereby lowering the free-energy barrier for membrane fusion. In such a model, immediately following the shedding of gp120, the N-heptad and C-heptad helices dissociate and melt into the host cell and viral membranes, respectively, pulling the destabilized membranes into juxtaposition, ready for fusion. Post-fusion, reaching the final 6-helix bundle (6 HB) conformation then involves competition between intermolecular interactions needed for formation of the symmetric 6 HB trimer and the membrane affinity of gp41's ectodomain, including its membrane-proximal regions. Our solution NMR study of the structural and dynamic properties of three constructs containing the ectodomain of gp41 with and without its membrane-proximal regions suggests that these segments do not form inter-helical interactions until the very late steps of the fusion process. Interactions between the polar termini of the heptad regions, which are not associating with the lipid surface, therefore may constitute the main driving force initiating formation of the final post-fusion states. The absence of significant intermolecular ectodomain interactions in the presence of dodecyl phosphocholine highlights the importance of trimerization of gp41's transmembrane helix to prevent complete dissociation of the trimer during the course of fusion.

Conformation of inhibitor-free HIV-1 protease derived from NMR spectroscopy in a weakly oriented solution.

Flexibility of the glycine-rich flaps is known to be essential for catalytic activity of the HIV-1 protease, but their exact conformations at the different stages of the enzymatic pathway remain subject to much debate. Although hundreds of crystal structures of protease-inhibitor complexes have been solved, only about a dozen inhibitor-free protease structures have been reported. These latter structures reveal a large diversity of flap conformations, ranging from closed to semi-open to wide open. To evaluate the average structure in solution, we measured residual dipolar couplings (RDCs) and compared these to values calculated for crystal structures representative of the closed, semi-open, and wide-open states. The RDC data clearly indicate that the inhibitor-free protease, on average, adopts a closed conformation in solution that is very similar to the inhibitor-bound state. By contrast, a highly drug-resistant protease mutant, PR20, adopts the wide-open flap conformation.

Cavities determine the pressure unfolding of proteins.

It has been known for nearly 100 years that pressure unfolds proteins, yet the physical basis of this effect is not understood. Unfolding by pressure implies that the molar volume of the unfolded state of a protein is smaller than that of the folded state. This decrease in volume has been proposed to arise from differences between the density of bulk water and water associated with the protein, from pressure-dependent changes in the structure of bulk water, from the loss of internal cavities in the folded states of proteins, or from some combination of these three factors. Here, using 10 cavity-containing variants of staphylococcal nuclease, we demonstrate that pressure unfolds proteins primarily as a result of cavities that are present in the folded state and absent in the unfolded one. High-pressure NMR spectroscopy and simulations constrained by the NMR data were used to describe structural and energetic details of the folding landscape of staphylococcal nuclease that are usually inaccessible with existing experimental approaches using harsher denaturants. Besides solving a 100-year-old conundrum concerning the detailed structural origins of pressure unfolding of proteins, these studies illustrate the promise of pressure perturbation as a unique tool for examining the roles of packing, conformational fluctuations, and water penetration as determinants of solution properties of proteins, and for detecting folding intermediates and other structural details of protein-folding landscapes that are invisible to standard experimental approaches.

The energetics of a three-state protein folding system probed by high-pressure relaxation dispersion NMR spectroscopy.

The energetic and volumetric properties of a three-state protein folding system, comprising a metastable triple mutant of the Fyn SH3 domain, have been investigated using pressure-dependent (15) N-relaxation dispersion NMR from 1 to 2500 bar. Changes in partial molar volumes (ΔV) and isothermal compressibilities (ΔκT ) between all the states along the folding pathway have been determined to reasonable accuracy. The partial volume and isothermal compressibility of the folded state are 100 mL mol(-1) and 40 µL mol(-1) bar(-1) , respectively, higher than those of the unfolded ensemble. Of particular interest are the findings related to the energetic and volumetric properties of the on-pathway folding intermediate. While the latter is energetically close to the unfolded state, its volumetric properties are similar to those of the folded protein. The compressibility of the intermediate is larger than that of the folded state reflecting the less rigid nature of the former relative to the latter.

Improved cross validation of a static ubiquitin structure derived from high precision residual dipolar couplings measured in a drug-based liquid crystalline phase.

The antibiotic squalamine forms a lyotropic liquid crystal at very low concentrations in water (0.3-3.5% w/v), which remains stable over a wide range of temperature (1-40 °C) and pH (4-8). Squalamine is positively charged, and comparison of the alignment of ubiquitin relative to 36 previously reported alignment conditions shows that it differs substantially from most of these, but is closest to liquid crystalline cetyl pyridinium bromide. High precision residual dipolar couplings (RDCs) measured for the backbone (1)H-(15)N, (15)N-(13)C', (1)H(α)-(13)C(α), and (13)C'-(13)C(α) one-bond interactions in the squalamine medium fit well to the static structural model previously derived from NMR data. Inclusion into the structure refinement procedure of these RDCs, together with (1)H-(15)N and (1)H(α)-(13)C(α) RDCs newly measured in Pf1, results in improved agreement between alignment-induced changes in (13)C' chemical shift, (3)JHNHα values, and (13)C(α)-(13)C(ß) RDCs and corresponding values predicted by the structure, thereby validating the high quality of the single-conformer structural model. This result indicates that fitting of a single model to experimental data provides a better description of the average conformation than does averaging over previously reported NMR-derived ensemble representations. The latter can capture dynamic aspects of a protein, thus making the two representations valuable complements to one another.

Elucidation of the Mechanisms of Inter-domain Coupling in the Monomeric State of Enzyme I by High-pressure NMR.

The catalytic cycle of Enzyme I (EI), a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate, is characterized by a series of local and global conformational rearrangements. This multistep process includes a monomer-to-dimer transition, followed by an open-to-closed rearrangement of the dimeric complex upon PEP binding. In the present study, we investigate the thermodynamics of EI dimerization using a range of high-pressure solution NMR techniques complemented by SAXS experiments. 1H-15N TROSY and 1H-13C methyl TROSY NMR spectra combined with 15N relaxation measurements revealed that a native-like engineered variant of full-length EI fully dissociates into stable monomeric state above 1.5 kbar. Conformational ensembles of EI monomeric state were generated via a recently developed protocol combining coarse-grained molecular simulations with experimental backbone residual dipolar coupling measurements. Analysis of the structural ensembles provided detailed insights into the molecular mechanisms driving formation of the catalytically competent dimeric state, and reveals that each step of EI catalytical cycle is associated with a significant reduction in either inter- or intra-domain conformational entropy. Altogether, this study completes a large body work conducted by our group on EI and establishes a comprehensive structural and dynamical description of the catalytic cycle of this prototypical multidomain, oligomeric enzyme.

APIS: an updated parentage assignment software managing triploids induced from diploid parents.

In aquaculture, sterile triploids are commonly used for production as sterility gives them potential gains in growth, yields, and quality. However, they cannot be reproduced, and DNA parentage assignment to their diploid or tetraploid parents is required to estimate breeding values for triploid phenotypes. No publicly available software has the ability to assign triploids to their parents. Here, we updated the R package APIS to support triploids induced from diploid parents. First, we created new exclusion and likelihood tables that account for the double allelic contribution of the dam and the recombination that can occur during female meiosis. As the effective recombination rate of each marker with the centromere is usually unknown, we set it at 0.5 and found that this value maximizes the assignment rate even for markers with high or low recombination rates. The number of markers needed for a high true assignment rate did not strongly depend on the proportion of missing parental genotypes. The assignment power was however affected by the quality of the markers (minor allele frequency, call rate). Altogether, 96-192 SNPs were required to have a high parentage assignment rate in a real rainbow trout dataset of 1,232 triploid progenies from 288 parents. The likelihood approach was more efficient than exclusion when the power of the marker set was limiting. When more markers were used, exclusion was more advantageous, with sensitivity reaching unity, very low false discovery rate (<0.01), and excellent specificity (0.96-0.99). Thus, APIS provides an efficient solution to assign triploids to their diploid parents.

Thermodynamic coupling of the tandem RRM domains of hnRNP A1 underlie its pleiotropic RNA binding functions.

The functional properties of RNA binding proteins (RBPs) require allosteric regulation through interdomain communication. Despite the importance of allostery to biological regulation, only a few studies have been conducted to describe the biophysical nature by which interdomain communication manifests in RBPs. Here, we show for hnRNP A1 that interdomain communication is vital for the unique stability of its amino-terminal domain, which consists of two RNA recognition motifs (RRMs). These RRMs exhibit drastically different stability under pressure. RRM2 unfolds as an individual domain but remains stable when appended to RRM1. Variants that disrupt interdomain communication between the tandem RRMs show a significant decrease in stability. Carrying these mutations over to the full-length protein for in vivo experiments revealed that the mutations affected the ability of the disordered carboxyl-terminal domain to engage in protein-protein interactions and influenced the protein's RNA binding capacity. Collectively, this work reveals that thermodynamic coupling between the tandem RRMs of hnRNP A1 accounts for its allosteric regulatory functions.

Enhanced stability of monomer fold correlates with extreme drug resistance of HIV-1 protease.

During treatment, mutations in HIV-1 protease (PR) are selected rapidly that confer resistance by decreasing affinity to clinical protease inhibitors (PIs). As these unique drug resistance mutations can compromise the fitness of the virus to replicate, mutations that restore conformational stability and activity while retaining drug resistance are selected on further evolution. Here we identify several compensating mechanisms by which an extreme drug-resistant mutant bearing 20 mutations (PR20) with >5-fold increased Kd and >4000-fold decreased affinity to the PI darunavir functions. (1) PR20 cleaves, albeit poorly, Gag polyprotein substrates essential for viral maturation. (2) PR20 dimer, which exhibits distinctly enhanced thermal stability, has highly attenuated autoproteolysis, thus likely prolonging its lifetime in vivo. (3) The enhanced stability of PR20 results from stabilization of the monomer fold. Both monomeric PR20(T26A) and dimeric PR20 exhibit Tm values 6-7.5 °C higher than those for their PR counterparts. Two specific mutations in PR20, L33F and L63P at sites of autoproteolysis, increase the Tm of monomeric PR(T26A) by ~8 °C, similar to PR20(T26A). However, without other compensatory mutations as seen in PR20, L33F and L63P substitutions, together, neither restrict autoproteolysis nor significantly reduce binding affinity to darunavir. To determine whether dimer stability contributes to binding affinity for inhibitors, we examined single-chain dimers of PR and PR(D25N) in which the corresponding identical monomer units were covalently linked by GGSSG sequence. Linking of the subunits did not appreciably change the ΔTm on inhibitor binding; thus stabilization by tethering appears to have little direct effect on enhancing inhibitor affinity.