Investigating Polypharmacology through Targeting Known Human Neutrophil Elastase Inhibitors to Proteinase 3.

Using a combination of multisite λ-dynamics (MSλD) together with in vitro IC50 assays, we evaluated the polypharmacological potential of a scaffold currently in clinical trials for inhibition of human neutrophil elastase (HNE), targeting cardiopulmonary disease, for efficacious inhibition of Proteinase 3 (PR3), a related neutrophil serine proteinase. The affinities we observe suggest that the dihydropyrimidinone scaffold can serve as a suitable starting point for the establishment of polypharmacologically targeting both enzymes and enhancing the potential for treatments addressing diseases like chronic obstructive pulmonary disease.

Norstictic Acid Is a Selective Allosteric Transcriptional Regulator.

Inhibitors of transcriptional protein-protein interactions (PPIs) have high value both as tools and for therapeutic applications. The PPI network mediated by the transcriptional coactivator Med25, for example, regulates stress-response and motility pathways, and dysregulation of the PPI networks contributes to oncogenesis and metastasis. The canonical transcription factor binding sites within Med25 are large (∼900 Å2) and have little topology, and thus, they do not present an array of attractive small-molecule binding sites for inhibitor discovery. Here we demonstrate that the depsidone natural product norstictic acid functions through an alternative binding site to block Med25-transcriptional activator PPIs in vitro and in cell culture. Norstictic acid targets a binding site comprising a highly dynamic loop flanking one canonical binding surface, and in doing so, it both orthosterically and allosterically alters Med25-driven transcription in a patient-derived model of triple-negative breast cancer. These results highlight the potential of Med25 as a therapeutic target as well as the inhibitor discovery opportunities presented by structurally dynamic loops within otherwise challenging proteins.

Modeling pH-Dependent NMR Chemical Shift Perturbations in Peptides.

Modeling the pH dependence of protein and peptide chemical shifts outside the range of physiological values (6.5-7) is key to understanding structure-function relationships of these systems. These capabilities are largely not available in current chemical shift prediction software. In this study, we utilize a combination of molecular dynamics and quantum mechanics to investigate the through-space and through-bond contributions of protonation-dependent chemical shift perturbations (CSPs) in model tripeptides. By altering the protonation state of the titratable group in the tripeptides, we observe a notable difference in the conformational ensembles and attendantly compute significant CSPs for all nuclei near the site of protonation. We thus demonstrate the ability to recapitulate experimental pH-dependent CSPs with good agreement (R = 0.85, 0.99, and 0.98 for 13C, 15N, and 1H, respectively). Broadly, we provide the groundwork for incorporating pH effects into empirical and semiempirical chemical shift predictors.

Investigating the Effect of Two-Point Surface Attachment on Enzyme Stability and Activity.

Immobilization on solid supports provides an effective way to improve enzyme stability and simplify downstream processing for biotechnological applications, which has been widely used in research and in applications. However, surface immobilization may disrupt enzyme structure due to interactions between the enzyme and the supporting substrate, leading to a loss of the enzyme catalytic efficiency and stability. Here, we use a model enzyme, nitroreductase (NfsB), to demonstrate that engineered variants with two strategically positioned surface-tethering sites exhibit improved enzyme stability when covalently immobilized onto a surface. Tethering sites were designed based on molecular dynamics (MD) simulations, and enzyme variants containing cysteinyl residues at these positions were expressed, purified, and immobilized on maleimide-terminated self-assembled monolayer (SAM) surfaces. Sum frequency generation (SFG) vibrational spectroscopy and attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy were used to deduce the NfsB enzyme orientations, which were found to be consistent with those predicted from the MD simulations. Thermal stability analyses demonstrated that NfsB variants immobilized through two tethering sites exhibited generally improved thermal stability compared with enzymes tethered at only one position. For example, NfsB enzyme chemically immobilized via positions 423 and 111 exhibits at least 60% stability increase compared to chemically immobilized NfsB mutant via a single site. This research develops a generally applicable and systematic approach using a combination of simulation and experimental methods to rationally select protein immobilization sites for the optimization of surface-immobilized enzyme activity and stability.

Conservation of coactivator engagement mechanism enables small-molecule allosteric modulators.

