Analysis of EGFR binding hotspots for design of new EGFR inhibitory biologics.

The epidermal growth factor (EGF) receptor (EGFR) is activated by the binding of one of seven EGF-like ligands to its ectodomain. Ligand binding results in EGFR dimerization and stabilization of the active receptor conformation subsequently leading to activation of downstream signaling. Aberrant activation of EGFR contributes to cancer progression through EGFR overexpression/amplification, modulation of its positive and negative regulators, and/or activating mutations within EGFR. EGFR targeted therapeutic antibodies prevent dimerization and interaction with endogenous ligands by binding the ectodomain of EGFR. However, these antibodies have had limited success in the clinic, partially due to EGFR ectodomain resistance mutations, and are only applicable to a subset of patients with EGFR-driven cancers. These limitations suggest that alternative EGFR targeted biologics need to be explored for EGFR-driven cancer therapy. To this end, we analyze the EGFR interfaces of known inhibitory biologics with determined structures in the context of endogenous ligands, using the Rosetta macromolecular modeling software to highlight the most important interactions on a per-residue basis. We use this analysis to identify the structural determinants of EGFR targeted biologics. We suggest that commonly observed binding motifs serve as the basis for rational design of new EGFR targeted biologics, such as peptides, antibodies, and nanobodies.

Dark nanodiscs for evaluating membrane protein thermostability by differential scanning fluorimetry.

Measuring protein thermostability provides valuable information on the biophysical rules that govern the structure-energy relationships of proteins. However, such measurements remain a challenge for membrane proteins. Here, we introduce a new experimental system to evaluate membrane protein thermostability. This system leverages a recently developed nonfluorescent membrane scaffold protein to reconstitute proteins into nanodiscs and is coupled with a nano-format of differential scanning fluorimetry (nanoDSF). This approach offers a label-free and direct measurement of the intrinsic tryptophan fluorescence of the membrane protein as it unfolds in solution without signal interference from the "dark" nanodisc. In this work, we demonstrate the application of this method using the disulfide bond formation protein B (DsbB) as a test membrane protein. NanoDSF measurements of DsbB reconstituted in dark nanodiscs loaded with 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dimyristoyl-sn-glycero-3-phosphorylglycerol (DMPG) lipids show a complex biphasic thermal unfolding pattern with a minor unfolding transition followed by a major transition. The inflection points of the thermal denaturation curve reveal two distinct unfolding midpoint melting temperatures (Tm) of 70.5°C and 77.5°C, consistent with a three-state unfolding model. Further, we show that the catalytically conserved disulfide bond between residues C41 and C130 drives the intermediate state of the unfolding pathway for DsbB in a DMPC and DMPG nanodisc. To extend the utility of this method, we evaluate and compare the thermostability of DsbB in different lipid environments. We introduce this method as a new tool that can be used to understand how compositionally and biophysically complex lipid environments drive membrane protein stability.

Dark nanodiscs as a model membrane for evaluating membrane protein thermostability by differential scanning fluorimetry.

Measuring protein thermostability provides valuable information on the biophysical rules that govern structure-energy relationships of proteins. However, such measurements remain a challenge for membrane proteins. Here, we introduce a new experimental system to evaluate membrane protein thermostability. This system leverages a recently-developed non-fluorescent membrane scaffold protein (MSP) to reconstitute proteins into nanodiscs and is coupled with a nano-format of differential scanning fluorimetry (nanoDSF). This approach offers a label-free and direct measurement of the intrinsic tryptophan fluorescence of the membrane protein as it unfolds in solution without signal interference from the "dark" nanodisc. In this work, we demonstrate the application of this method using the disulfide bond formation protein B (DsbB) as a test membrane protein. NanoDSF measurements of DsbB reconstituted in dark nanodiscs show a complex biphasic thermal unfolding pattern in the presence of lipids with a minor unfolding transition followed by a major transition. The inflection points of the thermal denaturation curve reveal two distinct unfolding midpoint melting temperatures (Tm) of 70.5 °C and 77.5 °C, consistent with a three-state unfolding model. Further, we show that the catalytically conserved disulfide bond between residues C41 and C130 drives the intermediate state of the unfolding pathway for DsbB in a nanodisc. We introduce this method as a new tool that can be used to understand how compositionally, and biophysically complex lipid environments drive membrane protein stability.

Computational modeling and prediction of deletion mutants.

