Gas-phase analytical techniques provide very valuable tools for tackling the structural complexity of macromolecular structures such as those encountered in biological systems. Conformational dynamics of polypeptides and polypeptide assemblies underlie most biological functionalities, yet great difficulties arise when investigating such phenomena with the well-established techniques of X-ray crystallography and NMR. In areas such as these ion mobility interfaced with mass spectrometry (IMMS) and molecular modelling can make a significant contribution. During an IMMS experiment analyte ions drift in a chamber filled with an inert gas; measurement of the transport properties of analyte ions under the influence of a weak electric field can lead to determination of the orientationally-averaged collision cross-section of all resolved ionic species. A comparison with cross-sections estimated for model molecular geometries can lead to structural assignments. Thus IMMS can be used effectively to separate gas-phase ions based on their conformation. The drift tube employed in the experiments described herein is thermally regulated, which also enables the determination of collision cross-sections over a range of temperatures, and can provide a view of temperature-dependent conformational dynamics over the experimental (low microsecond) timescale. Studies described herein employ IMMS and a gamut of other MS-based techniques, solution spectroscopy and – importantly – molecular mechanics simulations to assess a) conformational stability of isolated peptide ions, with a focus on small model peptides and proteins, especially the Trp cage miniprotein; and b) structural characteristics of oligomeric aggregates of an amyloidogenic peptide. The results obtained serve to clarify the factors which dominate the intrinsic stability of non-covalent structure in isolated peptides and peptide assemblies. Strong electrostatic interactions are found to play a pivotal role in determining the conformations of isolated proteins. Secondary structures held together by hydrogen bonding, such as helices, are stable in the absence of solvent, however gas-phase protein structures display loss of their hydrophobic cores. The absence of a polar solvent, “self-solvation” is by far the most potent force influencing the gas-phase configuration of these systems. Geometries that are more compact than the folded state observed in solution are routinely detected, indicating the existence of intrinsically stable compact non-native states in globular proteins, illuminating the nature of proteins’ ‘unfolded’ states.

No / D-amino acids are useful building blocks for de novo peptide design and they play a role in aging-related diseases associated with gradual protein racemization. For amino acids with achiral side chains, one should be able to presume that the conformational propensities of L- and D-amino acids are a reflection of one another due to the straightforward geometric inversion at the Cα atom. However, this presumption does not account for the directionality of the backbone dipole and the inverted propensities have never been definitively confirmed in this context. Furthermore, there is little known of how alternative side chain chirality affects the backbone conformations of isoleucine and threonine. Using a GGXGG host-guest pentapeptide system, we have completed exhaustive sampling of the conformational propensities of the D-amino acids, including D-allo-isoleucine and D-allo-threonine, using atomistic molecular dynamics simulations. Comparison of these simulations with the same systems hosting the cognate L-amino acids verifies that the intrinsic backbone conformational propensities of the D-amino acids are the inverse of their cognate L-enantiomers. Where amino acids have a chiral center in their side chain (Thr, Ile) the β-configuration affects the backbone sampling, which in turn can confer different biological properties. / NIH

Doctor of Philosophy / Department of Biochemistry and Molecular Biophysics / Jianhan Chen / Intrinsically disordered proteins (IDPs) are functional proteins that lack stable tertiary structures under physiological conditions. IDPs are key components of regulatory networks that dictate various aspects of cellular decision-making, and are over-represented in major disease pathways. For example, about 30% of eukaryotic proteins contain intrinsic disordered regions, and over 70% of cancer-associated proteins have been identified as IDPs. The highly heterogeneous nature of IDPs has presented significant challenge for experimental characterization using NMR, X-ray crystallography, or FRET. These challenges represent a unique opportunity for molecular mod- eling to make critical contributions. In this study, computer simulations at multiple scales were utilized to characterize the structural properties of unbound IDPs as well as to obtain a mechanistic understanding of IDP interactions. These studies of IDPs also reveal significant limitations in the current simulation methodology. In particular, successful simulations of biomolecules not only require accurate molecular models, but also depend on the ability to sufficiently sample the com- plex conformational space. By designing a realistic yet computationally tractable coarse-grained protein model, we demonstrated that the popular temperature replica exchange enhanced sampling is ineffective in driving faster reversible folding transitions for proteins. The second original contribution of this dissertation is the development of novel simulation methods for enhanced sampling of protein conformations, specifically, replica exchange with guided-annealing (RE-GA) method and multiscale enhanced sampling (MSES) method. We expect these methods to be highly useful in generating converged conformational ensembles.