Transcriptional coactivators are a molecular recognition marvel because a single domain within these proteins, the activator binding domain or ABD, interacts with multiple compositionally diverse transcriptional activators. Also remarkable is the structural diversity among ABDs, which range from conformationally dynamic helical motifs to those with a stable core such as a ß-barrel. A significant objective is to define conserved properties of ABDs that allow them to interact with disparate activator sequences. The ABD of the coactivator Med25 (activator interaction domain or AcID) is unique in that it contains secondary structural elements that are on both ends of the spectrum: helices and loops that display significant conformational mobility and a seven-stranded ß-barrel core that is structurally rigid. Using biophysical approaches, we build a mechanistic model of how AcID forms binary and ternary complexes with three distinct activators; despite its static core, Med25 forms short-lived, conformationally mobile, and structurally distinct complexes with each of the cognate partners. Further, ternary complex formation is facilitated by allosteric communication between binding surfaces on opposing faces of the ß-barrel. The model emerging suggests that the conformational shifts and cooperative binding is mediated by a flexible substructure comprised of two dynamic helices and flanking loops, indicating a conserved mechanistic model of activator engagement across ABDs. Targeting a region of this substructure with a small-molecule covalent cochaperone modulates ternary complex formation. Our data support a general strategy for the identification of allosteric small-molecule modulators of ABDs, which are key targets for mechanistic studies as well as therapeutic applications.

Effect of immobilization site on the orientation and activity of surface-tethered enzymes.

Tethering peptides and proteins to abiotic surfaces has the potential to create biomolecule-functionalized surfaces with useful properties. Commonly used methods of immobilization lack control over the orientation in which biological molecules are covalently or physically bound to the surface, leading to sub-optimal materials. Here we use an engineered beta-galactosidase that can be chemically immobilized on a surface with a well-defined orientation through unique surface-accessible cysteine residues. A combined study using sum frequency generation (SFG) vibrational spectroscopy and coarse grained molecular dynamics (MD) simulations was performed to determine the effects of enzyme immobilization site and abiotic surface chemistry on enzyme surface orientation, surface coverage, and catalytic activity. Two beta-galactosidase variants that were immobilized through cysteine introduced at positions 227 and 308 were studied. In both cases, when the abiotic surface was made more hydrophilic, the enzyme surface coverage decreased, but the activity increased. MD simulations indicated that this is due to the weakened interactions between the immobilized enzyme and the more hydrophilic surface. These studies provide improved understanding of how enzyme-surface interactions can be optimized to maximize the catalytic activity of surface tethered enzymes.

Mechanism of Vps4 hexamer function revealed by cryo-EM.

Vps4 is a member of AAA+ ATPase (adenosine triphosphatase associated with diverse cellular activities) that operates as an oligomer to disassemble ESCRT-III (endosomal sorting complex required for transport III) filaments, thereby catalyzing the final step in multiple ESCRT-dependent membrane remodeling events. We used electron cryo-microscopy to visualize oligomers of a hydrolysis-deficient Vps4 (vacuolar protein sorting-associated protein 4) mutant in the presence of adenosine 5'-triphosphate (ATP). We show that Vps4 subunits assemble into an asymmetric hexameric ring following an approximate helical path that sequentially stacks substrate-binding loops along the central pore. The hexamer is observed to adopt an open or closed ring configuration facilitated by major conformational changes in a single subunit. The structural transition of the mobile Vps4 subunit results in the repositioning of its substrate-binding loop from the top to the bottom of the central pore, with an associated translation of 33 Å. These structures, along with mutant-doping experiments and functional assays, provide evidence for a sequential and processive ATP hydrolysis mechanism by which Vps4 hexamers disassemble ESCRT-III filaments.

Molecular Mechanism for Isoform-Selective Inhibition of Acyl Protein Thioesterases 1 and 2 (APT1 and APT2).

