Protein flexibility directs DNA recognition by the papillomavirus E2 proteins

Molecular cloning, biophysical and in silico studies of Human papillomavirus 33 E2 DNA binding domain

Human papillomavirus 33, a high-risk HPV strain, is mainly responsible for HPV infection and cervical cancer in Asian countries. The E2 protein of HPV 33 is a DNA-binding protein that plays a crucial role in viral replication and transcription. We have cloned, overexpressed, and purified the DNA binding domain of the E2 protein. Size exclusion chromatography results suggested that the protein exists in a homodimeric state in the native form. Circular dichroism data showed that the protein has a higher content of β-sheet. The melting temperature obtained from differential scanning calorimetry is 52.59 °C, and the protein is stable at pH 8 and is in a dimeric form at basic pH. The protein is monomeric or unfolded at a very low pH. Chemical denaturation studies suggested that the protein denatured and dissociated simultaneously. The DNA binding activity of the protein was also confirmed and it showed binding affinity in the order of 106 M-1. The protein structure was modeled using homology modeling and other bioinformatic tools. The virtual screening and molecular dynamic simulation studies were performed to find compounds that can act as potent inhibitors against E2 DBD. This study expands the understanding of the conserved structural and binding properties of HPV33 E2 DBD and provides the first report on the characterization of the viral protein.Communicated by Ramaswamy H. Sarma.

The intrinsically disordered transcriptional activation domain of CIITA is functionally tuneable by single substitutions: An exception or a new paradigm?

During protein evolution, some amino acid substitutions modulate protein function ("tuneability"). In most proteins, the tuneable range is wide and can be sampled by a set of protein variants that each contains multiple amino acid substitutions. In other proteins, the full tuneable range can be accessed by a set of variants that each contains a single substitution. Indeed, in some globular proteins, the full tuneable range can be accessed by the set of site-saturating substitutions at an individual "rheostat" position. However, in proteins with intrinsically disordered regions (IDRs), most functional studies-which would also detect tuneability-used multiple substitutions or small deletions. In disordered transcriptional activation domains (ADs), studies with multiple substitutions led to the "acidic exposure" model, which does not anticipate the existence of rheostat positions. In the few studies that did assess effects of single substitutions on AD function, results were mixed: the ADs of two full-length transcription factors did not show tuneability, whereas a fragment of a third AD was tuneable by single substitutions. In this study, we tested tuneability in the AD of full-length human class II transactivator (CIITA). Sequence analyses and experiments showed that CIITA's AD is an IDR. Functional assays of singly-substituted AD variants showed that CIITA's function was highly tuneable, with outcomes not predicted by the acidic exposure model. Four tested positions showed rheostat behavior for transcriptional activation. Thus, tuneability of different IDRs can vary widely. Future studies are needed to illuminate the biophysical features that govern whether an IDR is tuneable by single substitutions.

Construction and Application of a High-Throughput In Vivo Screening Platform for the Evolution of Nitrile Metabolism-Related Enzymes Based on a Desensitized Repressive Biosensor

Transcription factor (TF)-based biosensors are expected to serve as powerful tools for the high-throughput screening of biocatalytic systems; however, most of them respond to ligands in a narrow concentration range, which limits their application. In this study, we constructed a heterogenous niacin biosensor using the repressive TF BsNadR and its target promoters from Bacillus subtilis. The fine-tunable output of the niacin biosensor was expanded to a wide range of niacin concentrations (0-50 mM) through desensitization engineering, which was suitable for the accurate identification of differences in enzyme activity. Structural mechanism analysis indicated that weakening the affinity of BsNadR with the ligand niacin and with DNA alters its regulatory properties. Based on the desensitized niacin biosensor, a high-throughput in vivo screening platform was developed for evolving nitrile metabolism-related enzymes. The evolved nitrilase, amidase, and nitrile hydratase with 6.6-, 2.1-, and 21.3-fold improvements in activity were achieved, respectively. In addition, these mutants also exhibited elevated activity toward other cognate substrates, indicating the broad applicability of the screening platform. This study not only provided a universal high-throughput screening platform for different nitrile metabolism-related enzymes but also demonstrated the advantages of repressive biosensors and the vital role of desensitization engineering of the TF in the development of high-throughput screening platforms for enzymes.

A Systematic Structure-Activity Study of a New Type of Small Peptidic Transfection Vector Reveals the Importance of a Special Oxo-Anion-Binding Motif for Gene Delivery

