Search

In situ structures of Platinum (Pt) nanoparticles (NPs) can be determined with graphene liquid cell transmission electron microscopy. Atomic-scale three-dimensional structural information about their physiochemical properties in solution is critical for understanding their chemical function. We here analyze eight atomic-resolution maps of small (<3 nm) colloidal Pt NPs. Their structures are composed of an ordered crystalline core surrounded by surface atoms with comparatively high mobility. 3D reconstructions calculated from cumulative doses of 8500 and 17,000 electrons/pixel, respectively, are characterized in terms of loss of atomic densities and atomic displacements. Less than 5% of the total number of atoms are lost due to dissolution or knock-on damage in five of the structures analyzed, whereas 10-16% are lost in the remaining three. Less than 5% of the atomic positions are displaced due to the increased electron irradiation in all structures. The surface dynamics will play a critical role in the diverse catalytic function of Pt NPs and must be considered in efforts to model Pt NP function computationally., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Determining the 3D atomic structures of multi-element nanoparticles in their native liquid environment is crucial to understanding their physicochemical properties. Graphene liquid cell (GLC) TEM offers a platform to directly investigate nanoparticles in their solution phase. Moreover, exploiting high-resolution TEM images of single rotating nanoparticles in GLCs, 3D atomic structures of nanoparticles are reconstructed by a method called "Brownian one-particle reconstruction". We here introduce a 3D atomic structure determination method for multi-element nanoparticle systems. The method, which is based on low-pass filtration and initial 3D model generation customized for different types of multi-element systems, enables reconstruction of high-resolution 3D Coulomb density maps for ordered and disordered multi-element systems and classification of the heteroatom type. Using high-resolution image datasets obtained from TEM simulations of PbSe, CdSe, and FePt nanoparticles that are structurally relaxed with first-principles calculations in the graphene liquid cell, we show that the types and positions of the constituent atoms are precisely determined with root mean square displacement values less than 24 pm. Our study suggests that it is possible to investigate the 3D atomic structures of synthesized multi-element nanoparticles in liquid phase., (© 2023. The Author(s).)

Analysis of the three-dimensional (3D) structures of nanocrystals with solution-phase transmission electron microscopy is beginning to reveal their unique physiochemical properties. We developed a "one-particle Brownian 3D reconstruction method" based on imaging of ensembles of colloidal nanocrystals using graphene liquid cell electron microscopy. Projection images of differently rotated nanocrystals are acquired using a direct electron detector with high temporal (<2.5 ms) resolution and analyzed to obtain an ensemble of 3D reconstructions. Here, we introduce computational methods required for successful atomic-resolution 3D reconstruction: (i) tracking of the individual particles throughout the time series, (ii) subtraction of the interfering background of the graphene liquid cell, (iii) identification and rejection of low-quality images, and (iv) tailored strategies for 2D/3D alignment and averaging that differ from those used in biological cryo-electron microscopy. Our developments are made available through the open-source software package SINGLE., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

We here introduce the third major release of the SIMPLE (Single-particle IMage Processing Linux Engine) open-source software package for analysis of cryogenic transmission electron microscopy (cryo-EM) movies of single-particles (Single-Particle Analysis, SPA). Development of SIMPLE 3.0 has been focused on real-time data processing using minimal CPU computing resources to allow easy and cost-efficient scaling of processing as data rates escalate. Our stream SPA tool implements the steps of anisotropic motion correction and CTF estimation, rapid template-based particle identification and 2D clustering with automatic class rejection. SIMPLE 3.0 additionally features an easy-to-use web-based graphical user interface (GUI) that can be run on any device (workstation, laptop, tablet or phone) and supports a remote multi-user environment over the network. The new project-based execution model automatically records the executed workflow and represents it as a flow diagram in the GUI. This facilitates meta-data handling and greatly simplifies usage. Using SIMPLE 3.0, it is possible to automatically obtain a clean SP data set amenable to high-resolution 3D reconstruction directly upon completion of the data acquisition, without the need for extensive image processing post collection. Only minimal standard CPU computing resources are required to keep up with a rate of ∼300 Gatan K3 direct electron detector movies per hour. SIMPLE 3.0 is available for download from simplecryoem.com., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Authors.)

