CONFORMATIONAL DYNAMICS OF K-RAS AND H-RAS PROTEINS: IS THERE FUNCTIONAL SPECIFICITY AT THE CATALYTIC DOMAIN?

Ras proteins serve as crucial signaling modulators in cell proliferation through their ability to hydrolyze GTP and exist in a GTP “on” state and GTP “off” state. There are three different human Ras isoforms: H-ras, N-ras and K-ras (4A and 4B). Although their sequence identity is very high at the catalytic domain, these isoforms differ in their ability to activate different effectors and hence different signaling pathways. Much of the previous work on this topic has attributed this difference to the hyper variable region of Ras proteins, which contains most of the sequence variance among the isoforms and encodes specificity for differential distribution in the membrane. However, we hypothesize that sequence variation on lobe II of Ras catalytic domain alters dynamics and leads to differential preference for different effectors or modulators. In this work, we used all atom molecular dynamics to analyze the dynamics in the catalytic domain of H-ras and K-ras. We have also analyzed the dynamics of a transforming mutant of H-ras and K-ras and further studied the dynamics of an effectorselective mutant of H-ras. Collectively we have determined that wild type K-ras is more dynamic than H-ras and that the structure of the effector binding loop more closely resembles that of the T35S Raf-selective mutant, possibly giving us a new view and insight into the v mode of effector specificity. Furthermore we have determined that specific mutations at the same location perturb the conformational equilibrium differently in H-ras and K-ras and that an enhanced oncogenic potential may arise from different structural perturbations for each point mutation of a specific isoform.

Conformational selection of inhibitors and substrates by proteolytic enzymes: Implications for drug design and polypeptide processing

Inhibitors of proteolytic enzymes (proteases) are emerging as prospective treatments for diseases such as AIDS and viral infections, cancers, inflammatory disorders, and Alzheimer's disease. Generic approaches to the design of protease inhibitors are limited by the unpredictability of interactions between, and structural changes to, inhibitor and protease during binding. A computer analysis of superimposed crystal structures for 266 small molecule inhibitors bound to 48 proteases (16 aspartic, 17 serine, 8 cysteine, and 7 metallo) provides the first conclusive proof that inhibitors, including substrate analogues, commonly bind in an extended beta-strand conformation at the active sites of all these proteases. Representative superimposed structures are shown for (a) multiple inhibitors bound to a protease of each class, (b) single inhibitors each bound to multiple proteases, and (c) conformationally constrained inhibitors bound to proteases. Thus inhibitor/substrate conformation, rather than sequence/composition alone, influences protease recognition, and this has profound implications for inhibitor design. This conclusion is supported by NMR, CD, and binding studies for HIV-1 protease inhibitors/ substrates which, when preorganized in an extended conformation, have significantly higher protease affinity. Recognition is dependent upon conformational equilibria since helical and turn peptide conformations are not processed by proteases. Conformational selection explains the resistance of folded/structured regions of proteins to proteolytic degradation, the susceptibility of denatured proteins to processing, and the higher affinity of conformationally constrained 'extended' inhibitors/substrates for proteases. Other approaches to extended inhibitor conformations should similarly lead to high-affinity binding to a protease.

Role of phosphorylation in the conformation of tau peptides implicated in Alzheimer's disease

A series of peptides corresponding to isolated regions of Tau (tau) protein have been synthesized and their conformations determined by H-1 NMR spectroscopy. Immunodominant peptides corresponding to tau(224-240) and a bisphosphorylated derivative in which a single Thr and a single Ser are phosphorylated at positions 231 and 235 respectively, and which are recognized by an Alzheimer's disease-specific monoclonal antibody, were the main focus of the study. The nonphosphorylated peptide adopts essentially a random coil conformation in aqueous solution, but becomes slightly more ordered into P-type structure as the hydrophobicity of the solvent is increased by adding up to 50% trifluoroethanol (TFE). Similar trends are observed for the bisphosphorylated peptide, with a somewhat stronger tendency to form an extended structure, There is tentative NMR evidence for a small population of species containing a turn at residues 229-231 in the phosphorylated peptide, and this is strongly supported by CD spectroscopy. A proposal that the selection of a bioactive conformation from a disordered solution ensemble may be an important step (in either tubulin binding or in the formation of PHF) is supported by kinetic data on Pro isomerization. A recent study showed that Thr231 phosphorylation affected the rate of prolyl isomerization and abolished tubulin binding. This binding was restored by the action of the prolyl isomerase Pin1. In the current study, we find evidence for the existence of both trans and cis forms of tau peptides in solution but no difference in the equilibrium distribution of cis-trans isomers upon phosphorylation. Increasing hydrophobicity decreases the prevalence of cis forms and increases the major trans conformation of each of the prolines present in these molecules. We also synthesized mutant peptides containing Tyr substitutions preceding the Pro residues and found that phosphorylation of Tyr appears to have an effect on the equilibrium ratio of cis-trans isomerization and decreases the cis content.

