Visualizing the active site oxyanion loop transition upon ensitrelvir binding and transient dimerization of SARS-CoV-2 main protease

N-terminal autoprocessing from its polyprotein precursor enables creating the mature-like stable dimer interface of SARS-CoV-2 main protease (MPro), concomitant with the active site oxyanion loop equilibrium transitioning to the active conformation (E*) and onset of catalytic activity. Through mutagenesis of critical interface residues and evaluating noncovalent inhibitor (ensitrelvir, ESV) facilitated dimerization through its weak binding to MPro, we demonstrate that residues extending from Ser1 through Glu14 are critical for dimerization. Combined mutations G11A, E290A and R298A (MProTM) restrict dimerization even upon binding of ESV to monomeric MProTM with an inhibitor dissociation constant of 7.4 ± 1.6 µM. Contrasting the covalent inhibitor NMV or GC373 binding to monomeric MPro, ESV binding enabled capturing the transition of the oxyanion loop conformations in the absence of a reactive warhead and independent of dimerization. Characterization of complexes by room-temperature X-ray crystallography reveals ESV bound to the E* state of monomeric MPro as well as an intermediate approaching the inactive state (E). It appears that the E* to E equilibrium shift occurs initially from G138-F140 residues, leading to the unwinding of the loop and formation of the 310-helix. Finally, we describe a transient dimer structure of the MPro precursor held together through interactions of residues A5-G11 with distinct states of the active sites, E and E*, likely representing an intermediate in the autoprocessing pathway.

Insights into the mechanism of SARS-CoV-2 main protease autocatalytic maturation from model precursors

A critical step for SARS-CoV-2 assembly and maturation involves the autoactivation of the main protease (MProWT) from precursor polyproteins. Upon expression, a model precursor of MProWT mediates its own release at its termini rapidly to yield a mature dimer. A construct with an E290A mutation within MPro exhibits time dependent autoprocessing of the accumulated precursor at the N-terminal nsp4/nsp5 site followed by the C-terminal nsp5/nsp6 cleavage. In contrast, a precursor containing E290A and R298A mutations (MProM) displays cleavage only at the nsp4/nsp5 site to yield an intermediate monomeric product, which is cleaved at the nsp5/nsp6 site only by MProWT. MProM and the catalytic domain (MPro1-199) fused to the truncated nsp4 region also show time-dependent conversion in vitro to produce MProM and MPro1-199, respectively. The reactions follow first-order kinetics indicating that the nsp4/nsp5 cleavage occurs via an intramolecular mechanism. These results support a mechanism involving an N-terminal intramolecular cleavage leading to an increase in the dimer population and followed by an intermolecular cleavage at the C-terminus. Thus, targeting the predominantly monomeric MPro precursor for inhibition may lead to the identification of potent drugs for treatment.

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MPro^WT) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373) and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ∼400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decreases in free energy (ΔG) and thermal stability (T_m), respectively, relative to MPro^WT. The H41A mutation results in a 20-fold increase in the dissociation constant (K_d), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and T_m, respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MPro^H41A, whereas no significant change is observed in binding to MPro^WT. Structures of the four inhibitor complexes with MPro^1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MPro^WT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MPro^WT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MPro^C145A also results in a significant, but smaller, increase in K_d and decrease in ΔG and T_m, relative to NMV.

Structures of brain-derived 42-residue amyloid-β fibril polymorphs with unusual molecular conformations and intermolecular interactions

Fibrils formed by the 42-residue amyloid-β peptide (Aβ42), a main component of amyloid deposits in Alzheimer's disease (AD), are known to be polymorphic, i.e., to contain multiple possible molecular structures. Previous studies of Aβ42 fibrils, including fibrils prepared entirely in vitro or extracted from brain tissue and using solid-state NMR (ssNMR) or cryogenic electron microscopy (cryo-EM) methods, have found polymorphs with differences in amino acid sidechain orientations, lengths of structurally ordered segments, and contacts between cross-β subunit pairs within a single filament. Despite these differences, Aβ42 molecules adopt a common S-shaped conformation in all previously described high-resolution Aβ42 fibril structures. Here we report two cryo-EM-based structures of Aβ42 fibrils that are qualitatively different, in samples derived from AD brain tissue by seeded growth. In type A fibrils, residues 12 to 42 adopt a ν-shaped conformation, with both intra-subunit and intersubunit hydrophobic contacts to form a compact core. In type B fibrils, residues 2 to 42 adopt an υ-shaped conformation, with only intersubunit contacts and internal pores. Type A and type B fibrils have opposite helical handedness. Cryo-EM density maps and molecular dynamics simulations indicate intersubunit K16-A42 salt bridges in type B fibrils and partially occupied K28-A42 salt bridges in type A fibrils. The coexistence of two predominant polymorphs, with differences in N-terminal dynamics, is supported by ssNMR data, as is faithful propagation of structures from first-generation to second-generation brain-seeded Aβ42 fibril samples. These results demonstrate that Aβ42 fibrils can exhibit a greater range of structural variations than seen in previous studies.

Autoprocessing and oxyanion loop reorganization upon GC373 and nirmatrelvir binding of monomeric SARS-CoV-2 main protease catalytic domain