Precise three-dimensional (3D) atomic structure determination of individual nanocrystals is a prerequisite for understanding and predicting their physical properties. Nanocrystals from the same synthesis batch display what are often presumed to be small but possibly important differences in size, lattice distortions, and defects, which can only be understood by structural characterization with high spatial 3D resolution. We solved the structures of individual colloidal platinum nanocrystals by developing atomic-resolution 3D liquid-cell electron microscopy to reveal critical intrinsic heterogeneity of ligand-protected platinum nanocrystals in solution, including structural degeneracies, lattice parameter deviations, internal defects, and strain. These differences in structure lead to substantial contributions to free energies, consequential enough that they must be considered in any discussion of fundamental nanocrystal properties or applications., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Motivation: No rigorous statistical tests for detecting point-group symmetry in three-dimensional (3D) charge density maps obtained by electron microscopy (EM) and related techniques have been developed., Results: We propose a method for determining the point-group symmetry of 3D charge density maps obtained by EM and related techniques. Our ab initio algorithm does not depend on atomic coordinates but utilizes the density map directly. We validate the approach for a range of publicly available single-particle cryo-EM datasets. In straightforward cases, our method enables fully automated single-particle 3D reconstruction without having to input an arbitrarily selected point-group symmetry. When pseudo-symmetry is present, our method provides statistics quantifying the degree to which the 3D density agrees with the different point-groups tested., Availability and Implementation: The software is freely available at https://github.com/hael/SIMPLE3.0., (© The Author(s) 2019. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Cryogenic electron microscopy (cryo-EM) and single-particle analysis enables determination of near-atomic resolution structures of biological molecules. However, large computational requirements limit throughput and rapid testing of new image processing tools. We developed PRIME, an algorithm part of the SIMPLE software suite, for determination of the relative 3D orientations of single-particle projection images. PRIME has primarily found use for generation of an initial ab initio 3D reconstruction. Here we show that the strategy behind PRIME, iterative estimation of per-particle orientation distributions with stochastic hill climbing, provides a competitive approach to near-atomic resolution single-particle 3D reconstruction. A number of mathematical techniques for accelerating the convergence rate are introduced, leading to a speedup of nearly two orders of magnitude. We benchmarked our developments on numerous publicly available data sets and conclude that near-atomic resolution ab initio 3D reconstructions can be obtained with SIMPLE in a matter of hours, using standard over-the-counter CPU workstations., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Protein dynamics manifested through structural flexibility play a central role in the function of biological molecules. Here we explore the substrate-mediated change in protein flexibility of an antibiotic target enzyme, Clostridium botulinum dihydrodipicolinate synthase. We demonstrate that the substrate, pyruvate, stabilizes the more active dimer-of-dimers or tetrameric form. Surprisingly, there is little difference between the crystal structures of apo and substrate-bound enzyme, suggesting protein dynamics may be important. Neutron and small-angle X-ray scattering experiments were used to probe substrate-induced dynamics on the sub-second timescale, but no significant changes were observed. We therefore developed a simple technique, coined protein dynamics-mass spectrometry (ProD-MS), which enables measurement of time-dependent alkylation of cysteine residues. ProD-MS together with X-ray crystallography and analytical ultracentrifugation analyses indicates that pyruvate locks the conformation of the dimer that promotes docking to the more active tetrameric form, offering insight into ligand-mediated stabilization of multimeric enzymes., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Cryogenic electron microscopy (cryo-EM) and single-particle analysis now enables the determination of high-resolution structures of macromolecular assemblies that have resisted X-ray crystallography and other approaches. We developed the SIMPLE open-source image-processing suite for analysing cryo-EM images of single-particles. A core component of SIMPLE is the probabilistic PRIME algorithm for identifying clusters of images in 2D and determine relative orientations of single-particle projections in 3D. Here, we extend our previous work on PRIME and introduce new stochastic optimization algorithms that improve the robustness of the approach. Our refined method for identification of homogeneous subsets of images in accurate register substantially improves the resolution of the cluster centers and of the ab initio 3D reconstructions derived from them. We now obtain maps with a resolution better than 10 Å by exclusively processing cluster centers. Excellent parallel code performance on over-the-counter laptops and CPU workstations is demonstrated., (© 2017 The Protein Society.)

