Thermodynamic mechanism for the evasion of antibody neutralization in flaviviruses

Conformational changes in Protein Kinase A along its activation cycle are rooted in the folding energetics of cyclic-nucleotide binding domains

Cyclic-nucleotide binding (CNB) domains are structurally and evolutionarily conserved signaling modules that regulate proteins with diverse folds and functions. Despite a wealth of structural information, the mechanisms by which CNB domains couple cyclic-nucleotide binding to conformational changes involved in signal transduction remain unknown. Here we combined single-molecule and computational approaches to investigate the conformation and folding energetics of the two CNB domains of the regulatory subunit of protein kinase A (PKA). We found that the CNB domains exhibit different conformational and folding signatures in the apo state, when bound to cAMP, or when bound to the PKA catalytic subunit, underscoring their ability to adapt to different binding partners. Moreover, we show while the two CNB domains have near-identical structures, their thermodynamic coupling signatures are divergent, leading to distinct cAMP responses and differential mutational effects. Specifically, we demonstrate the mutation W260A exerts local and allosteric effects that impact multiple steps of the PKA activation cycle. Taken together, these results highlight the complex interplay between folding energetics, conformational dynamics, and thermodynamic signatures that underlies structurally conserved signaling modules in response to ligand binding and mutational effects.

Signal Transmission in Escherichia coli Cyclic AMP Receptor Protein for Survival in Extreme Acidic Conditions

During the life cycle of enteric bacterium Escherichia coli, it encounters a wide spectrum of pH changes. The asymmetric dimer of the cAMP receptor protein, CRP, plays a key role in regulating the expressions of genes and the survival of E. coli. To elucidate the pH effects on the mechanism of signal transmission, we present a combination of results derived from ITC, crystallography, and computation. CRP responds to a pH change by inducing a differential effect on the affinity for the binding events to the two cAMP molecules, ensuing in a reversible conversion between positive and negative cooperativity at high and low pH, respectively. The structures of four crystals at pH ranging from 7.8 to 6.5 show that CRP responds by inducing a differential effect on the structures of the two subunits, particularly in the DNA binding domain. Employing the COREX/BEST algorithm, computational analysis shows the change in the stability of residues at each pH. The change in residue stability alters the connectivity between residues including those in cAMP and DNA binding sites. Consequently, the differential impact on the topology of the connectivity surface among residues in adjacent subunits is the main reason for differential change in affinity; that is, the pH-induced differential change in residue stability is the biothermodynamic basis for the change in allosteric behavior. Furthermore, the structural asymmetry of this homodimer amplifies the differential impact of any perturbations. Hence, these results demonstrate that the combination of these approaches can provide insights into the underlying mechanism of an apparent complex allostery signal and transmission in CRP.

Differential modulation of energy landscapes of cyclic AMP receptor protein (CRP) as a regulatory mechanism for class II CRP-dependent promoters

The Escherichia coli cAMP receptor protein, CRP, is a homodimeric global transcription activator that employs multiple mechanisms to modulate the expression of hundreds of genes. These mechanisms require different interfacial interactions among CRP, RNA, and DNA of varying sequences. The involvement of such a multiplicity of interfaces requires a tight control to ensure the desired phenotype. CRP-dependent promoters can be grouped into three classes. For decades scientists in the field have been puzzled over the differences in mechanisms between class I and II promoters. Using a new crystal structure, IR spectroscopy, and computational analysis, we defined the energy landscapes of WT and 14 mutated CRPs to determine how a homodimeric protein can distinguish nonpalindromic DNA sequences and facilitate communication between residues located in three different activation regions (AR) in CRP that are ∼30 Å apart. We showed that each mutation imparts differential effects on stability among the subunits and domains in CRP. Consequently, the energetic landscapes of subunits and domains are different, and CRP is asymmetric. Hence, the same mutation can exert different effects on ARs in class I or II promoters. The effect of a mutation is transmitted through a network by long-distance communication not necessarily relying on physical contacts between adjacent residues. The mechanism is simply the sum of the consequences of modulating the synchrony of dynamic motions of residues at a distance, leading to differential effects on ARs in different subunits. The computational analysis is applicable to any system and potentially with predictive capability.

Small molecule grp94 inhibitors block dengue and Zika virus replication

Two major flaviviruses, dengue virus (DENV) and Zika virus (ZIKV), cause severe health and economic burdens worldwide. Recently, genome-wide screenings have uncovered the importance of regulators of the Hrd1 ubiquitin ligase-mediated endoplasmic reticulum (ER)-associated degradation (ERAD) pathway for flavivirus replication in host cells. Here we report the identification of the compound Bardoxolone methyl (CDDO-me) as a potent inhibitor of the Hrd1 ubiquitin ligase-mediated ERAD, which possesses a broad-spectrum activity against both DENV and ZIKV. Cellular thermal shift assay (CETSA) suggested that CDDO-me binds to grp94, a key component of the Hrd1 pathway, at a low nanomolar concentration, whereas interaction was not detected with its paralog Hsp90. CDDO-me and the grp94 inhibitor PU-WS13 substantially suppressed DENV2 replication and the cytopathic effects caused by DENV and ZIKV infection. The antiviral activities of both compounds were demonstrated for all four DENV serotypes and four ZIKV strains in multiple human cell lines. This study defines grp94 as a crucial host factor for flavivirus replication and identified CDDO-me as a potent small molecule inhibitor of flavivirus infection. Inhibition of grp94 may contribute to the antiviral activity of CDDO-me. Further investigation of grp94 inhibitors may lead to a new class of broad-spectrum anti-flaviviral medications.

