Protein Folding Simulations: Confinement, External Fields and Sequence Design

The present Thesis looks at the problem of protein folding using Monte Carlo and Langevin simulations, three topics in protein folding have been studied: 1) the effect of confining potential barriers, 2) the effect of a static external field and 3) the design of amino acid sequences which fold in a short time and which have a stable native state (global minimum). Regarding the first topic, we studied the confinement of a small protein of 16 amino acids known as 1NJ0 (PDB code) which has a beta-sheet structure as a native state. The confinement of proteins occurs frequently in the cell environment. Some molecules called Chaperones, present in the cytoplasm, capture the unfolded proteins in their interior and avoid the formation of aggregates and misfolded proteins. This mechanism of confinement mediated by Chaperones is not yet well understood. In the present work we considered two kinds of potential barriers which try to mimic the confinement induced by a Chaperon molecule. The first kind of potential was a purely repulsive barrier whose only effect is to create a cavity where the protein folds up correctly. The second kind of potential was a barrier which includes both attractive and repulsive effects. We performed Wang-Landau simulations to calculate the thermodynamical properties of 1NJ0. From the free energy landscape plot we found that 1NJ0 has two intermediate states in the bulk (without confinement) which are clearly separated from the native and the unfolded states. For the case of the purely repulsive barrier we found that the intermediate states get closer to each other in the free energy landscape plot and eventually they collapse into a single intermediate state. The unfolded state is more compact, compared to that in the bulk, as the size of the barrier decreases. For an attractive barrier modifications of the states (native, unfolded and intermediates) are observed depending on the degree of attraction between the protein and the walls of the barrier. The strength of the attraction is measured by the parameter $\epsilon$. A purely repulsive barrier is obtained for $\epsilon=0$ and a purely attractive barrier for $\epsilon=1$. The states are changed slightly for magnitudes of the attraction up to $\epsilon=0.4$. The disappearance of the intermediate states of 1NJ0 is already observed for $\epsilon =0.6$. A very high attractive barrier ($\epsilon \sim 1.0$) produces a completely denatured state. In the second topic of this Thesis we dealt with the interaction of a protein with an external electric field. We demonstrated by means of computer simulations, specifically by using the Wang-Landau algorithm, that the folded, unfolded, and intermediate states can be modified by means of a field. We have found that an external field can induce several modifications in the thermodynamics of these states: for relatively low magnitudes of the field ($<2.06 \times 10^8$ V/m) no major changes in the states are observed. However, for higher magnitudes than ($6.19 \times 10^8$ V/m) one observes the appearance of a new native state which exhibits a helix-like structure. In contrast, the original native state is a $\beta$-sheet structure. In the new native state all the dipoles in the backbone structure are aligned parallel to the field. The design of amino acid sequences constitutes the third topic of the present work. We have tested the Rate of Convergence criterion proposed by D. Gridnev and M. Garcia ({\it work unpublished}). We applied it to the study of off-lattice models. The Rate of Convergence criterion is used to decide if a certain sequence will fold up correctly within a relatively short time. Before the present work, the common way to decide if a certain sequence was a good/bad folder was by performing the whole dynamics until the sequence got its native state (if it existed), or by studying the curvature of the potential energy surface. There are some difficulties in the last two approaches. In the first approach, performing the complete dynamics for hundreds of sequences is a rather challenging task because of the CPU time needed. In the second approach, calculating the curvature of the potential energy surface is possible only for very smooth surfaces. The Rate of Convergence criterion seems to avoid the previous difficulties. With this criterion one does not need to perform the complete dynamics to find the good and bad sequences. Also, the criterion does not depend on the kind of force field used and therefore it can be used even for very rugged energy surfaces.

Loss of Anticodon Wobble Uridine Modifications Affects tRNALys Function and Protein Levels in Saccharomyces cerevisiae