A critical step in the analysis of novel cryogenic electron microscopy (cryo-EM) single-particle datasets is the identification of homogeneous subsets of images. Methods for solving this problem are important for data quality assessment, ab initio 3D reconstruction, and analysis of population diversity due to the heterogeneous nature of macromolecules. Here we formulate a stochastic algorithm for identification of homogeneous subsets of images. The purpose of the method is to generate improved 2D class averages that can be used to produce a reliable 3D starting model in a rapid and unbiased fashion. We show that our method overcomes inherent limitations of widely used clustering approaches and proceed to test the approach on six publicly available experimental cryo-EM datasets. We conclude that, in each instance, ab initio 3D reconstructions of quality suitable for initialization of high-resolution refinement are produced from the cluster centers., (Crown Copyright © 2016. Published by Elsevier Ltd. All rights reserved.)

Pore Forming Toxins (PFTs) represent a key mechanism for permitting the passage of proteins and small molecules across the lipid membrane. These proteins are typically produced as soluble monomers that self-assemble into ring-like oligomeric structures on the membrane surface. Following such assembly PFTs undergo a remarkable conformational change to insert into the lipid membrane. While many different protein families have independently evolved such ability, members of the Membrane Attack Complex PerForin/Cholesterol Dependent Cytolysin (MACPF/CDC) superfamily form distinctive giant β-barrel pores comprised of up to 50 monomers and up to 300Å in diameter. In this review we focus on recent advances in understanding the structure of these giant MACPF/CDC pores as well as the underlying molecular mechanisms leading to their formation. Commonalities and evolved variations of the pore forming mechanism across the superfamily are discussed. This article is part of a Special Issue entitled: Pore-Forming Toxins edited by Mauro Dalla Serra and Franco Gambale., (Copyright © 2015 Elsevier B.V. All rights reserved.)

The membrane attack complex (MAC)/perforin-like protein complement component 9 (C9) is the major component of the MAC, a multi-protein complex that forms pores in the membrane of target pathogens. In contrast to homologous proteins such as perforin and the cholesterol-dependent cytolysins (CDCs), all of which require the membrane for oligomerisation, C9 assembles directly onto the nascent MAC from solution. However, the molecular mechanism of MAC assembly remains to be understood. Here we present the 8 Å cryo-EM structure of a soluble form of the poly-C9 component of the MAC. These data reveal a 22-fold symmetrical arrangement of C9 molecules that yield an 88-strand pore-forming β-barrel. The N-terminal thrombospondin-1 (TSP1) domain forms an unexpectedly extensive part of the oligomerisation interface, thus likely facilitating solution-based assembly. These TSP1 interactions may also explain how additional C9 subunits can be recruited to the growing MAC subsequent to membrane insertion.

The lethal factor in stonefish venom is stonustoxin (SNTX), a heterodimeric cytolytic protein that induces cardiovascular collapse in humans and native predators. Here, using X-ray crystallography, we make the unexpected finding that SNTX is a pore-forming member of an ancient branch of the Membrane Attack Complex-Perforin/Cholesterol-Dependent Cytolysin (MACPF/CDC) superfamily. SNTX comprises two homologous subunits (α and β), each of which comprises an N-terminal pore-forming MACPF/CDC domain, a central focal adhesion-targeting domain, a thioredoxin domain, and a C-terminal tripartite motif family-like PRY SPla and the RYanodine Receptor immune recognition domain. Crucially, the structure reveals that the two MACPF domains are in complex with one another and arranged into a stable early prepore-like assembly. These data provide long sought after near-atomic resolution insights into how MACPF/CDC proteins assemble into prepores on the surface of membranes. Furthermore, our analyses reveal that SNTX-like MACPF/CDCs are distributed throughout eukaryotic life and play a broader, possibly immune-related function outside venom.

