In silico protein fragmentation reveals the importance of critical nuclei on domain reassembly

Codon Selection Affects Recruitment of Ribosome-Associating Factors during Translation

An intriguing aspect of protein synthesis is how cotranslational events are managed inside the cell. In this study, we developed an in vivo bimolecular fluorescence complementation assay coupled to SecM stalling (BiFC-SecM) to study how codon usage influences the interactions of ribosome-associating factors that occur cotranslationally. We profiled ribosomal associations of a number of proteins, and observed differential association of chaperone proteins TF, DnaK, GroEL, and translocation factor Ffh as a result of introducing synonymous codon substitutions that change the affinity of the translating sequence to the ribosomal anti-Shine-Dalgarno (aSD) sequence. The use of pausing sequences within proteins regulates their transit within the translating ribosome. Our results indicate that the dynamics between cellular factors and the new polypeptide chain are affected by how codon composition is designed. Furthermore, associating factors may play a role in processes including protein quality control (folding and degradation) and cellular respiration.

Oriented covalent immobilization of antibodies for measurement of intermolecular binding forces between zipper-like contact surfaces of split inteins

In order to measure the intermolecular binding forces between two halves (or partners) of naturally split protein splicing elements called inteins, a novel thiol-hydrazide linker was designed and used to orient immobilized antibodies specific for each partner. Activation of the surfaces was achieved in one step, allowing direct intermolecular force measurement of the binding of the two partners of the split intein (called protein trans-splicing). Through this binding process, a whole functional intein is formed resulting in subsequent splicing. Atomic force microscopy (AFM) was used to directly measure the split intein partner binding at 1 μm/s between native (wild-type) and mixed pairs of C- and N-terminal partners of naturally occurring split inteins from three cyanobacteria. Native and mixed pairs exhibit similar binding forces within the error of the measurement technique (~52 pN). Bioinformatic sequence analysis and computational structural analysis discovered a zipper-like contact between the two partners with electrostatic and nonpolar attraction between multiple aligned ion pairs and hydrophobic residues. Also, we tested the Jarzynski's equality and demonstrated, as expected, that nonequilibrium dissipative measurements obtained here gave larger energies of interaction as compared with those for equilibrium. Hence, AFM coupled with our immobilization strategy and computational studies provides a useful analytical tool for the direct measurement of intermolecular association of split inteins and could be extended to any interacting protein pair.

Mesophilic and hyperthermophilic adenylate kinases differ in their tolerance to random fragmentation

The extent to which thermostability influences the location of protein fragmentation sites that allow retention of function is not known. To evaluate this, we used a novel transposase-based approach to create libraries of vectors that express structurally-related fragments of Bacillus subtilis adenylate kinase (BsAK) and Thermotoga neapolitana adenylate kinase (TnAK) with identical modifications at their termini, and we selected for variants in each library that complement the growth of Escherichia coli with a temperature-sensitive adenylate kinase (AK). Mutants created using the hyperthermophilic TnAK were found to support growth with a higher frequency (44%) than those generated from the mesophilic BsAK (6%), and selected TnAK mutants complemented E. coli growth more strongly than homologous BsAK variants. Sequencing of functional clones from each library also identified a greater dispersion of fragmentation sites within TnAK. Nondisruptive fission sites were observed within the AMP binding and core domains of both AK homologs. However, only TnAK contained sites within the lid domain, which undergoes dynamic fluctuations that are critical for catalysis. These findings implicate the flexible lid domain as having an increased sensitivity to fission events at physiological temperatures. In addition, they provide evidence that comparisons of nondisruptive fission sites in homologous proteins could be useful for finding dynamic regions whose conformational fluctuations are important for function, and they show that the discovery of protein fragments that cooperatively function in mesophiles can be aided by the use of thermophilic enzymes as starting points for protein design.

Kinetics and reaction coordinates of the reassembly of protein fragments via forward flux sampling