We discovered a new class of artificial peptidic transfection vectors based on an artificial anion-binding motif, the guanidiniocarbonylpyrrole (GCP) cation. This new type of vector is surprisingly smaller than traditional systems, and our previous work suggested that the GCP group was important for promoting critical endosomal escape. We now present here a systematic comparison of similar DNA ligands featuring our GCP oxo-anion-binding motif with DNA ligands only consisting of naturally occurring amino acids. Structure-activity studies showed that the artificial binding motif clearly outperformed natural amino acids such as histidine, lysine, and arginine. It improved the ability to shuttle foreign genetic material into cells, yet successfully mediated endosomal escape. Also, plasmids that were complexed by our artificial ligands were stabilized against cytosolic degradation to some extent. This resulted in the successful expression of plasmid information (comparable to gold standards such as polyethyleneimine). Hence, our study clearly demonstrates the importance of the tailor-made GCP anion-binding site for efficient gene transfection.

Microscopic understanding of the conformational features of a protein–DNA complex

Protein–DNA interactions play crucial roles in different stages of genetic activities, such as replication of genome, initiation of transcription,etc.

Importance of β2-β3 Loop Motion for the Increased Binding and Decreased Selectivity of the ΔLL Mutant of the Human Papillomavirus Type 6 E2 Protein

The binding affinity of the human papillomavirus type 6 E2 protein is strongly mediated by the sequence of the DNA linker region, with high affinity for the AATT linker and low affinity for the CCGG linker. When two terminal leucine residues are removed from the protein, the level of binding to both strands increases, but unequally, resulting in a significant decrease in selectivity for the AATT linker strand. To rationalize this behavior, we performed molecular dynamics simulations of the wild-type and mutant protein in the apo state and bound to DNA with high-affinity AATT and low-affinity CCGG linker strands. While no stable contacts were made between the β2-β3 loop and DNA in the wild type, this loop was repositioned in the mutant complexes and formed electrostatic contacts with the DNA backbone. More contacts were formed when the mutant was bound to the CCGG linker strand than to the AATT linker strand, resulting in a more favorable change in interaction energy for the CCGG strand. In addition, significant differences in correlated motions were found, which further explained the differences in binding. The simulations suggest that β2-β3 loop motions are responsible for the increased affinity and decreased selectivity of the mutant protein.

Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species

Background

The identification of the mechanisms of adaptation of protein structures to extreme environmental conditions is a challenging task of structural biology. We performed molecular dynamics (MD) simulations of the Nip7 protein involved in RNA processing from the shallow-water (P. furiosus) and the deep-water (P. abyssi) marine hyperthermophylic archaea at different temperatures (300 and 373 K) and pressures (0.1, 50 and 100 MPa). The aim was to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures.

Results

The current results demonstrate that the 3D models of the two proteins at all the examined values of pressures and temperatures are compact, stable and similar to the known crystal structure of the P. abyssi Nip7. The structural deviations and fluctuations in the polypeptide chain during the MD simulations were the most pronounced in the loop regions, their magnitude being larger for the C-terminal domain in both proteins. A number of highly mobile segments the protein globule presumably involved in protein-protein interactions were identified. Regions of the polypeptide chain with significant difference in conformational dynamics between the deep- and shallow-water proteins were identified.

Conclusions

The results of our analysis demonstrated that in the examined ranges of temperatures and pressures, increase in temperature has a stronger effect on change in the dynamic properties of the protein globule than the increase in pressure. The conformational changes of both the deep- and shallow-sea protein models under increasing temperature and pressure are non-uniform. Our current results indicate that amino acid substitutions between shallow- and deep-water proteins only slightly affect overall stability of two proteins. Rather, they may affect the interactions of the Nip7 protein with its protein or RNA partners.

Steric mechanism of auto-inhibitory regulation of specific and non-specific DNA binding by the ETS transcriptional repressor ETV6

DNA binding by the ETS transcriptional repressor ETV6 (or TEL) is auto-inhibited ~50-fold due to an α-helix that sterically blocks its ETS domain binding interface. Using NMR spectroscopy, we demonstrate that this marginally stable helix is unfolded, and not displaced to a non-inhibitory position, when ETV6 is bound to DNA containing a consensus (5')GGAA(3') recognition site. Although significantly lower in affinity, binding to non-specific DNA is auto-inhibited ~5-fold and is also accompanied by helix unfolding. Based on NMR chemical shift perturbations, both specific and non-specific DNA are bound via the same canonical ETS domain interface. However, spectral perturbations are smaller for the non-specific complex, suggesting weaker and less well-defined interactions than in the specific complex. In parallel, the crystal structure of ETV6 bound to a specific DNA duplex was determined. The structure of this complex reveals that a non-conserved histidine residue in the ETS domain recognition helix helps establish the specificity of ETV6 for DNA-binding sites containing (5')GGAA(3')versus(5')GGAT(3'). These studies provide a unified steric mechanism for attenuating ETV6 binding to both specific and non-specific DNA and expand the repertoire of characterized auto-inhibitory strategies utilized to regulate ETS factors.

DNA phosphate crowding correlates with protein cationic side chain density and helical curvature in protein/DNA crystal structures