The growing problem of antibiotic resistance underlies the critical need to develop new treatments to prevent and control resistant bacterial infection. Exogenous application of bacteriophage lysins results in rapid and specific destruction of Gram-positive bacteria and therefore lysins represent novel antibacterial agents. The PlyC phage lysin is the most potent lysin characterized to date and can rapidly lyse Group A, C and E streptococci. Previously, we have determined the X-ray crystal structure of PlyC, revealing a complicated and unique arrangement of nine proteins. The scaffold features a multimeric cell-wall docking assembly bound to two catalytic domains that communicate and work synergistically. However, the crystal structure appeared to be auto-inhibited and raised important questions as to the mechanism underlying its extreme potency. Here we use small angle X-ray scattering (SAXS) and reveal that the conformational ensemble of PlyC in solution is different to that in the crystal structure. We also investigated the flexibility of the enzyme using both normal mode (NM) analysis and molecular dynamics (MD) simulations. Consistent with our SAXS data, MD simulations show rotational dynamics of both catalytic domains, and implicate inter-domain communication in achieving a substrate-ready conformation required for enzyme function. Our studies therefore provide insights into how the domains in the PlyC holoenzyme may act together to achieve its extraordinary potency.

Membrane attack complex/perforin-like (MACPF) proteins comprise the largest superfamily of pore-forming proteins, playing crucial roles in immunity and pathogenesis. Soluble monomers assemble into large transmembrane pores via conformational transitions that remain to be structurally and mechanistically characterised. Here we present an 11 Å resolution cryo-electron microscopy (cryo-EM) structure of the two-part, fungal toxin Pleurotolysin (Ply), together with crystal structures of both components (the lipid binding PlyA protein and the pore-forming MACPF component PlyB). These data reveal a 13-fold pore 80 Å in diameter and 100 Å in height, with each subunit comprised of a PlyB molecule atop a membrane bound dimer of PlyA. The resolution of the EM map, together with biophysical and computational experiments, allowed confident assignment of subdomains in a MACPF pore assembly. The major conformational changes in PlyB are a ∼70° opening of the bent and distorted central β-sheet of the MACPF domain, accompanied by extrusion and refolding of two α-helical regions into transmembrane β-hairpins (TMH1 and TMH2). We determined the structures of three different disulphide bond-trapped prepore intermediates. Analysis of these data by molecular modelling and flexible fitting allows us to generate a potential trajectory of β-sheet unbending. The results suggest that MACPF conformational change is triggered through disruption of the interface between a conserved helix-turn-helix motif and the top of TMH2. Following their release we propose that the transmembrane regions assemble into β-hairpins via top down zippering of backbone hydrogen bonds to form the membrane-inserted β-barrel. The intermediate structures of the MACPF domain during refolding into the β-barrel pore establish a structural paradigm for the transition from soluble monomer to pore, which may be conserved across the whole superfamily. The TMH2 region is critical for the release of both TMH clusters, suggesting why this region is targeted by endogenous inhibitors of MACPF function.

