Circular codes revisited: a statistical approach

In 1996 Arquès and Michel [1996. A complementary circular code in the protein coding genes. J. Theor. Biol. 182, 45-58] discovered the existence of a common circular code in eukaryote and prokaryote genomes. Since then, circular code theory has provoked great interest and underwent a rapid development. In this paper we discuss some theoretical issues related to the synchronization properties of coding sequences and circular codes with particular emphasis on the problem of retrieval and maintenance of the reading frame. Motivated by the theoretical discussion, we adopt a rigorous statistical approach in order to try to answer different questions. First, we investigate the covering capability of the whole class of 216 self-complementary, C³ maximal codes with respect to a large set of coding sequences. The results indicate that, on average, the code proposed by Arquès and Michel has the best covering capability but, still, there exists a great variability among sequences. Second, we focus on such code and explore the role played by the proportion of the bases by means of a hierarchy of permutation tests. The results show the existence of a sort of optimization mechanism such that coding sequences are tailored as to maximize or minimize the coverage of circular codes on specific reading frames. Such optimization clearly relates the function of circular codes with reading frame synchronization.     

New advances in the study on the interaction of [Cr(phen)2(dppz)]3+ complex with biological models; association to transporting proteins

The present study reports a detailed investigation with the interaction of [Cr(phen)₂(dppz)]³⁺ with serum albumins, the key protein for the transport of drugs in the blood plasma, which allows us to understand further the role of [Cr(phen)₂(dppz)]³⁺ as sensitizer in Photodynamic Therapy (PDT).Chromium(III) complex [Cr(phen)₂(dppz)]³⁺, (dppz = dipyridophenazine and phen = 1,10-phenanthroline), where dppz is a planar bidentate ligand with an extended π system, has been found to bind strongly with bovine and human serum albumins (BSA and HSA) with an intrinsic binding constants, K_b, of (1.7 ± 0.3) × 10⁵ M^− 1 and (2.2 ± 0.3) × 10⁵ M^− 1 at 295 K, respectively. The interactions of serum albumins with [Cr(phen)₂(dppz)]³⁺ were assessed employing fluorescence spectroscopy, circular dichroism and UV-vis absorption spectroscopy. The serum albumins-[Cr(phen)₂(dppz)]³⁺ interactions caused conformational changes with the loss of helical stability of the protein and local perturbation in the domain IIA binding pocket. The relative fluorescence intensity of the albumin (BSA or HSA) bound to the Cr(III) complex decreased, suggesting that perturbation around the Trp 214 residue took place. The analysis of the thermodynamic parameters ΔG, ΔH, ΔS indicated that the hydrophobic interactions played a major role in both BSA-Cr(III) and HSA-Cr(III) association processes. The binding distances and transfer efficiencies for BSA-Cr(III) and HSA-Cr(III) binding reactions were calculated according to the Föster theory of non-radiation energy transfer. All these experimental results suggests that [Cr(phen)₂(dppz)]³⁺ binds to serum albumins, by which these proteins could act as carriers of this complex for further applications in PDT.     

Functional coupling between a distal interaction and the cleavage site in bacterial RNase-P-RNA-mediated cleavage

Bacterial RNase P consists of one protein and one RNA [RNase P RNA (RPR)]. RPR can process tRNA precursors correctly in the absence of the protein. Here we have used model hairpin loop substrates corresponding to the acceptor, T-stem, and T-loop of a precursor tRNA to study the importance of the T-loop structure in RPR-alone reaction. T-stem/loop (TSL) interacts with a region in RPR [TSL binding site (TBS)], forming TSL/TBS interaction. Altering the T-loop structure affects both cleavage site selection and rate of cleavage at the correct site + 1 and at the alternative site − 1. The magnitude of variation depended on the structures of the T-loop and the TBS region, with as much as a 150-fold reduction in the rate of cleavage at + 1. Interestingly, for one T-loop structure mutant, no difference in the rate at − 1 was detected compared to cleavage of the substrate with an unchanged T-loop, indicating that, in this case, the altered T-loop structure primarily influences events required for efficient cleavage at the correct site + 1. We also provide data supporting a functional link between a productive TSL/TBS interaction and events at the cleavage site. Collectively, our findings emphasize the interplay between separate regions upon formation of a productive RPR substrate that leads to efficient and accurate cleavage. These new data provide support for an induced-fit mechanism in bacterial RPR-mediated cleavage at the correct site + 1.     

