Transcription factors need to rapidly find one specific binding site among millions of nonspecific sites on the chromosomal DNA. In this thesis I use various aspects of reaction-diffusion theory to investigate the interaction between proteins and DNA and to explain the searching, finding and binding to specific operator sites. Using molecular dynamics methods we calculate the free energy profile for the model protein LacI as it leaves a nonspecific stretch of DNA and as it slides along DNA. Based on the free energy profiles we estimate the microscopic dissociation rate constant, kdmicro ~1.45×104s-1, and the 1D diffusion coefficient, D1 ~ 0.05-0.29 μm2s-1 (2-40μs to slide 1 basepair (bp)). At a non-atomistic level of detail we estimate the number of microscopic rebindings before a macroscopic dissociation occurs which leads to the  macroscopic residence time, τDmacro ~ 48±12ms resulting in a in vitro sliding length estimate of 135-345bp. When we fit the DNA interaction parameters for in vivo conditions to recent single molecule in vivo experiments we conclude that neither hopping nor intersegment transfer contribute to the target search for the LacI dimer, that it appears to bind the specific Osym operator site as soon as it slides into it, and that the sliding length is around 40bp in the cell. The estimated in vivo D1 ~ 0.025 μm2s-1 is higher than expected from estimates of D1 based on viscosity and the atomistic simulations. Surprisingly, we were also forced to conclude that the nonspecific association for the LacI dimer appeared reaction limited which is in conflict with the free energy profile. This inconsistency is resolved by allowing for steric effects. Using reaction-diffusion theory and simulations we show that an apparent reaction limited association can be diffusion limited if geometry and steric effects are taken into account. Furthermore, the simulations show that a protein binds ~2 times faster to a DNA molecule with a helical reactive patch than to a stripe patch running along the length of the DNA. This facilitated binding has a direct impact on the search time especially in the presence of other DNA binding proteins.

Orientador: Eudes Eterno Fileti / Dissertação (mestrado) - Universidade Federal do ABC. Programa de Pós-Graduação em Nanociências e Materiais Avançados, 2013

Orientador: Eudes Eterno Fileti / Tese (doutorado) - Universidade Federal do ABC. Programa de Pós-Graduação em Nanociências e Materiais Avançados, 2013

The current thesis presents work on the structure and dynamics of oligosaccharides and polysaccharides as well as the free energetics of carbohydrate-protein interactions. By applying various computational tools such as molecular dynamics simulation, our in-house fast sugar structure prediction software, replica exchange molecular dynamics, homology modeling, umbrella sampling, steered molecular dynamics as well as the thermodynamic integration formalism, we have been able to study the role of water on the surface of homopolysaccharides as well as complex oligosachharides, we have been able to produce a prediction of the bound structure of triantennary oligosaccride on the asialoglycoprotein receptor, we have been able to estimate the free energy of binding of ManΑ1→2Man to the HIV-1 inactivating protein, Cyanovirin-N as well as the relative binding free energies of mutants of Cyanovirin-N to the same ligand.

No description available.

Life is soft. From the fluid-like structure of lipid bilayers to the flexible folding of proteins, the realm of nanoscale soft matter is a complex and vibrant area of research. The lure of personalized medicine, advanced sensing technology, and understanding life at a fundamental level pushes research forward. This work considers to areas: (1) lipid bilayer dynamics in the presence of substrate defects and (2) the inverse temperature transition of elastic proteins. Molecular dynamics simulations as well as umbrella sampling were employed. The behavior of the bilayers discussed in the work provides evidence that small defects on confining surfaces can promote nucleation of lipid tethers. Results the second part of this work indicate elastin-like peptides experiencing inverse temperature transitions may be capable of performing amounts of work similar to RNA polymerase; additionally, resilin's inverse temperature transition may be closely linked to the molecule's ability to efficiently transmit energy through the similar coil-β secondary structure transition seen in both cases. These insights into the inverse transition temperature are relevant for the design of bio-inspired sensors and energy storage devices. / Master of Science

