Peptide-membrane interactions play an important role in a number of biological processes, such as antimicrobial defence mechanisms, viral translocation, membrane fusion and functions ofmembrane proteins. In particular, amphipathic α-helical peptides comprise a large family of membrane-active peptides that could exhibit a broad range of biological activities. A membrane, interacting with an amphipathic α-helical peptide, may experience a number of possible structural transitions, including stretching, reorganization of lipid molecules, formation of defects, transient and stable pores, formation of vesicles, endo- and pinocytosis and other phenomena. Naturally, theoretical and experimental studies of these interactions have been an intense on-going area of research. However, complete understanding of the relationship between the structure of the peptide and themechanismof interaction it induces, as well asmolecular details of this process, still remain elusive. Lack of this knowledge is a key challenge in our efforts to elucidate some of the biological functions of membrane active peptides or to design peptides with tailored functionalities that can be exploited in drug delivery or antimicrobial strategies. In principle,molecular dynamics is a powerful research tool to study peptide-membrane interactions, which can provide a detailed description of these processes on molecular level. However, a model operating on the appropriate time and length scale is imperative in this description. In this study, we adopt a coarse-grained approach where the accessible simulation time and length scales reach microseconds and tens of nanometers, respectively. Thus, the two key objectives of this study are to validate the applicability of the adopted coarse-grained approach to the study of peptide-membrane interactions and to provide a systematic description of these interactions as a function of peptide structure and surface chemistry. We applied the adopted strategy to a range of peptide systems, whose behaviour has been well established in either experiments or detailed atomistic simulations and outlined the scope and applicability of the coarse-grained model. We generated some useful insights on the relationship between the structure of the peptides and themechanism of peptide-membrane interactions. Particularly interesting results have been obtained for LS3, a membrane spanning peptide, with a propensity to self-assembly into ion-conducting channels. Firstly, we captured, for the first time, the complete process of self-assembly of LS3 into a hexameric ion-conducting channel and explored its properties. The channel has structure of a barrel-stave pore with peptides aligned along the lipid tails. However, we discovered that a shorter version of the peptide leads to a more disordered, less stable structure often classified as a toroidal pore. This link between two types of pores has been established for the first time and opens interesting opportunities in tuning peptide structures for a particular pore-inducing mechanism. We also established that different classes of peptides can be uniquely characterized by the distinct energy profile as they cross the membrane. Finally, we extended this investigation to the internalization mechanisms of more complex entities such as peptide complexes and nanoparticles. Coarse-grained steered molecular dynamics simulations of these model systems are performed and some preliminary results are presented in this thesis. To summarize, in this thesis, we demonstrate that coarse-grained models can be successfully used to underpin peptide interaction and self-assembly processes in the presence of membranes in their full complexity. We believe that these simulations can be used to guide the design of peptides with tailored functionalities for applications such as drug delivery vectors and antimicrobial systems. This study also suggests that coarse-grained simulations can be used as an efficient way to generate initial configurations for more detailed atomistic simulations. These multiscale simulation ideas will be a natural future extension of this work.

abstract: Improved knowledge connecting the chemistry, structure, and properties of polymers is necessary to develop advanced materials in a materials-by-design approach. Molecular dynamics (MD) simulations can provide tremendous insight into how the fine details of chemistry, molecular architecture, and microstructure affect many physical properties; however, they face well-known restrictions in their applicable temporal and spatial scales. These limitations have motivated the development of computationally-efficient, coarse-grained methods to investigate how microstructural details affect thermophysical properties. In this dissertation, I summarize my research work in structure-based coarse-graining methods to establish the link between molecular-scale structure and macroscopic properties of two different polymers. Systematically coarse-grained models were developed to study the viscoelastic stress response of polyurea, a copolymer that segregates into rigid and viscous phases, at time scales characteristic of blast and impact loading. With the application of appropriate scaling parameters, the coarse-grained models can predict viscoelastic properties with a speed up of 5-6 orders of magnitude relative to the atomistic MD models. Coarse-grained models of polyethylene were also created to investigate the thermomechanical material response under shock loading. As structure-based coarse-grained methods are generally not transferable to states different from which they were calibrated at, their applicability for modeling non-equilibrium processes such as shock and impact is highly limited. To address this problem, a new model is developed that incorporates many-body interactions and is calibrated across a range of different thermodynamic states using a least square minimization scheme. The new model is validated by comparing shock Hugoniot properties with atomistic and experimental data for polyethylene. Lastly, a high fidelity coarse-grained model of polyethylene was constructed that reproduces the joint-probability distributions of structural variables such as the distributions of bond lengths and bond angles between sequential coarse-grained sites along polymer chains. This new model accurately represents the structure of both the amorphous and crystal phases of polyethylene and enabling investigation of how polymer processing such as cold-drawing and bulk crystallization affect material structure at significantly larger time and length scales than traditional molecular simulations. / Dissertation/Thesis / Doctoral Dissertation Mechanical Engineering 2016