Post-translational S-palmitoylation directs the trafficking and membrane localization of hundreds of cellular proteins, often involving a coordinated palmitoylation cycle that requires both protein acyl transferases (PATs) and acyl protein thioesterases (APTs) to actively redistribute S-palmitoylated proteins toward different cellular membrane compartments. This process is necessary for the trafficking and oncogenic signaling of S-palmitoylated Ras isoforms, and potentially many peripheral membrane proteins. The depalmitoylating enzymes APT1 and APT2 are separately conserved in all vertebrates, suggesting unique functional roles for each enzyme. The recent discovery of the APT isoform-selective inhibitors ML348 and ML349 has opened new possibilities to probe the function of each enzyme, yet it remains unclear how each inhibitor achieves orthogonal inhibition. Herein, we report the high-resolution structure of human APT2 in complex with ML349 (1.64 Å), as well as the complementary structure of human APT1 bound to ML348 (1.55 Å). Although the overall peptide backbone structures are nearly identical, each inhibitor adopts a distinct conformation within each active site. In APT1, the trifluoromethyl group of ML348 is positioned above the catalytic triad, but in APT2, the sulfonyl group of ML349 forms hydrogen bonds with active site resident waters to indirectly engage the catalytic triad and oxyanion hole. Reciprocal mutagenesis and activity profiling revealed several differing residues surrounding the active site that serve as critical gatekeepers for isoform accessibility and dynamics. Structural and biochemical analysis suggests the inhibitors occupy a putative acyl-binding region, establishing the mechanism for isoform-specific inhibition, hydrolysis of acyl substrates, and structural orthogonality important for future probe development.

A Synthetic Loop Replacement Peptide That Blocks Canonical NF-κB Signaling.

Aberrant canonical NF-κB signaling is implicated in diseases from autoimmune disorders to cancer. A major therapeutic challenge is the need for selective inhibition of the canonical pathway without impacting the many non-canonical NF-κB functions. Here we show that a selective peptide-based inhibitor of canonical NF-κB signaling, in which a hydrogen bond in the NBD peptide is synthetically replaced by a non-labile bond, shows an about 10-fold increased potency relative to the original inhibitor. Not only is this molecule, NBD2, a powerful tool for dissection of canonical NF-κB signaling in disease models and healthy tissues, the success of the synthetic loop replacement suggests that the general strategy could be useful for discovering modulators of the many protein-protein interactions mediated by such structures.

Regulation of calreticulin-major histocompatibility complex (MHC) class I interactions by ATP.

The MHC class I peptide loading complex (PLC) facilitates the assembly of MHC class I molecules with peptides, but factors that regulate the stability and dynamics of the assembly complex are largely uncharacterized. Based on initial findings that ATP, in addition to MHC class I-specific peptide, is able to induce MHC class I dissociation from the PLC, we investigated the interaction of ATP with the chaperone calreticulin, an endoplasmic reticulum (ER) luminal, calcium-binding component of the PLC that is known to bind ATP. We combined computational and experimental measurements to identify residues within the globular domain of calreticulin, in proximity to the high-affinity calcium-binding site, that are important for high-affinity ATP binding and for ATPase activity. High-affinity calcium binding by calreticulin is required for optimal nucleotide binding, but both ATP and ADP destabilize enthalpy-driven high-affinity calcium binding to calreticulin. ATP also selectively destabilizes the interaction of calreticulin with cellular substrates, including MHC class I molecules. Calreticulin mutants that affect ATP or high-affinity calcium binding display prolonged associations with monoglucosylated forms of cellular MHC class I, delaying MHC class I dissociation from the PLC and their transit through the secretory pathway. These studies reveal central roles for ATP and calcium binding as regulators of calreticulin-substrate interactions and as key determinants of PLC dynamics.

Molecular-Level Insights into Orientation-Dependent Changes in the Thermal Stability of Enzymes Covalently Immobilized on Surfaces.