In eukaryotes, wobble uridines in the anticodons of tRNALysUUU, tRNAGluUUC and tRNAGlnUUG are modified to 5-methoxy-carbonyl-methyl-2-thio-uridine (mcm5s2U). While mutations in subunits of the Elongator complex (Elp1-Elp6), which disable mcm5 side chain formation, or removal of components of the thiolation pathway (Ncs2/Ncs6, Urm1, Uba4) are individually tolerated, the combination of both modification defects has been reported to have lethal effects on Saccharomyces cerevisiae. Contrary to such absolute requirement of mcm5s2U for viability, we demonstrate here that in the S. cerevisiae S288C-derived background, both pathways can be simultaneously inactivated, resulting in combined loss of tRNA anticodon modifications (mcm5U and s2U) without a lethal effect. However, an elp3 disruption strain displays synthetic sick interaction and synergistic temperature sensitivity when combined with either uba4 or urm1 mutations, suggesting major translational defects in the absence of mcm5s2U modifications. Consistent with this notion, we find cellular protein levels drastically decreased in an elp3uba4 double mutant and show that this effect as well as growth phenotypes can be partially rescued by excess of tRNALysUUU. These results may indicate a global translational or protein homeostasis defect in cells simultaneously lacking mcm5 and s2 wobble uridine modification that could account for growth impairment and mainly originates from tRNALysUUU hypomodification and malfunction.

"Genetic Diversity and Selection in Three Plasmodium vivax Merozoite Surface Protein 7 (Pvmsp-7) Genes in a Colombian Population"

A completely effective vaccine for malaria (one of the major infectious diseases worldwide) is not yet available; different membrane proteins involved in parasite-host interactions have been proposed as candidates for designing it. It has been found that proteins encoded by the merozoite surface protein (msp)-7 multigene family are antibody targets in natural infection; the nucleotide diversity of three Pvmsp-7 genes was thus analyzed in a Colombian parasite population. By contrast with P. falciparum msp-7 loci and ancestral P. vivax msp-7 genes, specie-specific duplicates of the latter specie display high genetic variability, generated by single nucleotide polymorphisms, repeat regions, and recombination. At least three major allele types are present in Pvmsp-7C, Pvmsp-7H and Pvmsp-7I and positive selection seems to be operating on the central region of these msp-7 genes. Although this region has high genetic polymorphism, the C-terminus (Pfam domain ID: PF12948) is conserved and could be an important candidate when designing a subunit-based antimalarial vaccine.

Genetic diversity and selection in three plasmodium vivax merozoite surface protein 7 (pvmsp-7) genes in a colombian population

A completely effective vaccine for malaria (one of the major infectious diseases worldwide) is not yet available; different membrane proteins involved in parasite-host interactions have been proposed as candidates for designing it. It has been found that proteins encoded by the merozoite surface protein (msp)-7 multigene family are antibody targets in natural infection; the nucleotide diversity of three Pvmsp-7 genes was thus analyzed in a Colombian parasite population. By contrast with P. falciparum msp-7 loci and ancestral P. vivax msp-7 genes, specie-specific duplicates of the latter specie display high genetic variability, generated by single nucleotide polymorphisms, repeat regions, and recombination. At least three major allele types are present in Pvmsp-7C, Pvmsp-7H and Pvmsp-7I and positive selection seems to be operating on the central region of these msp-7 genes. Although this region has high genetic polymorphism, the C-terminus (Pfam domain ID: PF12948) is conserved and could be an important candidate when designing a subunit-based antimalarial vaccine.

Identification, characterization and antigenicity of the Plasmodium vivax rhoptry neck protein 1 (PvRON1)