Structure and mechanical properties of DNA play a key role in its biological functioning. A lot of well-established conclusions about the DNA structure and its sequence-dependent variabil- ity came from various experimental and computational studies of the Dickerson-Drew dodecamer (DD), a prototypic B-DNA molecule of the sequence (5')CGCGAATTCGCG(3'). In this study we present a detailed analysis of structural and mechan- ical properties of DD based on extensive atomistic molecular dynamics (MD) simulations with explicit representation of wa- ter and ionic environment. We analyze three simulated systems covering different ionic conditions and water models. Two MD trajectories are reported for the first time, one of them being 2.4 µs long. An extensive comparsion with one recent NMR struc- ture and four recent X-ray structures is made. It is found that the end basepairs can adopt two different pairing motifs dur- ing the simulation: the canonical Watson-Crick pair or a non- canonical trans Watson-Crick/Sugar Edge pair. These states can significantly influence the structure of DD even at the third step from the end. A clear relationship is found between the BI/BII backbone substates and the basepair step conformation. A model of rigid bases is used to study mechanical properties of the DNA. The non-local...

Structural and mechanical properties of nucleic acids play a key role in a wide range of biological processes, as well as in the field of nucleic acid nanotechnology. The thesis presents results of several studies focused on modelling these properties. Extensive unrestrained atomic-resolution molecular dynamics (MD) simulations are used to investigate structural dynamics of nucleic acids, and to parametrize their mechanical models. The deformation energy is assumed to be a general quadratic function of suitably chosen internal coordinates. Two types of models are employed which differ in the level of coarse- graining. The first one is based on the description of conformation at the level of individual bases and the second, coarser one is used to study global bending and twisting flexibility. The models are applied to explain mechanical properties of A-tracts in the context of DNA looping and nucleosome positioning, to characterize twist-stretch cou- pled deformations in DNA and RNA, and to predict changes in the properties of damaged DNA that are likely to be relevant for damage recognition and repair. Besides that, we propose a general model of DNA allostery, applied to study the effect of minor groove binding of small ligands and the allosteric coupling between proteins mediated by the DNA. A careful...

The presented work is the application of recent methodologies on modeling and simulation of single stranded RNA viruses. We first present the methods of modeling RNA molecules using the coarse-grained modeling package, YUP. Coarse-grained models simplify complex structures such as viruses and let us study general behavior of the complex biological systems that otherwise cannot be studied with all-atom details. Second, we modeled the first all-atom T=3, icosahedral, single stranded RNA virus, Pariacoto virus (PaV). The x-ray structure of PaV shows only 35% of the total RNA genome and 88% of the capsid. We modeled both missing portions of RNA and protein. The final model of the PaV demonstrated that the positively charged protein N- terminus was located deep inside the RNA. We propose that the positively charged N- terminal tails make contact with the RNA genome and neutralize the negative charges in RNA and subsequently collapse the RNA/protein complex into an icosahedral virus. Third, we simulated T=1 empty capsids using a coarse-grained model of three capsid proteins as a wedge-shaped triangular capsid unit. We varied the edge angle and the potentials of the capsid units to perform empty capsid assembly simulations. The final model and the potential are further improved for the whole virus assembly simulations. Finally, we performed stability and assembly simulations of the whole virus using coarse-grained models. We tested various strengths of RNA-protein tail and capsid protein-capsid protein attractions in our stability simulations and narrowed our search for optimal potentials for assembly. The assembly simulations were carried out with two different protocols: co-transcriptional and post-transcriptional. The co-transcriptional assembly protocol mimics the assembly occurring during the replication of the new RNA. Proteins bind the partly transcribed RNA in this protocol. The post-transcriptional assembly protocol assumes that the RNA is completely transcribed in the absence of proteins. Proteins later bind to the fully transcribed RNA. We found that both protocols can assemble viruses, when the RNA structure is compact enough to yield a successful virus particle. The post-transcriptional protocol depends more on the compactness of the RNA structure compared to the co-transcriptional assembly protocol. Viruses can exploit both assembly protocols based on the location of RNA replication and the compactness of the final structure of the RNA.