Sequence-specific binding of proteins to their DNA targets involves a complex spectrum of processes that often induce DNA conformational variation in the bound complex. The forces imposed by protein binding that cause the helical deformations are intimately interrelated and difficult to parse or rank in importance. To investigate the role of electrostatics in helical deformation, we quantified the relationship between protein cationic residue density (Cpc) and DNA phosphate crowding (Cpp). The correlation between Cpc and Cpp was then calculated for a subset of 58 high resolution protein-DNA crystal structures. Those structures containing strong Cpc/Cpp correlation (>±0.25) were likely to contain DNA helical curvature. Further, the correlation factor sign predicted the direction of helical curvature with positive (16 structures) and negative (seven structures) correlation containing concave (DNA curved toward protein) and convex (DNA curved away from protein) curvature, respectively. Protein-DNA complexes without significant Cpc/Cpp (36 structures) correlation (-0.25<0<0.25) tended to contain DNA without significant curvature. Interestingly, concave and convex complexes also include more arginine and lysine phosphate contacts, respectively, whereas linear complexes included essentially equivalent numbers of Lys/Arg phosphate contacts. Together, these findings suggest an important role for electrostatic interactions in protein-DNA complexes involving helical curvature.

Molecular basis of the recognition of the ap65-1 gene transcription promoter elements by a Myb protein from the protozoan parasite Trichomonas vaginalis

Iron-inducible transcription of the ap65-1 gene in Trichomonas vaginalis involves at least three Myb-like transcriptional factors (tvMyb1, tvMyb2 and tvMyb3) that differentially bind to two closely spaced promoter sites, MRE-1/MRE-2r and MRE-2f. Here, we defined a fragment of tvMyb2 comprising residues 40-156 (tvMyb2₄₀₋₁₅₆) as the minimum structural unit that retains near full binding affinity with the promoter DNAs. Like c-Myb in vertebrates, the DNA-free tvMyb2₄₀₋₁₅₆ has a flexible and open conformation. Upon binding to the promoter DNA elements, tvMyb2₄₀₋₁₅₆ undergoes significant conformational re-arrangement and structure stabilization. Crystal structures of tvMyb2₄₀₋₁₅₆ in complex with promoter element-containing DNA oligomers showed that 5'-a/gACGAT-3' is the specific base sequence recognized by tvMyb2₄₀₋₁₅₆, which does not fully conform to that of the Myb binding site sequence. Furthermore, Lys⁴⁹, which is upstream of the R2 motif (amino acids 52-102) also participates in specific DNA sequence recognition. Intriguingly, tvMyb2₄₀₋₁₅₆ binds to the promoter elements in an orientation opposite to that proposed in the HADDOCK model of the tvMyb1₃₅₋₁₄₁/MRE-1-MRE-2r complex. These results shed new light on understanding the molecular mechanism of Myb-DNA recognition and provide a framework to study the molecular basis of transcriptional regulation of myriad Mybs in T. vaginalis.

Internal dynamics of the tryptophan repressor (TrpR) and two functionally distinct TrpR variants, L75F-TrpR and A77V-TrpR, in their l-Trp-bound forms

Backbone amide dynamics of the Escherichia coli tryptophan repressor protein (WT-TrpR) and two functionally distinct variants, L75F-TrpR and A77V-TrpR, in their holo (l-tryptophan corepressor-bound) form have been characterized using (15)N nuclear magnetic resonance (NMR) relaxation. The three proteins possess very similar structures, ruling out major conformational differences as the source of their functional differences, and suggest that changes in protein flexibility are at the origin of their distinct functional properties. Comparison of site specific (15)N-T(1), (15)N-T(2), (15)N-{(1)H} nuclear Overhauser effect, reduced spectral density, and generalized order (S(2)) parameters indicates that backbone dynamics in the three holo-repressors are overall very similar with a few notable and significant exceptions for backbone atoms residing within the proteins' DNA-binding domain. We find that flexibility is highly restricted for amides in core α-helices (i.e., helices A-C and F), and a comparable "stiffening" is observed for residues in the DNA recognition helix (helix E) of the helix D-turn-helix E (HTH) DNA-binding domain of the three holo-repressors. Unexpectedly, amides located in helix D and in adjacent turn regions remain flexible. These data support the concept that residual flexibility in TrpR is essential for repressor function, DNA binding, and molecular recognition of target operators. Comparison of the (15)N NMR relaxation parameters of the holo-TrpRs with those of the apo-TrpRs indicates that the single-point amino acid substitutions, L75F and A77V, perturb the flexibility of backbone amides of TrpR in very different ways and are most pronounced in the apo forms of the three repressors. Finally, we present these findings in the context of other DNA-binding proteins and the role of protein flexibility in molecular recognition.