Surface-immobilized enzymes are important for a wide range of technological applications, including industrial catalysis, drug delivery, medical diagnosis, and biosensors; however, our understanding of how enzymes and proteins interact with abiological surfaces on the molecular level remains extremely limited. We have compared the structure, activity, and thermal stability of two variants of a ß-galactosidase attached to a chemically well-defined maleimide-terminated self-assembled monolayer surface through a unique cysteinyl residue. In one case the enzyme is attached through an α helix and in the other case through an adjacent loop. Both enzymes exhibit similar specific activities and adopt similar orientations with respect to the surface normal, as determined by sum-frequency generation and attenuated total reflectance FT-IR spectroscopies. Surprisingly, however, the loop-tethered enzyme exhibits a thermal stability 10 °C lower than the helix-tethered enzyme and 13 °C lower than the enzyme in free solution. Using coarse-grain models, molecular dynamics simulations of the thermal unfolding of the surface-tethered enzymes were able to reproduce these differences in stability. Thus, revealing that tethering through the more flexible loop position provides more opportunity for surface residues on the protein to interact with the surface and undergo surface-induced unfolding. These observations point to the importance of the location of the attachment point in determining the performance of surface-supported biocatalysts and suggest strategies for optimizing their activity and thermal stability through molecular simulations.

Membrane environment modulates the pKa values of transmembrane helices.

In this work, we apply the recently developed constant pH molecular dynamics technique to study protonation equilibria of titratable side chains in the context of simple transmembrane (TM) helices and explore the effect of pH on their configurations in membrane bilayers. We observe that, despite a significant shift toward neutral states, considerable population of different side chains stay in the charged state that give rise to pKa values around 9.6 for Asp and Glu and 4.5 to 6 for His and Lys side chains, respectively. These charged states are highly stabilized by favorable interactions between head groups, water molecules, and the charged side chains that are facilitated by substantial changes in the configuration of the peptides. The pH dependent configurations and the measured pKa values are in good agreement with relatively recent solid state NMR measurements. Our results presented here demonstrate that all-atom constant pH molecular dynamics can be applied to membrane proteins and peptides to obtain reliable pKa values and pH dependent behavior for these systems.

Residue-level resolution of alphavirus envelope protein interactions in pH-dependent fusion.

Alphavirus envelope proteins, organized as trimers of E2-E1 heterodimers on the surface of the pathogenic alphavirus, mediate the low pH-triggered fusion of viral and endosomal membranes in human cells. The lack of specific treatment for alphaviral infections motivates our exploration of potential antiviral approaches by inhibiting one or more fusion steps in the common endocytic viral entry pathway. In this work, we performed constant pH molecular dynamics based on an atomic model of the alphavirus envelope with icosahedral symmetry. We have identified pH-sensitive residues that cause the largest shifts in thermodynamic driving forces under neutral and acidic pH conditions for various fusion steps. A series of conserved interdomain His residues is identified to be responsible for the pH-dependent conformational changes in the fusion process, and ligand binding sites in their vicinity are anticipated to be potential drug targets aimed at inhibiting viral infections.

Environmental effect on surface immobilized biological molecules.

Our recent sum frequency generation (SFG) vibrational spectroscopic experiment ( J. Phys. Chem. B 2014 , 118 , 2904 - 2912 ) showed that immobilized antimicrobial peptide cecropin P1 (cCP1) on a self-assembled monolayer (SAM) surface via N-terminus exhibited significantly different conformational and/or orientational behaviors when exposed to pure water vs a 50% (v/v) 2,2,2-trifluoroethanol (TFE)/water mixture. Meanwhile, our recent molecular dynamics (MD) simulations ( J. Phys. Chem. B 2014 , 118 , 5670 - 5680 ) further revealed that the immobilized cCP1 via N-terminus in pure water largely adopts an overall bent structure lying down on the SAM surface, consistent with the SFG observation. Here, MD simulations were performed on the immobilized cCP1 on a SAM surface via N-terminus while in contact with a 50% (v/v) TFE/water mixture to further investigate the effects of environment (water vs TFE/water mixture) on the interfacial structure and orientation of immobilized peptide. The simulation results demonstrated that the immobilized cCP1 on the SAM surface via the N-terminus with two different starting states with different orientations and conformations, when exposed to a 50% (v/v) TFE/water mixture, was eventually able to maintain a linear α-helical structure, standing upright on the SAM surface. Taken with the corresponding SFG observation, our simulation results indicate that the conformational behavior of the immobilized peptide is mediated by the local hydrophobic environments resulting from the TFE aggregation around the peptide. Such knowledge can be used to regulate the surface conformation and functionality of immobilized peptides via changing surrounding chemical environments (e.g., TFE cosolvent), which is important for the microbial detection and killing based on surface-immobilized antimicrobial peptides.