Selection of patients with germline MLH1 mutated Lynch syndrome by determination of MLH1 methylation and BRAF mutation.

Lynch syndrome is one of the most common hereditary colorectal cancer (CRC) syndrome and is caused by germline mutations of MLH1, MSH2 and more rarely MSH6, PMS2, MLH3 genes. Whereas the absence of MSH2 protein is predictive of Lynch syndrome, it is not the case for the absence of MLH1 protein. The purpose of this study was to develop a sensitive and cost effective algorithm to select Lynch syndrome cases among patients with MLH1 immunohistochemical silencing. Eleven sporadic CRC and 16 Lynch syndrome cases with MLH1 protein abnormalities were selected. The BRAF c.1799T> A mutation (p.Val600Glu) was analyzed by direct sequencing after PCR amplification of exon 15. Methylation of MLH1 promoter was determined by Methylation-Sensitive Single-Strand Conformation Analysis. In patients with Lynch syndrome, there was no BRAF mutation and only one case showed MLH1 methylation (6%). In sporadic CRC, all cases were MLH1 methylated (100%) and 8 out of 11 cases carried the above BRAF mutation (73%) whereas only 3 cases were BRAF wild type (27%). We propose the following algorithm: (1) no further molecular analysis should be performed for CRC exhibiting MLH1 methylation and BRAF mutation, and these cases should be considered as sporadic CRC; (2) CRC with unmethylated MLH1 and negative for BRAF mutation should be considered as Lynch syndrome; and (3) only a small fraction of CRC with MLH1 promoter methylation but negative for BRAF mutation should be true Lynch syndrome patients. These potentially Lynch syndrome patients should be offered genetic counselling before searching for MLH1 gene mutations.

Local selection rules that can determine specific pathways of DNA unknotting by type II DNA topoisomerases.

We performed numerical simulations of DNA chains to understand how local geometry of juxtaposed segments in knotted DNA molecules can guide type II DNA topoisomerases to perform very efficient relaxation of DNA knots. We investigated how the various parameters defining the geometry of inter-segmental juxtapositions at sites of inter-segmental passage reactions mediated by type II DNA topoisomerases can affect the topological consequences of these reactions. We confirmed the hypothesis that by recognizing specific geometry of juxtaposed DNA segments in knotted DNA molecules, type II DNA topoisomerases can maintain the steady-state knotting level below the topological equilibrium. In addition, we revealed that a preference for a particular geometry of juxtaposed segments as sites of strand-passage reaction enables type II DNA topoisomerases to select the most efficient pathway of relaxation of complex DNA knots. The analysis of the best selection criteria for efficient relaxation of complex knots revealed that local structures in random configurations of a given knot type statistically behave as analogous local structures in ideal geometric configurations of the corresponding knot type.

Genetic associations of sexual precocity with growth traits and visual scores of conformation, finishing, and muscling in Nelore cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic associations of conformation, finishing precocity and muscling visual scores with mature weight in Nelore cattle

Were estimate (co)variance and genetic associations between conformation, finishing precocity and muscling visual scores measured at weaning (SCW, SFW and SMW, respectively) and yearling (SCY. SFY and SMY, respectively) with mature weight (MW) in Nelore cows, in order to predict the possible changes that inclusion of visual scores in beef cattle selection indices would bring to female mature weight. The data set contained records of 36,757 females, born between 1993 and 2006, belonging to the Jacarezinho cattle raising farm. (Co)variance components were estimated by bivariate animal models using Bayesian inference method through Gibbs sampling, assuming a linear model for MW and a nonlinear (threshold) model for conformation, finishing precocity and muscling visual scores. The first 10,000 rounds were considered as the burn-in period and discarded. The posterior means of direct heritability distributions were: 0.16 +/- 0.02 (SCW); 0.20 +/- 0.02 (SFW); 0.19 +/- 0.02 (SMW); 0.24 +/- 0.02 (SCY); 0.31 +/- 0.02 (SFY); 0.32 +/- 0.02 (SMY) and 0.46 +/- 0.04 (MW). Estimates of genetic correlations between visual scores and MW were positive and moderate, ranging from 0.27 +/- 0.06 to 0.36 +/- 0.04. Visual scores and MW should respond favorably to direct selection. Mature weight can be used in Nelore breeding programs designed to monitor the cows' size. Selection of animals with higher conformation, finishing precocity and muscling scores, especially at yearling, should promote an increase in cows' mature weight. (c) 2010 Elsevier B.V. All rights reserved.