We studied the mechanism of the reassembly and folding process of two fragments of a split lattice protein by using forward flux sampling (FFS). Our results confirmed previous thermodynamics and kinetics analyses that suggested that the disruption of the critical core (of an unsplit protein that folds by a nucleation mechanism) plays a key role in the reassembly mechanism of the split system. For several split systems derived from a parent 48-mer model, we estimated the reaction coordinates in terms of collective variables by using the FFS least-square estimation method and found that the reassembly transition is best described by a combination of the total number of native contacts, the number of interchain native contacts, and the total conformational energy of the split system. We also analyzed the transition path ensemble obtained from FFS simulations using the estimated reaction coordinates as order parameters to identify the microscopic features that differentiate the reassembly of the different split systems studied. We found that in the fastest folding split system, a balanced distribution of the original-core amino acids (of the unsplit system) between protein fragments propitiates interchain interactions at early stages of the folding process. Only this system exhibits a different reassembly mechanism from that of the unsplit protein, involving the formation of a different folding nucleus. In the slowest folding system, the concentration of the folding nucleus in one fragment causes its early prefolding, whereas the second fragment tends to remain as a detached random coil. We also show that the reassembly rate can be either increased or decreased by tuning interchain cooperativeness via the introduction of a single point mutation that either strengthens or weakens one of the native interchain contacts (prevalent in the transition state ensemble).

Transition path sampling and forward flux sampling. Applications to biological systems

The last decade has seen a rapid growth in the number of simulation methods and applications dealing with the sampling of transition pathways of rare nanoscale events. Such studies are crucial, for example, for understanding the mechanism and kinetics of conformational transitions and enzymatic events associated with the function of biomolecules. In this review, a broad account of transition path sampling approaches is provided, starting from the general concepts, progressing to the specific principles that underlie some of the most important methods, and eventually singling out the so-called forward flux sampling method for a more detailed description. This is done because forward flux sampling, despite its appealing simplicity and potential efficiency, has thus far received limited attention from practitioners. While path sampling methods have a widespread application to many types of rare transitional events, here only recent applications involving biomolecules are reviewed, including isomerization, protein folding, and enzyme catalysis.

Prediction of the protein folding core: application to the immunoglobulin fold

We propose an algorithm that allows predicting residues important for the formation of the structure of globular proteins. It relies on a simulation that detects the amino acids presenting a maximum number of neighbours during the early steps of the folding process. They have been called MIR (Most Interacting Residues). Independently, description of the protein structures in fragments with closed ends shows the correlation between these extremities and the core of the globules. These fragments are of rather constant length, typically between 20 and 25 amino acids, and we have previously shown that their extremities are preferentially occupied by MIR. Introduction of rules derived from this fragment analysis of tertiary structures allows to smooth the distribution of MIR, for a better match between TEF ends and MIR. In order to assess this prediction of the folding core, a large family of structures has been used, with sequences as different as possible. A dataset of 56 immunoglobulin structures of various functions but common fold has been used in this study. This fold was chosen because it is one of the most populated with a large amount of data available on its nucleus. In the immunoglobulin domain, "functional and structural load is clearly separated: loops are responsible for binding and recognition while interactions between several residues of the buried core provide stability and fast folding"[1]. We then determined the positions susceptible of high importance for the folding process to occur and compared them to published data, either to High Throw Out Order (HTOO), Conservatism of Conservatism (CoC) or Phi value experiments. It results a reasonable agreement between the positions that we predict and experimental data. Besides, our prediction goes beyond the simple use of a null solvent accessibility of amino acids as a criterion to predict the core. We find the same quality of our prediction on the flavodoxin like superfamily.

Interaction studies and alanine scanning analysis of a semi-synthetic split intein reveal thiazoline ring formation from an intermediate of the protein splicing reaction

We recently reported an artificially split intein based on the Ssp DnaB mini-intein that consists of a synthetic N-terminal intein fragment (Int(N)) and a recombinant C-terminal part (Int(C)), which are 11 and 143 amino acids in length, respectively. This intein holds great promise for the preparation of semi-synthetic proteins by protein trans-splicing. In this work we synthesized a set of Int(N) peptide variants to investigate their structure-function relationship with regard to fragment association and promotion of protein trans-splicing. A further truncation of the Int(N) sequence below 11 amino acids resulted in loss of activity, whereas C-terminal extensions were tolerated. Alanine scanning analysis identified three essential hydrophobic residues, whereas substitutions at other positions were tolerated. We developed assays to monitor association of Int(N) with an Int(C) mutant blocked in protein splicing by native PAGE and fluorescence anisotropy. The kinetic parameters of intein complex formation were K(d) = 1.1 mum, k(on) = 16.8 m(-1) s(-1), and k(off) = 1.8 x 10(-5) s(-1) for the native Int(N11) sequence. Intriguingly, a G(-1)A substitution, previously known to significantly impair protein splicing, was revealed to result in thiazoline ring formation involving the catalytic Cys-1, likely by aberrant dehydration of a oxythiazolidine intermediate. This finding provides experimental evidence for the postulated intermediate during the initial N/S acyl shift and underlines the delicate spatial and temporal alignment required in the intein active site to prevent side reactions of the protein-splicing pathway.