Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations.

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

Basic Residues at Position 11 of α-Conotoxin LvIA Influence Subtype Selectivity between α3ß2 and α3ß4 Nicotinic Receptors via an Electrostatic Mechanism.

Understanding the determinants of α-conotoxin (α-CTX) selectivity for different nicotinic acetylcholine receptor (nAChR) subtypes is a prerequisite for the design of tool compounds to study nAChRs. However, selectivity optimization of these small, disulfide-rich peptides is difficult not only because of an absence of α-CTX/nAChR co-structures but also because it is challenging to predict how a mutation to an α-CTX will alter its potency and selectivity. As a prototypical system to investigate selectivity, we employed the α-CTX LvIA that is 25-fold selective for the α3ß2 nAChR over the related α3ß4 nAChR subtype, which is a target for nicotine addiction. Using two-electrode voltage clamp electrophysiology, we identified LvIA[D11R] that is 2-fold selective for the α3ß4 nAChR, reversing the subtype preference. This effect is specifically due to the change in charge and not shape of LvIA[D11R], as substitution of D11 with citrulline retains selectivity for the α3ß2 nAChR. Furthermore, LvIA[D11K] shows a stronger reversal, with 4-fold selectivity for the α3ß4 nAChR. Motivated by these findings, using site-directed mutagenesis, we found that ß2[K79A] (I79 on ß4), but not ß2[K78A] (N78 on ß4), largely restores the potency of basic mutants at position 11. Finally, to understand the structural basis of this effect, we used AlphaFold2 to generate models of LvIA in complex with both nAChR subtypes. Both models confirm the plausibility of an electrostatic mechanism to explain the data and also reproduce a broad range of potency and selectivity structure-activity relationships for LvIA mutants, as measured using free energy perturbation simulations. Our work highlights how electrostatic interactions can drive α-CTX selectivity and may serve as a strategy for optimizing the selectivity of LvIA and other α-CTXs.

Free Energy Perturbation Calculations of Mutation Effects on SARS-CoV-2 RBD::ACE2 Binding Affinity.

The strength of binding between human angiotensin converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of viral spike protein plays a role in the transmissibility of the SARS-CoV-2 virus. In this study we focus on a subset of RBD mutations that have been frequently observed in infected individuals and probe binding affinity changes to ACE2 using surface plasmon resonance (SPR) measurements and free energy perturbation (FEP) calculations. Our SPR results are largely in accord with previous studies but discrepancies do arise due to differences in experimental methods and to protocol differences even when a single method is used. Overall, we find that FEP performance is superior to that of other computational approaches examined as determined by agreement with experiment and, in particular, by its ability to identify stabilizing mutations. Moreover, the calculations successfully predict the observed cooperative stabilization of binding by the Q498R N501Y double mutant present in Omicron variants and offer a physical explanation for the underlying mechanism. Overall, our results suggest that despite the significant computational cost, FEP calculations may offer an effective strategy to understand the effects of interfacial mutations on protein-protein binding affinities and, hence, in a variety of practical applications such as the optimization of neutralizing antibodies.

Cryo-EM Structures of SARS-CoV-2 Spike without and with ACE2 Reveal a pH-Dependent Switch to Mediate Endosomal Positioning of Receptor-Binding Domains.

The SARS-CoV-2 spike employs mobile receptor-binding domains (RBDs) to engage the human ACE2 receptor and to facilitate virus entry, which can occur through low-pH-endosomal pathways. To understand how ACE2 binding and low pH affect spike conformation, we determined cryo-electron microscopy structures-at serological and endosomal pH-delineating spike recognition of up to three ACE2 molecules. RBDs freely adopted "up" conformations required for ACE2 interaction, primarily through RBD movement combined with smaller alterations in neighboring domains. In the absence of ACE2, single-RBD-up conformations dominated at pH 5.5, resolving into a solitary all-down conformation at lower pH. Notably, a pH-dependent refolding region (residues 824-858) at the spike-interdomain interface displayed dramatic structural rearrangements and mediated RBD positioning through coordinated movements of the entire trimer apex. These structures provide a foundation for understanding prefusion-spike mechanics governing endosomal entry; we suggest that the low pH all-down conformation potentially facilitates immune evasion from RBD-up binding antibody.

Cryo-EM Structures Delineate a pH-Dependent Switch that Mediates Endosomal Positioning of SARS-CoV-2 Spike Receptor-Binding Domains.