Cholesterol Dependent Cytolysins (CDCs) are important bacterial virulence factors that form large (200-300 Å) membrane embedded pores in target cells. Currently, insights from X-ray crystallography, biophysical and single particle cryo-Electron Microscopy (cryo-EM) experiments suggest that soluble monomers first interact with the membrane surface via a C-terminal Immunoglobulin-like domain (Ig; Domain 4). Membrane bound oligomers then assemble into a prepore oligomeric form, following which the prepore assembly collapses towards the membrane surface, with concomitant release and insertion of the membrane spanning subunits. During this rearrangement it is proposed that Domain 2, a region comprising three β-strands that links the pore forming region (Domains 1 and 3) and the Ig domain, must undergo a significant yet currently undetermined, conformational change. Here we address this problem through a systematic molecular modeling and structural bioinformatics approach. Our work shows that simple rigid body rotations may account for the observed collapse of the prepore towards the membrane surface. Support for this idea comes from analysis of published cryo-EM maps of the pneumolysin pore, available crystal structures and molecular dynamics simulations. The latter data in particular reveal that Domains 1, 2 and 4 are able to undergo significant rotational movements with respect to each other. Together, our data provide new and testable insights into the mechanism of pore formation by CDCs.

The human neuroendocrine enzyme glutamate decarboxylase (GAD) catalyses the synthesis of the inhibitory neurotransmitter gamma-aminobutyric acid (GABA) using pyridoxal 5'-phosphate as a cofactor. GAD exists as two isoforms named according to their respective molecular weights: GAD65 and GAD67. Although cytosolic GAD67 is typically saturated with the cofactor (holoGAD67) and constitutively active to produce basal levels of GABA, the membrane-associated GAD65 exists mainly as the inactive apo form. GAD65, but not GAD67, is a prevalent autoantigen, with autoantibodies to GAD65 being detected at high frequency in patients with autoimmune (type 1) diabetes and certain other autoimmune disorders. The significance of GAD65 autoinactivation into the apo form for regulation of neurotransmitter levels and autoantibody reactivity is not understood. We have used computational and experimental approaches to decipher the nature of the holo → apo conversion in GAD65 and thus, its mechanism of autoinactivation. Molecular dynamics simulations of GAD65 reveal coupling between the C-terminal domain, catalytic loop, and pyridoxal 5'-phosphate-binding domain that drives structural rearrangement, dimer opening, and autoinactivation, consistent with limited proteolysis fragmentation patterns. Together with small-angle X-ray scattering and fluorescence spectroscopy data, our findings are consistent with apoGAD65 existing as an ensemble of conformations. Antibody-binding kinetics suggest a mechanism of mutually induced conformational changes, implicating the flexibility of apoGAD65 in its autoantigenicity. Although conformational diversity may provide a mechanism for cofactor-controlled regulation of neurotransmitter biosynthesis, it may also come at a cost of insufficient development of immune self-tolerance that favors the production of GAD65 autoantibodies.

Bacteriophages deploy lysins that degrade the bacterial cell wall and facilitate virus egress from the host. When applied exogenously, these enzymes destroy susceptible microbes and, accordingly, have potential as therapeutic agents. The most potent lysin identified to date is PlyC, an enzyme assembled from two components (PlyCA and PlyCB) that is specific for streptococcal species. Here the structure of the PlyC holoenzyme reveals that a single PlyCA moiety is tethered to a ring-shaped assembly of eight PlyCB molecules. Structure-guided mutagenesis reveals that the bacterial cell wall binding is achieved through a cleft on PlyCB. Unexpectedly, our structural data reveal that PlyCA contains a glycoside hydrolase domain in addition to the previously recognized cysteine, histidine-dependent amidohydrolases/peptidases catalytic domain. The presence of eight cell wall-binding domains together with two catalytic domains may explain the extraordinary potency of the PlyC holoenyzme toward target bacteria.