In-frame deletion mutations can result in disease. The impact of these mutations on protein structure and subsequent functional changes remain understudied, partially due to the lack of comprehensive datasets including a structural readout. In addition, the recent breakthrough in structure prediction through deep learning demands an update of computational deletion mutation prediction. In this study, we deleted individually every residue of a small α-helical sterile alpha motif domain and investigated the structural and thermodynamic changes using 2D NMR spectroscopy and differential scanning fluorimetry. Then, we tested computational protocols to model and classify observed deletion mutants. We show a method using AlphaFold2 followed by RosettaRelax performs the best overall. In addition, a metric containing pLDDT values and Rosetta ΔΔG is most reliable in classifying tolerated deletion mutations. We further test this method on other datasets and show they hold for proteins known to harbor disease-causing deletion mutations.

Sparse pseudocontact shift NMR data obtained from a non-canonical amino acid-linked lanthanide tag improves integral membrane protein structure prediction.

A single experimental method alone often fails to provide the resolution, accuracy, and coverage needed to model integral membrane proteins (IMPs). Integrating computation with experimental data is a powerful approach to supplement missing structural information with atomic detail. We combine RosettaNMR with experimentally-derived paramagnetic NMR restraints to guide membrane protein structure prediction. We demonstrate this approach using the disulfide bond formation protein B (DsbB), an α-helical IMP. Here, we attached a cyclen-based paramagnetic lanthanide tag to an engineered non-canonical amino acid (ncAA) using a copper-catalyzed azide-alkyne cycloaddition (CuAAC) click chemistry reaction. Using this tagging strategy, we collected 203 backbone HN pseudocontact shifts (PCSs) for three different labeling sites and used these as input to guide de novo membrane protein structure prediction protocols in Rosetta. We find that this sparse PCS dataset combined with 44 long-range NOEs as restraints in our calculations improves structure prediction of DsbB by enhancements in model accuracy, sampling, and scoring. The inclusion of this PCS dataset improved the Cα-RMSD transmembrane segment values of the best-scoring and best-RMSD models from 9.57 Å and 3.06 Å (no NMR data) to 5.73 Å and 2.18 Å, respectively.

Epitope-focused immunogen design based on the ebolavirus glycoprotein HR2-MPER region.

The three human pathogenic ebolaviruses: Zaire (EBOV), Bundibugyo (BDBV), and Sudan (SUDV) virus, cause severe disease with high fatality rates. Epitopes of ebolavirus glycoprotein (GP) recognized by antibodies with binding breadth for all three ebolaviruses are of major interest for rational vaccine design. In particular, the heptad repeat 2 -membrane-proximal external region (HR2-MPER) epitope is relatively conserved between EBOV, BDBV, and SUDV GP and targeted by human broadly-neutralizing antibodies. To study whether this epitope can serve as an immunogen for the elicitation of broadly-reactive antibody responses, protein design in Rosetta was employed to transplant the HR2-MPER epitope identified from a co-crystal structure with the known broadly-reactive monoclonal antibody (mAb) BDBV223 onto smaller scaffold proteins. From computational analysis, selected immunogen designs were produced as recombinant proteins and functionally validated, leading to the identification of a sterile alpha motif (SAM) domain displaying the BDBV-HR2-MPER epitope near its C terminus as a promising candidate. The immunogen was fused to one component of a self-assembling, two-component nanoparticle and tested for immunogenicity in rabbits. Robust titers of cross-reactive serum antibodies to BDBV and EBOV GPs and moderate titers to SUDV GP were induced following immunization. To confirm the structural composition of the immunogens, solution NMR studies were conducted and revealed structural flexibility in the C-terminal residues of the epitope. Overall, our study represents the first report on an epitope-focused immunogen design based on the structurally challenging BDBV-HR2-MPER epitope.

Confirmation of association of TGFBI p.Ser591Phe mutation with variant lattice corneal dystrophy.

