Thesis (Ph. D.)--University of Washington, 1999. / Vita. Includes bibliographical references (leaves [188]-211).

Thesis (Ph. D.)--University of Texas at Austin, 2000. / Vita. Includes bibliographical references (leaves 94-99). Available also in a digital version from Dissertation Abstracts.

No description available.

Membrane proteins perform a variety of functions essential for the viability of the cell, including transport and signalling across the membrane. Most membrane proteins are formed from bundles of transmembrane helices. In this thesis molecular dynamics simulations have been used to investigate helix insertion into bilayers and helix association within bilayers. The potentials of mean force for the insertion of helices derived from the cystic fibrosis transmembrane conductance regulator into lipid bilayers were calculated using coarse-grained molecular dynamics simulations. The results showed that the insertion free energy increased with helix length and bilayer hydrophobic width. The insertion free energies obtained were significantly larger than comparable quantities obtained from translocon- mediated insertion experiments, consistent with a variety of previous studies. The implications of this observation for the interpretation of in vivo translocon-mediated insertion experiments, and the function of the translocon, are discussed. Coarse-grained and atomistic molecular dynamics simulations of the transmembrane region of the receptor tyrosine kinase EphA1 suggested that the transmembrane helix dimer was most stable when interacting via the glycine zipper motif, in agreement with a structure obtained by NMR spectroscopy. Coarse-grained simulations of the transmembrane region of EphA2 suggested that the dimer has two stable orientations, interacting via a glycine zipper or a heptad motif. Both structures showed right-handed dimers, although an NMR structure of the transmembrane region of EphA2 shows a left-handed dimer interacting via the heptad motif. Both structures obtained from coarse-grained simulations proved unstable when simulated at an atomistic level of detail. The potentials of mean force for dissociating the EphA1 and EphA2 dimers were calcu- lated using coarse-grained molecular dynamics calculations. Convergence of the detailed structure of the profiles was not conclusively shown, although association free energies cal- culated from the profiles were consistent over a variety of simulation times. The association free energies were slightly larger than experimental values obtained for comparable sys- tems, but consistent with similar computational calculations previously reported. However, direct comparisons are difficult owing to the influence of environmental factors on reported association free energies. The potential of mean force profiles showed that the interaction via the glycine zipper motif for EphA1 was significantly more stable than any other confor- mation. For EphA2 the potential of mean force profiles suggested that interaction via the glycine zipper and heptad motifs both provided stable or metastable conformations, with the interaction via the glycine zipper motif probably at least as stable as that via the heptad motif.

Thesis (Ph. D.)--University of Washington, 2002. / Vita. Includes bibliographical references (leaves 141-161).

