Phylogenetic diversification of immunoglobulin genes and the antibody repertoire

Immunoglobulins are encoded by a large multigene system that undergoes somatic rearrangement and additional genetic change during the development of immunoglobulin-producing cells. Inducible antibody and antibody-like responses are found in all vertebrates. However, immunoglobulin possessing disulfide-bonded heavy and light chains and domain-type organization has been described only in representatives of the jawed vertebrates. High degrees of nucleotide and predicted amino acid sequence identity are evident when the segmental elements that constitute the immunoglobulin gene loci in phylogenetically divergent vertebrates are compared. However, the organization of gene loci and the manner in which the independent elements recombine (and diversify) vary markedly among different taxa. One striking pattern of gene organization is the "cluster type" that appears to be restricted to the chondrichthyes (cartilaginous fishes) and limits segmental rearrangement to closely linked elements. This type of gene organization is associated with both heavy- and light-chain gene loci. In some cases, the clusters are "joined" or "partially joined" in the germ line, in effect predetermining or partially predetermining, respectively, the encoded specificities (the assumption being that these are expressed) of the individual loci. By relating the sequences of transcribed gene products to their respective germ-line genes, it is evident that, in some cases, joined-type genes are expressed. This raises a question about the existence and/or nature of allelic exclusion in these species. The extensive variation in gene organization found throughout the vertebrate species may relate directly to the role of intersegmental (V<==>D<==>J) distances in the commitment of the individual antibody-producing cell to a particular genetic specificity. Thus, the evolution of this locus, perhaps more so than that of others, may reflect the interrelationships between genetic organization and function.      

The distribution of residues in a polypeptide sequence is a determinant of aggregation optimized by evolution

It has been shown that the propensity of a protein to form amyloid-like fibrils can be predicted with high accuracy from the knowledge of its amino acid sequence. It has also been suggested, however, that some regions of the sequences are more important than others in determining the aggregation process. Here, we have addressed this issue by constructing a set of “sequence scrambled” variants of the first 29 residues of horse heart apomyoglobin (apoMb_1-29), in which the sequence was modified while maintaining the same amino acid composition. The clustering of the most amyloidogenic residues in one region of the sequence was found to cause a marked increase of the elongation rate (k_agg) and a remarkable shortening of the lag phase (t_lag) of the fibril growth, as determined by far-UV circular dichroism and thioflavin T fluorescence. We also show that taking explicitly into consideration the presence of aggregation-promoting regions in the predictive methods results in a quantitative agreement between the theoretical and observed k_agg and t_lag values of the apoMb_1-29 variants. These results, together with a comparison between homologous segments from the family of globins, indicate the existence of a negative selection against the clustering of highly amyloidogenic residues in one or few regions of polypeptide sequences.     

A superposition free method for protein conformational ensemble analyses and local clustering based on a differential geometry representation of backbone

Here a differential geometry (DG) representation of protein backbone is explored on the analyses of protein conformational ensembles. The protein backbone is described by curvature, κ, and torsion, τ, values per residue and we propose 1) a new dissimilarity and protein flexibility measurement and 2) a local conformational clustering method. The methods were applied to Ubiquitin and c-Myb-KIX protein conformational ensembles and results show that κ\τ metric space allows to properly judge protein flexibility by avoiding the superposition problem. The d_max measurement presents equally good or superior results when compared to RMSF, especially for the intrinsically unstructured protein. The clustering method is unique as it relates protein global to local dynamics by providing a global clustering solutions per residue. The methods proposed can be especially useful to the analyses of highly flexible proteins. The software written for the analyses presented here is available at https://github.com/AMarinhoSN/FleXgeo for academic usage only.     

Calculation of mutational free energy changes in transition states for protein folding

Recent advances in experimental and computational methods have made it possible to determine with considerable accuracy the structures whose formation is rate limiting for the folding of some small proteins-the transition state ensemble, or TSE. We present a method to analyze and validate all-atom models of such structures. The method is based on the comparison of experimental data with the computation of the change in free energy of the TSE resulting from specific mutations. Each mutation is modeled individually in all members of an ensemble of transition state structures using a method originally developed to predict mutational changes in the stability of native proteins. We first apply this method to six proteins for which we have determined the TSEs with a technique that uses experimental mutational data (Phi-values) as restraints in the structure determination and find a highly significant correlation between the calculated free energy changes and those derived from experimental kinetic data. We then use the procedure to analyze transition state structures determined by molecular dynamics simulations of unfolding, again finding a high correlation. Finally, we use the method to estimate changes in folding rates of several hydrophobic core mutants of Fyn SH3. Taken together, these results show that the procedure developed here is a tool of general validity for analyzing, assessing, and improving the quality of the structures of transition states for protein folding.     

