Leveraging crosslinking mass spectrometry in structural and cell biology

1. Download : Download high-res image (269KB)
2. Download : Download full-size image

     

Shape shifting: The multiple conformational substates of the PTEN N‐terminal PIP2 ‐binding domain

The Phosphatase and TENsin homolog deleted on chromosome 10 (PTEN) is a chief regulator of a variety of cellular processes including cell proliferation, migration, growth, and death. It is also a major tumor suppressor gene that is frequently mutated or lost under cancerous conditions. PTEN encodes a dual‐specificity (lipid and protein) phosphatase that negatively regulates the PI3K/AKT/mTOR signaling pathway where the PIP₂‐binding domain (PBD) regulates the lipid phosphatase function. Unfortunately, despite two decades of research, a full‐length structure of PTEN remains elusive, leaving open questions regarding PTEN''s disordered regions that mediate protein stability, post‐translational modifications, protein–protein interactions, while also hindering the design of small molecules that can regulate PTEN''s function. Here, we utilized a combination of crosslinking mass spectrometry, in silico predicted structural modeling (including AlphaFold2), molecular docking, molecular dynamics simulations, and residue interaction network modeling to obtain structural details and molecular insight into the behavior of the PBD of PTEN. Our study shows that the PBD exists in multiple conformations which suggests its ability to regulate PTEN''s variety of functions. Studying how these specific conformational substates contribute to PTEN function is imperative to defining its function in disease pathogenesis, and to delineate ways to modulate its tumor suppressor activity.     

Cross-linking and mass spectrometry as a tool for studying the structural biology of ribonucleoproteins

1. Download : Download high-res image (138KB)
2. Download : Download full-size image

     

Tightening the Crosslinking Distance Restraints for Better Resolution of Protein Structure and Dynamics

1. Download : Download high-res image (221KB)
2. Download : Download full-size image

     

Conformational dynamics in mixed alpha/beta-oligonucleotides containing polarity reversals: a molecular dynamics study using time-averaged restraints

Nucleic acid duplexes featuring a single alpha-anomeric thymidine inserted into each DNA strand via 3-3 and 5-5 phosphodiester linkages exhibit local conformational dynamics that are not adequately depicted by conventional restrained molecular dynamics (rMD) methods. We have used molecular dynamics with time-averaged NMR restraints (MDtar) to explore its applicability to describing the conformational dynamics of two -containing duplexes – d(GCGAAT-3-3-T-5-5-CGC)₂ and d(ATGG-3-3-T-5-5-GCTC)r(gagcaccau). In contrast to rMD, enforcing NOE-based distance restraints over a period of time in MDtar rather than instantaneously results in better agreement with the experimental NOE and J-data. This conclusion is based on the dramatic decreases in average distance and coupling constant violations (d _av, J _rms, and J _av) and improvements in sixth-root R-factors (R ^x). In both duplexes, the deoxyribose ring puckering behavior predicted independently by pseudorotation analysis is portrayed remarkably well using this approach compared to rMD. This indicates that the local dynamic behavior is encoded within the NOE data, although this is not obvious from the local R ^x values. In both systems, the backbone torsion angles comprising the 3-3 linkage as well as the (high S-) sugars of the -nucleotide and preceding residue (–1) are relatively static, while the conformations of the 5-5 linkage and the sugar in the neighboring -nucleotide (+1) show enhanced flexibility. To reduce the large ensembles generated by MDtar to more manageable clusters we utilized the PDQPRO program. The resulting PDQPRO clusters (in both cases, 13 structures and associated probabilities extracted from a pool of 300 structures) adequately represent the structural and dynamic characteristics predicted by the experimental data.     

One short cysteine‐rich sequence pattern – two different disulfide‐bonded structures – a molecular dynamics simulation study

The nematocyst walls of Hydra are formed by proteins containing small cysteine‐rich domains (CRDs) of ~25 amino acids. The first CRD of nematocyst outer all antigen (NW1) and the C‐terminal CRD of minicollagen‐1 (Mcol1C) contain six cysteines at identical sequence positions, however adopt different disulfide bonded structures. NW1 shows the disulfide connectivities C2‐C14/C6‐C19/C10‐C18 and Mcol1C C2‐C18/C6‐C14/C10‐C19. To analyze if both show structural preferences in the open, non‐disulfide bonded form, which explain the formation of either disulfide connectivity pattern, molecular dynamics (MD) simulations at different temperatures were performed. NW1 maintained in the 100‐ns MD simulations at 283 K a rather compact fold that is stabilized by specific hydrogen bonds. The Mcol1C structure fluctuated overall more, however stayed most of the time also rather compact. The analysis of the backbone Φ/ψ angles indicated different turn propensities for NW1 and Mcol1C, which mostly can be explained based on published data about the influence of different amino acid side chains on the local backbone conformation. Whereas a folded precursor mechanism may be considered for NW1, Mcol1C may fold according to the quasi‐stochastic folding model involving disulfide bond reshuffling and conformational changes, locking the native disulfide conformations. The study further demonstrates the power of MD simulations to detect local structural preferences in rather dynamic systems such as the open, non‐disulfide bonded forms of NW1 and Mcol1C, which complement published information from NMR backbone residual dipolar couplings. Because the backbone structural preferences encoded by the amino acid sequence embedding the cysteines influence which disulfide connectivities are formed, the data are generally interesting for a better understanding of oxidative folding and the design of disulfide stabilized therapeutics. Copyright © 2015 European Peptide Society and John Wiley & Sons, Ltd.     