Proteins are one of the most studied biological molecules of the last decades. A great amount of experimental techniques provide to researchers direct or indirect informations on proteins structure and function. In silico simulations can be used as a “computational microscope” giving the possibility to observe protein dynamic properties at atomistic resolution. In this work, various applications of computational methods to biological systems are presented. In particular, all-atom Molecular Dynamics (MD) simulations were employed to investigate the behaviour of proteins at atomstic resolution. The term “Molecular Dynamics” is usually referred to computational methods used for the simulation of classical many-body systems. These techniques are applied to microscopic systems and they represent a powerful approach for the study of physical processes, providing a tool for their interpretation. They have been widely used in the past decades to elucidate a large variety of molecular processes in different fields such as solid state physics, material science, chemistry, biochemistry and biophysics. Here, all-atom MD simulations were employed to observe equilibrium properties of several biologically relevant proteins. This allowed us to direct perform a comparison of molecular mechanisms occurring at the atomistic level as obtained from in silico studies with experimental data, which usually describe processes at larger length and time scales. These MD simulations were also meant as a starting point for the construction of simplified models, as they were processed through coarse-graining procedures to extrapolate crucial systems features, such as informative protein sites, on the basis of information theory approaches. Specifically we studied the dynamics of pembrolizumab, a humanized immunoglobulin of type G4 (IgG4) used as a therapeutic antibody. It is employed for the treatment of lung cancer, melanoma, stomach and head cancer and Hodgkin’s lymphoma. This antibody interacts with the programmed cell death protein 1 (PD-1) receptor, blocking the suppression of the immune response during cancer development. The studied systems are three: the apo state of pembrolizumab, the holo state (i.e. pembrolizumab bound to PD-1) and the glycosylated apo configuration. Each configuration was simulated for 2μs, for a total of 6μs. The analysis of the trajectories was carried out by combining standard structural analysis techniques and information theory-based measures of correlation. From MD trajectories we could extract valuable informations on the connectivity that exists among the structural domains that compose the antibody structure. Moreover, it was possible to infer which regions are involved in the structural rearrangement in the case of the antigen binding. We could observe that the presence of the antigen reduces the conformational variability of the molecule giving a greater stability to it. The second studied system is the P53 protein complex. In this case we focused on the tetramerization domain (TD) region that is composed by 2 identical dimers and has the function of bringing together the four monomers of the p53 complex. Starting from the observation that in case of the mutation of residue R337 several pathologies are developed in humans, we constructed computational models to reproduce the dynamics of the mutants and investigate their behaviour in silico. We performed simulations for a total of 16 μs divided in 8 different cases. In the first part of the study the wild type (WT) protein was compared to the R337C and the R337H mutant in three different protonation states: delta protonated Histidine, epsilon protonated Histidine ad double protonated Histidine. In the second part of the study we highlighted the differences between the WT configuration and three rationally designed mutants: R337D-352D, 337R-D352R, R337D-D352R. In this part of the investigation, the importance of the electrostatic interaction between residues R337 and D352 in the stability of the tetramerization do- main was discussed. Furthermore, we matched the obtained computational results of p53 tetramerization domain with functional experiments in yeasts (performed in collaboration with the CIBIO department) of all the simulated forms. The third simulated protein is the zinc sensing transcriptional repressor (CzrA), an homodimeric protein that binds DNA in Staphylococcus aureus. All-atom MD simulations of two different configurations were performed for a total of 4μs, the first one is the WT apo protein while the second is the WT holo system, where the protein is complexed with two Zn ions. In this case, in addition to standard analysis techniques, we applied the mapping entropy minimization protocol to highlight the most informative protein regions, from the perspective of information theory. Finally, our in silico results were compared to available NMR data of the protein itself.

No description available.