Cette thèse a pour but de développer le pendant in silico des étapes clés du Fragment-Based Drug Design (FBDD), et ce dans le cadre plus général du développement de l'outil S4MPLE. Le FBDD génère des ligands drug-like à partir de petites molécules (fragments). Après une étape de validation de S4MPLE et de sa fonction d’énergie, un recentrage autour du FBDD est réalisé, à travers le docking puis l'optimisation virtuelle de fragments par growing ou linking (G/L). Cette stratégie reposesur 1) la création d’une chimiothèque focalisée en connectant un ou deux fragment(s) avec des linkers pré-générés, et 2) l’échantillonnage avec S4MPLE des composés chimères dans le site avec des contraintes. Des simulations de G/L plus ou moins ambitieuses (site flexible, ajout de H2O libres) permettent de valider cette approche avec des études rétrospectives basées sur des données expérimentales. La dernière phase de la thèse a consisté à appliquer ce protocole in silico à un projet de l’entreprise. / This work aims to develop in silico methods targeting the key stages of Fragment-Based Drug Design (FBDD), participating to the development of the molecular modeling tool S4MPLE. Briefly, FBDD generates ıdrug-likeı ligands from small organic molecules called fragments. After a validation step of S4MPLE and its energy function, the work focused on FBDD: molecular docking of fragments and their subsequent virtual optimization. The latter mimics standard evolution strategies in FBDD(growing and linking). This in silico approach involves among other two key stages 1) building of a focused library by plugging in pre-generated linkers into reference fragments using rules and 2) sampling of these new compounds under atomic and binding site constraints. Validation simulations, relying on known experimental data, included ıclassicalı growing / linking and more challenging ones (site flexibility, free waters). Finally, this strategy is applied to one project of the company.

No description available.

No description available.

L'enjeu des travaux effectués au cours de cette thèse est l'extraction in silico de molécules potentiellement intéressantes dans le processus d'inhibition du récepteur tyrosine kinase c-Met. La faculté de cette protéine à interagir dans les phénomènes d'embryogenèse et de réparation tissulaires rendent son inhibition cruciale dans les traitements contre les développements tumoraux où c-Met se trouve impliquée. Pour cela, la stratégie employée implique l'utilisation de méthodes in silico de conception rationnelle de médicaments. Nous avons utilisé comme support les multiples structures cristallographiques publiées sur la ProteinData Base. Un travail de modélisation par homologie fut tout d'abord nécessaire pour combler les lacunes des structures cristallographiques collectées. Afin d'échantillonner au mieux l'espace conformationnel de la kinase c-Met et de caractériser sa flexibilité, une longue campagne de simulation de Dynamique Moléculaire fut menée. Pour compléter ces simulations, nous avons également utilisé la méthode des modes normaux de vibration. De ces 2 approches, nous avons extrait un ensemble de 10 conformères considérés comme les plus représentatifs de l'espace conformationnel simulé pour la kinase c-Met et avons proposé un mode de fonctionnement de ce récepteur. Utilisant les conformations représentatives, nous avons ensuite mené une importante campagne de criblage virtuel sur plusieurs chimiothèques constituant environ 70.000 composés. L'analyse des résultats de l'arrimage moléculaire nous a conduits à la sélection de composés intéressants possédant théoriquement une bonne affinité pour la kinase c-Met. Ces molécules ont été soumises aux tests expérimentaux. / The challenge of this PhD work is the in silico identification of potentially interesting molecules concerning the inhibitory process of tyrosine kinase receptor c-Met. The faculty of this protein to interact in embryogenesis and tissue repair phenomena makes its inhibition crucial for treatments against tumor development in which c-Met is involved. For that purpose, the employed strategy involves the use of several in silico methods for rational drug design. As the basement of this work, we used the multiple crystal structures published in the ProteinData Base (PDB). A preliminary homology modeling work was needed to fill gaps in the crystal structures. To sample at best the c-Met kinase conformational space and to characterize its flexibility, a long Molecular Dynamics (MD) simulation campaign was carried out both on apo and holo forms of available crystal structures. To complete these simulations, part of this work consisted to use normal modes of vibration (NM) method. From these two approaches (DM and NM), we extracted a set of 10 conformers considered as the most representative of the kinase simulated conformational space and we suggested a mode of operation of this kinase. Using extracted conformations from the conformational sampling has enabled us to conduct an extensive campaign on several virtual screening libraries constituting a total of approximately 70,000 compounds. Analysis of the molecular docking results has led us to the selection of several theoretically interesting molecules with good potential affinity for c-Met kinase. These molecules were submitted to experimental tests performed by the biologist team associated to our work.

No description available.