KRAS , the most frequently mutated oncogene in human cancer, produces two isoforms, KRAS4a and KRAS4b, through alternative splicing. These isoforms differ in exon 4, which encodes the final 15 residues of the G-domain and hypervariable regions (HVRs), vital for trafficking and membrane localization. While KRAS4b has been extensively studied, KRAS4a has been largely overlooked. Our multidisciplinary study compared the structural and functional characteristics of KRAS4a and KRAS4b, revealing distinct structural properties and thermal stability. Position 151 influences KRAS4a’s thermal stability, while position 153 affects binding to RAF1 CRD protein. Nuclear magnetic resonance analysis identified localized structural differences near sequence variations and provided a solution-state conformational ensemble. Notably, KRAS4a exhibits substantial transcript abundance in bile ducts, liver, and stomach, with transcript levels approaching KRAS4b in the colon and rectum. Functional disparities were observed in full-length KRAS variants, highlighting the impact of HVR variations on interaction with trafficking proteins and downstream effectors like RAF and PI3K within cells.

Fragment-based screening by ligand-observed 1D NMR and binding interface mapping by protein-observed 2D NMR are popular methods used in drug discovery. These methods allow researchers to detect compound binding over a wide range of affinities and offer a simultaneous assessment of solubility, purity, and chemical formula accuracy of the target compounds and the 15N-labeled protein when examined by 1D and 2D NMR, respectively. These methods can be applied for screening fragment binding to the active (GMPPNP-bound) and inactive (GDP-bound) states of oncogenic KRAS mutants.

The KRAS gene plays a pivotal role in numerous cancers by encoding a GTPase that upon association with the plasma membrane activates the MAPK pathway, promoting cellular proliferation. In our study, we investigated small molecules that disrupt KRAS's membrane interaction, hypothesizing that such disruption could in turn inhibit mutant RAS signaling. Native mass spectrometry screening of KRAS-FMe identified compounds with a preference for interacting with the hypervariable region (HVR), and surface plasmon resonance (SPR) further refined our selection to graveoline as a compound exhibiting preferential HVR binding. Subsequent nuclear magnetic resonance (NMR) analysis showed that graveoline's interaction with KRAS depends on C-terminal O-methylation. Moreover, our findings revealed multiple interaction sites, suggesting weak engagement with the KRAS G domain. Using nanodiscs as a membrane mimetic, further characterization through NMR and Förster resonance energy transfer (FRET) studies demonstrated graveoline's ability to perturb KRAS membrane interaction in a biochemical setting. Our biophysical approach sheds light on the intricate molecular mechanisms underlying KRAS-ligand interactions, providing valuable insights into understanding the KRAS-associated pathophysiology. These findings contribute to the translational aspect of our study, offering potential avenues for further research targeting KRAS membrane association with the potential to lead to a new class of RAS therapeutics.

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRASG13D, which suggests an exchange between two conformational states in solution. Methyl and 31P NMR spectra of active KRASG13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major Pγ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRASG13D and KRASG13D-RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRASG13D, showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

KRAS is the most frequently mutated RAS protein in cancer patients, and it is estimated that about 20% of the cancer patients in the United States carried mutant RAS proteins. To accelerate therapeutic development, structures and dynamics of RAS proteins had been extensively studied by various biophysical techniques for decades. Although 31P NMR studies revealed population equilibrium of the two major states in the active GMPPNP-bound form, more complex conformational dynamics in RAS proteins and oncogenic mutants subtly modulate the interactions with their downstream effectors. We established a set of customized NMR relaxation dispersion techniques to efficiently and systematically examine the ms-μs conformational dynamics of RAS proteins. This method allowed us to observe varying synchronized motions that connect the effector and allosteric lobes in KRAS. We demonstrated the role of conformational dynamics of KRAS in controlling its interaction with the Ras-binding domain of the downstream effector RAF1, the first kinase in the MAPK pathway. This allows one to explain, as well as to predict, the altered binding affinities of various KRAS mutants, which was neither previously reported nor apparent from the structural perspective.