Background: Plasmodium vivax malaria remains a major health problem in tropical and sub-tropical regions worldwide. Several rhoptry proteins which are important for interaction with and/or invasion of red blood cells, such as PfRONs, Pf92, Pf38, Pf12 and Pf34, have been described during the last few years and are being considered as potential anti-malarial vaccine candidates. This study describes the identification and characterization of the P. vivax rhoptry neck protein 1 (PvRON1) and examine its antigenicity in natural P. vivax infections. Methods: The PvRON1 encoding gene, which is homologous to that encoding the P. falciparum apical sushi protein (ASP) according to the plasmoDB database, was selected as our study target. The pvron1 gene transcription was evaluated by RT-PCR using RNA obtained from the P. vivax VCG-1 strain. Two peptides derived from the deduced P. vivax Sal-I PvRON1 sequence were synthesized and inoculated in rabbits for obtaining anti-PvRON1 antibodies which were used to confirm the protein expression in VCG-1 strain schizonts along with its association with detergent-resistant microdomains (DRMs) by Western blot, and its localization by immunofluorescence assays. The antigenicity of the PvRON1 protein was assessed using human sera from individuals previously exposed to P. vivax malaria by ELISA. Results: In the P. vivax VCG-1 strain, RON1 is a 764 amino acid-long protein. In silico analysis has revealed that PvRON1 shares essential characteristics with different antigens involved in invasion, such as the presence of a secretory signal, a GPI-anchor sequence and a putative sushi domain. The PvRON1 protein is expressed in parasite's schizont stage, localized in rhoptry necks and it is associated with DRMs. Recombinant protein recognition by human sera indicates that this antigen can trigger an immune response during a natural infection with P. vivax. Conclusions: This study shows the identification and characterization of the P. vivax rhoptry neck protein 1 in the VCG-1 strain. Taking into account that PvRON1 shares several important characteristics with other Plasmodium antigens that play a functional role during RBC invasion and, as shown here, it is antigenic, it could be considered as a good vaccine candidate. Further studies aimed at assessing its immunogenicity and protection-inducing ability in the Aotus monkey model are thus recommended.

The four major N- and C-terminal splice variants of the excitatory amino acid transporter GLT-1 form cell surface homomeric and heteromeric assemblies

The L-glutamate transporter GLT-1 is an abundant CNS membrane protein of the excitatory amino acid transporter (EAAT) family which controls extracellular L-glutamate levels and is important in limiting excitotoxic neuronal death. Using RT-PCR, we have determined that four mRNAs encoding GLT-1 exist in mouse brain, with the potential to encode four GLT-1 isoforms that differ in their N- and C-termini. We expressed all four isoforms (termed MAST-KREK, MPK-KREK, MAST-DIETCI and MPK-DIETCI according to amino acid sequence) in a range of cell lines and primary astrocytes and show that each isoform can reach the cell surface. In transfected HEK-293 or COS-7 cells, all four isoforms support high-affinity sodium-dependent L-glutamate uptake with identical pharmacological and kinetic properties. Inserting a viral epitope (V5, HA or FLAG) into the second extracellular domain of each isoform allowed co-immunoprecipitation and tr-FRET studies using transfected HEK-293 cells. Here we show for the first time that each of the four isoforms are able to combine to form homomeric and heteromeric assemblies, each of which are expressed at the cell surface of primary astrocytes. After activation of protein kinase C by phorbol ester, V5-tagged GLT-1 is rapidly removed from the cell surface of HEK-293 cells and degraded. This study provides direct biochemical evidence for oligomeric assembly of GLT-1 and reports the development of novel tools to provide insight into the trafficking of GLT-1.

Iron-regulated surface determinant protein a mediates adhesion of staphylococcus aureus to human corneocyte envelope proteins