The monomeric catalytic domain (residues 1–199) of SARS-CoV-2 main protease (MPro^1-199) fused to 25 amino acids of its flanking nsp4 region mediates its autoprocessing at the nsp4-MPro^1-199 junction. We report the catalytic activity and the dissociation constants of MPro^1-199 and its analogs with the covalent inhibitors GC373 and nirmatrelvir (NMV), and the estimated monomer-dimer equilibrium constants of these complexes. Mass spectrometry indicates the presence of the accumulated adduct of NMV bound to MPro^WT and MPro^1-199 and not of GC373. A room temperature crystal structure reveals a native-like fold of the catalytic domain with an unwound oxyanion loop (E state). In contrast, the structure of a covalent complex of the catalytic domain-GC373 or NMV shows an oxyanion loop conformation (E* state) resembling the full-length mature dimer. These results suggest that the E-E* equilibrium modulates autoprocessing of the main protease when converting from a monomeric polyprotein precursor to the mature dimer.

Structural characterization and catalytic activity of SARS-CoV-2 main protease reveal minimal interface regions enabling dimer formation driven by inhibitor-induced conformational changes of the oxyanion loop.

Real-time Exchange of the Lipid-bound Intermediate and Post-fusion States of the HIV-1 gp41 Ectodomain

The envelope glycoprotein gp41 of the HIV-1 virus mediates its entry into the host cell. During this process, gp41 undergoes large conformational changes and the energy released in the remodeling events is utilized to overcome the barrier associated with fusing the viral and host membranes. Although the structural intermediates of this fusion process are attractive targets for drug development, no detailed high-resolution structural information or quantitative thermodynamic characterization are available. By measuring the dynamic equilibrium between the lipid-bound intermediate and the post-fusion six-helical bundle (6HB) states of the gp41 ectodomain in the presence of bilayer membrane mimetics, we derived both the reaction kinetics and energies associated with these two states by solution NMR spectroscopy. At equilibrium, an exchange time constant of about 12 seconds at 38 °C is observed, and the post-fusion conformation is energetically more stable than the lipid-bound state by 3.4 kcal mol-1. The temperature dependence of the kinetics indicates that the folding occurs through a high-energy transition state which may resemble a 5HB structure. The energetics and kinetics of gp41 folding in the context of membrane bilayers provide a molecular basis for an improved understanding of viral membrane fusion.

Global Dynamics of a Protein on the Surface of Anisotropic Lipid Nanoparticles Derived from Relaxation-Based NMR Spectroscopy

The global motions of ubiquitin, a model protein, on the surface of anisotropically tumbling 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (POPG):1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC) bicelles are described. The shapes of POPG:DHPC bicelles prepared with high molar ratios q of POPG to DHPC can be approximated by prolate ellipsoids, with the ratio of ellipsoid dimensions and dimensions themselves increasing with higher values of q. Adaptation of the nuclear magnetic resonance (NMR) relaxation-based approach that we previously developed for interactions of ubiquitin with spherical POPG liposomes (Ceccon, A. J. Am. Chem. Soc. 2016, 138, 5789-5792) allowed us to quantitatively analyze the variation in lifetime line broadening of NMR signals (ΔR₂) measured for ubiquitin in the presence of q = 2 POPG:DHPC bicelles and the associated transverse spin relaxation rates (R_2,B) of bicelle-bound ubiquitin. Ubiquitin, transiently bound to POPG:DHPC bicelles, undergoes internal rotation about an axis orthogonal to the surface of the bicelle and perpendicular to the principal axis of its rotational diffusion tensor on the low microsecond time scale (∼3 μs), while the rotation axis itself wobbles in a cone on a submicrosecond time scale (≤ 500 ns).

Modulation of the monomer-dimer equilibrium and catalytic activity of SARS-CoV-2 main protease by a transition-state analog inhibitor

The role of dimer formation for the onset of catalytic activity of SARS-CoV-2 main protease (MPro^WT) was assessed using a predominantly monomeric mutant (MPro^M). Rates of MPro^WT and MPro^M catalyzed hydrolyses display substrate saturation kinetics and second-order dependency on the protein concentration. The addition of the prodrug GC376, an inhibitor of MPro^WT, to MPro^M leads to an increase in the dimer population and catalytic activity with increasing inhibitor concentration. The activity reaches a maximum corresponding to a dimer population in which one active site is occupied by the inhibitor and the other is available for catalytic activity. This phase is followed by a decrease in catalytic activity due to the inhibitor competing with the substrate. Detailed kinetics and equilibrium analyses are presented and a modified Michaelis-Menten equation accounts for the results. These observations provide conclusive evidence that dimer formation is coupled to catalytic activity represented by two equivalent active sites.

Structural, Electronic, and Electrostatic Determinants for Inhibitor Binding to Subsites S1 and S2 in SARS-CoV-2 Main Protease

Creating small-molecule antivirals specific for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins is crucial to battle coronavirus disease 2019 (COVID-19). SARS-CoV-2 main protease (Mpro) is an established drug target for the design of protease inhibitors. We performed a structure-activity relationship (SAR) study of noncovalent compounds that bind in the enzyme's substrate-binding subsites S1 and S2, revealing structural, electronic, and electrostatic determinants of these sites. The study was guided by the X-ray/neutron structure of Mpro complexed with Mcule-5948770040 (compound 1), in which protonation states were directly visualized. Virtual reality-assisted structure analysis and small-molecule building were employed to generate analogues of 1. In vitro enzyme inhibition assays and room-temperature X-ray structures demonstrated the effect of chemical modifications on Mpro inhibition, showing that (1) maintaining correct geometry of an inhibitor's P1 group is essential to preserve the hydrogen bond with the protonated His163; (2) a positively charged linker is preferred; and (3) subsite S2 prefers nonbulky modestly electronegative groups.