The SARS-CoV-2 spike employs mobile receptor-binding domains (RBDs) to engage the ACE2 receptor and to facilitate virus entry. Antibodies can engage RBD but some, such as CR3022, fail to inhibit entry despite nanomolar spike affinity. Here we show the SARS-CoV-2 spike to have low unfolding enthalpy at serological pH and up to 10-times more unfolding enthalpy at endosomal pH, where we observe significantly reduced CR3022 affinity. Cryo-EM structures -at serological and endosomal pH- delineated spike recognition of up to three ACE2 molecules, revealing RBD to freely adopt the 'up' conformation. In the absence of ACE2, single-RBD-up conformations dominated at pH 5.5, resolving into a locked all-down conformation at lower pH. Notably, a pH-dependent refolding region (residues 824-858) at the spike-interdomain interface displayed dramatic structural rearrangements and mediated RBD positioning and spike shedding of antibodies like CR3022. An endosomal mechanism involving spike-conformational change can thus facilitate immune evasion from RBD-'up'-recognizing antibody.

Rationally Designed Immunogens Targeting HIV-1 gp120 V1V2 Induce Distinct Conformation-Specific Antibody Responses in Rabbits.

The V1V2 region of HIV-1 gp120 harbors a major vulnerable site targeted by a group of broadly neutralizing monoclonal antibodies (MAbs) such as PG9 through strand-strand recognition. However, this epitope region is structurally polymorphic as it can also form a helical conformation recognized by RV144 vaccine-induced MAb CH58. This structural polymorphism is a potential mechanism for masking the V1V2 vulnerable site. Designing immunogens that can induce conformation-specific antibody (Ab) responses may lead to vaccines targeting this vulnerable site. We designed a panel of immunogens engrafting the V1V2 domain into trimeric and pentameric scaffolds in structurally constrained conformations. We also fused V1V2 to an Fc fragment to mimic the unconstrained V1V2 conformation. We tested these V1V2-scaffold proteins for immunogenicity in rabbits and assessed the responses by enzyme-linked immunosorbent assay (ELISA) and competition assays. Our V1V2 immunogens induced distinct conformation-specific Ab responses. Abs induced by structurally unconstrained immunogens reacted preferentially with unconstrained V1V2 antigens, suggesting recognition of the helical configuration, while Abs induced by the structurally constrained immunogens reacted preferentially with constrained V1V2 antigens, suggesting recognition of the ß-strand conformation. The Ab responses induced by the structurally constrained immunogens were more broadly reactive and had higher titers than those induced by the structurally unconstrained immunogens. Our results demonstrate that immunogens presenting the different structural conformations of the gp120 V1V2 vulnerable site can be designed and that these immunogens induce distinct Ab responses with epitope conformation specificity. Therefore, these structurally constrained V1V2 immunogens are vaccine prototypes targeting the V1V2 domain of the HIV-1 envelope. IMPORTANCE: The correlates analysis of the RV144 HIV-1 vaccine trial suggested that the presence of antibodies to the V1V2 region of HIV-1 gp120 was responsible for the modest protection observed in the trial. In addition, V1V2 harbors one of the key vulnerable sites of HIV-1 Env recognized by a family of broadly neutralizing MAbs such as PG9. Thus, V1V2 is a key target for vaccine development. However, this vulnerable site is structurally polymorphic, and designing immunogens that present different conformations is crucial for targeting this site. We show here that such immunogens can be designed and that they induced conformation-specific antibody responses in rabbits. Our immunogens are therefore prototypes of vaccine candidates targeting the V1V2 region of HIV-1 Env.

Structural Bases of Noncompetitive Inhibition of AMPA-Subtype Ionotropic Glutamate Receptors by Antiepileptic Drugs.

Excitatory neurotransmission plays a key role in epileptogenesis. Correspondingly, AMPA-subtype ionotropic glutamate receptors, which mediate the majority of excitatory neurotransmission and contribute to seizure generation and spread, have emerged as promising targets for epilepsy therapy. The most potent and well-tolerated AMPA receptor inhibitors act via a noncompetitive mechanism, but many of them produce adverse side effects. The design of better drugs is hampered by the lack of a structural understanding of noncompetitive inhibition. Here, we report crystal structures of the rat AMPA-subtype GluA2 receptor in complex with three noncompetitive inhibitors. The inhibitors bind to a novel binding site, completely conserved between rat and human, at the interface between the ion channel and linkers connecting it to the ligand-binding domains. We propose that the inhibitors stabilize the AMPA receptor closed state by acting as wedges between the transmembrane segments, thereby preventing gating rearrangements that are necessary for ion channel opening.