The bacterial pathogen Vibrio cholerae use a type III secretion system to inject effector proteins into a host cell. Recently, a putative Toxic GTPase Activating Protein (ToxGAP) called Vibrio outer protein E (VopE) was identified as a T3SS substrate and virulence factor that affected host mitochondrial dynamics and immune response. However, biophysical and structural characterization has been absent. Here, we describe solution NMR structure of the putative GTPase-activating protein (GAP) domain (73-204) of VopE. Using size exclusion chromatography coupled with small-angle x-ray scattering and residual dipolar coupling data, we restrained the MD process to efficiently determine the overall fold and improve the quality of the output calculated structures. Comparing the structure of VopE with other ToxGAP's revealed a similar overall fold with several features unique to VopE. Specifically, the "Bulge 1," α1 helix, and noteworthy "backside linker" elements on the N-terminus are dissimilar to the other ToxGAP's. By using NMR relaxation dispersion experiments, we demonstrate that these regions undergo motions on a > 6 s-1 timescale. Based on the disposition of these mobile regions relative to the putative catalytic arginine residue, we hypothesize that the protein may undergo structural changes to bind cognate GTPases.

The ATPase Family, AAA domain-containing protein 2 (ATAD2) bromodomain (BRD) has a canonical bromodomain structure consisting of four α-helices. ATAD2 functions as a co-activator of the androgen and estrogen receptors as well as the MYC and E2F transcription factors. ATAD2 also functions during DNA replication, recognizing newly synthesized histones. In addition, ATAD2 is shown to be up regulated in multiple forms of cancer including breast, lung, gastric, endometrial, renal, and prostate. Furthermore, up-regulation of ATAD2 is strongly correlated with poor prognosis in many types of cancer, making the ATAD2 bromodomain an innovative target for cancer therapeutics. In this study, we describe the recognition of histone acetyllysine modifications by the ATAD2 bromodomain. Residue-specific information on the complex formed between the histone tail and the ATAD2 bromodomain, obtained through nuclear magnetic resonance spectroscopy (NMR) and X-ray crystallography, illustrates key residues lining the binding pocket, which are involved in coordination of di-acetylated histone tails. Analytical ultracentrifugation, NMR relaxation data, and isothermal titration calorimetry further confirm the monomeric state of the functionally active ATAD2 bromodomain in complex with di-acetylated histone ligands. Overall, we describe histone tail recognition by ATAD2 BRD and illustrate that one acetyllysine group is primarily engaged by the conserved asparagine (N1064), the “RVF” shelf residues, and the flexible ZA loop. Coordination of a second acetyllysine group also occurs within the same binding pocket, but is essentially governed by unique hydrophobic and electrostatic interactions making the di-acetyllysine histone coordination more specific than previously presumed.

Arabidopsis thaliana contains 36 known members of a class of small proteins known as rapid alkalinization factors (RALF). We succeeded in overexpressing one of these proteins, RALF8, which we found to inhibit plant root growth, and labeled it uniformly with 13C and 15N for NMR analysis. RALF8, which contains four cysteine residues, proved to be a challenging target for structure determination because much of the protein is discorded. We used non-uniform sampling (NUS) for 3D NMR collection followed by spectral reconstruction. In determining the solution structure of RALF8, we used the latest version of the Integrative NMR platform, which includes the software tools described below. Backbone and sidechain signals were automatically detected and assigned by tools in NMRFAM-SPARKY. Tools in the I-PINE web server yielded structural information including the presence of two disulfide bridges. The disulfide pairing was determined by hetNOE analysis and confirmed by data from an optimized 1H,13C-HSQC-NOESY-1H,13C-HSQC pulse program. The 3D structure was determined automatically by AUDANA algorithm in the PONDEROSA-C/S suite, which automatically assigns NOE cross peaks and generates distance and angle restraints. After refining the restraints by the visual validation module in PONDEROSA-Analyzer, the Final Step with Explicit Water Refinement option was used in calculating 100 conformers and selecting the 20 conformers with lowest energy, which were used to represent the solution structure. The root mean square deviation (RMSD) between the backbone heavy atoms of all residues in the superimposed structures was 12.8 Å, consistent with overall dynamic disorder. Two loops between the disulfides (C21–29 and C41–47) were well-ordered with 0.6 and 0.7 Å RMSDs respectively. All assigned chemical shifts and atomic coordinates are available under PDB 6NU4 and BMRB 30565.

