In the late 1990s, a series of landmark publications described RNA interference (RNAi) and related RNA silencing phenomena in nematodes, plants and fungi. By manipulating RNA silencing, biologists have been able to create tools for specifically inactivating genes. In organisms from trypanosomes to insects, RNA silencing is now indispensible for studying gene function. RNA silencing has been used in a project aimed at systematically knocking out all genes in the model plant Arabidopsis thaliana. RNA silencing has a natural role in defending eukaryotic cells against virus replication. By assembling virus DNA sequences in a form that triggers RNA silencing, biologists have created plants resistant to specific viruses. In this study, we set out to test if a similar approach would protect plants against infection by the agriculturally important Brassica pathogen, Plasmodiophora brassicae. P. brassicae is an obligate intracellular biotroph, from the little studied eukaryotic supergroup, the Rhizaria. To identify the gene sequences that would be starting material for P. brassicae RNA silencing, new P. brassicae genes were gathered by cDNA cloning or genomic PCR-walking. Using suppression subtractive hybridisation (SSH) and oligo-capping cloning of full-length cDNAs, 76 new gene sequences were identified. A large proportion of the cDNAs were predicted to contain signal peptides for ER translocation. In addition to the new cDNA identified here, partial sequences for the P. brassicae actin and TPS genes were published by other researchers close to the beginning of this study. Using PCR-walking, full-length genomic DNA sequences from both genes were obtained. Later, genomic DNA sequences spanning or flanking a total of 24 P. brassicae genes were obtained. The P. brassicae genes were rich in typical eukaryotic spliceosomal introns. Transcription of P. brassicae genes also appears likely to begin from initiator elements rather than TATA-box-containing promoters. A segment of the P. brassicae actin gene was assembled in hairpin format and transformed into Arabidopsis thaliana. Observation of simultaneous knockdown of the GUS marker gene as well as detection of siRNAs indicated that the hpRNA sequences induced RNA silencing. However, inoculation of these plants with P. brassicae resulted in heavy club root infection. We were unable to detect decreases in actin gene expression in the infecting P. brassicae, at either early or late stages of infection. We conclude that, within the limits of the techniques used here, there is no evidence for induction of RNA silencing in P. brassicae by in planta produced siRNAs.

In the late 1990s, a series of landmark publications described RNA interference (RNAi) and related RNA silencing phenomena in nematodes, plants and fungi. By manipulating RNA silencing, biologists have been able to create tools for specifically inactivating genes. In organisms from trypanosomes to insects, RNA silencing is now indispensible for studying gene function. RNA silencing has been used in a project aimed at systematically knocking out all genes in the model plant Arabidopsis thaliana. RNA silencing has a natural role in defending eukaryotic cells against virus replication. By assembling virus DNA sequences in a form that triggers RNA silencing, biologists have created plants resistant to specific viruses. In this study, we set out to test if a similar approach would protect plants against infection by the agriculturally important Brassica pathogen, Plasmodiophora brassicae. P. brassicae is an obligate intracellular biotroph, from the little studied eukaryotic supergroup, the Rhizaria. To identify the gene sequences that would be starting material for P. brassicae RNA silencing, new P. brassicae genes were gathered by cDNA cloning or genomic PCR-walking. Using suppression subtractive hybridisation (SSH) and oligo-capping cloning of full-length cDNAs, 76 new gene sequences were identified. A large proportion of the cDNAs were predicted to contain signal peptides for ER translocation. In addition to the new cDNA identified here, partial sequences for the P. brassicae actin and TPS genes were published by other researchers close to the beginning of this study. Using PCR-walking, full-length genomic DNA sequences from both genes were obtained. Later, genomic DNA sequences spanning or flanking a total of 24 P. brassicae genes were obtained. The P. brassicae genes were rich in typical eukaryotic spliceosomal introns. Transcription of P. brassicae genes also appears likely to begin from initiator elements rather than TATA-box-containing promoters. A segment of the P. brassicae actin gene was assembled in hairpin format and transformed into Arabidopsis thaliana. Observation of simultaneous knockdown of the GUS marker gene as well as detection of siRNAs indicated that the hpRNA sequences induced RNA silencing. However, inoculation of these plants with P. brassicae resulted in heavy club root infection. We were unable to detect decreases in actin gene expression in the infecting P. brassicae, at either early or late stages of infection. We conclude that, within the limits of the techniques used here, there is no evidence for induction of RNA silencing in P. brassicae by in planta produced siRNAs.