Evaluation of mature cow weight: Genetic correlations with traits used in selection indices, correlated responses, and genetic trends in Nelore cattle

Genetic correlations of selection indices and the traits considered in these indices with mature weight (MW) of Nelore females and correlated responses were estimated to determine whether current selection practices will result in an undesired correlated response in MW. Genetic trends for weaning and yearling indices and MW were also estimated. Data from 612,244 Nelore animals born between 1984 and 2010, belonging to different beef cattle evaluation programs from Brazil and Paraguay, were used. The following traits were studied: weaning conformation (WC), weaning precocity (WP), weaning muscling (WM), yearling conformation (YC), yearling precocity (YP), yearling muscling (YM), weaning and yearling indices, BW gain from birth to weaning (BWG), postweaning BW gain (PWG), scrotal circumference (SC), and MW. The variance and covariance components were estimated by Bayesian inference in a multitrait analysis, including all traits in the same analysis, using a nonlinear (threshold) animal model for visual scores and a linear animal model for the other traits. The mean direct heritabilities were 0.21 ± 0.007 (WC), 0.22 ± 0.007 (WP), 0.20 ± 0.007 (WM), 0.43 ± 0.005 (YC), 0.40 ± 0.005 (YP), 0.40 ± 0.005 (YM), 0.17 ± 0.003 (BWG), 0.21 ± 0.004 (PWG), 0.32 ± 0.001 (SC), and 0.44 ± 0.018 (MW). The genetic correlations between MW and weaning and yearling indices were positive and of medium magnitude (0.30 ± 0.01 and 0.31 ± 0.01, respectively). The genetic changes in weaning index, yearling index, and MW, expressed as units of genetic SD per year, were 0.26, 0.27, and 0.01, respectively. The genetic trend for MW was nonsignificant, suggesting no negative correlated response. The selection practice based on the use of sires with high final index giving preference for those better ranked for yearling precocity and muscling than for conformation generates only a minimal correlated response in MW. © 2013 American Society of Animal Science. All rights reserved.

Body condition score of Nellore beef cows: a heritable measure to improve the selection of reproductive and maternal traits

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Analysis of pedigree and conformation data to explain genetic variability of the horse breed Franches-Montagnes

Franches-Montagnes is the only native horse breed in Switzerland, therefore special efforts should be made for ensuring its survival. The objectives of this study were to characterize the structure of this population as well as genetic variability with pedigree data, conformation traits and molecular markers. Studies were focused to clarify if this population is composed of a heavy- and a light-type subpopulation. Extended pedigree records of 3-year-old stallions (n = 68) and mares (n = 108) were available. Evaluations of body conformation traits as well as pedigree data and molecular markers did not support the two-subpopulation hypothesis. The generation interval ranged from 7.8 to 9.3 years. The complete generation equivalent was high (>12). The number of effective ancestors varied between 18.9 and 20.1, whereof 50% of the genetic variability was attributed to seven of them. Genetic contribution of Warmblood horses ranged from 36% to 42% and that of Coldblood horses from 4% to 6%. The average inbreeding coefficient reached 6%. Inbreeding effective population size was 114.5 when the average increase of the inbreeding coefficient per year since 1910 was taken. Our results suggest that bottleneck situations occurred because of selection of a small number of sire lines. Promotion of planned matings between parents that are less related is recommended in order to avoid a reduction of the genetic diversity.

Computational selection of inhibitors of Abeta aggregation and neuronal toxicity.

Alzheimer's disease (AD) is characterized by the cerebral accumulation of misfolded and aggregated amyloid-beta protein (Abeta). Disease symptoms can be alleviated, in vitro and in vivo, by 'beta-sheet breaker' pentapeptides that reduce plaque load. However the peptide nature of these compounds, made them biologically unstable and unable to penetrate membranes with high efficiency. The main goal of this study was to use computational methods to identify small molecule mimetics with better drug-like properties. For this purpose, the docked conformations of the active peptides were used to identify compounds with similar activities. A series of related beta-sheet breaker peptides were docked to solid state NMR structures of a fibrillar form of Abeta. The lowest energy conformations of the active peptides were used to design three dimensional (3D)-pharmacophores, suitable for screening the NCI database with Unity. Small molecular weight compounds with physicochemical features and a conformation similar to the active peptides were selected, ranked by docking and biochemical parameters. Of 16 diverse compounds selected for experimental screening, 2 prevented and reversed Abeta aggregation at 2-3microM concentration, as measured by Thioflavin T (ThT) fluorescence and ELISA assays. They also prevented the toxic effects of aggregated Abeta on neuroblastoma cells. Their low molecular weight and aqueous solubility makes them promising lead compounds for treating AD.