Structure-based functional characterization of repressor of toxin (Rot), a central regulator of Staphylococcus aureus virulence.

Staphylococcus aureus is responsible for a large number of diverse infections worldwide. In order to support its pathogenic lifestyle, S. aureus has to regulate the expression of virulence factors in a coordinated fashion. One of the central regulators of the S. aureus virulence regulatory networks is the transcription factor repressor of toxin (Rot). Rot plays a key role in regulating S. aureus virulence through activation or repression of promoters that control expression of a large number of critical virulence factors. However, the mechanism by which Rot mediates gene regulation has remained elusive. Here, we have determined the crystal structure of Rot and used this information to probe the contribution made by specific residues to Rot function. Rot was found to form a dimer, with each monomer harboring a winged helix-turn-helix (WHTH) DNA-binding motif. Despite an overall acidic pI, the asymmetric electrostatic charge profile suggests that Rot can orient the WHTH domain to bind DNA. Structure-based site-directed mutagenesis studies demonstrated that R(91), at the tip of the wing, plays an important role in DNA binding, likely through interaction with the minor groove. We also found that Y(66), predicted to bind within the major groove, contributes to Rot interaction with target promoters. Evaluation of Rot binding to different activated and repressed promoters revealed that certain mutations on Rot exhibit promoter-specific effects, suggesting for the first time that Rot differentially interacts with target promoters. This work provides insight into a precise mechanism by which Rot controls virulence factor regulation in S. aureus.

Rabbit anti-HIV-1 monoclonal antibodies raised by immunization can mimic the antigen-binding modes of antibodies derived from HIV-1-infected humans.

The rabbit is a commonly used animal model in studying antibody responses in HIV/AIDS vaccine development. However, no rabbit monoclonal antibodies (MAbs) have been developed previously to study the epitope-specific antibody responses against HIV-1 envelope (Env) glycoproteins, and little is known about how the rabbit immune system can mimic the human immune system in eliciting such antibodies. Here we present structural analyses of two rabbit MAbs, R56 and R20, against the third variable region (V3) of HIV-1 gp120. R56 recognizes the well-studied immunogenic region in the V3 crown, while R20 targets a less-studied region at the C terminus of V3. By comparison of the Fab/epitope complex structures of these two antibodies raised by immunization with that of the corresponding human antibodies derived from patients chronically infected with HIV-1, we found that rabbit antibodies can recognize immunogenic regions of gp120 and mimic the binding modes of human antibodies. This result can provide new insight into the use of the rabbit as an animal model in AIDS vaccine development.

Skin tattooing as a novel approach for DNA vaccine delivery.

Nucleic acid-based vaccination is a topic of growing interest, especially plasmid DNA (pDNA) encoding immunologically important antigens. After the engineered pDNA is administered to the vaccines, it is transcribed and translated into immunogen proteins that can elicit responses from the immune system. Many ways of delivering DNA vaccines have been investigated; however each delivery route has its own advantages and pitfalls. Skin tattooing is a novel technique that is safe, cost-effective, and convenient. In addition, the punctures inflicted by the needle could also serve as a potent adjuvant. Here, we a) demonstrate the intradermal delivery of plasmid DNA encoding enhanced green fluorescent protein (pCX-EGFP) in a mouse model using a tattooing device and b) confirm the effective expression of EGFP in the skin cells using confocal microscopy.

Crystal structure of a preacylation complex of the ß-lactamase inhibitor sulbactam bound to a sulfenamide bond-containing thiol-ß-lactamase.

The rise of inhibitor-resistant and other ß-lactamase variants is generating an interest in developing new ß-lactamase inhibitors to complement currently available antibiotics. To gain insight into the chemistry of inhibitor recognition, we determined the crystal structure of the inhibitor preacylation complex of sulbactam, a clinical ß-lactamase inhibitor, bound in the active site of the S70C variant of SHV-1 ß-lactamase, a resistance enzyme that is normally present in Klebsiella pneumoniae. The S70C mutation was designed to affect the reactivity of that catalytic residue to allow for capture of the preacylation complex. Unexpectedly, the 1.45 Å resolution inhibitor complex structure revealed that residue C70 is involved in a sulfenamide bond with K73. Such a covalent bond is not present in the wild-type SHV-1 or in an apo S70C structure also determined in this study. This bond likely contributed significantly to obtaining the preacylation complex with sulbactam due to further decreased reactivity toward substrates. The intact sulbactam is positioned in the active site such that its carboxyl moiety interacts with R244, S130, and T235 and its carbonyl moiety is situated in the oxyanion hole. To our knowledge, in addition to being the first preacylation inhibitor ß-lactamase complex, this is also the first observation of a sulfenamide bond between a cysteine and lysine in an active site. Not only could our results aid, therefore, structure-based inhibitor design efforts in class A ß-lactamases, but the sulfenamide-bond forming approach to yield preacylation complexes could also be applied to other classes of ß-lactamases and penicillin-binding proteins with the SXXK motif.