Computer simulations provide a potentially powerful complement to conventional experimental techniques in elucidating the structures, dynamics and interactions of macromolecules. In this thesis, I present three applications of computer simulations to investigate important biomolecules with sizes ranging from two-residue peptides, to proteins, and to whole chromosome structures. First, I describe the results of 441 independent explicit-solvent molecular dynamics (MD) simulations of all possible two-residue peptides that contain the 20 standard amino acids with neutral and protonated histidine. 3JHNHα coupling constants and δHα chemical shifts calculated from the MD simulations correlated quite well with recently published experimental measurements for a corresponding set of two-residue peptides. Neighboring residue effects (NREs) on the average 3JHNHα and δHα values of adjacent residues were also reasonably well reproduced. The intrinsic conformational preferences of each residue, and their NREs on the conformational preferences of adjacent residues, were analyzed. Finally, these NREs were compared with corresponding effects observed in a coil library and the average β-turn preferences of all residue types were determined. Second, I compare the abilities of three derivatives of the Amber ff99SB force field to reproduce a recent report of 3JHNHα scalar coupling constants for hundreds of two-residue peptides. All-atom MD simulations of 256 two-residue peptides were performed and the results showed that a recently-developed force field (RSFF2) produced a dramatic improvement in the agreement with experimental 3JHNHα coupling constants. I further show that RSFF2 also improved modestly agreement with experimental 3JHNHα coupling constants of five model proteins. However, an analysis of NREs on the 3JHNHα coupling constants of the two-residue peptides indicated little difference between the force fields’ abilities to reproduce experimental NREs. I speculate that this might indicate limitations in the force fields’ descriptions of nonbonded interactions between adjacent side chains or with terminal capping groups. Finally, coarse-grained (CG) models and multi-scale modeling methods are used to develop structural models of entire E. coli chromosomes confined within the experimentally-determined volume of the nucleoid. The final resolution of the chromosome structures built here was one-nucleotide-per-bead (1 NTB), which represents a significant increase in resolution relative to previously published CG chromosome models, in which one bead corresponds to hundreds or even thousands of basepairs. Based on the high-resolution final 1 NTB structures, important physical properties such as major and minor groove widths, distributions of local DNA bending angles, and topological parameters (Linking Number (Lk), Twist (Tw) and Writhe (Wr)) were accurately computed and compared with experimental measurements or predictions from a worm-like chain (WLC) model. All these analyses indicated that the chromosome models built in this study are reasonable at a microscopic level. This chromosome model provides a significant step toward the goal of building a whole-cell model of a bacterial cell.

<p>A thermodynamically consistent coarsed-grained phase field model was developed to find the conditions under which a heterogeneous porous electrode can be treated as homogeneous in the description of Lithium-ions in rechargeable batteries. Four regimes of behavior under which the transport phenomena can be homogenized to describe porous LIBs were identied: regime (a), where the model is inaccurate, for physically accessible particle packings of aspect ratios smaller than c/a = 0.5 and electrode porosities between 0.34 to 0.45; regime (b), where the model is valid, for particles of aspect ratios greater than c/a = 0.7 and electrode porosities greater than 0.35; regime (c), where the model is valid, but the microstructures are physically inaccessible, and correspond to particles with aspect ratios greater than c/a = 0.7 and electrode porosities smaller than 0.34; and regime (d), where the model is invalid and the porous microstructures are physically inaccessible, and correspond to particles with aspect ratios smaller than c/a = 1 and electrode porosities smaller than 0.34.</p> <p>The developed formulation was applied to the graphite | LixNi1/3Mn1/3Co1/3O2 system to analyze the effect of microstructure and coarsed-grained long-range chemomechanical effects on the electrochemical behavior. Specically, quantiable lithium distribution populations in the cathode, as a result of long range interactions of the diffuse interface, charge effects and mechanical stresses were identified: i) diffusion limited population due to negligible composition gradients, ii) stress-induced population as a result of chemically induced stresses, and iii) lithiation-induced population, as a consequence of the electrochemical potential gradients.</p>

