PubMed Health⌕ Search

Biomedical subjects

Hadas Leonov

Publications and source records attributed to Hadas Leonov.

5 recordsLinked to original sources

Genetic algorithm-based optimization of hydrophobicity tables.

SUMMARY: The genomic abundance and pharmacological importance of membrane proteins have fueled efforts to identify them based solely on sequence information. Previous methods based on the physicochemical principle of a sliding window of hydrophobicity (hydropathy analysis) have been replaced by approaches based on hidden Markov models or neural networks which prevail due to their probabilistic orientation. In the current study, an optimization of the hydrophobicity tables used in hydropathy analysis is performed using a genetic algorithm. As such, the approach can be viewed as a synthesis between the physicochemically and statistically based methods. The resulting hydrophobicity tables lead to significant improvement in the prediction accuracy of hydropathy analysis. Furthermore, since hydropathy analysis is less dependent on the basis set of membrane proteins is used to hone the statistically based methods, as well as being faster, it may be valuable in the analysis of new genomes. Finally, the values obtained for each of the amino acids in the new hydrophobicity tables are discussed.

Algorithms↗

A periodicity analysis of transmembrane helices.

Transmembrane helices and the helical bundles which they form are the major building blocks of membrane proteins. Since helices are characterized by a given periodicity, it is possible to search for patterns of traits which typify one side of the helix and not the other (e.g. amphipathic helices contain a polar and apolar sides). Using Fourier transformation we have analyzed solved membrane protein structures as well as sequences of membrane proteins from the Swiss-Prot database. The traits searched included aromaticity, volume and ionization. While a number of motifs were already recognized in the literature, many were not. One particular example involved helix VII of lactose permease which contains seven aromatic residues on six helical turns. Similarly six glycine residues in four consecutive helical turns were identified as forming a motif in the chloride channel. A tabulation of all the findings is presented as well as a possible rationalization of the function of the motif.

Algorithms↗

Monte Carlo estimation of the number of possible protein folds: effects of sampling bias and folds distributions.

The estimation of the number of protein folds in nature is a matter of considerable interest. In this study, a Monte Carlo method employing the broken stick model is used to assign a given number of proteins into a given number of folds. Subsequently, random, integer, non-repeating numbers are generated in order to simulate the process of fold discovery. With this conceptual framework at hand, the effects of two factors upon the fold identification process were investigated: (1) the nature of folds distributions and (2) preferential sampling bias of previously identified folds. Depending on the type of distribution, dividing 100,000 proteins into 1,000 folds resulted in 10-30% of the folds having 10 proteins or less per fold, approximately 10% of the folds having 10-20 proteins per fold, 31-45% having 20-100 proteins per fold, and >30% of the folds having more than 100 proteins per fold. After randomly sampling one tenth of the proteins, 68-96% of the folds were identified. These percentages depend both on folds distribution and biased/non-biased sampling. Only upon increasing the sampling bias for previously identified folds to 1,000, did the model result in a reduction of the number of proteins identified by an order of magnitude (approximately 9%). Thus, assuming the structures of one tenth of the population of proteins in nature have been solved, the results of the Monte Carlo simulation are more consistent with recent lower estimates of the number of folds, <or=1,000. Any deviation from this estimate would reflect significant bias in the experimental sampling of protein structure, and/or substantially nonuniform folds distribution, manifested in a large number of single-fold proteins.

Models, Theoretical↗

Syntaxin 1A modulates the voltage-gated L-type calcium channel (Ca(v)1.2) in a cooperative manner.

Syntaxin 1A (Sx1A) modifies the activity of voltage-gated Ca2+ channels acting via the cytosolic and the two vicinal cysteines (271 and 272) at the transmembrane domain. Here we show that Sx1A modulates the Lc-type Ca2+ channel, Cav1.2, in a cooperative manner, and we explore whether channel clustering or the Sx1A homodimer is responsible for this activity. Sx1A formed homodimers but, when mutated at the two vicinal transmembrane domain cysteines, was unable to either dimerize or modify the channel activity suggesting disulfide bond formation. Moreover, applying global molecular dynamic search established a theoretical prospect of generating a disulfide bond between two Sx1A transmembrane helices. Nevertheless, Sx1A activity was not correlated with Sx1A homodimer. Application of a vicinal thiol reagent, phenylarsine oxide, abolished Sx1A action indicating the accessibility of Cys-271,272 thiols. Sx1A inhibition of channel activity was restored by phenylarsine oxide antidote, 2,3-dimercaptopropanol, consistent with thiol interaction of Sx1A. In addition, the supralinear mode of channel inhibition was correlated to the monomeric form of Sx1A and was apparent only when the three channel subunits alpha11.2/alpha2delta1/beta2a were present. This functional demonstration of cooperativity suggests that the three-subunit channel responds as a cluster, and Sx1A monomers associate with a dimer (or more) of a three-subunit Ca2+ channel. Consistent with channel cluster linked to Sx1A, a conformational change driven by membrane depolarization and Ca2+ entry would rapidly be transduced to the exocytotic machinery. As shown herein, the supralinear relationship between Sx1A and the voltage-gated Ca2+ channel within the cluster could convey the cooperativity that distinguishes the process of neurotransmitter release.

Animals↗

Modeling the structure of the respiratory syncytial virus small hydrophobic protein by silent-mutation analysis of global searching molecular dynamics.

Human respiratory syncytial virus (RSV) encodes a small hydrophobic (SH) protein, whose function in the life cycle of the virus is unknown. Recent channel activity measurements of the protein suggest that like other viroporins, SH may assemble into a homo-oligomeric ion channel. To further our understanding of this potentially important protein, a new strategy was implemented in order to model the transmembrane oligomeric bundle of the protein. Global searching molecular dynamic simulations of SH proteins from eight different viral strains, each at different oligomeric states, as well as different lengths of the putative transmembrane domain, were undertaken. Taken together, a total of 45 different global molecular dynamic simulations pointed to a single pentameric structure for the protein that was found in all of the different variants. The model of the structure obtained is a channel-like homopentamer whose minimal transmembrane pore diameter is 1.46 A.

Amino Acid Sequence↗