The importance of the trans-enamine intermediate as a ß-lactamase inhibition strategy probed in inhibitor-resistant SHV ß-lactamase variants.

The ability of bacteria to express inhibitor-resistant (IR) ß-lactamases is stimulating the development of novel inhibitors of these enzymes. The 2'ß-glutaroxypenicillinate sulfone, SA2-13, was previously designed to enhance the stabilization of the deacylation-refractory, trans-enamine inhibitory intermediate. To test whether this mode of inhibition can overcome different IR mutations, we determined the binding mode of SA2-13 through X-ray crystallography, obtaining co-crystals of the inhibitor-protein complex by soaking crystals of the IR sulfhydryl variable (SHV) ß-lactamase variants S130G and M69V with the inhibitor. The 1.45 Å crystal structure of the S130G SHV:SA2-13 complex reveals that SA2-13 is still able to form the stable trans-enamine intermediate similar to the wild-type complex structure, yet with its carboxyl linker shifted deeper into the active site in the space vacated by the S130G mutation. In contrast, data from crystals of the M69V SHV:SA2-13 complex at 1.3 Å did not reveal clear inhibitor density indicating that this IR variant disfavors the trans-enamine conformation, likely due to a subtle shift in A237.

Functional and immunochemical cross-reactivity of V2-specific monoclonal antibodies from HIV-1-infected individuals.

The recent analysis of the first successful RV144 vaccine trial revealed that a high titer of plasma anti-V2 antibodies (Abs) correlated with a decreased risk of HIV-1 infection in vaccine recipients. To understand the mechanism of immune correlates, we studied seven anti-V2 monoclonal Abs (mAbs) developed from HIV-1 infected individuals. The V2 mAbs target conserved epitopes, including the binding site for α4ß7 integrin, and are broadly cross-reactive with various gp120 proteins. Preferential usage of the VH1-69 gene by V2 mAbs may depend on selection by the same antigenic structure. Six of seven V2 mAbs weakly neutralized four to eight of the 41 pseudoviruses tested and resistance to neutralization was correlated with longer V2 domains. The data suggest the presence of shared, conserved structural elements in the V2 loop, and these can be used in the design of vaccine immunogens inducing broadly reactive Abs with anti-viral activities.

Structural analysis of human and macaque mAbs 2909 and 2.5B: implications for the configuration of the quaternary neutralizing epitope of HIV-1 gp120.

The quaternary neutralizing epitope (QNE) of HIV-1 gp120 is preferentially expressed on the trimeric envelope spikes of intact HIV virions, and QNE-specific monoclonal antibodies (mAbs) potently neutralize HIV-1. Here, we present the crystal structures of the Fabs of human mAb 2909 and macaque mAb 2.5B. Both mAbs have long beta hairpin CDR H3 regions >20 Å in length that are each situated at the center of their respective antigen-binding sites. Computational analysis showed that the paratopes include the whole CDR H3, while additional CDR residues form shallow binding pockets. Structural modeling suggests a way to understand the configuration of QNEs and the antigen-antibody interaction for QNE mAbs. Our data will be useful in designing immunogens that may elicit potent neutralizing QNE Abs.

Ligand-dependent disorder of the Omega loop observed in extended-spectrum SHV-type beta-lactamase.