Les Récepteurs Couplés aux Protéines-G forment la plus importante famille de protéines membranaires chez l’homme et sont impliqués dans de nombreux processus de signalisation cellulaire. Aussi, ils forment un vivier très important de cibles thérapeutiques, déjà identifiées ou potentielles. L’activation d’un RCPG est amorcée par la liaison d’un ligand dans sa partie extra-cellulaire, modifiant ainsi ses propriétés dynamiques intrinsèques. Ces changements structuraux vont alors se répercuter le long des domaines trans-membranaires et promouvoir la dissociation de la Protéine-G hétéro-trimérique, de l’autre côté de la membrane, propageant ainsi le signal au compartiment intra-cellulaire. Ce processus peut être modulé par la liaison de nombreux autres partenaires des RCPGs. Malgré de nombreuses données structurales existantes, ces mécanismes restent encore mal connus à l’échelle moléculaire. Ainsi, la dynamique moléculaire s’est révélée être un outil formidable pour mieux comprendre ces mécanismes. Toutefois, les échelles de taille et de temps requises pour discuter de la dynamique de ces systèmes membranaires limitent ces études aux laboratoires ayant accès à une très grande puissance de calcul. L’objectif des travaux présentés dans ce manuscrit a été de prédire et de mieux comprendre la dynamique d’interaction de différents récepteurs de cette famille avec leurs partenaires, en développant un protocole de dynamique moléculaire, peu coûteux en ressources de calcul, combinant le champ de forces gros-grains MARTINI à un protocole de dynamique moléculaire « Replica-Exchange ».Dans un premier temps, nous présentons la validation de notre protocole pour la prédiction de la liaison de peptides à leur récepteur avec l’étude des peptides Neurotensine, agoniste du Récepteur de la Neurotensine-1, et CVX15, antagoniste du Récepteur Chemokine C-X-C de type-4. Nous montrons également que notre protocole est capable de prédire la sélectivité de plusieurs peptides dérivés de la Neurotensine envers plusieurs récepteurs sauvages et mutés, ne présentant qu’un résidu de différence.Dans un second temps, nous nous sommes intéressés à la dynamique de formation d’un hétéro-dimère de RCPGs impliquant le Récepteur de la Ghréline et le récepteur de la Dopamine D2, couplés aux protéines Gq et Gi respectivement. Ce modèle validé au laboratoire par des mesures LRET montre une interface impliquant une forte complémentarité entre les protéines-G. En se basant sur notre modèle, nous avons conçu et synthétisé des peptides inhibiteurs de la formation de cet hétéro-dimère de protéines-G.Enfin, nous présentons d’autres exemples d’applications de notre protocole et comment il peut être utilisé de concert avec l’expérience avec : la prédiction de la liaison de toxines de serpents aux Récepteurs de la Vasopressine-1a et V2 ; la prédiction de la liaison des peptides Ghréline et Leap2 au Récepteur GHSR-1a et la prédiction de la sélectivité de couplage de différents récepteurs aux peptides C-terminaux de la sous-unité α des protéines-G. / G-Protein Coupled Receptors form the largest family of human membrane proteins and are involved in many cellular signaling processes. Thus, they constitute a pool of already identified or potential pharmacological targets. The activation of a GPCR starts with the binding of a ligand in its extra-cellular part, further modifying its intrinsic dynamical properties. These structural rearrangements are then transmitted along the transmembrane domains and promote the dissociation of the G-protein on the other side of the bilayer, thus propagating the signal into the intra-cellular compartment. This activation process can be modulated by the binding of many other partners of GPCRs. Despite many structural data now available, these mechanisms are still badly known at the molecular scale. In agreement, molecular dynamics simulations appear to be a method of choice to get a better description of these mechanisms. Nevertheless, the size and the time scales required for the simulation of these membrane systems limit such studies to laboratories having access to large computational facilities.The objective of this work was to predict and get a dynamical view of the interactions of several GPCRs with their partners, by developing an affordable molecular dynamics protocol that combines the coarse-grained MARTINI force field to Replica-Exchange MD simulations.In a first step, we validated our protocol by showing its ability to predict the dynamical binding of peptides to their receptors, through the study of Neurotensin, an agonist of the Neurotensin-1 receptor and CVX15, an antagonist of the CXCR4 chemokine receptor. We also show that the same protocol is able to predict the selectivity of several Neurotensin derived peptides against several wild-type/mutated receptors differing by a single residue.In a second step, we were concerned by the dynamical assembly of a GPCR heterodimer involving the Ghrelin and the Dopamine D2 receptors, respectively coupled to Gq and Gi proteins. Our model was validated by LRET measurements confirming a large protein:protein interface and a high complementarity between G-proteins. Based on this model, we designed and synthesized some peptides able to inhibit the assembly of this G-proteins heterodimer.Finally, we describe other applications of our protocol and how it can be employed and confronted to experiments to : predict the dynamical binding of toxins from snake’s venom to the Vasopressin-1a and Vasopressin-2 receptors ; predict the binding of the Ghrelin and Leap2 peptides to their GHSR-1a receptor and predict the coupling selectivity of several receptors to peptides mimicking the C-terminus of the α subunit of G-proteins.