The ability of Staphylococcus aureus to colonize the human nares is a crucial prerequisite for disease. IsdA is a major S. aureus surface protein that is expressed during human infection and required for nasal colonization and survival on human skin. In this work, we show that IsdA binds to involucrin, loricrin, and cytokeratin K10, proteins that are present in the cornified envelope of human desquamated epithelial cells. To measure the forces and dynamics of the interaction between IsdA and loricrin (the most abundant protein of the cornified envelope), single-molecule force spectroscopy was used, demonstrating high-specificity binding. IsdA acts as a cellular adhesin to the human ligands, promoting whole-cell binding to immobilized proteins, even in the absence of other S. aureus components (as shown by heterologous expression in Lactococcus lactis). Inhibition experiments revealed the binding of the human ligands to the same IsdA region. This region was mapped to the NEAT domain of IsdA. The NEAT domain also was found to be required for S. aureus whole-cell binding to the ligands as well as to human nasal cells. Thus, IsdA is an important adhesin to human ligands, which predominate in its primary ecological niche.

Contribution of redox status to hepatitis C virus E2 envelope protein function and antigenicity

Disulfide bonding contributes to the function and antigenicity of many viral envelope glycoproteins. We assessed here its significance for the hepatitis C virus E2 envelope protein and a counterpart deleted for hypervariable region-1 (HVR1). All 18 cysteine residues of the antigens were involved in disulfides. Chemical reduction of up to half of these disulfides was compatible with anti-E2 monoclonal antibody reaction, CD81 receptor binding, and viral entry, whereas complete reduction abrogated these properties. The addition of 5,5'-dithiobis-2-nitrobenzoic acid had no effect on viral entry. Thus, E2 function is only weakly dependent on its redox status, and cell entry does not require redox catalysts, in contrast to a number of enveloped viruses. Because E2 is a major neutralizing antibody target, we examined the effect of disulfide bonding on E2 antigenicity. We show that reduction of three disulfides, as well as deletion of HVR1, improved antibody binding for half of the patient sera tested, whereas it had no effect on the remainder. Small scale immunization of mice with reduced E2 antigens greatly improved serum reactivity with reduced forms of E2 when compared with immunization using native E2, whereas deletion of HVR1 only marginally affected the ability of the serum to bind the redox intermediates. Immunization with reduced E2 also showed an improved neutralizing antibody response, suggesting that potential epitopes are masked on the disulfide-bonded antigen and that mild reduction may increase the breadth of the antibody response. Although E2 function is surprisingly independent of its redox status, its disulfide bonds mask antigenic domains. E2 redox manipulation may contribute to improved vaccine design.

Fas ligand-induced apoptosis is regulated by nitric oxide through the inhibition of fas receptor clustering and the nitrosylation of protein kinase Cepsilon

Apoptosis induced by the death-inducing ligand FasL (CD95L) is a major mechanism of cell death. Trophoblast cells express the Fas receptor yet survive in an environment that is rich in the ligand. We report that basal nitric oxide (NO) production is responsible for the resistance of trophoblasts to FasL-induced apoptosis. In this study we demonstrate that basal NO production resulted in the inhibition of receptor clustering following ligand binding. In addition NO also protected cells through the selective nitrosylation, and inhibition, of protein kinase Cepsilon (PKCepsilon) but not PKCalpha. In the absence of NO production PKCepsilon interacted with, and phosphorylated, the anti-apoptotic protein cFLIP. The interaction is predominantly with the short form of cFLIP and its phosphorylation reduces its recruitment to the death-inducing signaling complex (DISC) that is formed following binding of a death-inducing ligand to its receptor. Inhibition of cFLIP recruitment to the DISC leads to increased activation of caspase 8 and subsequently to apoptosis. Inhibition of PKCepsilon using siRNA significantly reversed the sensitivity to apoptosis induced by inhibition of NO synthesis suggesting that NO-mediated inhibition of PKCepsilon plays an important role in the regulation of Fas-induced apoptosis.

Analysis of protein-protein interactions in the feline calicivirus replication complex