Michaelis-like complex of SARS-CoV-2 main protease visualized by room-temperature X-ray crystallography

SARS-CoV-2 emerged at the end of 2019 to cause an unprecedented pandemic of the deadly respiratory disease COVID-19 that continues to date. The viral main protease (Mpro) is essential for SARS-CoV-2 replication and is therefore an important drug target. Understanding the catalytic mechanism of Mpro, a cysteine protease with a catalytic site comprising the noncanonical Cys145-His41 dyad, can help in guiding drug design. Here, a 2.0 Å resolution room-temperature X-ray crystal structure is reported of a Michaelis-like complex of Mpro harboring a single inactivating mutation C145A bound to the octapeptide Ac-SAVLQSGF-CONH2 corresponding to the nsp4/nsp5 autocleavage site. The peptide substrate is unambiguously defined in subsites S5 to S3' by strong electron density. Superposition of the Michaelis-like complex with the neutron structure of substrate-free Mpro demonstrates that the catalytic site is inherently pre-organized for catalysis prior to substrate binding. Induced fit to the substrate is driven by P1 Gln binding in the predetermined subsite S1 and rearrangement of subsite S2 to accommodate P2 Leu. The Michaelis-like complex structure is ideal for in silico modeling of the SARS-CoV-2 Mpro catalytic mechanism.

Transient lipid-bound states of spike protein heptad repeats provide insights into SARS-CoV-2 membrane fusion

Entry of SARS-CoV-2 into a host cell is mediated by spike, a class I viral fusion protein responsible for merging the viral and host cell membranes. Recent studies have revealed atomic-resolution models for both the postfusion 6-helix bundle (6HB) and the prefusion state of spike. However, a mechanistic understanding of the molecular basis for the intervening structural transition, important for the design of fusion inhibitors, has remained elusive. Using nuclear magnetic resonance spectroscopy and other biophysical methods, we demonstrate the presence of α-helical, membrane-bound, intermediate states of spike’s heptad repeat (HR1 and HR2) domains that are embedded at the lipid-water interface while in a slow dynamic equilibrium with the postfusion 6HB state. These results support a model where the HR domains lower the large energy barrier associated with membrane fusion by destabilizing the host and viral membranes, while 6HB formation actively drives their fusion by forcing physical proximity.

Constraints on the Structure of Fibrils Formed by a Racemic Mixture of Amyloid-β Peptides from Solid-State NMR, Electron Microscopy, and Theory

Previous studies have shown that racemic mixtures of 40- and 42-residue amyloid-β peptides (d,l-Aβ40 and d,l-Aβ42) form amyloid fibrils with accelerated kinetics and enhanced stability relative to their homochiral counterparts (l-Aβ40 and l-Aβ42), suggesting a "chiral inactivation" approach to abrogating the neurotoxicity of Aβ oligomers (Aβ-CI). Here we report a structural study of d,l-Aβ40 fibrils, using electron microscopy, solid-state nuclear magnetic resonance (NMR), and density functional theory (DFT) calculations. Two- and three-dimensional solid-state NMR spectra indicate molecular conformations in d,l-Aβ40 fibrils that resemble those in known l-Aβ40 fibril structures. However, quantitative measurements of 13C-13C and 15N-13C distances in selectively labeled d,l-Aβ40 fibril samples indicate a qualitatively different supramolecular structure. While cross-β structures in mature l-Aβ40 fibrils are comprised of in-register, parallel β-sheets, our data indicate antiparallel β-sheets in d,l-Aβ40 fibrils, with alternation of d and l molecules along the fibril growth direction, i.e., antiparallel "rippled sheet" structures. The solid-state NMR data suggest the coexistence of d,l-Aβ40 fibril polymorphs with three different registries of intermolecular hydrogen bonds within the antiparallel rippled sheets. DFT calculations support an energetic preference for antiparallel alignments of the β-strand segments identified by solid-state NMR. These results provide insight into the structural basis for Aβ-CI and establish the importance of rippled sheets in self-assembly of full-length, naturally occurring amyloidogenic peptides.

MWC allosteric model explains unusual hemoglobin-oxygen binding curves from sickle cell drug binding

An oxygen-affinity-modifying drug, voxelotor, has very recently been approved by the FDA for treatment of sickle cell disease. The proposed mechanism of action is by preferential binding of the drug to the R quaternary conformation, which cannot copolymerize with the T conformation to form sickle fibers. Here, we report widely different oxygen dissociation and oxygen association curves for normal blood in the presence of voxelotor and interpret the results in terms of the allosteric model of Monod, Wyman, and Changeux with the addition of drug binding. The model does remarkably well in quantitatively explaining a complex data set with just the addition of drug binding and dissociation rates for the R and T conformations. Whereas slow dissociation of the drug from R results in time-independent dissociation curves, the changing association curves result from slow dissociation of the drug from T, as well as extremely slow binding of the drug to T. By calculating true equilibrium curves from the model parameters, we show that there would be a smaller decrease in oxygen delivery from the left shift in the dissociation curve caused by drug binding if drug binding and dissociation for both R and T were rapid. Our application of the Monod, Wyman, and Changeux model demonstrates once more its enormous power in explaining many different kinds of experimental results for hemoglobin. It should also be helpful in analyzing oxygen binding and in vivo delivery in future investigations of oxygen-affinity-modifying drugs for sickle cell disease.