Bromodomain-containing proteins are often part of chromatin-modifying complexes, and their activity can lead to altered expression of genes that drive cancer, inflammation and neurological disorders in humans. Bromodomain-PHD finger protein 1 (BRPF1) is part of the MOZ (monocytic leukemic zinc-finger protein) HAT (histone acetyltransferase) complex, which is associated with chromosomal translocations known to contribute to the development of acute myeloid leukemia (AML). BRPF1 contains a unique combination of chromatin reader domains including two plant homeodomain (PHD) fingers separated by a zinc knuckle (PZP domain), a bromodomain, and a proline-tryptophan-tryptophan-proline (PWWP) domain. BRPF1 is known to recruit the MOZ HAT complex to chromatin by recognizing acetylated lysine residues on the N-terminal histone tail region through its bromodomain. However, histone proteins can contain several acetylation modifications on their N-terminus, and it is unknown how additional marks influence bromodomain recruitment to chromatin. Here, we identify the BRPF1 bromodomain as a selective reader of di-acetyllysine modifications on histone H4. We used ITC assays to characterize the binding of di-acetylated histone ligands to the BRPF1 bromodomain and found that the domain binds preferentially to histone peptides H4K5acK8ac and H4K5acK12ac. Analytical ultracentrifugation (AUC) experiments revealed that the monomeric state of the BRPF1 bromodomain coordinates di-acetylated histone ligands. NMR chemical shift perturbation studies, along with binding and mutational analyses, revealed non-canonical regions of the bromodomain-binding pocket that are important for histone tail recognition. Together, our findings provide critical information on how the combinatorial action of post-translational modifications can modulate BRPF1 bromodomain binding and specificity.

Bromodomains exhibit preferences for specific patterns of post-translational modifications on core and variant histone proteins. We examined the ligand specificity of the ATAD2B bromodomain and compared it to its closely related paralogue in ATAD2. We show that the ATAD2B bromodomain recognizes mono- and diacetyllysine modifications on histones H4 and H2A. A structure-function approach was used to identify key residues in the acetyllysine-binding pocket that dictate the molecular recognition process, and we examined the binding of an ATAD2 bromodomain inhibitor by ATAD2B. Our analysis demonstrated that critical contacts required for bromodomain inhibitor coordination are conserved between the ATAD2/B bromodomains, with many residues playing a dual role in acetyllysine recognition. We further characterized an alternative splice variant of ATAD2B that results in a loss of function. Our results outline the structural and functional features of the ATAD2B bromodomain and identify a novel mechanism regulating the interaction of the ATAD2B protein with chromatin.

Streptococcus pneumoniae is an opportunistic human pathogen that utilizes the competence regulon, a quorum-sensing circuitry, to acquire antibiotic resistance genes and initiate its attack on the human host. Interception of the competence regulon can therefore be utilized to study S. pneumoniae cell-cell communication and behavioral changes, as well as attenuate S. pneumoniae infectivity. Herein we report the design and synthesis of cyclic dominant negative competence-stimulating peptide (dnCSP) analogs capable of intercepting the competence regulon in both S. pneumoniae specificity groups with activities at the low nanomolar range. Structural analysis of lead analogs provided important insights as to the molecular mechanism that drives CSP receptor binding and revealed that the pan-group cyclic CSPs exhibit a chimeric hydrophobic patch conformation that resembles the hydrophobic patches required for both ComD1 and ComD2 binding. Moreover, the lead cyclic dnCSP, CSP1-E1A-cyc(Dap6E10), was found to possess superior pharmacological properties, including improved resistance to enzymatic degradation, while remaining nontoxic. Lastly, CSP1-E1A-cyc(Dap6E10) was capable of attenuating mouse mortality during acute pneumonia caused by both group 1 and group 2 S. pneumoniae strains. This cyclic pan-group dnCSP is therefore a promising drug lead scaffold against S. pneumoniae infections that could be administered individually or utilized in combination therapy to augment the effects of current antimicrobial agents.

The J-domain protein Zuotin is a multi-domain eukaryotic Hsp70 co-chaperone. Though it is primarily ribosome-associated, positioned at the exit of the 60S subunit tunnel where it promotes folding of nascent polypeptide chains, Zuotin also has off-ribosome functions. Domains of Zuotin needed for 60S association and interaction with Hsp70 are conserved in eukaryotes. However, whether the 4-helix bundle (4HB) domain is conserved remains an open question. We undertook evolutionary and structural approaches to clarify this issue. We found that the 4HB segment of human Zuotin also forms a bundle of 4 helices. The positive charge of Helix I, which in Saccharomyces cerevisiae is responsible for interaction with the 40S subunit, is particularly conserved. However, the C-termini of fungal and human 4HBs are not similar. In fungi the C-terminal segment forms a plug that folds back into the bundle; in S. cerevisiae it plays an important role in bundle stability and, off the ribosome, in transcriptional activation. In human, C-terminal helix IV of the 4HB is extended, protruding from the bundle. This extension serves as a linker to the regulatory SANT domains, which are present in animals, plants and protists, but not fungi. Further analysis of Zuotin sequences revealed that the plug likely arose as a result of genomic rearrangement upon SANT domain loss early in the fungal lineage. In the lineage leading to S. cerevisiae, the 4HB was subjected to positive selection with the plug becoming increasingly hydrophobic. Eventually, these hydrophobic plug residues were coopted for a novel regulatory function—activation of a recently emerged transcription factor, Pdr1. Our data suggests that Zuotin evolved off-ribosome functions twice—once involving SANT domains, then later in fungi, after SANT domain loss, by coopting the hydrophobic plug. Zuotin serves as an example of complex intertwining of molecular chaperone function and cell regulation.