SHAMS: Combining chemical modification of RNA with mass spectrometry to examine polypurine tract-containing RNA/DNA hybrids

Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) has gained popularity as a facile method of examining RNA structure both in vitro and in vivo, exploiting accessibility of the ribose 2′-OH to acylation by N-methylisatoic anhydride (NMIA) in unpaired or flexible configurations. Subsequent primer extension terminates at the site of chemical modification, and these products are fractionated by high-resolution gel electrophoresis. When applying SHAPE to investigate structural features associated with the wild-type and analog-substituted polypurine tract (PPT)–containing RNA/DNA hybrids, their size (20–25 base pairs) rendered primer extension impractical. As an alternative method of detection, we reasoned that chemical modification could be combined with tandem mass spectrometry, relying on the mass increment of RNA fragments containing the NMIA adduct (M_r = 133 Da). Using this approach, we demonstrate both specific modification of the HIV-1 PPT RNA primer and variations in its acylation pattern induced by replacing template nucleotides with a non-hydrogen-bonding thymine isostere. Our selective 2′-hydroxyl acylation analyzed by mass spectrometry strategy (SHAMS) should find utility when examining the structure of small RNA fragments or RNA/DNA hybrids where primer extension cannot be performed.     

Understanding the acylation mechanisms of active-site serine penicillin-recognizing proteins: a molecular dynamics simulation study

Beta-lactam antibiotics inhibit enzymes involved in the last step of peptidoglycan synthesis. These enzymes, also identified as penicillin-binding proteins (PBPs), form a long-lived acyl-enzyme complex with beta-lactams. Antibiotic resistance is mainly due to the production of beta-lactamases, which are enzymes that hydrolyze the antibiotics and so prevent them reaching and inactivating their targets, and to mutations of the PBPs that decrease their affinity for the antibiotics. In this study, we present a theoretical study of several penicillin-recognizing proteins complexed with various beta-lactam antibiotics. Hybrid quantum mechanical/molecular mechanical potentials in conjunction with molecular dynamics simulations have been performed to understand the role of several residues, and pK(a) calculations have also been done to determine their protonation state. We analyze the differences between the beta-lactamase TEM-1, the membrane-bound PBP2x of Streptococcus pneumoniae, and the soluble DD-transpeptidase of Streptomyces K15.     

Exploring the selectivity of a ligand complex with CDK2/CDK1: a molecular dynamics simulation approach

Cyclin‐dependent kinases (CDKs) are core components of the cell cycle machinery that govern the transition between phases during cell cycle progression. Abnormalities in CDKs activity and regulation are common features of cancer, making CDK family members attractive targets for the development of anticancer drugs. Their inhibitors have entered in clinical trials to treat cancer. Very recently, Heathcote et al. (J. Med. Chem. 2010, 53:8508–8522) have found a ligand BS194 that has a high affinity with CDK2 (IC₅₀ = 3 nm ) but shows low affinity with CDK1 (IC₅₀ = 30 nm ). To understand the selectivity, we used homology modeling, molecular docking, molecular dynamics, and free‐energy calculation to analyze the interactions. A rational three‐dimensional model of the CDK1/BS194 complex is built. We found that Leu83 is a key residue that recognizes BS194 more effectively with CDK2 with good binding free energies rather than CDK1. Energetic analysis reveals that van der Waals interaction and non‐polar contributions to solvent are favorable in the formation of complexes and amine group of the ligand, which plays a crucial role for binding selectivity between CDK2 and CDK1. Copyright © 2012 John Wiley & Sons, Ltd.     

Restrained molecular dynamics simulations of HIV-1 protease: the first step in validating a new target for drug design