A weakened interface in the P182L variant of HSP27 associated with severe Charcot-Marie-Tooth neuropathy causes aberrant binding to interacting proteins

HSP27 is a human molecular chaperone that forms large, dynamic oligomers and functions in many aspects of cellular homeostasis. Mutations in HSP27 cause Charcot‐Marie‐Tooth (CMT) disease, the most common inherited disorder of the peripheral nervous system. A particularly severe form of CMT disease is triggered by the P182L mutation in the highly conserved IxI/V motif of the disordered C‐terminal region, which interacts weakly with the structured core domain of HSP27. Here, we observed that the P182L mutation disrupts the chaperone activity and significantly increases the size of HSP27 oligomers formed in vivo, including in motor neurons differentiated from CMT patient‐derived stem cells. Using NMR spectroscopy, we determined that the P182L mutation decreases the affinity of the HSP27 IxI/V motif for its own core domain, leaving this binding site more accessible for other IxI/V‐containing proteins. We identified multiple IxI/V‐bearing proteins that bind with higher affinity to the P182L variant due to the increased availability of the IxI/V‐binding site. Our results provide a mechanistic basis for the impact of the P182L mutation on HSP27 and suggest that the IxI/V motif plays an important, regulatory role in modulating protein–protein interactions.

Probing the Interaction between HIV-1 Protease and the Homodimeric p66/p66' Reverse Transcriptase Precursor by Double Electron-Electron Resonance EPR Spectroscopy

Following excision from the Gag-Pol polyprotein, HIV-1 reverse transcriptase is released as an asymmetric homodimer comprising two p66 subunits that are structurally dissimilar but identical in amino acid sequence. Subsequent cleavage of the RNase H domain from only one of the subunits, denoted p66', results in the formation of the mature p66/p51 enzyme in which catalytic activity resides in the p66 subunit, and the p51 subunit (derived from p66') provides a supporting structural scaffold. Here, we probe the interaction of the p66/p66' asymmetric reverse transcriptase precursor with HIV-1 protease by pulsed Q-band double electron-electron resonance EPR spectroscopy to measure distances between nitroxide labels introduced at surface-engineered cysteine residues. The data suggest that the flexible, exposed linker between the RNaseH and connection domains in the open state of the p66' subunit binds to the active site of protease in a configuration that is similar to that of extended peptide substrates.

Concentration-Dependent Structural Transition of the HIV-1 gp41 MPER Peptide into α-Helical Trimers

The membrane proximal external region (MPER) of HIV-1 gp41 contains epitopes for at least four broadly neutralizing antibodies. Depending on solution conditions and construct design, different structures have been reported for this segment. We show that in aqueous solution the MPER fragment (gp160^660-674 ) exists in a monomer-trimer equilibrium with an association constant in the micromolar range. Thermodynamic analysis reveals that the association is exothermic, more favorable in D₂ O than H₂ O, and increases with ionic strength, indicating hydrophobically driven intermolecular interactions. Circular dichroism, ¹³ C^α chemical shifts, NOE, and hydrogen exchange rates reveal that MPER undergoes a structural transition from predominately unfolded monomer at low concentrations to an α-helical trimer at high concentrations. This result has implications for antibody recognition of MPER prior to and during the process where gp41 switches from a pre-hairpin intermediate to its post-fusion 6-helical bundle state.

Visualizing Tetrahedral Oxyanion Bound in HIV-1 Protease Using Neutrons: Implications for the Catalytic Mechanism and Drug Design

HIV-1 protease is indispensable for virus propagation and an important therapeutic target for antiviral inhibitors to treat AIDS. As such inhibitors are transition-state mimics, a detailed understanding of the enzyme mechanism is crucial for the development of better anti-HIV drugs. Here, we used room-temperature joint X-ray/neutron crystallography to directly visualize hydrogen atoms and map hydrogen bonding interactions in a protease complex with peptidomimetic inhibitor KVS-1 containing a reactive nonhydrolyzable ketomethylene isostere, which, upon reacting with the catalytic water molecule, is converted into a tetrahedral intermediate state, KVS-1_TI. We unambiguously determined that the resulting tetrahedral intermediate is an oxyanion, rather than the gem-diol, and both catalytic aspartic acid residues are protonated. The oxyanion tetrahedral intermediate appears to be unstable, even though the negative charge on the oxyanion is delocalized through a strong n → π* hyperconjugative interaction into the nearby peptidic carbonyl group of the inhibitor. To better understand the influence of the ketomethylene isostere as a protease inhibitor, we have also examined the protease structure and binding affinity with keto-darunavir (keto-DRV), which similar to KVS-1 includes the ketomethylene isostere. We show that keto-DRV is a significantly less potent protease inhibitor than DRV. These findings shed light on the reaction mechanism of peptide hydrolysis catalyzed by HIV-1 protease and provide valuable insights into further improvements in the design of protease inhibitors.