Motivation: The β-barrel is a ubiquitous fold that is deployed to accomplish a wide variety of biological functions including membrane-embedded pores. Key influences of β-barrel lumen diameter include the number of β-strands (n) and the degree of shear (S), the latter value measuring the extent to which the β-sheet is tilted within the β-barrel. Notably, it has previously been reported that the shear value for small antiparallel β-barrels (n≤24) typically ranges between n and 2n. Conversely, it has been suggested that the β-strands in giant antiparallel β-barrels, such as those formed by pore forming cholesterol-dependent cytolysins (CDC), are parallel relative to the axis of the β-barrel, i.e. S=0. The S=0 arrangement, however, has never been observed in crystal structures of small β-barrels. Therefore, the structural basis for how CDCs form a β-barrel and span a membrane remains to be understood., Results: Through comparison of molecular models with experimental data, we are able to identify how giant CDC β-barrels utilize a 'near parallel' arrangement of β-strands where S=n/2. Furthermore, we show how side-chain packing within the β-barrel lumen is an important limiting factor with respect to the possible shear values for small β-barrels (n≤24 β-strands). In contrast, our models reveal no such limitation restricts the shear value of giant β-barrels (n>24 β-strands). Giant β-barrels can thus access a different architecture compared with smaller β-barrels.

Dihydrodipicolinate synthase (DHDPS) is an essential enzyme involved in the lysine biosynthesis pathway. DHDPS from E. coli is a homotetramer consisting of a 'dimer of dimers', with the catalytic residues found at the tight-dimer interface. Crystallographic and biophysical evidence suggest that the dimers associate to stabilise the active site configuration, and mutation of a central dimer-dimer interface residue destabilises the tetramer, thus increasing the flexibility and reducing catalytic efficiency and substrate specificity. This has led to the hypothesis that the tetramer evolved to optimise the dynamics within the tight-dimer. In order to gain insights into DHDPS flexibility and its relationship to quaternary structure and function, we performed comparative Molecular Dynamics simulation studies of native tetrameric and dimeric forms of DHDPS from E. coli and also the native dimeric form from methicillin-resistant Staphylococcus aureus (MRSA). These reveal a striking contrast between the dynamics of tetrameric and dimeric forms. Whereas the E. coli DHDPS tetramer is relatively rigid, both the E. coli and MRSA DHDPS dimers display high flexibility, resulting in monomer reorientation within the dimer and increased flexibility at the tight-dimer interface. The mutant E. coli DHDPS dimer exhibits disorder within its active site with deformation of critical catalytic residues and removal of key hydrogen bonds that render it inactive, whereas the similarly flexible MRSA DHDPS dimer maintains its catalytic geometry and is thus fully functional. Our data support the hypothesis that in both bacterial species optimal activity is achieved by fine tuning protein dynamics in different ways: E. coli DHDPS buttresses together two dimers, whereas MRSA dampens the motion using an extended tight-dimer interface.

The crystal structures of unliganded and liganded pMHC molecules provide a structural basis for TCR recognition yet they represent 'snapshots' and offer limited insight into dynamics that may be important for interaction and T cell activation. MHC molecules HLA-B*3501 and HLA-B*3508 both bind a 13 mer viral peptide (LPEP) yet only HLA-B*3508-LPEP induces a CTL response characterised by the dominant TCR clonetype SB27. HLA-B*3508-LPEP forms a tight and long-lived complex with SB27, but the relatively weak interaction between HLA-B*3501-LPEP and SB27 fails to trigger an immune response. HLA-B*3501 and HLA-B*3508 differ by only one amino acid (L/R156) located on α2-helix, but this does not alter the MHC or peptide structure nor does this polymorphic residue interact with the peptide or SB27. In the absence of a structural rationalisation for the differences in TCR engagement we performed a molecular dynamics study of both pMHC complexes and HLA-B*3508-LPEP in complex with SB27. This reveals that the high flexibility of the peptide in HLA-B*3501 compared to HLA-B*3508, which was not apparent in the crystal structure alone, may have an under-appreciated role in SB27 recognition. The TCR pivots atop peptide residues 6-9 and makes transient MHC contacts that extend those observed in the crystal structure. Thus MD offers an insight into 'scanning' mechanism of SB27 that extends the role of the germline encoded CDR2α and CDR2β loops. Our data are consistent with the vast body of experimental observations for the pMHC-LPEP-SB27 interaction and provide additional insights not accessible using crystallography.