We report the solution-phase structures of native signal peptides and related analogs capable of either strongly agonizing or antagonizing the AgrC quorum sensing (QS) receptor in the emerging pathogen Staphylococcus epidermidis. Chronic S. epidermidis infections are often recalcitrant to traditional therapies due to antibiotic resistance and formation of robust biofilms. The accessory gene regulator (agr) QS system plays an important role in biofilm formation in this opportunistic pathogen, and the binding of an autoinducing peptide (AIP) signal to its cognate transmembrane receptor (AgrC) is responsible for controlling agr. Small molecules or peptides capable of modulating this binding event are of significant interest as probes to investigate both the agr system and QS as a potential antivirulence target. We used NMR spectroscopy to characterize the structures of the three native S. epidermidis AIP signals and five non-native analogs with distinct activity profiles in the AgrC-I receptor from S. epidermidis. These studies revealed a suite of structural motifs critical for ligand activity. Interestingly, a unique beta-turn was present in the macrocycles of the two most potent AgrC-I modulators, in both an agonist and an antagonist, which was distinct from the macrocycle conformation in the less-potent AgrC-I modulators and in the native AIP-I itself. This previously unknown beta-turn provides a structural rationale for these ligands' respective biological activity profiles. Development of analogs to reinforce the beta-turn resulted in our first antagonist with subnanomolar potency in AgrC-I, while analogs designed to contain a disrupted beta-turn were dramatically less potent relative to their parent compounds. Collectively, these studies provide new insights into the AIP:AgrC interactions crucial for QS activation in S. epidermidis and advance the understanding of QS at the molecular level.

We report the recombinant preparation from Escherichia coli cells of samples of two closely related, small, secreted cysteine-rich plant peptides: rapid alkalinization factor 1 (RALF1) and rapid alkalinization factor 8 (RALF8). Purified samples of the native sequence of RALF8 exhibited well-resolved nuclear magnetic resonance (NMR) spectra and also biological activity through interaction with a plant receptor kinase, cytoplasmic calcium mobilization, and in vivo root growth suppression. By contrast, RALF1 could only be isolated from inclusion bodies as a construct containing an N-terminal His-tag; its poorly resolved NMR spectrum was indicative of aggregation. We prepared samples of the RALF8 peptide labeled with 15 N and 13 C for NMR analysis and obtained near complete 1 H, 13 C, and 15 N NMR assignments; determined the disulfide pairing of its four cysteine residues; and examined its solution structure. RALF8 is mostly disordered except for the two loops spanned by each of its two disulfide bridges.