PURPOSE: To provide the initial confirmation of the c.1772C>T (p.Ser591Phe) mutation in the transforming growth factor-ß-induced (TGFBI) gene as being associated with variant lattice corneal dystrophy (LCD). METHODS: Ophthalmologic examination of the proband was performed with slit lamp biomicroscopy. Saliva was collected as a source of DNA for screening all 17 exons of TGFBI, after which three family members were selectively screened for variants in exon 13. Rosetta-based structure prediction was used to calculate changes in TGFBI protein (TGFBIp) stability secondary to the c.1772C>T (p.Ser591Phe) missense mutation. RESULTS: Slit lamp examination of the 38-year-old proband revealed a clear cornea right eye and unilateral, discrete, and branching lattice lines in the anterior and mid-stroma of the central cornea left eye. Screening of TGFBI in the proband revealed a heterozygous missense mutation in exon 13 (c.1772C>T (p.Ser591Phe)) that was also identified in her affected mother but not in her brother or maternal grandmother. Calculated energy change in Rosetta (ΔΔG) for the TGFBIp variant p.Ser591Phe was 23.5, indicating a thermodynamic destabilization resulting from energetic frustration. CONCLUSIONS: The p.Ser591Phe mutation in TGFBI is associated with an unilateral variant of LCD. Rosetta-predicted stability changes indicate that the p.Ser591Phe variant is destabilizing, which is consistent with other observations for LCD-causing mutations.

A Conformationally Gated Model of Methadone and Loperamide Transport by P-Glycoprotein.

P-glycoprotein (Pgp) is a multidrug resistance transporter that limits the penetration of a wide range of neurotherapeutics into the brain including opioids. The diphenylpropylamine opioids methadone and loperamide are structurally similar, but loperamide has about a 4-fold higher Pgp-mediated transport rate. In addition to these differences, they showed significant differences in their effects on Pgp-mediated adenosine triphosphate (ATP) hydrolysis. The activation of Pgp-mediated ATP hydrolysis by methadone was monophasic, whereas loperamide activation of ATP hydrolysis was biphasic implying methadone has a single binding site and loperamide has 2 binding sites on Pgp. Quenching of tryptophan fluorescence with these drugs and digoxin showed competition between the opioids and that loperamide does not compete for the digoxin-binding site. Acrylamide quenching of tryptophan fluorescence to probe Pgp conformational changes revealed that methadone- and loperamide-induced conformational changes were distinct. These results were used to develop a model for Pgp-mediated transport of methadone and loperamide where opioid binding and conformational changes are used to explain the differences in the opioid transport rates between methadone and loperamide.

Cardiovascular Ion Channel Inhibitor Drug-Drug Interactions with P-glycoprotein.

P-glycoprotein (Pgp) is an ATP-binding cassette (ABC) transporter that plays a major role in cardiovascular drug disposition by effluxing a chemically and structurally diverse range of cardiovascular therapeutics. Unfortunately, drug-drug interactions (DDIs) with the transporter have become a major roadblock to effective cardiovascular drug administration because they can cause adverse drug reactions (ADRs) or reduce the efficacy of drugs. Cardiovascular ion channel inhibitors are particularly susceptible to DDIs and ADRs with Pgp because they often have low therapeutic indexes and are commonly coadministered with other drugs that are also Pgp substrates. DDIs from cardiovascular ion channel inhibitors with the transporter occur because of inhibition or induction of the transporter and the transporter's tissue and cellular localization. Inhibiting Pgp can increase absorption and reduce excretion of drugs, leading to elevated drug plasma concentrations and drug toxicity. In contrast, inducing Pgp can have the opposite effect by reducing the drug plasma concentration and its efficacy. A number of in vitro and in vivo studies have already demonstrated DDIs from several cardiovascular ion channel inhibitors with human Pgp and its animal analogs, including verapamil, digoxin, and amiodarone. In this review, Pgp-mediated DDIs and their effects on pharmacokinetics for different categories of cardiovascular ion channel inhibitors are discussed. This information is essential for improving pharmacokinetic predictions of cardiovascular therapeutics, for safer cardiovascular drug administration and for mitigating ADRs emanating from Pgp.

Cooperativity between verapamil and ATP bound to the efflux transporter P-glycoprotein.