Caliciviruses are a major cause of gastroenteritis in humans and cause a wide variety of other diseases in animals. Here, the characterization of protein-protein interactions between the individual proteins of Feline calicivirus (FCV), a model system for other members of the family Caliciviridae, is reported. Using the yeast two-hybrid system combined with a number of other approaches, it is demonstrated that the p32 protein (the picornavirus 2B analogue) of FCV interacts with p39 (2C), p30 (3A) and p76 (3CD). The FCV protease/RNA polymerase (ProPol) p76 was found to form homo-oligomers, as well as to interact with VPg and ORF2, the region encoding the major capsid protein VP1. A weak interaction was also observed between p76 and the minor capsid protein encoded by ORF3 (VP2). ORF2 protein was found to interact with VPg, p76 and VP2. The potential roles of the interactions in calicivirus replication are discussed.

Influence of mutations affecting gonadotropin production or responsiveness on expression of inhibin subunit mRNA and protein in the mouse ovary

Measurement of inhibins A and B in the serum of normal cyclic rodents has implicated FSH in the regulation of these peptides within the ovary. To extend these observations we have used a panel of mutant mice carrying mutations which affect either the production of, or the ability to respond to, FSH and LH. As a consequence, the females are infertile and show different degrees of follicular development. The aim of this study was to measure inhibin gene transcription in the ovaries of these mutant females together with inhibin protein levels in ovaries and serum and to relate these to follicular development within the ovary. Comparison was made with a pool of normal/heterozygous females. In hpg females where lack of GnRH production results in the absence of gonadotropin synthesis, in FSHbeta knockout (FSHbetaKO) females where disruption of the gene encoding FSHbeta results in the absence of FSH production, and in FSH receptor knockout (FSHRKO) females which are unable to respond to circulating FSH, follicular development remains at the pre-antral stage in these three mutants. Only in the hpg females were common inhibin alpha subunit mRNA levels significantly lower than normal. In these three mutants, however, mRNA levels for both the betaA and betaB subunits were extremely low compared with normal mice. At the protein level, neither inhibin A nor B was detected in the serum of these three mutants; however inhibin B, albeit at very low levels, was detectable within the ovaries. These observations confirm a major role for FSH in the control of transcription of the RA and betaB genes but suggest that the constitutive transcription of the alpha subunit is less dependent on FSH. In contrast, in LH receptor knockout (LuRKO) female mice inhibin betaA subunit mRNA levels were similar to those measured in normal/heterozygous females but levels of inhibin alpha and betaB subunit mRNAs were significantly higher than in the normal group. This was reflected in significantly higher inhibin B protein levels in ovaries and serum. An inability to respond to LH combined with high circulating levels of FSH leads to a high proportion of antral follicles in LuRKO females, with granulosa cells constituting the major cell type within the ovary. The high percentage of antral granulosa cells is likely to account for the significantly higher levels of inhibin B production in these ovaries.

Cleavage and serum reactivity of the severe acute respiratory syndrome coronavirus spike protein

Severe acute respiratory syndrome (SARS) coronavirus (SCoV) spike (S) protein is the major surface antigen of the virus and is responsible for receptor binding and the generation of neutralizing antibody. To investigate SCoV S protein, full-length and individual domains of S protein were expressed on the surface of insect cells and were characterized for cleavability and reactivity with serum samples obtained from patients during the convalescent phase of SARS. S protein could be cleaved by exogenous trypsin but not by coexpressed furin, suggesting that the protein is not normally processed during infection. Reactivity was evident by both flow cytometry and Western blot assays, but the pattern of reactivity varied according to assay and sequence of the antigen. The antibody response to SCoV S protein involves antibodies to both linear and conformational epitopes, with linear epitopes associated with the carboxyl domain and conformational epitopes associated with the amino terminal domain. Recombinant SCoV S protein appears to be a suitable antigen for the development of an efficient and sensitive diagnostic test for SARS, but our data suggest that assay format and choice of S antigen are important considerations.

Stress relaxation behavior of wheat dough, gluten, and gluten protein fractions