To test the anticorrelated relationship that was recently displayed in conventional molecular dynamics (MD) simulations, several different restrained MD simulations on a wild type and on the V82F/I84V drug-resistant mutant of HIV-1 protease were performed. This anticorrelated relationship refers to the observation that compression of the peripheral ear-to-cheek region of HIV protease (i.e., the elbow of the flap to the fulcrum and the cantilever) occurred as the active site flaps were opening, and, conversely, expansion of that ear-to-cheek region occurred as both flaps were closing. An additional examination of this anticorrelated relationship was necessary to determine whether it can be harnessed in a useful manner. Consequently, six different MD experiments were performed that incorporated pairwise distance restraints in that ear-to-cheek region (i.e., the distance between the alpha-carbons of Gly40 and Gln61 was restrained to either 7.7 or 10.5 A, in both monomers). Pushing the backbones of the ear and the cheek regions away from each other slightly did force the flaps that guard the active site to remain closed in both the wild type and the mutant systems-even though there were no ligands in the active sites. Thus, these restrained MD simulations provided evidence that the anticorrelated relationship can be exploited to affect the dynamic behavior of the flaps that guard the active site of HIV-1 protease. These simulations supported our hypothesis of the mechanism governing flap motion, and they are the first step towards validating that peripheral surface as a new target for drug design.     

Using hydrogen deuterium exchange mass spectrometry to engineer optimized constructs for crystallization of protein complexes: Case study of PI4KIIIβ with Rab11

The ability of proteins to bind and interact with protein partners plays fundamental roles in many cellular contexts. X‐ray crystallography has been a powerful approach to understand protein‐protein interactions; however, a challenge in the crystallization of proteins and their complexes is the presence of intrinsically disordered regions. In this article, we describe an application of hydrogen deuterium exchange mass spectrometry (HDX‐MS) to identify dynamic regions within type III phosphatidylinositol 4 kinase beta (PI4KIIIβ) in complex with the GTPase Rab11. This information was then used to design deletions that allowed for the production of diffraction quality crystals. Importantly, we also used HDX‐MS to verify that the new construct was properly folded, consistent with it being catalytically and functionally active. Structures of PI4KIIIβ in an Apo state and bound to the potent inhibitor BQR695 in complex with both GTPγS and GDP loaded Rab11 were determined. This hybrid HDX‐MS/crystallographic strategy revealed novel aspects of the PI4KIIIβ‐Rab11 complex, as well as the molecular mechanism of potency of a PI4K specific inhibitor (BQR695). This approach is widely applicable to protein‐protein complexes, and is an excellent strategy to optimize constructs for high‐resolution structural approaches.     

Identification of a conserved 8 aa insert in the PIP5K protein in the Saccharomycetaceae family of fungi and the molecular dynamics simulations and structural analysis to investigate its potential functional role

Homologs of the phosphatidylinositol‐4‐phosphate‐5‐kinase (PIP5K), which controls a multitude of essential cellular functions, contain a 8 aa insert in a conserved region that is specific for the Saccharomycetaceae family of fungi. Using structures of human PIP4K proteins as templates, structural models were generated of the Saccharomyces cerevisiae and human PIP5K proteins. In the modeled S. cerevisiae PIP5K, the 8 aa insert forms a surface exposed loop, present on the same face of the protein as the activation loop of the kinase domain. Electrostatic potential analysis indicates that the residues from 8 aa conserved loop form a highly positively charged surface patch, which through electrostatic interaction with the anionic portions of phospholipid head groups, is expected to play a role in the membrane interaction of the yeast PIP5K. To unravel this prediction, molecular dynamics (MD) simulations were carried out to examine the binding interaction of PIP5K, either containing or lacking the conserved signature insert, with two different membrane lipid bilayers. The results from MD studies provide insights concerning the mechanistic of interaction of PIP5K with lipid bilayer, and support the contention that the identified 8 aa conserved insert in fungal PIP5K plays an important role in the binding of this protein with membrane surface. Proteins 2017; 85:1454–1467. © 2017 Wiley Periodicals, Inc.     

Intramolecular hydrogen bonding in articaine can be related to superior bone tissue penetration: a molecular dynamics study

Local anesthetics (LAs) are drugs that cause reversible loss of nociception during surgical procedures. Articaine is a commonly used LA in dentistry that has proven to be exceptionally effective in penetrating bone tissue and induce anesthesia on posterior teeth in maxilla and mandibula. In the present study, our aim was to gain a deeper understanding of the penetration of articaine through biological membranes by studying the interactions of articaine with a phospholipid membrane. Our approach involves Langmuir monolayer experiments combined with molecular dynamics simulations. Membrane permeability of LAs can be modulated by pH due to a titratable amine group with a pKa value close to physiological pH. A change in protonation state is thus known to act as a lipophilicity switch in LAs. Our study shows that articaine has an additional unique lipophilicity switch in its ability to form an intramolecular hydrogen bond. We suggest this intramolecular hydrogen bond as a novel and additional solvent-dependent mechanism for modulation of lipophilicity of articaine which may enhance its diffusion through membranes and connective tissue.