Among Gram-negative bacteria, resistance to ß-lactams is mediated primarily by ß-lactamases (EC 3.2.6.5), periplasmic enzymes that inactivate ß-lactam antibiotics. Substitutions at critical amino acid positions in the class A ß-lactamase families result in enzymes that can hydrolyze extended-spectrum cephalosporins, thus demonstrating an "extended-spectrum" ß-lactamase (ESBL) phenotype. Using SHV ESBLs with substitutions in the Ω loop (R164H and R164S) as target enzymes to understand this enhanced biochemical capability and to serve as a basis for novel ß-lactamase inhibitor development, we determined the spectra of activity and crystal structures of these variants. We also studied the inactivation of the R164H and R164S mutants with tazobactam and SA2-13, a unique ß-lactamase inhibitor that undergoes a distinctive reaction chemistry in the active site. We noted that the reduced Ki values for the R164H and R164S mutants with SA2-13 are comparable to those with tazobactam (submicromolar). The apo enzyme crystal structures of the R164H and R164S SHV variants revealed an ordered Ω loop architecture that became disordered when SA2-13 was bound. Important structural alterations that result from the binding of SA2-13 explain the enhanced susceptibility of these ESBL enzymes to this inhibitor and highlight ligand-dependent Ω loop flexibility as a mechanism for accommodating and hydrolyzing ß-lactam substrates.

Novel insights into the mode of inhibition of class A SHV-1 beta-lactamases revealed by boronic acid transition state inhibitors.

Boronic acid transition state inhibitors (BATSIs) are potent class A and C ß-lactamase inactivators and are of particular interest due to their reversible nature mimicking the transition state. Here, we present structural and kinetic data describing the inhibition of the SHV-1 ß-lactamase, a clinically important enzyme found in Klebsiella pneumoniae, by BATSI compounds possessing the R1 side chains of ceftazidime and cefoperazone and designed variants of the latter, compounds 1 and 2. The ceftazidime and cefoperazone BATSI compounds inhibit the SHV-1 ß-lactamase with micromolar affinity that is considerably weaker than their inhibition of other ß-lactamases. The solved crystal structures of these two BATSIs in complex with SHV-1 reveal a possible reason for SHV-1's relative resistance to inhibition, as the BATSIs adopt a deacylation transition state conformation compared to the usual acylation transition state conformation when complexed to other ß-lactamases. Active-site comparison suggests that these conformational differences might be attributed to a subtle shift of residue A237 in SHV-1. The ceftazidime BATSI structure revealed that the carboxyl-dimethyl moiety is positioned in SHV-1's carboxyl binding pocket. In contrast, the cefoperazone BATSI has its R1 group pointing away from the active site such that its phenol moiety moves residue Y105 from the active site via end-on stacking interactions. To work toward improving the affinity of the cefoperazone BATSI, we synthesized two variants in which either one or two extra carbons were added to the phenol linker. Both variants yielded improved affinity against SHV-1, possibly as a consequence of releasing the strain of its interaction with the unusual Y105 conformation.

Design, synthesis, and crystal structures of 6-alkylidene-2'-substituted penicillanic acid sulfones as potent inhibitors of Acinetobacter baumannii OXA-24 carbapenemase.

Class D ß-lactamases represent a growing and diverse class of penicillin-inactivating enzymes that are usually resistant to commercial ß-lactamase inhibitors. As many such enzymes are found in multi-drug resistant (MDR) Acinetobacter baumannii and Pseudomonas aeruginosa, novel ß-lactamase inhibitors are urgently needed. Five unique 6-alkylidene-2'-substituted penicillanic acid sulfones (1-5) were synthesized and tested against OXA-24, a clinically important ß-lactamase that inactivates carbapenems and is found in A. baumannii. Based upon the roles Tyr112 and Met223 play in the OXA-24 ß-lactamase, we also engineered two variants (Tyr112Ala and Tyr112Ala,Met223Ala) to test the hypothesis that the hydrophobic tunnel formed by these residues influences inhibitor recognition. IC(50) values against OXA-24 and two OXA-24 ß-lactamase variants ranged from 10 ± 1 (4 vs WT) to 338 ± 20 nM (5 vs Tyr112Ala, Met223Ala). Compound 4 possessed the lowest K(i) (500 ± 80 nM vs WT), and 1 possessed the highest inactivation efficiency (k(inact)/K(i) = 0.21 ± 0.02 µM(-1) s(-1)). Electrospray ionization mass spectrometry revealed a single covalent adduct, suggesting the formation of an acyl-enzyme intermediate. X-ray structures of OXA-24 complexed to four inhibitors (2.0-2.6 Å) reveal the formation of stable bicyclic aromatic intermediates with their carbonyl oxygen in the oxyanion hole. These data provide the first structural evidence that 6-alkylidene-2'-substituted penicillin sulfones are effective mechanism-based inactivators of class D ß-lactamases. Their unique chemistry makes them developmental candidates. Mechanisms for class D hydrolysis and inhibition are discussed, and a pathway for the evolution of the BlaR1 sensor of Staphylococcus aureus to the class D ß-lactamases is proposed.