The P-glycoprotein (Pgp) transporter plays a central role in drug disposition by effluxing a chemically diverse range of drugs from cells through conformational changes and ATP hydrolysis. A number of drugs are known to activate ATP hydrolysis of Pgp, but coupling between ATP and drug binding is not well understood. The cardiovascular drug verapamil is one of the most widely studied Pgp substrates and therefore, represents an ideal drug to investigate the drug-induced ATPase activation of Pgp. As previously noted, verapamil-induced Pgp-mediated ATP hydrolysis kinetics was biphasic at saturating ATP concentrations. However, at subsaturating ATP concentrations, verapamil-induced ATPase activation kinetics became monophasic. To further understand this switch in kinetic behavior, the Pgp-coupled ATPase activity kinetics was checked with a panel of verapamil and ATP concentrations and fit with the substrate inhibition equation and the kinetic fitting software COPASI. The fits suggested that cooperativity between ATP and verapamil switched between low and high verapamil concentration. Fluorescence spectroscopy of Pgp revealed that cooperativity between verapamil and a non-hydrolyzable ATP analog leads to distinct global conformational changes of Pgp. NMR of Pgp reconstituted in liposomes showed that cooperativity between verapamil and the non-hydrolyzable ATP analog modulate each other's interactions. This information was used to produce a conformationally-gated model of drug-induced activation of Pgp-mediated ATP hydrolysis.

Unravelling the complex drug-drug interactions of the cardiovascular drugs, verapamil and digoxin, with P-glycoprotein.

Drug-drug interactions (DDIs) and associated toxicity from cardiovascular drugs represents a major problem for effective co-administration of cardiovascular therapeutics. A significant amount of drug toxicity from DDIs occurs because of drug interactions and multiple cardiovascular drug binding to the efflux transporter P-glycoprotein (Pgp), which is particularly problematic for cardiovascular drugs because of their relatively low therapeutic indexes. The calcium channel antagonist, verapamil and the cardiac glycoside, digoxin, exhibit DDIs with Pgp through non-competitive inhibition of digoxin transport, which leads to elevated digoxin plasma concentrations and digoxin toxicity. In the present study, verapamil-induced ATPase activation kinetics were biphasic implying at least two verapamil-binding sites on Pgp, whereas monophasic digoxin activation of Pgp-coupled ATPase kinetics suggested a single digoxin-binding site. Using intrinsic protein fluorescence and the saturation transfer double difference (STDD) NMR techniques to probe drug-Pgp interactions, verapamil was found to have little effect on digoxin-Pgp interactions at low concentrations of verapamil, which is consistent with simultaneous binding of the drugs and non-competitive inhibition. Higher concentrations of verapamil caused significant disruption of digoxin-Pgp interactions that suggested overlapping and competing drug-binding sites. These interactions correlated to drug-induced conformational changes deduced from acrylamide quenching of Pgp tryptophan fluorescence. Also, Pgp-coupled ATPase activity kinetics measured with a range of verapamil and digoxin concentrations fit well to a DDI model encompassing non-competitive and competitive inhibition of digoxin by verapamil. The results and previous transport studies were combined into a comprehensive model of verapamil-digoxin DDIs encompassing drug binding, ATP hydrolysis, transport and conformational changes.

The role of cytochrome P450 BM3 phenylalanine-87 and threonine-268 in binding organic hydroperoxides.

BACKGROUND: Cytochrome P450 (P450) BM3, from Bacillus megaterium, catalyzes a wide range of chemical reactions and is routinely used as a model system to study mammalian P450 reactions and structure. METHODS: The metabolism of 2,6-di-tert-butyl-4-hydroperoxy-4-methyl-2,5-cyclohexadienone (BHTOOH) and 2-tert-butyl-4-hydroperoxy-4-methyl-2,5-cyclohexadien-1-one (BMPOOH) was examined with P450 BM3 and with the conserved T268 and F87 residues mutated to investigate their effects on organic hydroperoxide metabolism. To determine the effects of the mutations on the active site volume and architecture, the X-ray crystal structure of the F87A/T268A P450 BM3 heme domain (BMP) was determined and compared to previous structures. To investigate the interactions of the substrates with the F87 and T268 residues, BHTOOH and BMPOOH were docked into the BMP X-ray crystal structures. RESULTS: Lower metabolism of BHTOOH and BMPOOH was observed in the WT P450 BM3 and the T268A P450 BM3 mutant than in the F87A and F87A/T268A P450 BM3 mutants. Large differences were found in the F-G loop regions and active site cavity volumes for the F87A mutated structures. CONCLUSIONS: Analysis of the metabolism, X-ray crystal structures, and molecular docking simulations suggests that P450 BM3 activity toward BHTOOH and BMPOOH is mediated through substrate recognition by T268 and F87, and the active site cavity volume. Based on this information, a simplified representation is presented with the relative orientation of organic hydroperoxides in the P450 BM3 active site. GENERAL SIGNIFICANCE: The metabolism results and structural analysis of this model P450 allowed us to rationalize the structural factors that influence organic hydroperoxide metabolism.