Determinants of Unfolding Cooperativity and Binding Are Decoupled in a DNA Binding Domain.

The relative magnitudes of noncovalent stabilization energies or the coupling free energies in folded proteins are anisotropically distributed, uniquely influencing folding and functional behaviors. In this regard, the fructose repressor (FruR) DBD belonging to the LacR repressor family harbors a three-residue insertion─KQY─between the canonical second and third helices. This sequence insertion promotes a strong Tyr-Tyr stacking interaction that is not observed in related homologues. Combining experiments with simulations, we show that the Tyr-Tyr stacking contributes to a decoupled unfolding due to the localization of a large part of the stabilization energy in this specific structural region. This leads to melting temperatures from different probes spanning nearly 10 K, while concomitantly stabilizing a partially structured intermediate state. Disruption of the aromatic stacking interaction via an alanine mutation promotes a molten-globular state whose native ensemble is replete with non-native interactions while displaying enhanced thermodynamic fluctuations and minimal calorimetric cooperativity. Surprisingly, the molten-globular variant of FruR DBD binds to the operator site on DNA with an affinity similar to that of the wild-type but with altered secondary-structure characteristics in the bound state, underscoring the chaperone-like role of DNA through its large negative electrostatic potential. FruR DBD thus appears to be at the verge of disorder as expected of an entropically destabilizing three-residue insertion but is rescued by the aromatic stacking interaction that distinctly dictates the finer details of stability, cooperativity, and binding.

Structural basis of DNA recognition by BEN domain proteins reveals a role for oligomerization in unmethylated DNA selection by BANP.

The BEN domain is a newly discovered type of DNA-binding domain that exists in a variety of species. There are nine BEN domain-containing proteins in humans, and most have been shown to have chromatin-related functions. NACC1 preferentially binds to CATG motif-containing sequences and functions primarily as a transcriptional coregulator. BANP and BEND3 preferentially bind DNA bearing unmethylated CpG motifs, and they function as CpG island-binding proteins. To date, the DNA recognition mechanism of quite a few of these proteins remains to be determined. In this study, we solved the crystal structures of the BEN domains of NACC1 and BANP in complex with their cognate DNA substrates. We revealed the details of DNA binding by these BEN domain proteins and unexpectedly revealed that oligomerization is required for BANP to select unmethylated CGCG motif-containing DNA substrates. Our study clarifies the controversies surrounding DNA recognition by BANP and demonstrates a new mechanism by which BANP selects unmethylated CpG motifs and functions as a CpG island-binding protein. This understanding will facilitate further exploration of the physiological functions of the BEN domain proteins in the future.

Structural insights into BirA from Haemophilus influenzae, a bifunctional protein as a biotin protein ligase and a transcriptional repressor.