Antibody microarrays constitute a next-generation sensing platform that has the potential to revolutionize the way that molecular detection is conducted in many scientific fields. Unfortunately, current technologies have not found mainstream use because of reliability problems that undermine trust in their results. Although several factors are involved, it is believed that undesirable protein interactions with the array surface are a fundamental source of problems where little detail about the molecular-level biophysics are known. A better understanding of antibody stability and antibody-antigen binding on the array surface is needed to improve microarray technology. Despite the availability of many laboratory methods for studying protein stability and binding, these methods either do not work when the protein is attached to a surface or they do not provide the atomistic structural information that is needed to better understand protein behavior on the surface. As a result, molecular simulation has emerged as the primary method for studying proteins on surfaces because it can provide metrics and views of atomistic structures and molecular motion. Using an advanced, coarse-grain, protein-surface model this study investigated how antibodies react to and function on different types of surfaces. Three topics were addressed: (1) the stability of individual antibodies on surfaces, (2) antibody binding to small antigens while on a surface, and (3) antibody binding to large antigens while on a surface. The results indicate that immobilizing antibodies or antibody fragments in an upright orientation on a hydrophilic surface can provide the molecules with thermal stability similar to their native aqueous stability, enhance antigen binding strength, and minimize the entropic cost of binding. Furthermore, the results indicate that it is more difficult for large antigens to approach the surface than small antigens, that multiple binding sites can aid antigen binding, and that antigen flexiblity simultaneously helps and hinders the binding process as it approaches the surface. The results provide hope that next-generation microarrays and other devices decorated with proteins can be improved through rational design.

abstract: Conformational changes in biomolecules often take place on longer timescales than are easily accessible with unbiased molecular dynamics simulations, necessitating the use of enhanced sampling techniques, such as adaptive umbrella sampling. In this technique, the conformational free energy is calculated in terms of a designated set of reaction coordinates. At the same time, estimates of this free energy are subtracted from the potential energy in order to remove free energy barriers and cause conformational changes to take place more rapidly. This dissertation presents applications of adaptive umbrella sampling to a variety of biomolecular systems. The first study investigated the effects of glycosylation in GalNAc2-MM1, an analog of glycosylated macrophage activating factor. It was found that glycosylation destabilizes the protein by increasing the solvent exposure of hydrophobic residues. The second study examined the role of bound calcium ions in promoting the isomerization of a cis peptide bond in the collagen-binding domain of Clostridium histolyticum collagenase. This study determined that the bound calcium ions reduced the barrier to the isomerization of this peptide bond as well as stabilizing the cis conformation thermodynamically, and identified some of the reasons for this. The third study represents the application of GAMUS (Gaussian mixture adaptive umbrella sampling) to on the conformational dynamics of the fluorescent dye Cy3 attached to the 5' end of DNA, and made predictions concerning the affinity of Cy3 for different base pairs, which were subsequently verified experimentally. Finally, the adaptive umbrella sampling method is extended to make use of the roll angle between adjacent base pairs as a reaction coordinate in order to examine the bending both of free DNA and of DNA bound to the archaeal protein Sac7d. It is found that when DNA bends significantly, cations from the surrounding solution congregate on the concave side, which increases the flexibility of the DNA by screening the repulsion between phosphate backbones. The flexibility of DNA on short length scales is compared to the worm-like chain model, and the contribution of cooperativity in DNA bending to protein-DNA binding is assessed. / Dissertation/Thesis / Ph.D. Chemistry 2011