Validity of Gō models: comparison with a solvent-shielded empirical energy decomposition

Do Gō-type model potentials provide a valid approach for studying protein folding? They have been widely used for this purpose because of their simplicity and the speed of simulations based on their use. The essential assumption in such models is that only contact interactions existing in the native state determine the energy surface of a polypeptide chain, even for non-native configurations sampled along folding trajectories. Here we use an all-atom molecular mechanics energy function to investigate the adequacy of Gō-type potentials. We show that, although the contact approximation is accurate, non-native contributions to the energy can be significant. The assumed relation between residue-residue interaction energies and the number of contacts between them is found to be only approximate. By contrast, individual residue energies correlate very well with the number of contacts. The results demonstrate that models based on the latter should give meaningful results (e.g., as used to interpret phi values), whereas those that depend on the former are only qualitative, at best.     

Toward an energy function for the contact map representation of proteins

We analyzed several energy functions for predicting the native state of proteins from an energy minimization procedure. We derived the parameters of a given energy function by imposing the basic requirement that the energy of the native conformation of a protein is lower than that of any conformation chosen from a set of decoys. Our work is motivated by a recent result which proved that the simple pairwise contact approximation of the energy is insufficient to satisfy simultaneously such a basic requirement for all the proteins in a database. Here, we investigate the reasons of such negative results and show how to improve the predictive power of methods based on energy minimization. We generated decoys by gapless threading, and we derive energy parameters by perceptron learning. We first considered hydrophobic contributions to the energy, defined in several ways, and showed that the additional hydrophobic terms enlarge slightly the number of proteins that can be stabilized together. Next, we performed various modifications of the pairwise energy term. We introduced (1) a distinction between inter-residue contacts on the surface and in the core of a protein and (2) a simple distance-dependent pairwise interaction in which a two-tier definition of contact replaces the original (single-tier) one. Our results suggest that a detailed treatment of the pairwise potential is likely to be more relevant than the consideration of other forces.     

A new fourier transform approach for protein coding measure based on the format of the Z curve

MOTIVATION: At the core of most protein gene-finding algorithms are the coding measures used to make a decision on coding/non-coding. Of the protein coding measures, the Fourier measure is one of the most important. However, due to the limited length of the windows usually used, the accuracy of the measure is not satisfactory. This paper is devoted to improving the accuracy by lengthening the sequence to amplify the periodicity of 3 in the coding regions. RESULTS: A new algorithm is presented called the lengthen-shuffle Fourier transform algorithm. For the same window length, the percentage accuracy of the new algorithm is 6-7% higher than that of the ordinary Fourier transform algorithm. The resulting percentage accuracy (average of specificity and sensitivity) of the new measure is 84.9% for the window length 162 bp. AVAILABILITY: The program is available on request fromC.- T. Zhang. Contact: ctzhang@tju.edu.cn      

Physicochemical and mechanical characterization of galactomannan from Mimosa scabrella: Effect of drying method

The galactomannan from Mimosa scabrella Bentham was extracted on a pilot plant scale and dried either by vacuum oven (GVO) or by spray-drier (GSD) to evaluate the effect of the drying technique on the powder quality and its applicability as excipient in solid dosage forms. The analysis by high performance size exclusion chromatography (HPSEC) coupled to multiangle laser light scattering (MALLS) suggests that both products behave as semi-flexible polymers, although the GSD showed more aggregation at molecular level (～10%) and higher chain stiffness (Lp 9.1 nm). TG and DSC analysis showed weight loss event with peak at 299.7 and 311.9 °C to GVO and GSD, respectively, as well and higher ash content for GSD sample, in both inert and oxidant atmosphere. The X-ray diffraction confirmed the amorphous nature of both galactomannans although GVO showed high crystallinity. The GSD showed lower density (1.009 g/cm³), higher cohesiveness (repose angle 35.5°, compressibility 32.2% and absence of flow), smaller and more spherical particles than GVO sample, both with high polydispersion. As vacuum oven-drying resulted in a like fibrous material, spray-drying appears as an alternative method easy to extrapolate in industry, requiring a glidant incorporation to improve the powder flowability.