Proton transfer and drug binding details revealed in neutron diffraction studies of wild-type and drug resistant HIV-1 protease

HIV-1 protease is an essential therapeutic target for the design and development of antiviral inhibitors to treat AIDS. We used room temperature neutron crystallography to accurately determine hydrogen atom positions in several protease complexes with clinical drugs, amprenavir and darunavir. Hydrogen bonding interactions were carefully mapped to provide an unprecedented picture of drug binding to the protease target. We demonstrate that hydrogen atom positions within the enzyme catalytic site can be altered by introducing drug resistant mutations and by protonating surface residues that trigger proton transfer reactions between the catalytic Asp residues and the hydroxyl group of darunavir. When protein perdeuteration is not feasible, we validate the use of initial H/D exchange with unfolded protein and partial deuteration in pure D₂O with hydrogenous glycerol to maximize deuterium incorporation into the protein, with no detrimental effects on the growth of quality crystals suitable for neutron diffraction experiments.

Diverse Folding Pathways of HIV-1 Protease Monomer on a Rugged Energy Landscape

The modern energy landscape theory of protein folding predicts multiple folding pathways connecting a myriad of unfolded conformations and a well-defined folded state. However, direct experimental observation of heterogeneous folding pathways is difficult. Naturally evolved proteins typically exhibit a smooth folding energy landscape for fast and efficient folding by avoiding unfavorable kinetic traps. In this case, rapid fluctuations between unfolded conformations result in apparent two-state behavior and make different pathways indistinguishable. However, the landscape roughness can be different, depending on the selection pressures during evolution. Here, we characterize the unusually rugged folding energy landscape of human immunodeficiency virus-1 protease monomer using single-molecule Förster resonance energy transfer spectroscopy. Our data show that fluctuations between unfolded conformations are slow, which enables the experimental observation of heterogeneous folding pathways as predicted by the landscape theory. Although the landscape ruggedness is sensitive to the mutations and fluorophore locations, the folding rate is similar for various protease constructs. The natural evolution of the protease to have a rugged energy landscape likely results from intrinsic pressures to maintain robust folding when human immunodeficiency virus-1 mutates frequently, which is essential for its survival.

Importance of time-ordered non-uniform sampling of multi-dimensional NMR spectra of Aβ1-42 peptide under aggregating conditions

Although the order of the time steps in which the non-uniform sampling (NUS) schedule is implemented when acquiring multi-dimensional NMR spectra is of limited importance when sample conditions remain unchanged over the course of the experiment, it is shown to have major impact when samples are unstable. In the latter case, time-ordering of the NUS data points by the normalized radial length yields a reduction of sampling artifacts, regardless of the spectral reconstruction algorithm. The disadvantage of time-ordered NUS sampling is that halting the experiment prior to its completion will result in lower spectral resolution, rather than a sparser data matrix. Alternatively, digitally correcting for sample decay prior to reconstruction of randomly ordered NUS data points can mitigate reconstruction artifacts, at the cost of somewhat lower sensitivity. Application of these sampling schemes to the Alzheimer's amyloid beta (Aβ1-42) peptide at an elevated concentration, low temperature, and 3 kbar of pressure, where approximately 75% of the peptide reverts to an NMR-invisible state during the collection of a 3D 15N-separated NOESY spectrum, highlights the improvement in artifact suppression and reveals weak medium-range NOE contacts in several regions, including the C-terminal region of the peptide.

Inhibition of HIV Maturation via Selective Unfolding and Cross-Linking of Gag Polyprotein by a Mercaptobenzamide Acetylator

For HIV to become infectious, any new virion produced from an infected cell must undergo a maturation process that involves the assembly of viral polyproteins Gag and Gag-Pol at the membrane surface. The self-assembly of these viral proteins drives formation of a new viral particle as well as the activation of HIV protease, which is needed to cleave the polyproteins so that the final core structure of the virus will properly form. Molecules that interfere with HIV maturation will prevent any new virions from infecting additional cells. In this manuscript, we characterize the unique mechanism by which a mercaptobenzamide thioester small molecule (SAMT-247) interferes with HIV maturation via a series of selective acetylations at highly conserved cysteine and lysine residues in Gag and Gag-Pol polyproteins. The results provide the first insights into how acetylation can be utilized to perturb the process of HIV maturation and reveal a new strategy to limit the infectivity of HIV.

Three-Color Single-Molecule FRET and Fluorescence Lifetime Analysis of Fast Protein Folding