La flexibilité est une caractéristique intrinsèque des protéines qui doivent, dès le mo- ment de leur synthèse, passer d’un état de chaîne linéaire à un état de structure tridimen- sionnelle repliée et enzymatiquement active. Certaines protéines restent flexibles une fois repliées et subissent des changements de conformation de grande amplitude lors de leur cycle enzymatique. D’autres contiennent des segments si flexibles que leur structure ne peut être résolue par des méthodes expérimentales. Dans cette thèse, nous présentons notre application de méthodes in silico d’analyse de la flexibilité des protéines : • À l’aide des méthodes de dynamique moléculaire dirigée et d’échantillonnage pa- rapluie, nous avons caractérisé les trajectoires de liaison de l’inhibiteur Z-pro- prolinal à la protéine Prolyl oligopeptidase et identifié la trajectoire la plus pro- bable. Nos simulations ont aussi identifié un mode probable de recrutement des ligands utilisant une boucle flexible de 19 acides aminés à l’interface des deux domaines de la protéine. • En utilisant les méthodes de dynamique moléculaire traditionnelle et dirigée, nous avons examiné la stabilité de la protéine SAV1866 dans sa forme fermée insérée dans une membrane lipidique et étudié un des modes d’ouverture possibles par la séparation de ses domaines liant le nucléotide. • Nous avons adapté auproblème de la prédiction de la structure des longues boucles flexibles la méthode d’activation et de relaxation ART-nouveau précédemment uti- lisée dans l’étude du repliement et de l’agrégation de protéines. Appliqué au replie- ment de boucles de 8 à 20 acides aminés, la méthode démontre une dépendance quadratique du temps d’exécution sur la longueur des boucles, rendant possible l’étude de boucles encore plus longues. / Flexibility is an intrinsic characteristic of proteins who from the moment of synthesis into a linear chain of amino acids, have to adopt an enzymatically active tridimensionnel structure. Some proteins stay flexible once folded and display large amplitude confor- mational changes during their enzymatic cycles. Others contain parts that are so flexible that their structure can’t be resolved using experimental methods. In this thesis, we present our application of in silico methods to the study of protein flexibility. • Using steered molecular dynamics and umbrella sampling, we characterized the binding trajectories of the Z-pro-prolinal inhibiter to the Prolyl oligopeptidase pro- tein and we identified the most probable trajectory. Our simulations also found a possible ligand recrutement mechanism that involves a 19 amino acids flexible loop at the interface of the two domains of the protein. • Using traditional and steered molecular dynamics, we examined the stability of the SAV1866 protein in its closed conformation in a lipid membrane and we studied one of its proposed opening modes by separating its nucleotide binding domains. • We also adapted the activation-relaxation technique ART-nouveau which was pre- viously used to study protein folding and aggregation to the problem of structure prediction of large flexible loops. When tested on loops of 8 to 20 amino acids, the method demonstrate a quadratic execution time dependance on the loop length, which makes it possible to use the method on even larger loops.

La flexibilité est une caractéristique intrinsèque des protéines qui doivent, dès le mo- ment de leur synthèse, passer d’un état de chaîne linéaire à un état de structure tridimen- sionnelle repliée et enzymatiquement active. Certaines protéines restent flexibles une fois repliées et subissent des changements de conformation de grande amplitude lors de leur cycle enzymatique. D’autres contiennent des segments si flexibles que leur structure ne peut être résolue par des méthodes expérimentales. Dans cette thèse, nous présentons notre application de méthodes in silico d’analyse de la flexibilité des protéines : • À l’aide des méthodes de dynamique moléculaire dirigée et d’échantillonnage pa- rapluie, nous avons caractérisé les trajectoires de liaison de l’inhibiteur Z-pro- prolinal à la protéine Prolyl oligopeptidase et identifié la trajectoire la plus pro- bable. Nos simulations ont aussi identifié un mode probable de recrutement des ligands utilisant une boucle flexible de 19 acides aminés à l’interface des deux domaines de la protéine. • En utilisant les méthodes de dynamique moléculaire traditionnelle et dirigée, nous avons examiné la stabilité de la protéine SAV1866 dans sa forme fermée insérée dans une membrane lipidique et étudié un des modes d’ouverture possibles par la séparation de ses domaines liant le nucléotide. • Nous avons adapté auproblème de la prédiction de la structure des longues boucles flexibles la méthode d’activation et de relaxation ART-nouveau précédemment uti- lisée dans l’étude du repliement et de l’agrégation de protéines. Appliqué au replie- ment de boucles de 8 à 20 acides aminés, la méthode démontre une dépendance quadratique du temps d’exécution sur la longueur des boucles, rendant possible l’étude de boucles encore plus longues. / Flexibility is an intrinsic characteristic of proteins who from the moment of synthesis into a linear chain of amino acids, have to adopt an enzymatically active tridimensionnel structure. Some proteins stay flexible once folded and display large amplitude confor- mational changes during their enzymatic cycles. Others contain parts that are so flexible that their structure can’t be resolved using experimental methods. In this thesis, we present our application of in silico methods to the study of protein flexibility. • Using steered molecular dynamics and umbrella sampling, we characterized the binding trajectories of the Z-pro-prolinal inhibiter to the Prolyl oligopeptidase pro- tein and we identified the most probable trajectory. Our simulations also found a possible ligand recrutement mechanism that involves a 19 amino acids flexible loop at the interface of the two domains of the protein. • Using traditional and steered molecular dynamics, we examined the stability of the SAV1866 protein in its closed conformation in a lipid membrane and we studied one of its proposed opening modes by separating its nucleotide binding domains. • We also adapted the activation-relaxation technique ART-nouveau which was pre- viously used to study protein folding and aggregation to the problem of structure prediction of large flexible loops. When tested on loops of 8 to 20 amino acids, the method demonstrate a quadratic execution time dependance on the loop length, which makes it possible to use the method on even larger loops.