Biological processes are so complicated that to understand the mechanisms underlying the functioning of biomolecules it is inevitable to study them from various perspectives and with a wide range of tools. Understanding the function at the molecular level obviously requires the knowledge of the three dimensional structure of the biomolecules. Experimentally this can be obtained by techniques such as X‐ray crystallography and NMR studies. Computational biology has also played an important role in elucidating the structure function relationship in biomolecules. Computationally one can obtain the temporal as well as ensemble behavior of biomolecules at atomic level under conditions that are experimentally not accessible. Molecular dynamics(MD) study is a technique that can be used to obtain information of the dynamic behavior of the biomolecules. Dynamics of large systems like proteins can be investigated by classical force fields. However, the changes at the level of covalent bond involve the reorganization of electron density distribution which can be addressed only at Quantum mechanical level. In the present thesis, some of the biological systems have been characterized both at the classical and quantum mechanical level. The systems investigated by MD simulations and the insights brought from these studies are presented in Chapters 3 and 4. The unusual bonds such as pyrophosphate linkage in ATP and short strong hydrogen bonds in proteins, investigated through high level quantum chemical methods, are presented in Chapters 5, 6 and 7. Part of this thesis is aimed to address some important issues related to the dynamics of Tryptophanyl tRNA synthetase (TrpRS) which belongs to classic of aminoacyl‐tRNA synthetases (aaRS). aaRSs are extremely important class of enzymes involved in the translation of genetic code. These enzymes catalyze the aminoacylation of tRNAs to relate the cognate amino acids to the anticodon trinucleotide sequences. aaRSs are modular enzymes with distinct domains on which extensive kinetic and mutational experiments as well as structural analyses have been carried out, highlighting the role of inter‐domain communication (Alexander and Schimmel, 2001). The overall architecture of tRNA synthetases consists of primarily two domains. The active site domain is responsible for the activation of an amino acid with ATP in synthesizing an enzyme‐bound aminoacyl‐adenylate, and transfer of the aminoacyl‐adenylate intermediate to the 3’end of tRNA. The second domain is responsible for selection and binding of the cognate tRNA. aaRSs are allosteric proteins in which the binding of tRNA at the anticodon domain influences the activity at the catalytic region. These two binding sites are separated by a large distance. One of the aims of this thesis is to characterize such long distance communication (allosteric communication) at atomic level in Tryptophanyl tRNA synthetase. This is achieved by generating ensembles of conformations by MD simulations and analyzing the trajectories by novel graph theoretic approach. Graph and network based approaches are well established in the field of protein structure analysis for analyzing protein structure, stability and function (Kannan and Vishveshwara, 1999; Brinda and Vishveshwara, 2005). The parameters such as clusters, hubs and shortest paths provide valuable information on the structure and dynamics of the proteins. In this thesis, network parameters are used for the analysis of molecular dynamics MD) simulation data, to represent the global dynamic behavior of protein in a more elegant way. MD simulations are performed on some available (and modeled) structures of TrpRS bound to a variety of ligands, and the protein structure networks( PSN) of non‐covalent interactions are characterized in dynamical equilibrium. The ligand induced conformational changes are investigated through structure networks. These networks are used to understand the mode of communication between the anticodon domain and the active site. The interface dynamics is crucial for the function of TrpRS (since it is a functional dimer) and it is investigated through interface clusters. The matter embodied in the thesis is presented as 9 chapters. Chapter 1 lays the suitable background and foundation for the study, surveying relevant literature from different fields .Chapter 2 describes in detail the various materials, methods and techniques employed in the different analyses and studies presented in this thesis. A brief description of well‐known methods of molecular dynamics simulations, essential dynamics calculations, cross correlation maps, conformational clustering etc.is presented. The methods for constructing protein structure graphs and networks, developed in our lab, are described in detail. The use of network parameters for the analysis of MD simulation data to address the problem of communication between the two distal sites is also presented. Some descriptions of the ab initio quantum mechanical methods, which are used to investigate the unusual bonds in biomolecules, are also presented in this chapter. Chapter 3 is devoted in discussing the results from several normal as well as high temperature MD simulations of ligand‐free and ligand bound Bacillus stearothermophilus Tryptophanyl‐tRNA synthetase (bsTrpRS). The essential modes of the protein in the presence of different ligands are captured by essential dynamics calculations. Different conformations of the protein associated with the catalysis process of TrpRS, as captured through experiments, are discussed in the context of conformational sampling. High temperature simulations are carried out to explore the larger conformational space. Chapter 4 is focused on the results obtained from the MD simulation of human Tryptophanyl‐tRNA synthetase (hTrpRS). The structure of human TrpRS bound to the activated ligand (TrpAMP) and the cognate tRNA(tRNATRP) is modeled since no structure in the presence of both TrpAMP and tRNATRP is available. MD simulations on these modeled as well as other complexes of hTrpRS are performed to capture the dynamical process of ligand induced conformational changes (Hansiaetal., communicated). Both the local and the global changes in the protein conformation from the protein structure network (PSN) of MD snapshots are analyzed. Several important information such as the ligand induced correlation between different residues of the protein, asymmetric binding of the ligands to the two subunits of the protein, and the path of communication between the anticodon region and the aminoacylation site are obtained. Also, the role of the dimmer interface, from a dynamic perspective, is obtained for the first time. The interface dynamics which stabilize different quaternary structures of lectins (with high sequence and structure similarity) were investigated in a collaborative work (Hansiaetal.,2007). The lectin peanut agglutinin (PNA) is a tetramer with three different types of interfaces. The interface dynamics of this protein in the presence and in the absence of metal ions was investigated and the paper reporting the results from this study is included as appendix in this thesis. Chapter 5 deals with high level ab initio quantum chemical calculations on tri‐ and diphosphate fragments of adenosine triphosphate (ATP). Pyrophosphate prototypes such as methyl triphosphate and methyl diphosphate molecules in their different protonation states have been investigated at high levels of calculations (Hansiaetal., 2006a). The optimized geometries, the thermochemistry of the hydrolysis and the molecular orbitals contributing to the high energy of these compounds have been analyzed. These investigations provide insights into the‘‘highenergy’’character of ATP molecule. Further, the dependence of vibrational frequencies on the number of phosphate groups and the charged states has also been presented. These results aid in the interpretation of spectra obtained by experiments on complexes containing pyrophosphate prototypes. Hydrogen bonding is fundamental in understanding the structure and properties of molecules of biological interest including proteins. A recent analysis carried out in our lab showed that a significant number of short hydrogen bonds (SHB) are present in proteins (Rajagopal and Vishveshwara, 2005). Chapters 6 and 7 elucidate the results obtained from ab initio quantum chemical calculations on some of these SHBs to get aquantitative estimation of their geometry and strength. In chapter 6, asystematic analysis of the geometries and the energetics of possible SHB systems, which are frequently encountered in proteins, are presented at different levels of theory (HF,DFTandMP2). It is found that the SHBs involving both charged residues in the proteins are intrinsic in nature. However, two neutral residues form a SHB in the protein crystal structures either due to geometric constraints or due to the environment of these residues. This analysis enables one to distinguish SHBs which are formed because of geometric constraints from those which are formed because of the inherent property of the chemical groups involved in the hydrogen bonding. These results are useful in refining protein structures determined by crystallographic or NMR methods. In addition, sulfur atom of methionine and cysteinein proteins also participate in SHBs, which are not so well characterized. Chapter 7 presents the similar analysis carried out on short hydrogen bonds in proteins involving sulfur atom. A detailed analysis of SHBs of sulfur containing groups in a data set of proteins has been carried out. Some of the residue pairs from this analysis were considered for ab initio calculations. However, the optimization of these examples resulted in breaking of the hydrogen bonds involving sulfur atoms and formation of new hydrogen bonds with oxygen and/or nitrogen atoms. Hence model systems, which mimic the real examples, were designed to carry out ab initio studies and to investigate the short hydrogen bonds involving sulfur atoms. Another study on the protein‐water interaction, which does not fall under the realm of the main objective of the thesis, is discussed in Chapter 8. Protein–water interaction is crucial for accomplishing many biological functions of proteins. In the recent past, natural probe tryptophan, located at the protein surfaces, has been extensively investigated using femtosecond spectroscopy experiments to understand salvation dynamics (Peonetal.,2002). In this chapter a method is described to follow up the molecular events of the protein–water interactions in detail. Tryptophan–water interaction in the protein Monellin is investigated in order to get the atomic level insights into the hydration dynamics, by carrying out MD simulations on Monellin (Hansiaetal.,2006b). The results are compared with those obtained from femtosecond resolved fluorescence spectroscopy. The time constants of the survival correlation function match well with the reported experimental values.This validates the procedure, adapted here for Monellin, to investigate the hydration dynamics in general. The last chapter (Chapter9) summarizes the results obtained from various studies and discusses the future directions. First part of this thesis aims to present the analysis by carrying out MD simulations on monomeric and dimeric TrpRS protein in order to understand the two steps of the aminoacylation reaction: activation of the aminoacid Trp in the first step and the transfer of the activated amino acid in the next step. In the second part, quantitative estimation of the geometry and the strength of pyrophosphate bond and short hydrogen bonds in proteins are reported in detail by subjecting the systems to high levels of quantum mechanical calculations(QM). The use of ab initio QM/MM calculations by combining the quantum mechanics(QM) with the molecular mechanics(MM) in order to study the enzymatic reactions is discussed as the future