We describe the theory, experiment, and analysis of three-color Förster resonance energy transfer (FRET) spectroscopy for probing conformational dynamics of a fast-folding protein, α₃D. In three-color FRET, site-specific labeling of fluorophores is required to avoid ambiguity resulting from various species with different combinations of labeling positions. To this end, we first attached two dyes to a cysteine residue and an unnatural amino acid and then appended a cysteine residue to the C-terminus of the protein by the sortase-mediated ligation for attaching the third dye. To determine all three FRET efficiencies, we used alternating excitation of the donor and acceptor 1 with two picosecond-pulsed lasers. Since the folded and unfolded states are not distinguishable in binned fluorescence trajectories due to fast-folding on a millisecond time scale, we used a maximum likelihood method that analyzes photon trajectories without binning the data. The extracted kinetic parameters agree very well with the previously measured parameters for the same protein with two-color FRET, suggesting that the addition of the third fluorophore does not affect the folding dynamics of the protein. From the extracted fractions of acceptor photon counts, the FRET efficiencies for all three dye pairs were calculated after various corrections. They were compared with the FRET efficiencies obtained from the global analysis of two-color segments collected in the same experiment. The FRET efficiencies of the folded state from the three-color segments agree with those from the two-color segments, whereas the three-color and two-color FRET efficiencies of the unfolded state are different. This happens because fluctuations of all three interdye distances contribute to the FRET efficiency measured in three-color FRET. We show that this difference can be accounted for by using the Gaussian chain model for the unfolded state with the parameters obtained from the analysis of two-color segments. This result shows that three-color FRET provides additional information on the flexibility of molecules that cannot be obtained from a combination of two-color FRET experiments with three dye pairs. Using the delay times of photons from the laser pulse, fluorescence lifetimes were determined using the maximum likelihood analysis. The correlation between FRET efficiencies and lifetimes of the donor, acceptor 1, and acceptor 2 was visualized in two-dimensional FRET efficiency-lifetime histograms. These histograms can be used to demonstrate the presence of conformational dynamics in a protein.

Co-Evolutionary Fitness Landscapes for Sequence Design

Efficient and accurate models to predict the fitness of a sequence would be extremely valuable in protein design. We have explored the use of statistical potentials for the coevolutionary fitness landscape, extracted from known protein sequences, in conjunction with Monte Carlo simulations, as a tool for design. As proof of principle, we created a series of predicted high-fitness sequences for three different protein folds, representative of different structural classes: the GA (all-α) and GB (α/β) binding domains of streptococcal protein G, and an SH3 (all-β) domain. We found that most of the designed proteins can fold stably to the target structure, and a structure for a representative of each for GA, GB and SH3 was determined. Several of our designed proteins were also able to bind to native ligands, in some cases with higher affinity than wild-type. Thus, a search using a statistical fitness landscape is a remarkably effective tool for finding novel stable protein sequences.

Tilted, Uninterrupted, Monomeric HIV-1 gp41 Transmembrane Helix from Residual Dipolar Couplings

Cryo-electron microscopy and X-ray crystallography have shown that the pre- and postfusion states of the HIV-1 gp41 viral coat protein, although very different from one another, each adopt C₃ symmetric structures. A stable homotrimeric structure for the transmembrane domain (TM) also was modeled and supported by experimental data. For a C₃ symmetric structure, alignment in an anisotropic medium must be axially symmetric, with the unique axis of the alignment tensor coinciding with the C₃ axis. However, NMR residual dipolar couplings (RDCs) measured under three different alignment conditions were found to be incompatible with C₃ symmetry. Subsequent measurements by paramagnetic relaxation enhancement, analytical ultracentrifugation, and DEER EPR, indicate that the transmembrane domain is monomeric. ¹⁵N NMR relaxation data and RDCs show that TM is highly ordered and uninterrupted for a total length of 32 residues, extending well into the membrane proximal external region.

Room Temperature Neutron Crystallography of Drug Resistant HIV-1 Protease Uncovers Limitations of X-ray Structural Analysis at 100 K

HIV-1 protease inhibitors are crucial for treatment of HIV-1/AIDS, but their effectiveness is thwarted by rapid emergence of drug resistance. To better understand binding of clinical inhibitors to resistant HIV-1 protease, we used room-temperature joint X-ray/neutron (XN) crystallography to obtain an atomic-resolution structure of the protease triple mutant (V32I/I47V/V82I) in complex with amprenavir. The XN structure reveals a D⁺ ion located midway between the inner Oδ1 oxygen atoms of the catalytic aspartic acid residues. Comparison of the current XN structure with our previous XN structure of the wild-type HIV-1 protease-amprenavir complex suggests that the three mutations do not significantly alter the drug-enzyme interactions. This is in contrast to the observations in previous 100 K X-ray structures of these complexes that indicated loss of interactions by the drug with the triple mutant protease. These findings, thus, uncover limitations of structural analysis of drug binding using X-ray structures obtained at 100 K.

Structural Studies of a Rationally Selected Multi-Drug Resistant HIV-1 Protease Reveal Synergistic Effect of Distal Mutations on Flap Dynamics

We report structural analysis of HIV protease variant PR^S17 which was rationally selected by machine learning to represent wide classes of highly drug-resistant variants. Crystal structures were solved of PR^S17 in the inhibitor-free form and in complex with antiviral inhibitor, darunavir. Despite its 17 mutations, PR^S17 has only one mutation (V82S) in the inhibitor/substrate binding cavity, yet exhibits high resistance to all clinical inhibitors. PR^S17 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 PR^S17 shows an open flap conformation with a curled tip correlating with G48V flap mutation. NMR studies on inactive PR^S17_D25N unambiguously confirm that the flaps adopt mainly an open conformation in solution very similar to that in the inhibitor-free crystal structure. In PR^S17, 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 PR^S17/DRV complex alter the Phe53 conformation by steric hindrance between the side chains. Unlike the L10F mutation in PR20, L10I in PR^S17 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 PR^S17. PR^S17 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.