Relaxation behavior was measured for dough, gluten and gluten protein fractions obtained from the U.K. biscuitmaking flour, Riband, and the U.K. breadmaking flour, Hereward. The relaxation spectrum, in which relaxation times (tau) are related to polymer molecular size, for dough showed a broad molecular size distribution, with two relaxation processes: a major peak at short times and a second peak at times longer than 10 sec, which is thought to correspond to network structure, and which may be attributed to entanglements and physical cross-links of polymers. Relaxation spectra of glutens were similar to those for the corresponding doughs from both flours. Hereward gluten clearly showed a much more pronounced second peak in relaxation spectrum and higher relaxation modulus than Riband gluten at the same water content. In the gluten protein fractions, gliadin and acetic acid soluble glutenin only showed the first relaxation process, but gel protein clearly showed both the first and second relaxation processes. The results show that the relaxation properties of dough depend on its gluten protein and that gel protein is responsible for the network structure for dough and gluten.

Selective separation of the major whey proteins using ion exchange membranes

Synthetic microporous membranes with functional groups covalently attached were used to selectively separate beta-lactoglobulin, BSA, and alpha-lactalbumin from rennet whey. The selectivity and membrane performance of strong (quaternary ammonium) and weak (diethylamine) ion-exchange membranes were studied using breakthrough curves, measurement of binding capacity, and protein composition of the elution fraction to determine the binding behavior of each membrane. When the weak and strong anion exchange membranes were saturated with whey, they were both selective primarily for beta-lactoglobulin with less than 1% of the eluate consisting of alpha-lactalbumin or BSA. The binding capacity of a pure alpha-lactoglobulin solution was in excess of 1.5 mg/cm(2) of membrane. This binding capacity was reduced to approximately 1.2 mg/cm(2) when using a rennet whey solution (pH 6.4). This reduction in protein binding capacity can be explained by both the competitive effects of other whey proteins and the effect of ions present in whey. Using binary solution breakthrough curves and rennet whey breakthrough curves, it was shown that alpha-lactalbumin and BSA were displaced from the strong and weak anion exchange membranes by beta-lactoglobulin. Finally, the effect of ionic strength on the binding capacity of individual proteins for each membrane was determined by comparing model protein solutions in milk permeate (pH 6.4) and a 10 mM sodium phosphate buffer (pH 6.4). Binding capacities of beta-lactoglobulin, alpha-lactalbumin, and BSA in milk permeate were reduced by as much as 50%. This reduction in capacity coupled with the low binding capacity of current ion exchange membranes are 2 serious considerations for selectively separating complex and concentrated protein solutions.

Localization to the nucleolus is a common feature of coronavirus nucleoproteins and the protein may disrupt host cell division

The subcellular localization of transmissible gastroenteritis virus (TGEV) and mouse hepatitis virus (MHV) (group I and group II coronaviruses, respectively) nucleoproteins (N proteins) were examined by confocal microscopy. The proteins were shown to localize either to the cytoplasm alone or to the cytoplasm and a structure in the nucleus. This feature was confirmed to be the nucleolus by using specific antibodies to nucleolin, a major component of the nucleolus, and by confocal microscopy to image sections through a cell expressing N protein. These findings are consistent with our previous report for infectious bronchitis virus (group III coronavirus) (J. A. Hiscox et al., J. Virol. 75:506-512, 2001), indicating that nucleolar localization of the N protein is a common feature of the coronavirus family and is possibly of functional significance. Nucleolar localization signals were identified in the domain III region of the N protein from all three coronavirus groups, and this suggested that transport of N protein to the nucleus might be an active process. In addition, our results suggest that the N protein might function to disrupt cell division. Thus, we observed that approximately 30% of cells transfected with the N protein appeared to be undergoing cell division. The most likely explanation for this is that the N protein induced a cell cycle delay or arrest, most likely in the G2/M phase. In a fraction of transfected cells expressing coronavirus N proteins, we observed multinucleate cells and dividing cells with nucleoli (which are only present during interphase). These findings are consistent with the possible inhibition of cytokinesis in these cells.