The selection of the forcefield is a crucial issue in any MD related work and there is no clear indication as to which of the many available forcefields is the best for protein analysis. Many recent literature surveys indicate that MD work may be hindered by two limitations, namely conformational sampling and forcefields used (inaccuracies in the potential energy function may bias the simulation toward incorrect conformations). However, the advances in computing infrastructures, theoretical and computing aspects of MD have paved the way to carry out a sampling on a sufficiently longtime scale, putting a need for the accuracies in the forcefield. Because there are established differences in MD results when using forcefields, we have sought to ask how we could assess common mobility segments from a protein by analysis of trajectories using three forcefields in a similar environment. This is important because, disparate fluctuations appear to be more at flexible regions compared to stiff regions; in particular, flexible regions are more relevant to functional activities of the protein molecule. Therefore, we have tried to assess the similarity in the dynamics using three well-known forcefields ENCAD, CHARMM27 and AMBERFF99SB for 61 monomeric proteins and identify the properties of dynamic residues, which may be important for function. The comparison of popular forcefields with different parameterization philosophy may give hints to improve some of the currently existing agnostics in forcefields and characterization of mobile regions based on dynamics of proteins with diverse folds. These may also give some signature on the proteins at the level of dynamics in relation to function, which can be used in protein engineering studies. Nanosecond level MD simulation(30ns) on 61 monomeric proteins were carried out using CHARMM and AMBER forcefields and the trajectories with ENCAD forcefield obtained from Dynameomics database. The trajectories were first analyzed to check whether structural and dynamic properties from the three forcefields similar choosing few parameters in each case. The gross dynamic properties calculated (root mean square deviation (RMSD), TM-score derived RMSD, radius of gyration and accessible surface area) indicated similarity in many proteins. Flexibility index analysis on 17 proteins, which showed a notable difference in the flexibility, indicated that tertiary interactions (fraction of nonnative stable hydrogen bonds and salt bridges) might be responsible for the difference in the flexibility index. The normalized subspace overlap and shape overlap score taken based on the covariance matrices derived from trajectories indicated that majority of the proteins show a range between 0.3-0.5 indicating that the first principal components from these proteins in different combinations may not match well. These results indicate that although dynamic properties in general are similar in many proteins. However, flexibility index and normalized subspace overlap score indicate that subspaces on the first principal component in many proteins may not match completely. The number of proteins showing a better correlation is higher in CHARMM-AMBER combinations than the other two. The structural features from trajectories have been computed in terms of fraction of secondary structure, hydrogen bonds, salt bridges and native contacts. Although secondary structures and native contacts are well preserved during the simulations, the tertiary interactions (hydrogen bonds) are lost in many proteins and may be responsible for the difference in the some of properties among forcefields. Comparison of simulation results to experimental structures in terms of Root mean square fluctuations, Accessible surface area and radius of gyration indicates that the simulations results are on par with the ones derived from experimental structures. We have tried to assess the flexibility in the proteins using normalized Root mean square fluctuations (nRMSF), which for a residue is the ratio of RMSF from simulation to that of crystal structure. We have selected a threshold for this nRMSF to indicate the mobile regions in a protein based on secondary structure analysis. Based on the threshold of nRMSF and conformational properties (deviation in the dihedral angles), we have classified the residue and evaluated the properties of rigid hinge residues and corresponding mobile residues in terms of residue propensity, secondary structure preference and accessible surface area ranges. Since the rigid dynamic residues represent the inherent mobility, they might be important for function. Therefore, we have tried to assess the functional relevance considering the dynamic mobile residues from each protein from each forcefield simulation with the residues important for the function (taken from literature and databases). It is observed that some residues found to be mobile from the simulation are found to match with the experimental ones, although in many cases the number of these mobile residues is higher compared to the experimental ones. In summary, an analysis of protein simulation trajectories using three forcefields on a set of monomeric protein has shown that the gross structural properties and secondary structures from many proteins remain similar, but there are differences as may be seen from flexibility index. However correlation in parameters from CHARMM and AMBER force field is better compared to other two combinations. The differences seen in some of structural properties may arise mainly due to the loss of few tertiary interactions as indicated by the fraction of native hydrogen bonds and salt bridges. Based on the nRMSF, mobile segments obtained from the simulations were identified, and some of the mobile segments are found to match the functionally important residues from the experimental ones. Our work indicates that there are still some differences in the properties from the simulations, which indicates that care must be exercised when choosing a forcefield, especially assessing the functionally relevant residues from the simulations.