Prepaying the entropic cost for allosteric regulation in KIX.

The kinase-inducible domain interacting (KIX) domain of the CREB binding protein (CBP) is capable of simultaneously binding two intrinsically disordered transcription factors, such as the mixed-lineage leukemia (MLL) and c-Myb peptides, at isolated interaction sites. In vitro, the affinity for binding c-Myb is approximately doubled when KIX is in complex with MLL, which suggests a positive cooperative binding mechanism, and the affinity for MLL is also slightly increased when KIX is first bound by c-Myb. Expanding the scope of recent NMR and computational studies, we explore the allosteric mechanism at a detailed molecular level that directly connects the microscopic structural dynamics to the macroscopic shift in binding affinities. To this end, we have performed molecular dynamics simulations of free KIX, KIX-c-Myb, MLL-KIX, and MLL-KIX-c-Myb using a topology-based Go-like model. Our results capture an increase in affinity for the peptide in the allosteric site when KIX is prebound by a complementary effector and both peptides follow an effector-independent folding-and-binding mechanism. More importantly, we discover that MLL binding lowers the entropic cost for c-Myb binding, and vice versa, by stabilizing the L12-G2 loop and the C-terminal region of the α3 helix on KIX. This work demonstrates the importance of entropy in allosteric signaling between promiscuous molecular recognition sites and can inform the rational design of small molecule stabilizers to target important regions of conformationally dynamic proteins.

Molecular structures of C- and N-terminus cysteine modified cecropin P1 chemically immobilized onto maleimide-terminated self-assembled monolayers investigated by molecular dynamics simulation.

Biosensors using peptides or proteins chemically immobilized on surfaces have many advantages such as better sensitivity, improved stability, and longer shelf life compared to those prepared using physically adsorbed biomolecules. Chemical immobilization can better control the interfacial conformation and orientation of peptides and proteins, leading to better activity of these biomolecules. In this research, molecular dynamics (MD) simulations were employed to systematically investigate the structure and dynamics of surface-tethered antimicrobial peptide cecropin P1 (CP1) modified with a cysteine residue at the C- (CP1c) or N-terminus (cCP1). Such CP1c and cCP1 molecules were chemically immobilized onto a silane-EG4-maleimide self-assembled monolayer (SAM) surface by forming a thio-ether bond between the cysteine group in CP1c or cCP1 and the surface maleimide group. The simulation results showed that the immobilized cCP1 (via the N-terminus) tends to bend and gradually lie down onto the SAM surface, due to the large structural fluctuation of the C-terminus induced by unfavorable interactions between the hydrophobic C-terminal residues and water. Differently, the tethered CP1c (via the C-terminus) more or less stands up on the surface, only tilting slightly even after 60 ns. The simulation results can be well correlated to the recent experimental results obtained from sum frequency generation (SFG) vibrational spectroscopic study. The current simulation data provide more atomic level details on how the hydrophobicity difference in the C-terminus and N-terminus of the amphiphilic peptide can lead to different structures of the same peptide tethered to the surface via different termini. This knowledge can be used to rationally design chemically immobilized peptides to achieve desired structure and functionality.

Different interfacial behaviors of peptides chemically immobilized on surfaces with different linker lengths and via different termini.

Molecular structures such as conformation and orientation are crucial in determining the activity of peptides immobilized to solid supports. In this study, sum frequency generation (SFG) vibrational spectroscopy was applied to investigate such structures of peptides immobilized on self-assembled monolayers (SAMs). Here cysteine-modified antimicrobial peptide cecropin P1 (CP1) was chemically immobilized onto SAM with a maleimide terminal group. Two important characteristics, length of the poly(ethylene glycol) (PEG) segment in the SAM and location of the cysteine residue in the peptide, were examined using SFG spectroscopy to determine the effect of each on surface immobilization as well as peptide secondary structure and its orientation in the immobilized state. Results have shown that while each length of PEG chain studied promotes chemical immobilization of the target peptide and prevents nonspecific adsorption, CP1 immobilized on long-chain (PEG2k) maleimide SAMs shows random coil structure in water, whereas CP1 demonstrates α-helical structure when immobilized on short-chain (with four ethylene glycol units - (EG4)) maleimide SAMs. Placement of the cysteine residue at the C-terminus promotes the formation of α-helical structure of CP1 with a single orientation when tethered to EG4 maleimide SAM surfaces. In contrast, immobilization via the N-terminal cysteine of CP1 results in a random coil or lying-down helical structure. The bacteria capturing/killing capability was tested, showing that the surface-immobilized CP1 molecules via C- and N- terminal cysteine exhibit only slight difference, even though they have different secondary structures and orientations.