Evolution of the folding ability of proteins through functional selection

An evolutionary process is simulated with a simple spin-glass-like model of proteins to examine the origin of folding ability. At each generation, sequences are randomly mutated and subjected to a simulation of the folding process based on the model. According to the frequency of local configurations at the active sites, sequences are selected and passed to the next generation. After a few hundred generations, a sequence capable of folding globally into a native conformation emerges. Moreover, the selected sequence has a distinct energy minimum and an anisotropic funnel on the energy surface, which are the imperative features for fast folding of proteins. The proposed model reveals that the functional selection on the local configurations leads a sequence to fold globally into a conformation at a faster rate.

Conformation, energy, and folding ability of selected amino acid sequences.

Evolutionary selection of sequences is studied with a knowledge-based Hamiltonian to find the design principle for folding to a model protein structure. With sequences selected by naive energy minimization, the model structure tends to be unstable and the folding ability is low. Sequences with high folding ability have only the low-lying energy minimum but also an energy landscape which is similar to that found for the native sequence over a wide region of the conformation space. Though there is a large fluctuation in foldable sequences, the hydrophobicity pattern and the glycine locations are preserved among them. Implications of the design principle for the molecular mechanism of folding are discussed.

Directed selection of recombinant human monoclonal antibodies to herpes simplex virus glycoproteins from phage display libraries.

Human monoclonal antibodies have considerable potential in the prophylaxis and treatment of viral disease. However, only a few such antibodies suitable for clinical use have been produced to date. We have previously shown that large panels of human recombinant monoclonal antibodies against a plethora of infectious agents, including herpes simplex virus types 1 and 2, can be established from phage display libraries. Here we demonstrate that facile cloning of recombinant Fab fragments against specific viral proteins in their native conformation can be accomplished by panning phage display libraries against viral glycoproteins "captured" from infected cell extracts by specific monoclonal antibodies immobilized on ELISA plates. We have tested this strategy by isolating six neutralizing recombinant antibodies specific for herpes simplex glycoprotein gD or gB, some of which are against conformationally sensitive epitopes. By using defined monoclonal antibodies for the antigen-capture step, this method can be used for the isolation of antibodies to specific regions and epitopes within the target viral protein. For instance, monoclonal antibodies to a nonneutralizing epitope can be used in the capture step to clone antibodies to neutralizing epitopes, or antibodies to a neutralizing epitope can be used to clone antibodies to a different neutralizing epitope. Furthermore, by using capturing antibodies to more immunodominant epitopes, one can direct the cloning to less immunogenic ones. This method should be of value in generating antibodies to be used both in the prophylaxis and treatment of viral infections and in the characterization of the mechanisms of antibody protective actions at the molecular level.

CTL recognition of a bulged viral peptide involves biased TCR selection

MHC class I molecules generally present peptides of 8-10 aa long, forming an extended coil in the HLA cleft. Although longer peptides can also bind to class I molecules, they tend to bulge from the cleft and it is not known whether the TCR repertoire has sufficient plasticity to recognize these determinants during the antiviral CTL response. In this study, we show that unrelated individuals infected with EBV generate a significant CTL response directed toward an HLA-B*3501-restricted, 11-mer epitope from the BZLF1 Ag. The 11-mer determinant adopts a highly bulged conformation with seven of the peptide side chains being solvent-exposed and available for TCR interaction. Such a complex potentially creates a structural challenge for TCR corecognition of both HLA-B*3501 and the peptide Ag. Surprisingly, unrelated B*3501 donors recognizing the 11-mer use identical or closely related alpha beta TCR sequences that share particular CDR3 motifs. Within the small number of dominant CTL clonotypes observed, each has discrete fine specificity for the exposed-side chain residues of the peptide. The data show that bulged viral peptides are indeed immunogenic but suggest that the highly constrained TCR repertoire reflects a limit to TCR diversity when responding to some unusual MHC peptide ligands.