Dihydrodipicolinate synthase (DHDPS) catalyzes the rate limiting step in lysine biosynthesis in bacteria and plants. The structure of DHDPS has been determined from several bacterial species and shown in most cases to form a homotetramer or dimer of dimers. However, only one plant DHDPS structure has been determined to date from the wild tobacco species, Nicotiana sylvestris (Blickling et al. (1997) J. Mol. Biol. 274, 608-621). Whilst N. sylvestris DHDPS also forms a homotetramer, the plant enzyme adopts a 'back-to-back' dimer of dimers compared to the 'head-to-head' architecture observed for bacterial DHDPS tetramers. This raises the question of whether the alternative quaternary architecture observed for N. sylvestris DHDPS is common to all plant DHDPS enzymes. Here, we describe the structure of DHDPS from the grapevine plant, Vitis vinifera, and show using analytical ultracentrifugation, small-angle X-ray scattering and X-ray crystallography that V. vinifera DHDPS forms a 'back-to-back' homotetramer, consistent with N. sylvestris DHDPS. This study is the first to demonstrate using both crystal and solution state measurements that DHDPS from the grapevine plant adopts an alternative tetrameric architecture to the bacterial form, which is important for optimizing protein dynamics as suggested by molecular dynamics simulations reported in this study.

Background: Canonical serine protease inhibitors commonly bind to their targets through a rigid loop stabilised by an internal hydrogen bond network and disulfide bond(s). The smallest of these is sunflower trypsin inhibitor (SFTI-1), a potent and broad-range protease inhibitor. Recently, we re-engineered the contact β-sheet of SFTI-1 to produce a selective inhibitor of kallikrein-related peptidase 4 (KLK4), a protease associated with prostate cancer progression. However, modifications in the binding loop to achieve specificity may compromise structural rigidity and prevent re-engineered inhibitors from reaching optimal binding affinity., Methodology/principal Findings: In this study, the effect of amino acid substitutions on the internal hydrogen bonding network of SFTI were investigated using an in silico screen of inhibitor variants in complex with KLK4 or trypsin. Substitutions favouring internal hydrogen bond formation directly correlated with increased potency of inhibition in vitro. This produced a second generation inhibitor (SFTI-FCQR Asn(14)) which displayed both a 125-fold increased capacity to inhibit KLK4 (K(i) = 0.0386±0.0060 nM) and enhanced selectivity over off-target serine proteases. Further, SFTI-FCQR Asn(14) was stable in cell culture and bioavailable in mice when administered by intraperitoneal perfusion., Conclusion/significance: These findings highlight the importance of conserving structural rigidity of the binding loop in addition to optimising protease/inhibitor contacts when re-engineering canonical serine protease inhibitors.

Background: The C-terminal domain of MotB (MotB-C) shows high sequence similarity to outer membrane protein A and related peptidoglycan (PG)-binding proteins. It is believed to anchor the power-generating MotA/MotB stator unit of the bacterial flagellar motor to the peptidoglycan layer of the cell wall. We previously reported the first crystal structure of this domain and made a puzzling observation that all conserved residues that are thought to be essential for PG recognition are buried and inaccessible in the crystal structure. In this study, we tested a hypothesis that peptidoglycan binding is preceded by, or accompanied by, some structural reorganization that exposes the key conserved residues., Methodology/principal Findings: We determined the structure of a new crystalline form (Form B) of Helicobacter pylori MotB-C. Comparisons with the existing Form A revealed conformational variations in the petal-like loops around the carbohydrate binding site near one end of the β-sheet. These variations are thought to reflect natural flexibility at this site required for insertion into the peptidoglycan mesh. In order to understand the nature of this flexibility we have performed molecular dynamics simulations of the MotB-C dimer. The results are consistent with the crystallographic data and provide evidence that the three loops move in a concerted fashion, exposing conserved MotB residues that have previously been implicated in binding of the peptide moiety of peptidoglycan., Conclusion/significance: Our structural analysis provides a new insight into the mechanism by which MotB inserts into the peptidoglycan mesh, thus anchoring the power-generating complex to the cell wall.