Conformational dynamics of a regulator of G-protein signaling protein reveals a mechanism of allosteric inhibition by a small molecule.

Regulators of G protein signaling (RGS) proteins are key players in regulating signaling via G protein-coupled receptors. RGS proteins directly bind to the Gα-subunits of activated heterotrimeric G-proteins, and accelerate the rate of GTP hydrolysis, thereby rapidly deactivating G-proteins. Using atomistic simulations and NMR spectroscopy, we have studied in molecular detail the mechanism of action of CCG-50014, a potent small molecule inhibitor of RGS4 that covalently binds to cysteine residues on RGS4. We apply temperature-accelerated molecular dynamics (TAMD) to carry out enhanced conformational sampling of apo RGS4 structures, and consistently find that the α5-α6 helix pair of RGS4 can spontaneously span open-like conformations, allowing binding of CCG-50014 to the buried side-chain of Cys95. Both NMR experiments and MD simulations reveal chemical shift perturbations in residues in the vicinity of inhibitor binding site as well as in the RGS4-Gα binding interface. Consistent with a loss of G-protein binding, GAP activity, and allosteric mechanism of action of CCG-50014, our simulations of the RGS4-Gα complex in the presence of inhibitor suggest a relatively unstable protein-protein interaction. These results have potential implications for understanding how the conformational dynamics among RGS proteins may play a key role in the sensitivity of inhibitors.

Hexameric helicase deconstructed: interplay of conformational changes and substrate coupling.

Hexameric helicases are molecular motor proteins that utilize energy obtained from ATP hydrolysis to translocate along and/or unwind nucleic acids. In this study, we investigate the dynamic behavior of the Simian Virus 40 hexameric helicase bound to DNA by performing molecular dynamics simulations employing a coarse-grained model. Our results elucidate the two most important molecular features of the helicase motion. First, the attractive interactions between the DNA-binding domain of the helicase and the DNA backbone are essential for the helicase to exhibit a unidirectional motion along the DNA strand. Second, the sequence of ATP binding at multiple binding pockets affects the helicase motion. Specifically, concerted ATP binding does not generate a unidirectional motion of the helicase. It is only when the binding of ATP occurs sequentially from one pocket to the next that the helicase moves unidirectionally along the DNA. Interestingly, in the reverse order of sequential ATP binding, the helicase also moves unidirectionally but in the opposite direction. These observations suggest that in nature ATP molecules must distinguish between different available ATP binding pockets of the hexameric helicase in order to function efficiently. To this end, simulations reveal that the binding of ATP in one pocket induces an opening of the next ATP-binding pocket and such an asymmetric deformation may coordinate the sequential ATP binding in a unidirectional manner. Overall, these findings may provide clues toward understanding the mechanism of substrate translocation in other motor proteins.

Community-wide assessment of GPCR structure modelling and ligand docking: GPCR Dock 2008.

Recent breakthroughs in the determination of the crystal structures of G protein-coupled receptors (GPCRs) have provided new opportunities for structure-based drug design strategies targeting this protein family. With the aim of evaluating the current status of GPCR structure prediction and ligand docking, a community-wide, blind prediction assessment - GPCR Dock 2008 - was conducted in coordination with the publication of the crystal structure of the human adenosine A(2A) receptor bound to the ligand ZM241385. Twenty-nine groups submitted 206 structural models before the release of the experimental structure, which were evaluated for the accuracy of the ligand binding mode and the overall receptor model compared with the crystal structure. This analysis highlights important aspects for success and future development, such as accurate modelling of structurally divergent regions and use of additional biochemical insight such as disulphide bridges in the extracellular loops.