Shortened engineered human antibody CH2 domains: increased stability and binding to the human neonatal Fc receptor

The immunoglobulin (Ig) constant CH2 domain is critical for antibody effector functions. Isolated CH2 domains are promising scaffolds for construction of libraries containing diverse binders that could also confer some effector functions. We have shown previously that an isolated human CH2 domain is relatively unstable to thermally induced unfolding, but its stability can be improved by engineering an additional disulfide bond (Gong, R., Vu, B. K., Feng, Y., Prieto, D. A., Dyba, M. A., Walsh, J. D., Prabakaran, P., Veenstra, T. D., Tarasov, S. G., Ishima, R., and Dimitrov, D. S. (2009) J. Biol. Chem. 284, 14203-14210). We have hypothesized that the stability of this engineered antibody domain could be further increased by removing unstructured residues. To test our hypothesis, we removed the seven N-terminal residues that are in a random coil as suggested by our analysis of the isolated CH2 crystal structure and NMR data. The resulting shortened engineered CH2 (m01s) was highly soluble, monomeric, and remarkably stable, with a melting temperature (T(m)) of 82.6 °C, which is about 10 and 30 °C higher than those of the original stabilized CH2 (m01) and CH2, respectively. m01s and m01 were more resistant to protease digestion than CH2. A newly identified anti-CH2 antibody that recognizes a conformational epitope bound to m01s significantly better (>10-fold higher affinity) than to CH2 and slightly better than to m01. m01s bound to a recombinant soluble human neonatal Fc receptor at pH 6.0 more strongly than CH2. These data suggest that shortening the m01 N terminus significantly increases stability without disrupting its conformation and that our approach for increasing stability and decreasing size by removing unstructured regions may also apply to other proteins.     

NMR, biophysical, and biochemical studies reveal the minimal Calmodulin binding domain of the HIV-1 matrix protein

Subcellular distribution of Calmodulin (CaM) in human immunodeficiency virus type-1 (HIV-1)-infected cells is distinct from that observed in uninfected cells. CaM has been shown to interact and co-localize with the HIV-1 Gag protein in infected cells. However, the precise molecular mechanism of this interaction is not known. Binding of Gag to CaM is dependent on calcium and is mediated by the N-terminal-myristoylated matrix (myr(+)MA) domain. We have recently shown that CaM binding induces a conformational change in the MA protein, triggering exposure of the myristate group. To unravel the molecular mechanism of CaM-MA interaction and to identify the minimal CaM binding domain of MA, we devised multiple approaches utilizing NMR, biochemical, and biophysical methods. Short peptides derived from the MA protein have been examined. Our data revealed that whereas peptides spanning residues 11-28 (MA-(11-28)) and 31-46 (MA-(31-46)) appear to bind preferentially to the C-terminal lobe of CaM, a peptide comprising residues 11-46 (MA-(11-46)) appears to engage both domains of CaM. Limited proteolysis data conducted on the MA-CaM complex yielded a MA peptide (residues 8-43) that is protected by CaM and resistant to proteolysis. MA-(8-43) binds to CaM with a very high affinity (dissociation constant = 25 nm) and in a manner that is similar to that observed for the full-length MA protein. The present findings provide new insights on how MA interacts with CaM that may ultimately help in identification of the functional role of CaM-Gag interactions in the HIV replication cycle.     

The triazatruxene derivative azatrux binds to the parallel form of the human telomeric G-quadruplex under molecular crowding conditions: biophysical and molecular modeling studies

The present study has employed a combination of spectroscopic, calorimetric and computational methods to explore the binding of the three side-chained triazatruxene derivative, termed azatrux, to a human telomeric G-quadruplex sequence, under conditions of molecular crowding. The binding of azatrux to the tetramolecular parallel [d(TGGGGT)]₄ quadruplex in the presence and absence of crowding conditions, was also characterized. The data indicate that azatrux binds in an end-stacking mode to the parallel G-quadruplex scaffold and highlights the key structural elements involved in the binding. The selectivity of azatrux for the human telomeric G-quadruplex relative to another biologically relevant G-quadruplex (c-Kit87up) and to duplex DNA was also investigated under molecular crowding conditions, showing that azatrux has good selectivity for the human telomeric G-quadruplex over the other investigated DNA structures.     