Biotin is an essential coenzyme involved in various metabolic processes across all known organisms, with biotinylation being crucial for the activity of carboxylases. BirA from Haemophilus influenzae is a bifunctional protein that acts as a biotin protein ligase and a transcriptional repressor. This study reveals the crystal structures of Hin BirA in both its apo- and holo-(biotinyl-5'-AMP bound) forms. As a class II BirA, it consists of three domains: N-terminal DNA binding domain, central catalytic domain, and C-terminal SH3-like domain. The structural analysis shows that the biotin-binding loop forms an ordered structure upon biotinyl-5'-AMP binding. This facilitates its interaction with the ligand and promotes protein dimerization. Comparative studies with other BirA homologs from different organisms indicate that the residues responsible for binding biotinyl-5'-AMP are highly conserved. This study also utilized AlphaFold2 to model the potential heterodimeric interaction between Hin BirA and biotin carboxyl carrier protein, thereby providing insights into the structural basis for biotinylation. These findings enhance our understanding of the structural and functional characteristics of Hin BirA, highlighting its potential as a target for novel antibiotics that disrupt the bacterial biotin synthesis pathways.

Structure based computational RNA design towards MafA transcriptional repressor implicated in multiple myeloma.

Multiple myeloma is recognized as the second most common hematological cancer. MafA transcriptional repressor is an established mediator of myelomagenesis. While there are multitude of drugs available for targeting various effectors in multiple myeloma, current literature lacks a candidate RNA based MafA modulator. Thus, using the structure of MafA homodimer-consensus target DNA, a computational effort was implemented to design a novel RNA based chemical modulator against MafA. First, available MafA-consensus DNA structure was employed to generate an RNA library. This library was further subjected to global docking to select the most plausible RNA candidates, preferring to bind DNA binding region of MafA. Following global docking, MD-ready complexes that were prepared via local docking program, were subjected to 500 ns of MD simulations. First, each of these MD simulations were analyzed for relative binding free energy through MM-PBSA method, which pointed towards a strong RNA based MafA binder, RNA1. Second, through a detailed MD analysis, RNA1 was shown to prefer binding to a single monomer of the dimeric DNA binding domain of MafA using higher number of hydrophobic interactions compared with positive control MafA-DNA complex. At the final phase, a principal component analyses was conducted, which led us to identify the actual interaction region of RNA1 and MafA monomer. Overall, to our knowledge, this is the first computational study that presents an RNA molecule capable of potentially targeting MafA protein. Furthermore, limitations of our study together with possible future implications of RNA1 in multiple myeloma were also discussed.

Unveiling the potential of phytochemicals to inhibit nuclear receptor binding SET domain protein 2 for cancer: Pharmacophore screening, molecular docking, ADME properties, and molecular dynamics simulation investigations.

Nuclear receptor binding SET domain protein 2 (NSD2) significantly contributes to the development of cancer, making it a promising target for cancer drug discovery. This research explores natural compounds as potential selective inhibitors for NSD2 in cancer treatment. Employing a comprehensive in silico approach, the study utilized pharmacophore modeling, molecular docking, pharmacokinetic profiling, and molecular dynamics simulations. An e-pharmacophore model-based screening using the first selective and potent ligand bound to NSD2 identified 49,248 natural compounds from the SuperNatural 3.0 database (containing 449,008 molecules) with acceptable alignment with the developed pharmacophore hypotheses. Subsequently, molecular docking was executed to assess the standout compounds which led to the selection of ten candidates that surpassed the reference inhibitor in accordance w the binding affinity expressed as a G score. Ligand-residue interaction analyses of the top three hits (SN0450102, SN0410255, and SN0142336) revealed diverse crucial interactions with the NSD2 active site, including hydrogen bonds, pi-pi stacking, and hydrophobic contacts with key amino acid residues in the NSD2-PWWP1 domain. Pharmacokinetic profiling confirmed the drug-likability for the refined hits, indicating good cellular permeability and minimal blood-brain barrier penetration. Molecular dynamics simulations for 200 nanoseconds affirmed the stability of protein-ligand complexes, with minimal fluctuations in root mean square deviation and root mean square fluctuation analyses. Overall, this study identified promising natural compounds as potential pharmaceutical agents in the treatment of NSD2-associated cancers.

Structural Perspectives on Metal Dependent Roles of Ferric Uptake Regulator (Fur).

Iron is crucial for the metabolism and growth of most prokaryotic cells. The ferric uptake regulator (Fur) protein plays a central role in regulating iron homeostasis and metabolic processes in bacteria. It ensures the proper utilization of iron and the maintenance of cellular functions in response to environmental cues. Fur proteins are composed of an N-terminal DNA-binding domain (DBD) and a C-terminal dimerization domain (DD), typically existing as dimers in solution. Fur proteins have conserved metal-binding sites named S1, S2, and S3. Among them, site S2 serves as a regulatory site, and metal binding at S2 results in conformational changes. Additionally, as a transcriptional regulator, Fur specifically binds to a consensus DNA sequence called the Fur box. To elucidate the structural and functional properties of Fur proteins, various structures of metal- or DNA-bound Fur proteins or apo-Fur proteins have been determined. In this review, we focus on the structural properties of Fur proteins according to their ligand-bound state and the drug development strategies targeting Fur proteins. This information provides valuable insights for drug discovery.

Engineered transcription activator-like effector dimer proteins confer DNA loop-dependent gene repression comparable to Lac repressor.

Natural prokaryotic gene repression systems often exploit DNA looping to increase the local concentration of gene repressor proteins at a regulated promoter via contributions from repressor proteins bound at distant sites. Using principles from the Escherichia coli lac operon we design analogous repression systems based on target sequence-programmable Transcription Activator-Like Effector dimer (TALED) proteins. Such engineered switches may be valuable for synthetic biology and therapeutic applications. Previous TALEDs with inducible non-covalent dimerization showed detectable, but limited, DNA loop-based repression due to the repressor protein dimerization equilibrium. Here, we show robust DNA loop-dependent bacterial promoter repression by covalent TALEDs and verify that DNA looping dramatically enhances promoter repression in E. coli. We characterize repression using a thermodynamic model that quantitates this favorable contribution of DNA looping. This analysis unequivocally and quantitatively demonstrates that optimized TALED proteins can drive loop-dependent promoter repression in E. coli comparable to the natural LacI repressor system. This work elucidates key design principles that set the stage for wide application of TALED-dependent DNA loop-based repression of target genes.

The Curious Case of A31P, a Topology-Switching Mutant of the Repressor of Primer Protein: A Molecular Dynamics Study of Its Folding and Misfolding.

The effect of mutations on protein structures is usually rather localized and minor. Finding a mutation that can single-handedly change the fold and/or topology of a protein structure is a rare exception. The A31P mutant of the homodimeric Repressor of primer (Rop) protein is one such exception: This single mutation ─and as demonstrated by two independent crystal structure determinations─ can convert the canonical (left-handed/all-antiparallel) 4-α-helical bundle of Rop to a new form (right-handed/mixed parallel and antiparallel bundle) displaying a previously unobserved "bisecting U" topology. The main problem with understanding the dramatic effect of this mutation on the folding of Rop is to understand its very existence: Most computational methods appear to agree that the mutation should have had no appreciable effect, with the majority of energy minimization methods and protein structure prediction protocols indicating that this mutation is fully consistent with the native Rop structure, requiring only a local and minor change at the mutation site. Here we use two long (10 µs each) molecular dynamics simulations to compare the stability and dynamics of the native Rop versus a hypothetical structure that is identical with the native Rop but is carrying this single Alanine31 to Proline mutation. Comparative analysis of the two trajectories convincingly shows that, in contrast to the indications from energy minimization ─but in agreement with the experimental data─, this hypothetical native-like A31P structure is unstable, with its turn regions almost completely unfolding, even under the relatively mild 320 K NpT simulations that we have used for this study. We discuss the implication of these findings for the folding of the A31P mutant, especially with respect to the proposed model of a double-funneled energy landscape.

Dissecting the structural and functional consequences of the evolutionary proline-glycine deletion in the wing 1 region of the forkhead domain of human FoxP1.

The human FoxP transcription factors dimerize via three-dimensional domain swapping, a unique feature among the human Fox family, as result of evolutionary sequence adaptations in the forkhead domain. This is the case for the conserved glycine and proline residues in the wing 1 region, which are absent in FoxP proteins but present in most of the Fox family. In this work, we engineered both glycine (G) and proline-glycine (PG) insertion mutants to evaluate the deletion events in FoxP proteins in their dimerization, stability, flexibility, and DNA-binding ability. We show that the PG insertion only increases protein stability, whereas the single glycine insertion decreases the association rate and protein stability and promotes affinity to the DNA ligand.

In-silico study unveils potential phytocompounds in Andrographis paniculata against E6 protein of the high-risk HPV-16 subtype for cervical cancer therapy.

Despite therapeutic advancements, cervical cancer caused by high-risk subtypes of the human papillomavirus (HPV) remains a leading cause of cancer-related deaths among women worldwide. This study aimed to discover potential drug candidates from the Asian medicinal plant Andrographis paniculata, demonstrating efficacy against the E6 protein of high-risk HPV-16 subtype through an in-silico computational approach. The 3D structures of 32 compounds (selected from 42) derived from A. paniculata, exhibiting higher binding affinity, were obtained from the PubChem database. These structures underwent subsequent analysis and screening based on criteria including binding energy, molecular docking, drug likeness and toxicity prediction using computational techniques. Considering the spectrometry, pharmacokinetic properties, docking results, drug likeliness, and toxicological effects, five compounds-stigmasterol, 1H-Indole-3-carboxylic acid, 5-methoxy-, methyl ester (AP7), andrographolide, apigenin and wogonin-were selected as the potential inhibitors against the E6 protein of HPV-16. We also performed 200 ns molecular dynamics simulations of the compounds to analyze their stability and interactions as protein-ligand complexes using imiquimod (CID-57469) as a control. Screened compounds showed favorable characteristics, including stable root mean square deviation values, minimal root mean square fluctuations and consistent radius of gyration values. Intermolecular interactions, such as hydrogen bonds and hydrophobic contacts, were sustained throughout the simulations. The compounds displayed potential affinity, as indicated by negative binding free energy values. Overall, findings of this study suggest that the selected compounds have the potential to act as inhibitors against the E6 protein of HPV-16, offering promising prospects for the treatment and management of CC.

Targeting pericentric non-consecutive motifs for heterochromatin initiation.

Pericentric heterochromatin is a critical component of chromosomes marked by histone H3 K9 (H3K9) methylation1-3. However, what recruits H3K9-specific histone methyltransferases to pericentric regions in vertebrates remains unclear4, as does why pericentric regions in different species share the same H3K9 methylation mark despite lacking highly conserved DNA sequences2,5. Here we show that zinc-finger proteins ZNF512 and ZNF512B specifically localize at pericentric regions through direct DNA binding. Notably, both ZNF512 and ZNF512B are sufficient to initiate de novo heterochromatin formation at ectopically targeted repetitive regions and pericentric regions, as they directly recruit SUV39H1 and SUV39H2 (SUV39H) to catalyse H3K9 methylation. SUV39H2 makes a greater contribution to H3K9 trimethylation, whereas SUV39H1 seems to contribute more to silencing, probably owing to its preferential association with HP1 proteins. ZNF512 and ZNF512B from different species can specifically target pericentric regions of other vertebrates, because the atypical long linker residues between the zinc-fingers of ZNF512 and ZNF512B offer flexibility in recognition of non-consecutively organized three-nucleotide triplets targeted by each zinc-finger. This study addresses two long-standing questions: how constitutive heterochromatin is initiated and how seemingly variable pericentric sequences are targeted by the same set of conserved machinery in vertebrates.

Analysis of Fluorescent Proteins for Observing Single Gene Locus in a Live and Fixed Escherichia coli Cell.

Fluorescent proteins (FPs) are essential tools for advanced microscopy techniques such as super-resolution imaging, single-particle tracking, and quantitative single-molecule counting. Various FPs fused to DNA-binding proteins have been used to observe the subcellular location and movement of specific gene loci in living and fixed bacterial cells. However, quantitative assessments of the properties of FPs for gene locus measurements are still lacking. Here, we assessed various FPs to observe specific gene loci in live and fixed Escherichia coli cells using a fluorescent repressor-operator binding system (FROS), tet operator-Tet repressor proteins (TetR). Tsr-fused FPs were used to assess the intensity and photostability of various FPs (five red FPs: mCherry2, FusionRed, mRFP, mCrimson3, and dKatushka; and seven yellow FPs: SYFP2, Venus, mCitrine, YPet, mClover3, mTopaz, and EYFP) at the single-molecule level in living cells. These FPs were then used for gene locus measurements using FROS. Our results indicate that TetR-mCrimson3 (red) and TetR-EYFP (yellow) had better properties for visualizing gene loci than the other TetR-FPs. Furthermore, fixation procedures affected the clustering of diffusing TetR-FPs and altered the locations of the TetR-FP foci. Fixation with formaldehyde consistently disrupted proper DNA locus observations using TetR-FPs. Notably, the foci measured using TetR-mCrimson3 remained close to their original positions in live cells after glyoxal fixation. This in vivo study provides a cell-imaging guide for the use of FPs for gene-locus observation in E. coli and a scheme for evaluating the use of FPs for other cell-imaging purposes.

Oligomerization of protein arginine methyltransferase 1 and its effect on methyltransferase activity and substrate specificity.

Proper protein arginine methylation by protein arginine methyltransferase 1 (PRMT1) is critical for maintaining cellular health, while dysregulation is often associated with disease. How the activity of PRMT1 is regulated is therefore paramount, but is not clearly understood. Several studies have observed higher order oligomeric species of PRMT1, but it is unclear if these exist at physiological concentrations and there is confusion in the literature about how oligomerization affects activity. We therefore sought to determine which oligomeric species of PRMT1 are physiologically relevant, and quantitatively correlate activity with specific oligomer forms. Through quantitative western blotting, we determined that concentrations of PRMT1 available in a variety of human cell lines are in the sub-micromolar to low micromolar range. Isothermal spectral shift binding data were modeled to a monomer/dimer/tetramer equilibrium with an EC50 for tetramer dissociation of ~20 nM. A combination of sedimentation velocity and Native polyacrylamide gel electrophoresis experiments directly confirmed that the major oligomeric species of PRMT1 at physiological concentrations would be dimers and tetramers. Surprisingly, the methyltransferase activity of a dimeric PRMT1 variant is similar to wild type, tetrameric PRMT1 with some purified substrates, but dimer and tetramer forms of PRMT1 show differences in catalytic efficiencies and substrate specificity for other substrates. Our results define an oligomerization paradigm for PRMT1, show that the biophysical characteristics of PRMT1 are poised to support a monomer/dimer/tetramer equilibrium in vivo, and suggest that the oligomeric state of PRMT1 could be used to regulate substrate specificity.

Dynamic BTB-domain filaments promote clustering of ZBTB proteins.

The formation of dynamic protein filaments contributes to various biological functions by clustering individual molecules together and enhancing their binding to ligands. We report such a propensity for the BTB domains of certain proteins from the ZBTB family, a large eukaryotic transcription factor family implicated in differentiation and cancer. Working with Xenopus laevis and human proteins, we solved the crystal structures of filaments formed by dimers of the BTB domains of ZBTB8A and ZBTB18 and demonstrated concentration-dependent higher-order assemblies of these dimers in solution. In cells, the BTB-domain filamentation supports clustering of full-length human ZBTB8A and ZBTB18 into dynamic nuclear foci and contributes to the ZBTB18-mediated repression of a reporter gene. The BTB domains of up to 21 human ZBTB family members and two related proteins, NACC1 and NACC2, are predicted to behave in a similar manner. Our results suggest that filamentation is a more common feature of transcription factors than is currently appreciated.

Unveiling critical structural features for effective HDAC8 inhibition: a comprehensive study using quantitative read-across structure-activity relationship (q-RASAR) and pharmacophore modeling.

Histone deacetylases constitute a group of enzymes that participate in several biological processes. Notably, inhibiting HDAC8 has become a therapeutic strategy for various diseases. The current inhibitors for HDAC8 lack selectivity and target multiple HDACs. Consequently, there is a growing recognition of the need for selective HDAC8 inhibitors to enhance the effectiveness of therapeutic interventions. In our current study, we have utilized a multi-faceted approach, including Quantitative Structure-Activity Relationship (QSAR) combined with Quantitative Read-Across Structure-Activity Relationship (q-RASAR) modeling, pharmacophore mapping, molecular docking, and molecular dynamics (MD) simulations. The developed q-RASAR model has a high statistical significance and predictive ability (Q2F1:0.778, Q2F2:0.775). The contributions of important descriptors are discussed in detail to gain insight into the crucial structural features in HDAC8 inhibition. The best pharmacophore hypothesis exhibits a high regression coefficient (0.969) and a low root mean square deviation (0.944), highlighting the importance of correctly orienting hydrogen bond acceptor (HBA), ring aromatic (RA), and zinc-binding group (ZBG) features in designing potent HDAC8 inhibitors. To confirm the results of q-RASAR and pharmacophore mapping, molecular docking analysis of the five potent compounds (44, 54, 82, 102, and 118) was performed to gain further insights into these structural features crucial for interaction with the HDAC8 enzyme. Lastly, MD simulation studies of the most active compound (54, mapped correctly with the pharmacophore hypothesis) and the least active compound (34, mapped poorly with the pharmacophore hypothesis) were carried out to validate the observations of the studies above. This study not only refines our understanding of essential structural features for HDAC8 inhibition but also provides a robust framework for the rational design of novel selective HDAC8 inhibitors which may offer insights to medicinal chemists and researchers engaged in the development of HDAC8-targeted therapeutics.

Discovery of NSD2 non-histone substrates and design of a super-substrate.

The human protein lysine methyltransferase NSD2 catalyzes dimethylation at H3K36. It has very important roles in development and disease but many mechanistic features and its full spectrum of substrate proteins are unclear. Using peptide SPOT array methylation assays, we investigate the substrate sequence specificity of NSD2 and discover strong readout of residues between G33 (-3) and P38 (+2) on H3K36. Unexpectedly, we observe that amino acid residues different from natural ones in H3K36 are preferred at some positions. Combining four preferred residues led to the development of a super-substrate which is methylated much faster by NSD2 at peptide and protein level. Molecular dynamics simulations demonstrate that this activity increase is caused by distinct hyperactive conformations of the enzyme-peptide complex. To investigate the substrate spectrum of NSD2, we conducted a proteome wide search for nuclear proteins matching the specificity profile and discovered 22 peptide substrates of NSD2. In protein methylation studies, we identify K1033 of ATRX and K819 of FANCM as NSD2 methylation sites and also demonstrate their methylation in human cells. Both these proteins have important roles in DNA repair strengthening the connection of NSD2 and H3K36 methylation to DNA repair.

The Importance of Structural Water in HDAC8 for Correct Binding Pose Applied for Drug Design of Anticancer Molecules.

AIMS: Validating the docking procedure and maintaining the structural water molecules at HDAC8 catalytic site. BACKGROUND: Molecular docking simulations play a significant role in Computer-Aided Drug Design, contributing to the development of new molecules. To ensure the reliability of these simulations, a validation process called "self-docking or re-docking" is employed, focusing on the binding mode of a ligand co-crystallized with the protein of interest. OBJECTIVE: In this study, several molecular docking studies were conducted using five X-ray structures of HDAC8-ligand complexes from the PDB. METHODS: Ligands initially complexed with HDAC8 were removed and re-docked onto the free protein, revealing a poor reproduction of the expected binding mode. In response to this, we observed that most HDAC8-ligand complexes contained one to two water molecules in the catalytic site, which were crucial for maintaining the cocrystallized ligand. RESULTS: These water molecules enhance the binding mode of the co-crystallized ligand by stabilizing the proteinligand complex through hydrogen bond interactions between ligand and water molecules. Notably, these interactions are lost if water molecules are removed, as is often done in classical docking methodologies. Considering this, molecular docking simulations were repeated, both with and without one or two conserved water molecules near Zn+2 in the catalytic cavity. Simulations indicated that replicating the native binding pose of co-crystallized ligands on free HDAC8 without these water molecules was challenging, showing greater coordinate displacements (RMSD) compared to those including conserved water molecules from crystals. CONCLUSION: The study highlighted the importance of conserved water molecules within the active site, as their presence significantly influenced the successful reproduction of the ligands' native binding modes. The results suggest an optimal molecular docking procedure for validating methods suitable for filtering new HDAC8 inhibitors for future experimental assays.

Rapid Determination of Kinetic Constants for Slow-Binding Inhibitors and Inactivators of Human Histone Deacetylase 8.

The kinetics and mechanism of drug binding to its target are critical to pharmacological efficacy. A high throughput (HTS) screen often results in hundreds of hits, of which usually only simple IC50 values are determined during reconfirmation. However, kinetic parameters such as residence time for reversible inhibitors and the kinact/KI ratio, which is the critical measure for evaluating covalent inactivators, are early predictive measures to assess the chances of success of the hits in the clinic. Using the promising cancer target human histone deacetylase 8 as an example, we present a robust method that calculates concentration-dependent apparent rate constants for the inhibition or inactivation of HDAC8 from dose-response curves recorded after different pre-incubation times. With these data, hit compounds can be classified according to their mechanism of action, and the relevant kinetic parameters can be calculated in a highly parallel fashion. HDAC8 inhibitors with known modes of action were correctly assigned to their mechanism, and the binding mechanisms of some hits from an internal HDAC8 screening campaign were newly determined. The oxonitriles SVE04 and SVE27 were classified as fast reversible HDAC8 inhibitors with moderate time-constant IC50 values of 4.2 and 2.6 µM, respectively. The hit compound TJ-19-24 and SAH03 behave like slow two-step inactivators or reversible inhibitors, with a very low reverse isomerization rate.

A parameterized two-domain thermodynamic model explains diverse mutational effects on protein allostery.

New experimental findings continue to challenge our understanding of protein allostery. Recent deep mutational scanning study showed that allosteric hotspots in the tetracycline repressor (TetR) and its homologous transcriptional factors are broadly distributed rather than spanning well-defined structural pathways as often assumed. Moreover, hotspot mutation-induced allostery loss was rescued by distributed additional mutations in a degenerate fashion. Here, we develop a two-domain thermodynamic model for TetR, which readily rationalizes these intriguing observations. The model accurately captures the in vivo activities of various mutants with changes in physically transparent parameters, allowing the data-based quantification of mutational effects using statistical inference. Our analysis reveals the intrinsic connection of intra- and inter-domain properties for allosteric regulation and illustrate epistatic interactions that are consistent with structural features of the protein. The insights gained from this study into the nature of two-domain allostery are expected to have broader implications for other multi-domain allosteric proteins.

Synthesis and Assay of Vibrio Quorum Sensing Inhibitors.

Bacteria detect local population numbers using quorum sensing, a method of cell-cell communication broadly utilized to control bacterial behaviors. In Vibrio species, the master quorum sensing regulators LuxR/HapR control hundreds of quorum sensing genes, many of which influence virulence, metabolism, motility, and more. Thiophenesulfonamides are potent inhibitors of LuxR/HapR that bind the ligand pocket in these transcription factors and block downstream quorum sensing gene expression. This class of compounds served as the basis for the development of a set of simple, robust, and educational procedures for college students to assimilate their chemistry and biology skills using a CURE model: course-based undergraduate research experience. Optimized protocols are described that comprise three learning stages in an iterative and multi-disciplinary platform to engage students in a year-long CURE: (1) design and synthesize new small molecule inhibitors based on the thiophenesulfonamide core, (2) use structural modeling to predict binding affinity to the target, and (3) assay the compounds for efficacy in microbiological assays against specific Vibrio LuxR/HapR proteins. The described reporter assay performed in E. coli successfully predicts the efficacy of the compounds against target proteins in the native Vibrio species.