Currently, over a hundred high-resolution structures of serpins are available, exhibiting a wide range of conformations. However, our understanding of serpin dynamics and conformational change is still limited, mainly due to challenges of monitoring structural changes and characterizing transient conformations using experimental methods. Insight can be provided, however, by employing theoretical and computational approaches. In this chapter, we present an overview of such methods, focusing on molecular dynamics and simulation. As serpin conformational dynamics span a wide range of timescales, we discuss the relative merits of each method and suggest which method is suited to specific conformational phenomena., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Proteases are tightly regulated by specific inhibitors, such as serpins, which are able to undergo considerable and irreversible conformational changes in order to trap their targets. There has been a considerable effort to investigate serpin structure and functions in the past few decades; however, the specific interactions between proteases and serpins remain elusive. In this chapter, we describe detailed experimental protocols to determine and characterize the extended substrate specificity of proteases based on a substrate phage display technique. We also describe how to employ a bioinformatics system to analyze the substrate specificity data obtained from this technique and predict the potential inhibitory serpin partners of a protease (in this case, the immune protease, granzyme B) in a step-by-step manner. The method described here could also be applied to other proteases for more generalized substrate specificity analysis and substrate discovery., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Background: The crystallographic determination of protein structures can be computationally demanding and for difficult cases can benefit from user-friendly interfaces to high-performance computing resources. Molecular replacement (MR) is a popular protein crystallographic technique that exploits the structural similarity between proteins that share some sequence similarity. But the need to trial permutations of search models, space group symmetries and other parameters makes MR time- and labour-intensive. However, MR calculations are embarrassingly parallel and thus ideally suited to distributed computing. In order to address this problem we have developed MrGrid, web-based software that allows multiple MR calculations to be executed across a grid of networked computers, allowing high-throughput MR., Methodology/principal Findings: MrGrid is a portable web based application written in Java/JSP and Ruby, and taking advantage of Apple Xgrid technology. Designed to interface with a user defined Xgrid resource the package manages the distribution of multiple MR runs to the available nodes on the Xgrid. We evaluated MrGrid using 10 different protein test cases on a network of 13 computers, and achieved an average speed up factor of 5.69., Conclusions: MrGrid enables the user to retrieve and manage the results of tens to hundreds of MR calculations quickly and via a single web interface, as well as broadening the range of strategies that can be attempted. This high-throughput approach allows parameter sweeps to be performed in parallel, improving the chances of MR success.

Background: A central tenet of structural biology is that related proteins of common function share structural similarity. This has key practical consequences for the derivation and analysis of protein structures, and is exploited by the process of "molecular sieving" whereby a common core is progressively distilled from a comparison of two or more protein structures. This paper reports a novel web server for "sieving" of protein structures, based on the multiple structural alignment program MUSTANG., Methodology/principal Findings: "Sieved" models are generated from MUSTANG-generated multiple alignment and superpositions by iteratively filtering out noisy residue-residue correspondences, until the resultant correspondences in the models are optimally "superposable" under a threshold of RMSD. This residue-level sieving is also accompanied by iterative elimination of the poorly fitting structures from the input ensemble. Therefore, by varying the thresholds of RMSD and the cardinality of the ensemble, multiple sieved models are generated for a given multiple alignment and superposition from MUSTANG. To aid the identification of structurally conserved regions of functional importance in an ensemble of protein structures, Lesk-Hubbard graphs are generated, plotting the number of residue correspondences in a superposition as a function of its corresponding RMSD. The conserved "core" (or typically active site) shows a linear trend, which becomes exponential as divergent parts of the structure are included into the superposition., Conclusions: The application addresses two fundamental problems in structural biology: first, the identification of common substructures among structurally related proteins--an important problem in characterization and prediction of function; second, generation of sieved models with demonstrated uses in protein crystallographic structure determination using the technique of Molecular Replacement.