Insights into the Conformation of the Membrane Proximal Regions Critical to the Trimerization of the HIV-1 gp41 Ectodomain Bound to Dodecyl Phosphocholine Micelles

The transitioning of the ectodomain of gp41 from a pre-hairpin to a six-helix bundle conformation is a crucial aspect of virus-cell fusion. To gain insight into the intermediary steps of the fusion process we have studied the pH and dodecyl phosphocholine (DPC) micelle dependent trimer association of gp41 by systematic deletion analysis of an optimized construct termed 17-172 (residues 528 to 683 of Env) that spans the fusion peptide proximal region (FPPR) to the membrane proximal external region (MPER) of gp41, by sedimentation velocity and double electron-electron resonance (DEER) EPR spectroscopy. Trimerization at pH 7 requires the presence of both the FPPR and MPER regions. However, at pH 4, the protein completely dissociates to monomers. DEER measurements reveal a partial fraying of the C-terminal MPER residues in the 17-172 trimer while the other regions, including the FPPR, remain compact. In accordance, truncating nine C-terminal MPER residues (675-683) in the 17-172 construct does not shift the trimer-monomer equilibrium significantly. Thus, in the context of the gp41 ectodomain spanning residues 17-172, trimerization is clearly dependent on FPPR and MPER regions even when the terminal residues of MPER unravel. The antibody Z13e1, which spans both the 2F5 and 4E10 epitopes in MPER, binds to 17-172 with a Kd of 1 ± 0.12 μM. Accordingly, individual antibodies 2F5 and 4E10 also recognize the 17-172 trimer/DPC complex. We propose that binding of the C-terminal residues of MPER to the surface of the DPC micelles models a correct positioning of the trimeric transmembrane domain anchored in the viral membrane.

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.

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

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.

Long-Range Electrostatics-Induced Two-Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

Neutron crystallography was used to directly locate two protons before and after a pH-induced two-proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low-pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.

Analysis of Fluorescence Lifetime and Energy Transfer Efficiency in Single-Molecule Photon Trajectories of Fast-Folding Proteins

In single-molecule Förster resonance energy transfer (FRET) spectroscopy, the dynamics of molecular processes are usually determined by analyzing the fluorescence intensity of donor and acceptor dyes. Since FRET efficiency is related to fluorescence lifetimes, additional information can be extracted by analyzing fluorescence intensity and lifetime together. For fast processes where individual states are not well separated in a trajectory, it is not easy to obtain the lifetime information. Here, we present analysis methods to utilize fluorescence lifetime information from single-molecule FRET experiments, and apply these methods to three fast-folding, two-state proteins. By constructing 2D FRET efficiency-lifetime histograms, the correlation can be visualized between the FRET efficiency and fluorescence lifetimes in the presence of the submicrosecond to millisecond dynamics. We extend the previously developed method for analyzing delay times of donor photons to include acceptor delay times. To determine the kinetics and lifetime parameters accurately, we used a maximum likelihood method. We found that acceptor blinking can lead to inaccurate parameters in the donor delay time analysis. This problem can be solved by incorporating acceptor blinking into a model. While the analysis of acceptor delay times is not affected by acceptor blinking, it is more sensitive to the shape of the delay time distribution resulting from a broad conformational distribution in the unfolded state.

The C34 Peptide Fusion Inhibitor Binds to the Six-Helix Bundle Core Domain of HIV-1 gp41 by Displacement of the C-Terminal Helical Repeat Region

The conformational transition of the core domain of HIV-1 gp41 from a prehairpin intermediate to a six-helix bundle is responsible for virus-cell fusion. Several inhibitors which target the N-heptad repeat helical coiled-coil trimer that is fully accessible in the prehairpin intermediate have been designed. One such inhibitor is the peptide C34 derived from the C-heptad repeat of gp41 that forms the exterior of the six-helix bundle. Here, using a variety of biophysical techniques, including dye tagging, size-exclusion chromatography combined with multiangle light scattering, double electron-electron resonance EPR spectroscopy, and circular dichroism, we investigate the binding of C34 to two six-helix bundle mimetics comprising N- and C-heptad repeats either without (core(SP)) or with (core(S)) a short spacer connecting the two. In the case of core(SP), C34 directly exchanges with the C-heptad repeat. For core(S), up to two molecules of C34 bind the six-helix bundle via displacement of the C-heptad repeat. These results suggest that fusion inhibitors such as C34 can target a continuum of transitioning conformational states from the prehairpin intermediate to the six-helix bundle prior to the occurrence of irreversible fusion of viral and target cell membranes.

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.

Evidence of Distinct Channel Conformations and Substrate Binding Affinities for the Mitochondrial Outer Membrane Protein Translocase Pore Tom40