In the late 1990s, a series of landmark publications described RNA interference (RNAi) and related RNA silencing phenomena in nematodes, plants and fungi. By manipulating RNA silencing, biologists have been able to create tools for specifically inactivating genes. In organisms from trypanosomes to insects, RNA silencing is now indispensible for studying gene function. RNA silencing has been used in a project aimed at systematically knocking out all genes in the model plant Arabidopsis thaliana. RNA silencing has a natural role in defending eukaryotic cells against virus replication. By assembling virus DNA sequences in a form that triggers RNA silencing, biologists have created plants resistant to specific viruses. In this study, we set out to test if a similar approach would protect plants against infection by the agriculturally important Brassica pathogen, Plasmodiophora brassicae. P. brassicae is an obligate intracellular biotroph, from the little studied eukaryotic supergroup, the Rhizaria. To identify the gene sequences that would be starting material for P. brassicae RNA silencing, new P. brassicae genes were gathered by cDNA cloning or genomic PCR-walking. Using suppression subtractive hybridisation (SSH) and oligo-capping cloning of full-length cDNAs, 76 new gene sequences were identified. A large proportion of the cDNAs were predicted to contain signal peptides for ER translocation. In addition to the new cDNA identified here, partial sequences for the P. brassicae actin and TPS genes were published by other researchers close to the beginning of this study. Using PCR-walking, full-length genomic DNA sequences from both genes were obtained. Later, genomic DNA sequences spanning or flanking a total of 24 P. brassicae genes were obtained. The P. brassicae genes were rich in typical eukaryotic spliceosomal introns. Transcription of P. brassicae genes also appears likely to begin from initiator elements rather than TATA-box-containing promoters. A segment of the P. brassicae actin gene was assembled in hairpin format and transformed into Arabidopsis thaliana. Observation of simultaneous knockdown of the GUS marker gene as well as detection of siRNAs indicated that the hpRNA sequences induced RNA silencing. However, inoculation of these plants with P. brassicae resulted in heavy club root infection. We were unable to detect decreases in actin gene expression in the infecting P. brassicae, at either early or late stages of infection. We conclude that, within the limits of the techniques used here, there is no evidence for induction of RNA silencing in P. brassicae by in planta produced siRNAs.

Le virus des nervures jaunes et nécrotiques de la betterave (Beet necrotic yellow vein virus, BNYVV) est l’agent infectieux responsable de la rhizomanie de la betterave sucrière, une maladie caractérisée par une prolifération anarchique du chevelu racinaire. Le Beet soil-borne mosaic virus (BSBMV) appartient également au genre Benyvirus mais n’est retrouvé qu’en Amérique du Nord. Ce virus, identifié pour la première fois au Texas, est morphologiquement et génétiquement semblable au BNYVV mais sérologiquement éloigné. Compte tenu des différences moléculaires existant, le BSBMV et BNYVV correspondent à deux espèces virales distinctes. Mon projet de thèse a consisté à étudier les interactions moléculaires entre le BNYVV et le BSBMV et rechercher les mécanismes impliqués dans la pathogénicité de ces deux virus. Des clones complets cDNA infectieux du BNYVV étaient disponibles, tout comme ceux de BSBMV. Compte tenu de l’aspect versatile de l’obtention de transcrits infectieux de ces différents clones, j’ai entrepris de produire des clones cDNA de chacun des ARN viraux sous contrôle d’un promoteur constitutive végétal pour initier l’infection par agroinfiltration. Les plantes hôtes Chenopodium quinoa et Nicotiana benthamiana ont été inoculées par des transcrits et agroinfiltrées pour initier l’infection virale et étudier l’interaction entre les ARN génomiques 1 et 2 des deux virus et étudier les propriétés de constructions chimères. En parallèle à ce travail, j’ai réalisé la caractérisation du suppresseur de RNA silencing du BSBMV en le comparant à celui du BNYVV. / The genus Benyvirus includes the most important and widespread sugar beet viruses transmitted through the soil by the plasmodiophorid Polymyxa betae. In particular Beet necrotic yellow vein virus (BNYVV), the leading infectious agent that affects sugar beet, causes an abnormal rootlet proliferation known as rhizomania. Beet soil-borne mosaic virus (BSBMV) is widely distributed in the United States and, up to date has not been reported in others countries. My PhD project aims to investigate molecular interactions between BNYVV and BSBMV and the mechanisms involved in the pathogenesis of these viruses.BNYVV full-length infectious cDNA clones were available as well as full-length cDNA clones of BSBMV RNA-1, -2, -3 and -4. Handling of these cDNA clones in order to produce in vitro infectious transcripts need sensitive and expensive steps, so Ideveloped agroclones of BNYVV and BSBMV RNAs, as well as viral replicons allowing the expression of different proteins.Chenopodium quinoa and Nicotiana benthamiana plants have been infected with in vitro transcripts and agroclones to investigate the interaction between BNYVV and BSBMV RNA-1 and -2 and the behavior of artificial viral chimeras. Simultaneously I characterized BSBMV p14 and demonstrated that it is a suppressor of posttranscriptional gene silencing sharing common features with BNYVV p14.