Networked Communication between Polymerase and Exonuclease Active Sites in Human Mitochondrial DNA Polymerase

High fidelity human mitochondrial DNA polymerase (Pol γ) contains two active sites, a DNA polymerization site (pol) and a 3'-5' exonuclease site (exo) for proofreading. Although separated by 35 Å, coordination between the pol and exo sites is crucial to high fidelity replication. The biophysical mechanisms for this coordination are not completely understood. To understand the communication between the two active sites, we used a statistical-mechanical model of the protein ensemble to calculate the energetic landscape and local stability. We compared a series of structures of Pol γ, complexed with primer/template DNA, and either a nucleotide substrate or a series of nucleotide analogues, which are differentially incorporated and excised by pol and exo activity. Despite the nucleotide or its analogues being bound in the pol, Pol γ residue stability varied across the protein, particularly in the exo domain. This suggests that substrate presence in the pol can be "sensed" in the exo domain. Consistent with this hypothesis, in silico mutations made in one active site mutually perturbed the energetics of the other. To identify specific regions of the polymerase that contributed to this communication, we constructed an allosteric network connectivity map that further demonstrates specific pol-exo cooperativity. Thus, a cooperative network underlies energetic connectivity. We propose that Pol γ and other dual-function polymerases exploit an energetic coupling network that facilitates domain-domain communication to enhance discrimination between correct and incorrect nucleotides.

A comprehensive ensemble model for comparing the allosteric effect of ordered and disordered proteins

Intrinsically disordered proteins/regions (IDPs/IDRs) are prevalent in allosteric regulation. It was previously thought that intrinsic disorder is favorable for maximizing the allosteric coupling. Here, we propose a comprehensive ensemble model to compare the roles of both order-order transition and disorder-order transition in allosteric effect. It is revealed that the MWC pathway (order-order transition) has a higher probability than the EAM pathway (disorder-order transition) in allostery, suggesting a complicated role of IDPs/IDRs in regulatory proteins. In addition, an analytic formula for the maximal allosteric coupling response is obtained, which shows that too stable or too unstable state is unfavorable to endow allostery, and is thus helpful for rational design of allosteric drugs.

Allosteric effect is an important regulation mechanism in biological processes, where the binding of a ligand at one site of a protein influences the function of a distant site. Conventionally, allostery was thought to originate from structural transition. However, in recent years, intrinsically disordered proteins (IDPs) were found to be widely involved in allosteric regulation in despite of their lack of ordered structure under physiological condition. It is still a mystery why IDPs are prevalent in allosteric proteins and how they differ from ordered proteins in allostery. Here, we propose a comprehensive ensemble model which includes both ordered and disordered states of a two-domain protein, and investigate the role of various state combinations in allosteric effect. By sampling the parameter space, we conclude that disordered proteins are less competitive than ordered proteins in performing allostery from a thermodynamic point of view. The prevalence of IDPs in allosteric regulation is likely determined by all their advantage, but not only by their capacity in endowing allostery.

Long-Range Communication Network in the Type 1B Bone Morphogenetic Protein Receptor

Protein-protein interactions are recognized as a fundamental phenomenon that is intimately associated with biological functions and thus are ideal targets for developing modulators for regulating biological functions. A challenge is to identify a site that is situated away from but functionally connected to the protein-protein interface. We employed bone morphogenetic proteins (BMPs) and their receptors as a model system to develop a strategy for identifying such a network of communication. Accordingly, using computational analyses with the COREX/BEST algorithm, we uncovered an overall pattern connecting various regions of BMPR-1B ectodomain, including the four conserved residues in the protein-protein interface. In preparation for testing the long-range effects of mutations of distal residues for future studies, we examined the extent of measurable perturbation of the four conserved residues by determination of the conformation and relative affinities of these BMPR-1B mutants for ligands BMP-2, -6, and -7 and GDF-5. Results suggest no significant structural changes in the receptor but do suggest that the four residues play different roles in defining ligand affinity and both intra- and intermolecular interactions play a role in defining ligand affinity. Thus, these results established two primary but necessary goals: (1) the baseline knowledge of perturbation of conserved interfacial residues for future reference and (2) the ability of the computational approach to identify the distal residues connecting to the interfacial residues. The data presented here provide the foundation for future experiments to identify the effects of distal residues that affect the specificity and affinity of BMP recognition. Protein-protein interactions are integral reactions in essentially all biological activities such as gene regulation and age-related development. Often, diseases are consequences of the alteration of these intermacromolecular interactions, which are thus recognized as a legitimate target for developing modulators for regulating biological functions. One approach is to design ligands that bind to the protein-protein interface. Another is to identify an allosteric site, an advantage of which is bypassing the potential challenge in competing for high-affinity interfacial interactions or a specific interface in a superassembly of multiple macromolecules. However, a challenge of this approach is identifying a site that is situated away from but functionally connected to the protein-protein interface.

Broad spectrum assessment of the epitope fluctuation--Immunogenicity hypothesis

Prediction of immunogenicity is a substantial barrier in vaccine design. Here, a molecular dynamics approach to assessing the immunogenicity of nanoparticles based on structure is presented. Molecular properties of epitopes on nonenveloped viral particles are quantified via a set of metrics. One such metric, epitope fluctuation (and implied flexibility), is shown to be inversely correlated with immunogenicity for each of a broad spectrum of nonenveloped viruses. The molecular metrics and experimentally determined immunogenicities for these viruses are archived in the open-source vaccine computer-aided design database. Results indicate the promise of computer-aided vaccine design to bring greater efficiency to traditional lab-based vaccine discovery approaches.