Nearly all mitochondrial proteins are coded by the nuclear genome and must be transported into mitochondria by the translocase of the outer membrane complex. Tom40 is the central subunit of the translocase complex and forms a pore in the mitochondrial outer membrane. To date, the mechanism it utilizes for protein transport remains unclear. Tom40 is predicted to comprise a membrane-spanning β-barrel domain with conserved α-helical domains at both the N and C termini. To investigate Tom40 function, including the role of the N- and C-terminal domains, recombinant forms of the Tom40 protein from the yeast Candida glabrata, and truncated constructs lacking the N- and/or C-terminal domains, were functionally characterized in planar lipid membranes. Our results demonstrate that each of these Tom40 constructs exhibits at least four distinct conductive levels and that full-length and truncated Tom40 constructs specifically interact with a presequence peptide in a concentration- and voltage-dependent manner. Therefore, neither the first 51 amino acids of the N terminus nor the last 13 amino acids of the C terminus are required for Tom40 channel formation or for the interaction with a presequence peptide. Unexpectedly, substrate binding affinity was dependent upon the Tom40 state corresponding to a particular conductive level. A model where two Tom40 pores act in concert as a dimeric protein complex best accounts for the observed biochemical and electrophysiological data. These results provide the first evidence for structurally distinct Tom40 conformations playing a role in substrate recognition and therefore in transport function.

Mutations Proximal to Sites of Autoproteolysis and the α-Helix That Co-evolve under Drug Pressure Modulate the Autoprocessing and Vitality of HIV-1 Protease

N-Terminal self-cleavage (autoprocessing) of the HIV-1 protease precursor is crucial for liberating the active dimer. Under drug pressure, evolving mutations are predicted to modulate autoprocessing, and the reduced catalytic activity of the mature protease (PR) is likely compensated by enhanced conformational/dimer stability and reduced susceptibility to self-degradation (autoproteolysis). One such highly evolved, multidrug resistant protease, PR20, bears 19 mutations contiguous to sites of autoproteolysis in retroviral proteases, namely clusters 1-3 comprising residues 30-37, 60-67, and 88-95, respectively, accounting for 11 of the 19 mutations. By systematically replacing corresponding clusters in PR with those of PR20, and vice versa, we assess their influence on the properties mentioned above and observe no strict correlation. A 10-35-fold decrease in the cleavage efficiency of peptide substrates by PR20, relative to PR, is reflected by an only ∼4-fold decrease in the rate of Gag processing with no change in cleavage order. Importantly, optimal N-terminal autoprocessing requires all 19 PR20 mutations as evaluated in vitro using the model precursor TFR-PR20 in which PR is flanked by the transframe region. Substituting PR20 cluster 3 into TFR-PR (TFR-PR(PR20-3)) requires the presence of PR20 cluster 1 and/or 2 for autoprocessing. In accordance, substituting PR clusters 1 and 2 into TFR-PR20 affects the rate of autoprocessing more drastically (>300-fold) compared to that of TFR-PR(PR20-3) because of the cumulative effect of eight noncluster mutations present in TFR-PR20(PR-12). Overall, these studies imply that drug resistance involves a complex synchronized selection of mutations modulating all of the properties mentioned above governing PR regulation and function.

Testing Landscape Theory for Biomolecular Processes with Single Molecule Fluorescence Spectroscopy

Although Kramers' theory for diffusive barrier crossing on a 1D free energy profile plays a central role in landscape theory for complex biomolecular processes, it has not yet been rigorously tested by experiment. Here we test this 1D diffusion scenario with single molecule fluorescence measurements of DNA hairpin folding. We find an upper bound of 2.5  μs for the average transition path time, consistent with the predictions by theory with parameters determined from optical tweezer measurements.

Substituted Bis-THF Protease Inhibitors with Improved Potency against Highly Resistant Mature HIV-1 Protease PR20

An extremely drug resistant mutant of HIV-1 protease (PR) bearing 20 mutations (PR20) has been studied with two potent antiviral investigational inhibitors. GRL-5010A and GRL-4410A were designed to introduce hydrogen bond interactions with the flexible flaps of the PR by incorporating gem-difluorines and alkoxy, respectively, at the C4 position of the bis-THF of darunavir. PR20 provides an excellent model for high level resistance, since clinical inhibitors are >1000-fold less active on PR20 than on wild-type enzyme. GRL-5010A and GRL-4410A show inhibition constants of 4.3 ± 7.0 and 1.7 ± 1.8 nM, respectively, for PR20, compared to the binding affinity of 41 ± 1 nM measured for darunavir. Crystal structures of PR20 in complexes with the two inhibitors confirmed the new hydrogen bond interactions with Gly 48 in the flap of the enzyme. The two new compounds are more effective than darunavir in inhibiting mature PR20 and show promise for further development of antiviral agents targeting highly resistant PR mutants.

Dependence of distance distributions derived from double electron-electron resonance pulsed EPR spectroscopy on pulse-sequence time

Pulsed double electron-electron resonance (DEER) provides pairwise P(r) distance distributions in doubly spin labeled proteins. We report that in protonated proteins, P(r) is dependent on the length of the second echo period T owing to local environmental effects on the spin-label phase memory relaxation time Tm . For the protein ABD, this effect results in a 1.4 Å increase in the P(r) maximum from T=6 to 20 μs. Protein A has a bimodal P(r) distribution, and the relative height of the shorter distance peak at T=10 μs, the shortest value required to obtain a reliable P(r), is reduced by 40 % relative to that found by extrapolation to T=0. Our results indicate that data at a series of T values are essential for quantitative interpretation of DEER to determine the extent of the T dependence and to extrapolate the results to T=0. Complete deuteration (99 %) of the protein was accompanied by a significant increase in Tm and effectively abolished the P(r) dependence on T.