A presumed RNA cloverleaf (5'CL), located at the 5'-most end of the noncoding region of the enterovirus genome, is the primary established site for initiation of genomic replication. Stem-loop B (SLB) and stem-loop D (SLD), the two largest stem-loops within the 5'CL, serve as recognition sites for protein interactions that are essential for replication. Here we present the solution structure of rhinovirus serotype 14 5'CL using a combination of nuclear magnetic resonance spectroscopy and small-angle X-ray scattering. In the absence of magnesium, the structure adopts an open, somewhat extended conformation. In the presence of magnesium, the structure compacts, bringing SLB and SLD into close contact, a geometry that creates an extensive accessible major groove surface, and permits interaction between the proteins that target each stem-loop.

Computational molecular dynamics, energy minimization, and modeling of molecular interactions are widely used in studies involving natural products, metabolites, and drugs. Manually directed computational steps commonly utilize an evolving collection of experimental and computational data, to which new data sources are added or modified as needed. Several software packages capable of incorporating sources of data are available, but the process remains error prone owing to the complexities of preparing and maintaining a consistent set of input files and the proper post-processing of derived data. We have devised a methodology and implemented it using an extensible software pipeline called RUNER (for Robust and Unique Nomenclature for Enhanced Reproducibility) that creates a robust and standardized computational process. The pipeline combines a web service and a graphical user interface (GUI) to enable seamless modifications and verified maintenance of atom force field parameters. The GUI provides an implementation for the widely used molecular modeling software package Xplor-NIH. We describe the RUNER software and demonstrate the rationale for the pipeline through examples of structural studies of small molecules and natural products. The software, pipeline, force field parameters, and file verification data for more than 4,100 compounds (including FDA-approved drugs and natural products) are freely accessible from [ http://runer.nmrfam.wisc.edu ]. Author Summary We describe an automated and verifiable computational pipeline for calculating the force field parameters of small molecules. The pipeline integrates several software tools and guarantees reproducibility of the parameters by utilizing a standard nomenclature across multiple computational steps and by maintaining file verification identifiers. We demonstrate the application of this pipeline to (a) processing of more than 4,100 compounds in high-throughput mode, and (b) structural studies of natural products. The graphical user interface (GUI) associated with the pipeline facilitates the manually tedious steps of force field parameters adjustments and supports visualization of the process.

Streptococcus pneumoniae is an important pathogen that utilizes quorum sensing (QS) to regulate genetic transformation, virulence, and biofilm formation. The competence-stimulating peptide (CSP) is a 17-amino acid signal peptide that is used by S. pneumoniae to trigger QS. S. pneumoniae strains can be divided into two main specificity groups based on the CSP signal they produce (CSP1 or CSP2) and their compatible receptors (ComD1 or ComD2, respectively). Modulation of QS in S. pneumoniae can be achieved by targeting the CSP:ComD interaction using synthetic CSP analogues. However, to rationally design CSP-based QS modulators with enhanced activities, an in-depth understanding of the structural features that are required for receptor binding is needed. Herein, we report a comprehensive in-solution three-dimensional structural characterization of eight CSP1 and CSP2 analogues with varied biological activities using nuclear magnetic resonance spectroscopy. Analysis of these structures revealed two distinct hydrophobic patches required for effective ComD1 and ComD2 binding.

A polyether antibiotic, ecteinamycin (1), was isolated from a marine Actinomadura sp., cultivated from the ascidian Ecteinascidia turbinata. C-13 enrichment, high resolution NMR spectroscopy, and molecular modeling enabled elucidation of the structure of 1, which was validated on the basis of comparisons with its recently reported crystal structure. Importantly, ecteinamycin demonstrated potent activity against the toxigenic strain of Clostridium difficile NAP1/B1/027 (MIC = 59 ng/mu L), as well as other toxigenic and nontoxigenic C. difficile isolates both in vitro and in vivo. Additionally, chemical genomics studies using Escherichia coli barcoded deletion mutants led to the identification of sensitive mutants such as trkA and kdpD involved in potassium cation transport and homeostasis supporting a mechanistic proposal that ecteinamycin acts as an ionophore antibiotic. This is the first antibacterial agent whose mechanism of action has been studied using E. coli chemical genomics. On the basis of these data, we propose ecteinamycin as an ionophore antibiotic that causes C. difficile detoxification and cell death via potassium transport dysregulation.

By binding to a multitude of polypeptide substrates, Hsp70-based molecular chaperone systems perform a range of cellular functions. All J-protein co-chaperones play the essential role, via action of their J-domains, of stimulating the ATPase activity of Hsp70, thereby stabilizing its interaction with substrate. In addition, J-proteins drive the functional diversity of Hsp70 chaperone systems through action of regions outside their J-domains. Targeting to specific locations within a cellular compartment and binding of specific substrates for delivery to Hsp70 have been identified as modes of J-protein specialization. To better understand J-protein specialization, we concentrated on Saccharomyces cerevisiae SIS1, which encodes an essential J-protein of the cytosol/nucleus. We selected suppressors that allowed cells lacking SIS1 to form colonies. Substitutions changing single residues in Ydj1, a J-protein, which, like Sis1, partners with Hsp70 Ssa1, were isolated. These gain-of-function substitutions were located at the end of the J-domain, suggesting that suppression was connected to interaction with its partner Hsp70, rather than substrate binding or subcellular localization. Reasoning that, if YDJ1 suppressors affect Ssa1 function, substitutions in Hsp70 itself might also be able to overcome the cellular requirement for Sis1, we carried out a selection for SSA1 suppressor mutations. Suppressing substitutions were isolated that altered sites in Ssa1 affecting the cycle of substrate interaction. Together, our results point to a third, additional means by which J-proteins can drive Hsp70’s ability to function in a wide range of cellular processes—modulating the Hsp70-substrate interaction cycle.