Structural and functional analysis of substrate recognition by the 250s loop in amylomaltase from Thermus brockianus

Amylomaltase, or 4‐α‐glucanotransferase (EC 2.4.1.25), is involved in glycogen and maltooligosaccharide metabolism in microorganisms, catalyzing both the hydrolysis and transfer of an α‐1,4‐oligosacchraride to other sugar molecules. In this study, we determined the crystal structure of amylomaltase from Thermus brockianus at a resolution of 2.3 Å and conducted a biochemical study to understand the detailed mechanism for its activity. Careful comparison with previous amylomaltase structures showed a pattern of conformational flexibility in the 250s loop with higher B‐factor. Amylomaltase from T. brockianus exhibited a high transglycosylation factor for glucose and a lower value for maltose. Mutation of Gln256 resulted in increased K_m for maltotriose and a sharp decrease of the transglycosylation factor for maltose, suggesting the involvement of Gln 256 in substrate binding between subsites +1 and +2. Mutation of Phe251 resulted in significantly lower glucose production but increased maltose production from maltopentose substrates, showing an altered substrate‐binding affinity. The mutational data suggest the conformational flexibility of the loop may be involved in substrate binding in the GH77 family. Here, we present an action model of the 250s loop providing the molecular basis for the involvement of residues Phe251, Gln256, and Trp258 in the hydrolysis and transglycosylation activities in amylomaltase. Proteins 2011. © 2010 Wiley‐Liss, Inc.     

Crystal structure of the novel PaiA N-acetyltransferase from Thermoplasma acidophilum involved in the negative control of sporulation and degradative enzyme production

GCN5-related N-acetyltransferases (GNATs) are the most widely distributed acetyltransferase systems among all three domains of life. GNATs appear to be involved in several key processes, including microbial antibiotic resistance, compacting eukaryotic DNA, controlling gene expression, and protein synthesis. Here, we report the crystal structure of a putative GNAT Ta0374 from Thermoplasma acidophilum, a hyperacidophilic bacterium, that has been determined in an apo-form, in complex with its natural ligand (acetyl coenzyme A), and in complex with a product of reaction (coenzyme A) obtained by cocrystallization with spermidine. Sequence and structural analysis reveals that Ta0374 belongs to a novel protein family, PaiA, involved in the negative control of sporulation and degradative enzyme production. The crystal structure of Ta0374 confirms that it binds acetyl coenzyme A in a way similar to other GNATs and is capable of acetylating spermidine. Based on structural and docking analysis, it is expected that Glu53 and Tyr93 are key residues for recognizing spermidine. Additionally, we find that the purification His-Tag in the apo-form structure of Ta0374 prevents binding of acetyl coenzyme A in the crystal, though not in solution, and affects a chain-flip rotation of "motif A" which is the most conserved sequence among canonical acetyltransferases.     

Protein loop selection using orientation-dependent force fields derived by parameter optimization

We used the orientation‐dependent Optimized Side Chain Atomic eneRgy (OSCAR‐o), derived in an early study, for protein loop selection. The prediction accuracy of OSCAR‐o was better than that of physics‐based force fields or statistical potential energy functions for both the RAPPER decoy set and the Jacobson decoy set. The native conformer was frequently ranked as lowest energy among the decoys. Furthermore, strong correlation was observed between the OSCAR‐o score and the root mean square deviation (RMSD) from the native structure for energy‐minimized decoys. In practical use, we applied OSCAR‐o to rescore decoys generated by a widely used loop‐modeling program, LOOPY. As a result, the mean RMSD values of top‐ranked decoys were reduced by 0.3 Å for loop targets of seven to nine residues. We expect similar performance for OSCAR‐o with other loop‐modeling algorithms in the context of decoy rescoring. A loop selection program (OSCAR‐ls) based on OSCAR‐o is available at http://sysimm.ifrec.osaka‐u.ac.jp/OSCAR/ . Proteins 2011; © 2011 Wiley‐Liss, Inc.