El virus de la tristeza de los cítricos (Citrus tristeza virus; CTV) es el agente causal de unas de las enfermedades virales de los árboles cítricos más devastadoras en el mundo. CTV está restringido al floema en su huésped cítrico natural, y ha desarrollado tres proteínas supresoras de silenciamiento que actúan a nivel intra-(p23 y p20) e intercelular (p20 y p25) para superar la fuerte defensa antiviral del huésped. La interferencia de RNA, una aproximación basada en el uso de dsRNA para desencadenar el silenciamiento de RNA, ha sido utilizada ampliamente para generar plantas transgénicas resistentes a virus. Considerando el importante papel de p23, p20 y p25 en la patogénesis de CTV, hemos transformado plantas de lima Mexicana con un vector intrón-horquilla que porta la secuencia completa en versión no traducible de los genes p25, p20, p23 y el extremo 3¿-UTR de la cepa T36 de CTV, para intentar silenciar su expresión en células infectadas. Se ha observado resistencia completa a la infección viral en tres líneas transgénicas, manteniéndose todas sus propagaciones asintomáticas y libres de virus tras ser inoculadas mediante injerto con CTV-T36, tanto en el portainjertos no transgénico como directamente sobre la variedad transgénica. La acumulación de siRNA derivados del transgén fue necesaria pero no suficiente para lograr resistencia frente a CTV en las plantas. Al inocular propagaciones de las líneas transgénicas inmunes con una cepa de CTV divergente, la resistencia fue parcialmente superada, destacando la importancia de la identidad de secuencia en el mecanismo subyacente a la interferencia de RNA. Este trabajo es el primero en que se consigue resistencia completa a CTV en un huésped cítrico muy sensible, actuando simultáneamente sobre los tres supresores virales de silenciamiento mediante interferencia de RNA. La proteína p23 codificada por el virus es además un importante factor de patogenicidad. La expresión ectópica de p23 en plantas de cítricos induce aberraciones fenológicas semejantes a síntomas de CTV. Para estudiar en más detalle el papel de p23 en la patogénesis de CTV, se ha sobre-expresado en lima Mexicana el gen p23 de CTV T36 y tres versiones truncadas del mismo bajo el control del promotor 35S del virus del mosaico de la coliflor (Cauliflower mosaic virus). Solo la versión truncada, que expresa los aminoácidos del 1 al 157 (p23-¿157) indujo síntomas similares a los producidos por CTV, aunque más suaves que los inducidos por la expresión de la proteína p23 entera (209 aminoácidos), permitiendo delimitar la región responsable de la patogénesis de p23 en cítricos a un fragmento de 157 aminoácidos que incluye el dedo de zinc y los motivos básicos flanqueantes de la proteína. La actividad de p23 como supresor de silenciamiento de RNA en N. benthamiana se perdía en todos los mutantes de p23 probados, lo cual indica que la supresión de silenciamiento implica a la mayoría de las regiones de la proteína. Para profundizar más en el papel de p23 en la patogénesis, en un siguiente paso hemos restringido la expresión de transgenes derivados de p23 a células asociadas al floema de lima Mexicana mediante el uso del promotor especifico de floema del virus del moteado amarillo de la comelina (Commelina yellow mottle virus, CoYMV). Se transformó lima Mexicana con construcciones que portaban el gen p23 completo, ya sea de la cepa agresiva de CTV T36 o de la suave T317, o con un fragmento que comprende el dedo de zinc y los motivos básicos flanqueantes de la primera, todas ellas bajo el control bien del promotor de CoYMV o bien del promotor constitutivo 35S. La expresión de estas construcciones en el floema dio lugar a aberraciones semejantes a los síntomas específicos de CTV, pero no a los síntomas inespecíficos observados cuando se expresaba p23 de forma constitutiva. Por otra parte, la apariencia e intensidad de las aberraciones fenotípicas más notorias similares a síntomas inducidos por CTV generadas por la expresión específica en floema del gen p23 se relacionó positivamente con la agresividad de la cepa origen utilizada. Además, la expresión en tejidos floemáticos del fragmento de p23 que comprende el dominio de dedo de zinc y los motivos básicos flanqueantes fue suficiente para inducir síntomas semejantes a los producidos por la infección con CTV, confirmando así que la región N-terminal delimitada por los aminoácidos 1 y 157 podría determinar, al menos en parte, la patogénesis de CTV en lima Mexicana. / Citrus tristeza virus (CTV) is the causal agent of one of the most devastating viral diseases of citrus trees in the world. CTV is phloem-restricted in natural citrus hosts, and has evolved three silencing suppressor proteins acting at intra- (p23 and p20) and inter-cellular level (p20 and p25) to overcome strong host antiviral defense in citrus. RNA interference (RNAi), an approach based on using dsRNA to trigger RNA silencing, has been widely used for generating transgenic plants resistant against viruses. Considering the important role of p23, p20 and p25 in CTV pathogenesis, we have transformed Mexican lime plants with an intron-hairpin vector carrying full untranslatable versions of genes p25, p20, p23 and the 3¿-UTR from the CTV strain T36, to attempt silencing their expression in CTV-infected cells. Complete resistance to viral infection was observed in three transgenic lines, with all their propagations remaining symptomless and virus-free after graft-inoculation with CTV-T36, either in the non-transgenic rootstock or directly in the transgenic scion. Accumulation of transgene-derived siRNAs was necessary but not sufficient for CTV resistance. Challenging immune transformants with a divergent CTV strain resulted in partial breakage of the resistance, stressing the importance of sequence identity in the underlying RNAi mechanism. This is the first evidence that it is possible to achieve full resistance to CTV in a highly sensitive citrus host by targeting simultaneously its three viral silencing suppressors through RNAi. The p23 protein encoded by the virus is additionally an important pathogenicity factor. Ectopic expression of p23 in transgenic citrus plants induces developmental aberrations resembling CTV symptoms. To explore in more detail the role of p23 in CTV pathogenesis, the p23 gene from CTV T36 and three truncated versions thereof under the control of the Cauliflower mosaic virus 35S promoter were used to transform Mexican lime. Only the truncated version expressing amino acids 1 to 157 (p23¿158-209) elicited CTV-like symptoms, similar to, albeit milder than, those incited by expressing the whole p23 protein (209 amino acids), thus delimiting the region responsible for p23 pathogenesis in citrus to a 157 amino acid fragment including the Zn finger and flanking basic motifs of the protein. RNA silencing suppressor activity of p23 in N. benthamiana was abolished by all mutants tested, indicating that silencing suppression involves most p23 regions. To better define the role of p23 in CTV pathogenesis, we next restricted the expression of p23-derived transgenes to phloem-associated cells in Mexican lime plants by means of using the phloem-specific promoter from Commelina yellow mottle virus (CoYMV). Constructions carrying the complete gene p23 from either the severe T36 or the mild T317 CTV strains, or a fragment comprising the zinc-finger and flanking basic motifs from the former, either under the control of the CoYMV promoter or the constitutive 35S promoter were used for genetic transformation of Mexican lime. Expression of these constructs in the phloem incited aberrations resembling CTV-specific symptoms, but not the unspecific symptoms observed when p23 was constitutively expressed. Moreover, appearance and intensity of the most notorious CTV-like phenotypic aberrations induced by the phloem-specific expression of the p23 gene were positively related with the aggressiveness of the source CTV strain used. Additionally, expression in phloem-tissues of the p23 fragment comprising the zinc-finger domain and flanking basic motifs was sufficient to induce CTV-like symptoms, corroborating that the N-terminal region (delimited by amino acids 1 and 157) determines, at least in part, CTV pathogenesis in Mexican lime. / Soler Calvo, N. (2013). Transgenic resistance against Citrus tristeza virus (CTV) and analysis of the viral p23 protein as pathogenicity determinant in citrus [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/31631 / TESIS

Molecular modeling is a term referring to the study of proteins, nucleic acids, lipids, and other bio- or macro- or small molecules at the atomistic level using a combination of computational methods, physico-chemical principles, and mathematical functions. It can be generally sub-divided into two areas: molecular mechanics, which is the treatment of atoms and bonds as Newtonian particles and springs, and quantum mechanics, which models electronic behaviors using the Schrödinger equation and wavefunctions. Each technique is a powerful tool that, when used alone or in combination with wet lab experiments, can yield useful results, the products of which have broad applications in studying human disease models, oxidative damage, and other biomolecular processes that are otherwise not easily observed by experiment alone. Within this document, we study seven different such systems. This includes the mode of inhibitor binding to the enzyme monoamine oxidase B, the active site mechanism of that same enzyme, the dynamics of the unstructured p53 C-terminal domain in complex with globular, structured proteins, the process of the viral protein B2 unbinding from double-stranded RNA, and a focus on the dynamics of a variable loop in the antigenic peanut protein Ara h 2. In addition to those conventional molecular modeling studies, several of which were done in tandem with wet lab experiment, we also discuss the validation of charges and charge group parameters for small molecules used in molecular mechanics, and the development of software for the analysis of lipid bilayer systems in molecular mechanics simulations. As computational resources continue to evolve, and as more structural information becomes available, these methods are becoming an integral part of the study of biomolecules in the context of disease. / Ph. D.

Banana bunchy top virus (BBTV) causes one of the most devastating diseases of banana. Transgenic virus resistance is now considered one of the most promising strategies to control BBTV. Pathogen-derived resistance (PDR) strategies have been applied successfully to generate plants that are resistant to numerous different viruses, primarily against those viruses with RNA genomes. BBTV is a circular, single-stranded (css) DNA virus of the family Nanoviridae, which is closely related to the family Geminiviridae. Although there are some successful examples of PDR against geminiviruses, PDR against the nanoviruses has not been reported. Therefore, the aim of this thesis was to investigate the potential of BBTV genes to interfere with virus replication when used as transgenes for engineering banana plants resistance to BBTV. The replication initiation protein (Rep) of nanoviruses is the only viral protein essential for viral replication and represents an ideal target for PDR. Therefore, this thesis focused on the effect of wild-type or mutated Rep genes from BBTV satellite DNAs or the BBTV integral genome on the replication of BBTV in banana embryogenic cell suspensions. A new Rep-encoding satellite DNA, designated BBTV DNA-S4, was isolated from a Vietnamese BBTV isolate and characterised. When the effect of DNA-S4 on the replication of BBTV was examined, it was found that DNA-S4 enhanced the replication of BBTV. When the replicative capabilities of DNA-S4 and the previously characterised Rep-encoding BBTV satellite, DNA-S1, were compared, it was found that the amount of DNA-S4 accumulated to higher levels than DNA-S1. The interaction between BBTV and DNA-S1 was also examined. It was found that over-expression of the Rep encoded by DNA-S1 using ubi1 maize polyubiquitin promoter enhanced replication of BBTV. However, when the Rep-encoded by DNA-S1 was expressed by the native S1 promoter (in plasmid pBT1.1-S1), it suppressed the replication of BBTV. Based on this result, the use of DNA-S1 as a possible transgene to generate PDR against BBTV was investigated. The roles of the Rep-encoding and U5 genes of BBTV DNA-R, and the effects of over-expression of these two genes on BBTV replication were also investigated. Three mutants of BBTV DNA-R were constructed; plasmid pUbi-RepOnly-nos contained the ubi1 promoter driving Rep expression from DNA-R, plasmid pUbi-IntOnly-nos contained the ubi1 promoter driving expression of the DNA-R internal gene product (U5), while plasmid pUbi-R.ORF-nos contained the ubi1 promoter driving the expression of both Rep and the internal U5 gene product. The replication of BBTV was found to be significantly suppressed by pUbi-RepOnly-nos, weakly suppressed by pUbi-IntOnly-nos, but strongly enhanced by pUbi-R.ORF-nos. The effect of mutations in three conserved residues within the BBTV Rep on BBTV replication was also assessed. These mutations were all made in the regions in the ATPase motifs and resulted in changes from hydrophilic to hydrophobic residues (i.e. K187→M, D224→I and N268→L). None of these Rep mutants was able to initiate BBTV replication. However, over-expression of Reps containing the K187→M or N268→L mutations significantly suppressed the replication of BBTV. In summary, the Rep constructs that significantly suppressed replication of DNA-R and -C in banana embryogenic cell suspensions have the potential to confer resistance against BBTV by interfering with virus replication. It may be concluded that BBTV satellite DNAs are not ideal for conferring PDR because they did not suppress BBTV replication consistently. Wild-type Rep transcripts and mutated (i.e. K187→M and N248→L) Rep proteins of BBTV DNA-R, however, when over-expressed by a strong promoter, are all promising candidates for generating BBTV-resistant banana plants.

Knowledge of the dissociation constants of the ionizable protons of weak acids in aqueous media is of fundamental importance in many areas of chemistry and biochemistry. The pKa value, or equilibrium dissociation constant, of a molecule determines the relative concentration of its protonated and deprotonated forms at a specified pH and is therefore an important descriptor of its chemical reactivity. Considerable efforts have been devoted to the determination of pKa values by different experimental techniques. Although in most cases the determination of pKa values from experimental is straightforward, there are situations where interpretation is difficult and the results ambiguous. It is, therefore, not surprising that the capability to provide accurate estimates of the pKa value has been a central goal in theoretical chemistry and there has been a large effort in developing methodologies for predicting pKa values for a variety of chemical systems by differing quantum chemical techniques. A prediction accuracy within 0.5 pKa units of experiment is the desirable level of accuracy. This is a non-trivial exercise, for an error of 1 kcal/mol in estimates of the free energy value would result in an error of 0.74 pKa units. In this thesis ab initio Car-Parrinello molecular dynamics (CPMD) has been used for investigating the Brϕnsted acid-base chemistry of weak acids in aqueous solution. A key issue in any dissociation event is how the solvating water molecules arrange themselves spatially and dynamically around the neutral and dissociated acid molecule. Ab initio methods have the advantage that all solvent water molecules can, in principle, be con- sidered explicitly. One of the factors that has inhibited the widespread use of ab initio MD methods to study the dissociation reaction is that dissociation of weak acids are rare events that require extremely long simulation times before one is observed. The metady- namics formalism provides a solution to this conundrum by preventing the system from revisiting regions of configuration space where it has been in the past. The formalism allows the system to escape the free-energy minima by biasing the dynamics with a history dependent potential (or force) that acts on select degrees of freedom, referred to as collective variables. The bias potentials, modeled by repulsive inverted Gaussians that are dropped during propagation, drive the system out of any free-energy minima and allow it to explore the configurational space by a relatively quick and efficient sampling. The the- sis deals with a detailed investigation of the Brϕnsted acid-base chemistry of weak acids in aqueous solutions by the CPMD-metadynamics procedure. In Chapter 1, current approaches for the theoretical estimation of pKa values are summarized while in Chapter 2 the simulation methodology and the metadynamics sampling techniques used in this study are described. The potential of the CPMD-metadynamics procedure to provide estimates of the acid dissociation constant (pKa) is explored in Chapter 3, using acetic acid as a test sys- tem. Using the bond-distance dependent coordination number of protons bound to the dissociating carboxylic groups as the collective variable, the free-energy profile for the dissociation reaction of acetic acid in water was computed. Convergence of the free-energy profiles and barriers for the simulations parameters is demonstrated. The free-energy profiles exhibit two distinct minima corresponding to the dissociated and neutral states of the acid and the deterrence in their values provides the estimate for pKa. The estimated value of pKa for acetic acid from the simulations, 4.80, is in good agreement with the experiment at value of 4.76. It is shown that the good agreement with experiment is a consequence of the cancellation of errors, as the pKa values are computed as the divergence in the free energy values at the minima corresponding to the neutral and dissociated state. The chapter further explores the critical factors required for obtaining accurate estimates of the pKa values by the CPMD-metadynamics procedure. It is shown that having water molecules sufficient to complete three hydration shells as well as maintaining water density in the simulation cell as close to unity is important. In Chapter 4, the CPMD-metadynamics procedure described in Chapter-3 has been used to investigate the dissociation of a series of weak organic acids in aqueous solutions. The acids studied were chosen to highlight some of the major factors that influence the dissociation constant. These include the influence of the inductive effect, the stabilization of the dissociated anion by H-bonding as well as the presence of multiple ionizable groups. The acids investigated were aliphatic carboxylic acids, chlorine-substituted carboxylic acids, cis- and trans-butenedioic, the isomers of hydroxybenzoic acid and ophthalmic acids and its isomers. It was found that in each of these examples the CPMD-metadynamics procedure correctly estimates the pKa values, indicating that the formulism is capable of capturing these influences and equally importantly indicating that the cancellation of errors is indeed universal. Further, it is shown that the procedure can provide accurate estimates of the successive pKa values of polypro tic acids as well as the subtle difference in their values for different isomers of the acid molecule. Changes in protonation-deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. It is shown that CPMD simulations in conjunction with metadynamics calculations of the free energy profile of the protonation- deprotonation reaction can provide estimates of the multiple pKa values of the 20 canonical α-amino acids in aqueous solutions in good agreement with experiment (Chapter 5). The distance-dependent coordination number of the protons bound to the hydroxyl oxygen of the carboxylic and the amine groups is used as the collective variable to explore the free energy profiles of the Brϕnsted acid-base chemistry of amino acids in aqueous solutions. Water molecules, sufficient to complete three hydration shells surrounding the acid molecule were included explicitly in the computation procedure. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pKa values with a mean relative error with respect to experimental results, of 0.2 pKa units. The tripeptide Glutathione (GSH) is one of the most abundant peptides and the major repository for non-protein sulfur in both animal and plant cells. It plays a critical role in intracellular oxidative stress management by the reversible formation of glutathione disulfide with the thioldisulfide pair acting as a redox buffer. The state of charge of the ionizable groups of GSH can influences the redox couple and hence the pKa value of the cysteine residue of GSH is critical to its functioning. In Chapter 6, it has been reported that ab initio Car-Parrinello Molecular Dynamics simulations of glutathione solvated by 200 water molecules, all of which are considered in the simulation. It is shown that the free-energy landscape for the protonation - deprotonation reaction of the cysteine residue of GSH computed using metadynamics sampling provides accurate estimates of the pKa and correctly predicts the shift in the dissociation constant values as compared to the isolated cysteine amino acid. The dissociation constants of weak acids are commonly determined from pH-titration curves. For simple acids the determination of the pKa from the titration curves using the Henderson-Hasselbalch equation is relatively straightforward. There are situations, however, especially in polyprotic acids with closely spaced dissociation constants, where titration curves do not exhibit clear inflexion and equivalence stages and consequently the estimation of multiple pKa values from a single titration curve is no longer straightfor- ward resulting in uncertainties in the determined pKa values. In Chapter 7, the multiple dissociation constant of the hexapeptide glutathione disulfide (GSSG) with six ionizable groups and six associated dissociation constants has been investigated. The six pKa values of GSSG were estimated using the CPMD-metadynamics procedure from the free-energy profiles for each dissociation reaction computed using the appropriate collective variable. The six pKa values of GSSG were estimated and the theoretical pH-titration curve was then compared with the experimentally measured pH-titration curve and found to be in excellent agreement. The object of the exercise was to establish whether interpretation of pH-titration curves of complex molecules with multiple ionizable groups could be facilitated using results of ab initio molecular dynamics simulations.