New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
 Agricultural Sciences
 Anthropology
 Applied Biological Sciences
 Biochemistry
 Biophysics and Computational Biology
 Cell Biology
 Developmental Biology
 Ecology
 Environmental Sciences
 Evolution
 Genetics
 Immunology and Inflammation
 Medical Sciences
 Microbiology
 Neuroscience
 Pharmacology
 Physiology
 Plant Biology
 Population Biology
 Psychological and Cognitive Sciences
 Sustainability Science
 Systems Biology
Electrostatics of nanosystems: Application to microtubules and the ribosome

Communicated by Peter G. Wolynes, University of California at San Diego, La Jolla, CA (received for review May 8, 2001)
Abstract
Evaluation of the electrostatic properties of biomolecules has become a standard practice in molecular biophysics. Foremost among the models used to elucidate the electrostatic potential is the PoissonBoltzmann equation; however, existing methods for solving this equation have limited the scope of accurate electrostatic calculations to relatively small biomolecular systems. Here we present the application of numerical methods to enable the trivially parallel solution of the PoissonBoltzmann equation for supramolecular structures that are orders of magnitude larger in size. As a demonstration of this methodology, electrostatic potentials have been calculated for large microtubule and ribosome structures. The results point to the likely role of electrostatics in a variety of activities of these structures.
The importance of electrostatic modeling to biophysics is well established; electrostatics have been shown to influence various aspects of nearly all biochemical reactions. Advances in NMR, xray, and cryoelectron microscopy techniques for structure elucidation have drastically increased the size and number of biomolecules and molecular complexes for which coordinates are available. However, although the biophysical community continues to examine macromolecular systems of increasing scale, the computational evaluation of electrostatic properties for these systems is limited by methodology that can handle only relatively small systems, typically consisting of fewer than 50,000 atoms. Despite these limitations, such computational methods have been immensely useful in analyses of the stability, dynamics, and association of proteins, nucleic acids, and their ligands (1–3). Here we describe algorithms that open the way to similar analyses of much larger subcellular structures.
One of the most widespread models for the evaluation of electrostatic properties is the PoissonBoltzmann equation (PBE) (4, 5) 1 a secondorder nonlinear elliptic partial differential equation that relates the electrostatic potential (φ) to the dielectric properties of the solute and solvent (ɛ), the ionic strength of the solution and the accessibility of ions to the solute interior (κ̄^{2}), and the distribution of solute atomic partial charges (f). To expedite solution of the equation, this nonlinear PBE is often approximated by the linearized PBE (LPBE) by assuming sinhφ(x) ≈ φ(x). Several numerical techniques have been used to solve the nonlinear PBE and LPBE, including boundary element (6–8), finite element (9–11), and finite difference (12–14) algorithms. However, despite the variety of solution methods, none of these techniques has been satisfactorily applied to large molecular structures at the scales currently accessible to modern biophysical methods. To accommodate arbitrarily large biomolecules, algorithms for solving the PBE must be both efficient and amenable to implementation on a parallel platform in a scalable fashion, requirements that current methods have been unable to satisfy. Although boundary element LPBE solvers provide an efficient representation of the problem domain, they are not useful for the nonlinear problem and have not been applied to the PBE on parallel platforms. Similarly, adaptive finite element methods have shown some success in parallel evaluation of both the LPBE and nonlinear PBE (15), but limitations in current solver technology and difficulty with efficient representation of the biomolecular data prohibits their practical application to large biomolecular systems. Finally, unlike the boundary and finite element techniques, finite difference methods have the advantage of very efficient multilevel solvers (12, 16) and applicability to both the linear and nonlinear forms of the PBE; however, existing parallel finite difference algorithms often require costly interprocessor communication that limits both the nature and scale of their execution on parallel platforms (17–21) [see especially Van de Velde (19) for reviews of the various methods].
Multigrid Solution Through Parallel Focusing
Recently, a new algorithm (BankHolst) was described for the parallel adaptive finite element solution of elliptic partial differential equations with negligible interprocess communication (22). Unlike finite difference methods, which use a fixed resolution of the problem domain, adaptive finite element techniques generate very accurate solutions of these equations by locally enriching the basis set in regions of high error through refinement of the domain discretization. The first step in the BankHolst parallel finite element algorithm is the solution of the equation over the entire problem domain using a coarse resolution basis set by each processor. This solution is then used, in conjunction with an a posteriori error estimator, to partition the problem domain into P subdomains, which are assigned to the P processors of a parallel computer. The algorithm is load balanced by equidistributing a posteriori error estimates across the subdomains. Each processor then solves the partial differential equation over the global mesh, but confines its adaptive refinement to the local subdomain and a small surrounding overlap region. This procedure results in a very accurate representation of the solution over the local subdomain. After all processors have completed their local adaptive solution of the equation, a global solution is constructed by the piecewise assembly of the solutions in each subdomain. It can be rigorously shown that the piecewiseassembled solution is as accurate as the solution to the problem on a single global mesh (23). Because each processor performs all of its computational work independently, the BankHolst algorithm requires very little interprocess communication and exhibits excellent parallel scaling. Unfortunately, although this algorithm works well for adaptive techniques, such as finite elements, it is not directly applicable to the fixedresolution finite difference methods that currently offer the most efficient solution of the PBE for large molecular systems. However, by combining the BankHolst method with other techniques, we have been able to extend this algorithm to finite difference solvers.
Electrostatic “focusing” is a popular technique in finite difference methods for generating accurate solutions to the PBE in subsets of the problem domain, such as a binding or titratable sites within a protein (4, 14). Like the BankHolst method, the first step in electrostatic focusing is the calculation of a lowaccuracy solution on a coarse finite difference mesh spanning the entire problem domain. This coarse solution is then used to define the boundary conditions for a much more accurate calculation on a finer discretization of the desired subdomain. As noted previously (22), this focusing technique is superficially related to the BankHolst algorithm, where the local enrichment of a coarse, global solution has been replaced by a the solution of a fine, local multigrid problem using the solution from a coarse, global problem for boundary conditions.
We have combined standard focusing techniques and the BankHolst algorithm into a “parallel focusing” method for the solution of the PBE. Unlike previous parallel algorithms for solving the PBE, this method has excellent parallel complexity, permitting the treatment of very large biomolecular systems on massively parallel computational platforms. Furthermore, the finite difference discretization on a regular mesh allows for fast solution by certain highly efficient multigrid solvers (12). The algorithm is summarized below:
Given the problem data and P processors of a parallel machine:
(i) Each processor i = 1...,P (a) obtains a coarse solution of the PBE over the global domain, (b) heuristically subdivides the global domain into P subdomains (Ω_{1}...,Ω_{P}), which are assigned to processors 1...,P, (c) assigns boundary conditions to its subdomain Ω_{i} using the coarse global solution, and (d) solves the PBE on subdomain Ω_{i}.
(ii) A master processor collects observable data (free energy, etc.) from the other processors and controls I/O.
This parallel focusing algorithm begins with each processor independently solving a coarse global problem. The perprocessor subdomains are chosen in a heuristic fashion as outlined in Fig. 1. Like standard focusing calculations, only a subset of the global mesh surrounding the area of interest is used for the parallel calculations. This subset is partitioned into P approximately equal subdomains, which are distributed among the processors. Each processor then performs a finescale finite difference calculation over this subdomain and an overlap region that usually spans about 5–10% of the neighboring subdomains. The overlap regions are included to compensate for inaccuracies in the boundary conditions derived from the global coarse solution; however, these regions are not used to assemble the finescale global solution and do not contribute to calculations of observables such as forces and energies. After the finescale calculations are complete, a master processor accumulates the desired data from the other processes and controls I/O of the solution. Like the BankHolst algorithm, each process computes independently on both the global coarse problem and its subdomain. These independent computations result in an algorithm that requires negligible interprocess communication and offers excellent parallel performance.
Because of its low communication requirements, this parallel focusing algorithm can be used on parallel platforms ranging from lowbandwidth networks of workstations to highperformance supercomputers. Because step d of the parallel focusing algorithm typically dominates the overall computation, this algorithm provides excellent time complexity with respect to the number of processors. Most importantly for PoissonBoltzmann calculations on large molecules, the mesh resolution increases linearly with the number of processors.** This linear behavior has been observed on calculations involving hundreds of processors; Fig. 2 shows the scaling results for calculations on the biomolecular systems discussed below. The parallel focusing algorithm was tested on up to 700 processors and shows linear scaling over the entire range. The implications of the improved performance available through this parallel algorithm are discussed in more detail in the next section.
Parallel focusing has been implemented as an extension (23) of the apbs (Adaptive PoissonBoltzmann Solver) software package (9) by using the pmg multigrid library to assemble and solve the PBE algebraic systems (12, 16). The parallel routines have been implemented by using maloc, a portable hardware abstraction layer library (23), which allows interprocess communication through a variety of methods.
The following sections describe the application of the parallel focusing techniques to very large biomolecular systems whose electrostatic properties were unattainable by using traditional methods for solving the PBE. Specifically, we examine two important cellular components: a millionatom microtubule fragment and the large and small subunits of the ribosome.
Electrostatic Properties of Microtubules
Within every eukaryotic cell is a complex system of proteins and filaments called the cytoskeleton. The largest structures within the cytoskeleton are microtubules, cylindrical polymers formed from the protein tubulin. Microtubules perform many functions within the cell, from providing basic structure to cell division and transport, and many of these roles involve electrostatic interactions. The recent solution of the structure of tubulin (24) has allowed us to construct a working model of a microtubule containing 90 tubulin dimers. This structure contains more than 1.25 million atoms [the negatively charged C termini of both the α and β subunits were not resolved in the tubulin structure (24) and have not been rebuilt in our microtubule structure], and the structure has dimensions of ≈300 Å × 300 Å × 600 Å. The LPBE was solved at 0.54Å mesh spacing by using 686 processors of the National Partnership for Advanced Computational Infrastructure (NPACI) Blue Horizon supercomputer. Blue Horizon is a massively parallel IBM RS6000 supercomputer consisting of 1,152 375MHz Power3 processors distributed among 144 eightway symmetric multiprocessing nodes. These nodes are connected with the IBM Colony switch, which delivers 350 MB/sec peak messagepassing interface performance with 17ms latency. Following the parallel focusing algorithm, each processor solved the LPBE on a coarse 97^{3} global mesh, as well as a finer 97^{3} spanning a subset of the overall domain. The global solution was obtained in less than 1 h. This system provides an excellent example of the improved performance available through the parallel focusing method. A sequential PBE solver treating the same system at a similar resolution would discretize the problem on a 1,111 × 555 × 555 mesh, requiring more than 350 times more memory and time to solve.
Although more work is presently underway to elucidate the implications of electrostatics on microtubule structure and function, Fig. 3 shows some early results. The most obvious feature is the overall negative electrostatic potential of the microtubule with smaller regions of positive potential. This large negative potential likely plays a significant role in the interaction with microtubuleassociated proteins, motor proteins such as kinesin and dynein, and other microtubules. Fig. 3 b and c also exhibits dramatic differences between the + end (β tubulin monomer exposed) and − end (α tubulin monomer exposed) of the microtubule, which, as has been shown for actin filaments (3), may contribute to differences in the polymerization properties (25) and relative stabilities of the two ends. Fig. 4 shows a microtubule inner surface in crosssection. The binding regions of compounds such as vinblastine, colchicine, and taxol (26) are shown in Fig. 4a, while electrostatic isocontours are shown, for the same view, in Fig. 4b. Qualitative examination of the electrostatic data shows interesting variations in the drug binding regions; however, much more work needs to be done to quantitatively elucidate contribution of electrostatics to the interaction of these pharmaceuticals with microtubules.
The microtubule also provides a good example for examining the improved performance of this parallel algorithm with respect to sequential solvers.
Electrostatic Properties of the Ribosome
Ribosomes are macromolecular complexes responsible for the translation of the mRNA into protein. These complexes consist of two subunits: the large 50S subunit and the small 30S subunit, both of which are composed of RNA and protein constituents. During translation, the large and small subunits associate to form an active ribosome. Recently, highresolution structures of the 50S subunit from Haloarcula marismortui (27) and the 30S subunit from Thermus thermophilus (28, 29) were solved by xray crystallography. These structures present a unique opportunity to examine the electrostatic potential of large ribonucleoprotein complexes.
The protonated 30S structure consists of more than 88,000 atoms and roughly spans a 200Ålong cubic box. Using apbs on 343 processors of the NPACI Blue Horizon, LPBE was solved to give the electrostatic potential of the small ribosomal subunit at 0.41Å resolution. Calculations also were performed on the protonated 50S subunit structure, which consisted of 94,854 atoms‡‡ and has similar dimensions to the 30S subunit. The electrostatic potential of the large subunit was obtained at 0.45Å resolution by using apbs to solve the LPBE on 343 processors of the Blue Horizon. Both systems were solved by using the parallel focusing algorithm: each processor first solved a coarse problem defined by a 97^{3} mesh covering the entire problem domain, and the solved the LPBE on a 97^{3} mesh covering the particular subset of the global problem. Examination of the electrostatic potential mapped to the molecular surfaces of the 30S (Fig. 5) and 50S (Fig. 6) reveals large areas of negative potential with smaller regions of positive potential corresponding to some of the proteins implicated in binding of the 30S and 50S subunits to form the active ribosomal complex (27, 29–31). Not surprisingly, the potential surface maps between the two subunits exhibit qualitative electrostatic complementarity. Despite the fact that the 30S and 50S structures are from distinct species that thrive in very different environments, such comparisons of electrostatic potential are not unwarranted. Formation of active ribosomal complexes has been observed by using subunits from different species (32, 33), indicating some evolutionary conservation of the elements involved in subunit association.
Fig. 5 b and d shows the electrostatic potential mapped to the 30S subunit molecular surface. As evidenced by the 30S structure (Fig. 5 a and c), regions of positive potential typically correspond to protein components of the small subunit or cocrystallized counterions. Of particular interest is the active site of the 30S subunit, shown in Fig. 5 a and b. There are interesting variations in the electrostatic potential near regions of antibiotic binding (34) and the A, P, and E tRNA binding sites (29, 34). For example, the codon binding sites and “platform region” (protuberance on the right side of Fig. 5 a and b) are surrounded by regions of positive potential, which could help stabilize complexes of the ribosome with mRNA and tRNA. However, more quantitative calculations will be needed to fully explore these hypotheses.
The electrostatic potential of the 50S subunit is shown in Fig. 6 b and d. Like the small subunit, the electrostatic surface potential is largely negative, with scattered regions of positive potential typically associated with ribosomal proteins. Some of the most interesting aspects of this data set are the regions of positive potential on the 50S “crown” (see upper protuberance in Fig. 6 a and b), which correspond to proteins of the large subunit involved in tRNA binding. Specifically, these positive regions are due to proteins L44e, L5, and L10e (see Fig. 6a), which have been implicated in tRNA binding to the 50S subunit (35), and contribute large regions of positive potential to the molecular surface in the upper portions of the A, P, and E tRNA binding sites. Furthermore, the Ploop (shown as blue spheres in Fig. 6a), an important component of the 50S Psite (35, 36), shows significant positive surface potential due to a nearby Mg^{2+} ion. These areas of positive potential may provide stabilizing interactions with mRNA or tRNA bound to the ribosomal complex during translation. However, as with the 30S calculations, more detailed work is required to accurately resolve the role of electrostatics in tRNA and mRNA binding to the 50S subunit.
Conclusions
We have described the combination of standard finite difference focusing techniques and the BankHolst algorithm into a parallel focusing method to facilitate the solution of the PBE for nanoscale systems. Unlike previous multiprocessor algorithms for solving the PBE, this method has excellent parallel complexity that permits the solution of these problems on massively parallel computational platforms. Furthermore, the finite difference discretization on a regular mesh allows for fast solution by highly efficient multigrid solvers. This algorithm has been implemented in the apbs software by using the pmg multigrid solver and tested on the NPACI IBM Blue Horizon supercomputer.
Using these methods for the parallel multigrid solution of elliptic differential equations, the electrostatic properties of very large biomolecular assemblages are now amenable to computation. This technique relies on the efficient solution of the PBE combined with parallel focusing techniques to solve these large problems in a variety of distributed computational environments. Solution of the LPBE for the 1.2 millionatom microtubule system provided electrostatic potential data, which revealed interesting features near drug binding sites and provided possible insight into stability differences at the + and − ends of the microtubule. Such detailed electrostatic information will be central to future studies that examine the possible collective effects involved in the formation of structural defects and the stabilizing effects of taxol binding to the interior of microtubules. Likewise, application of this methodology to ribosome systems elucidated intriguing electrostatic properties of the ribosomal active site, which will provide the starting point for investigation of the stabilization of the tRNA and mRNAribosome complexes during translation and the rational design of novel antibiotics. Finally, the ability to determine the contribution of electrostatics to the forces and energies of nanoscale systems should extend the scale of implicit solvent dynamics methods to much larger macromolecular complexes.
Acknowledgments
N.A.B. thanks A. Elcock and C. Wong for helpful discussions on the PBE and A. Majumdar and G. Chukkapalli for assistance with technical issues on the Blue Horizon supercomputer. N.A.B. was supported by predoctoral fellowships from the Howard Hughes Medical Institute and the BurroughsWellcome La Jolla Interfaces in Science program. This work was supported, in part, by an IBM/American Chemical Society award to N.A.B., providing time on the Minnesota Supercomputing Institute IBM SP2; by grants to J.A.M. from the National Institutes of Health, National Science Foundation, and NPACI/San Diego Supercompter Center; by National Science Foundation CAREER Award 9875856 to M.J.H.; and by National Science Foundation Grant MCB0078322 to S.J. Additional support has been provided by the W. M. Keck Foundation and the National Biomedical Computation Resource.
Footnotes

↵‡ To whom reprint requests should be addressed at: University of California at San Diego, 9500 Gilman Drive, Mail Code 0365, La Jolla, CA 920930365. Email: nbaker{at}mccammon.ucsd.edu.

↵** For example, let v(P) = h_{x}(P)h_{y}(P)h_{z}(P) be the volume of a grid element in a Pprocessor calculation with 100f% overlap on a n_{x} × n_{y} × n_{z} mesh with spacings h_{x}(P),h_{y}(P),h_{z}(P) over a problem domain of volume V. One measure of the mesh resolution is the inverse element volume, which behaves as v(P)^{−1} ∼(1 − 2f)^{3}P(n_{x} − 1)(n_{y} − 1)(n_{x}− 1)/V, which is linear in P.

↵‡‡ The 50S Protein Data Bank coordinate file (1FFK) (27) contained only Cα coordinates for the protein constituents. Therefore, each protein residue was simply represented by its Cα atom, which was assigned the total charge of the residue and a radius of 4.0 Å.
Abbreviations
 PBE,
 PoissonBoltzmann equation;
 LPBE,
 linearized PBE;
 NPACI,
 National Partnership for Advanced Computational Infrastructure
 Received May 8, 2001.
 Accepted July 5, 2001.
 Copyright © 2001, The National Academy of Sciences
References
 ↵
 Elcock A H,
 Sept D,
 McCammon J A
 ↵
 ↵
 Davis M E,
 McCammon J A
 ↵
 Honig B,
 Nicholls A
 ↵
 ↵
 ↵
 ↵
 ↵

 Davis M E,
 McCammon J A
 ↵
 ↵
 Baker N A,
 Sept D,
 Holst M J,
 McCammon J A
 ↵
 Keyes D E,
 Xu J
 Holst M,
 Saied F
 ↵
 Schwarz H A

 Bjørstad P E,
 Widlund O B
 ↵
 Van de Velde E F

 Ilin A,
 Bagheri B,
 Scott L R,
 Briggs J M,
 McCammon J A
 ↵
 Tuminaro R S,
 Womble D E
 ↵
 Bank R,
 Holst M
 ↵
Holst, M. (2001) Adv. Comput. Math., in press.
 ↵
 ↵
 Walker R A,
 O'Brien E T,
 Pryer N K,
 Soboeiro M F,
 Voter W A,
 Erikson H P,
 Salmon E D
 ↵
 Giannakakou P,
 Sackett D L,
 Kang YK,
 Zhan Z R,
 Buters J T M,
 Fojo T,
 Poruchynsky M S
 ↵
 Ban N,
 Nissen P,
 Hansen J,
 Moore P B,
 Steitz T A
 ↵
 ↵

 Cate J H,
 Yuspov M M,
 Yusupova G Z,
 Earnest T N,
 Noller H F
 ↵
 Culver G M,
 Cate J H,
 Yusupova G Z,
 Yusupov M M,
 Noller H F
 ↵
 Klein H A,
 Ochoa S
 ↵
 ↵
 ↵
 Nissen P,
 Hansen J,
 Ban N,
 Moore P B,
 Steitz T A
 ↵
Citation Manager Formats
More Articles of This Classification
Physical Sciences
Chemistry
Biological Sciences
Related Content
 No related articles found.
Cited by...
 Modeling and Mutational Analysis of the Binding Mode for the Multimodal Antidepressant Drug Vortioxetine to the Human 5HT3A Receptor
 NadA3 Structures Reveal Undecad Coiled Coils and LOX1 Binding Regions Competed by Meningococcus B VaccineElicited Human Antibodies
 Variants with increased negative electrostatic potential in the Cx50 gap junction pore increased unitary channel conductance and magnesium modulation
 Molecular mechanism of activation of the immunoregulatory amidase NAAA
 Glycoside hydrolase family 18 and 20 enzymes are novel targets of the traditional medicine berberine
 Roles of glycosaminoglycans as regulators of ligand/receptor complexes
 Structurefunction analyses of the bacterial zinc metalloprotease effector protein GtgA uncover key residues required for deactivating NF{kappa}B
 Pseudomonas aeruginosa pyoverdine maturation enzyme PvdP has a noncanonical domain architecture and affords insight into a new subclass of tyrosinases
 Structure of the {gamma}{varepsilon} complex of cyanobacterial F1ATPase reveals a suppression mechanism of the {gamma} subunit on ATP hydrolysis in phototrophs
 Redoxcoupled quinone dynamics in the respiratory complex I
 Sucrose Nonfermenting 1Related Protein Kinase 1 Phosphorylates a Geminivirus Rep Protein to Impair Viral Replication and Infection
 Guarding the gateway to histidine biosynthesis in plants: Medicago truncatula ATPphosphoribosyltransferase in relaxed and tense states
 Structure of the EmrE multidrug transporter and its use for inhibitor peptide design
 Bioinformatics analysis of CYP1B1 mutation hotspots in Chinese primary congenital glaucoma patients
 Amino Acid Substitutions within HLAB*27Restricted T Cell Epitopes Prevent Recognition by Hepatitis Delta VirusSpecific CD8+ T Cells
 Allatom simulations disentangle the functional dynamics underlying gene maturation in the intron lariat spliceosome
 Two distinct sites of client protein interaction with the chaperone cpSRP43
 Gating modifier toxins isolated from spider venom: Modulation of voltagegated sodium channels and the role of lipid membranes
 A chargesensing region in the stromal interaction molecule 1 luminal domain confers stabilizationmediated inhibition of SOCE in response to Snitrosylation
 Insights into the binding behavior of native and nonnative cytochromes to photosystem I from Thermosynechococcus elongatus
 Structure of an EIIC sugar transporter trapped in an inwardfacing conformation
 The highaffinity calcium sensor synaptotagmin7 serves multiple roles in regulated exocytosis
 The functional principle of eukaryotic molybdenum insertases
 The AB loop and Dhelix in binding site III of human Oncostatin M (OSM) are required for OSM receptor activation
 Probing the quality control mechanism of the Escherichia coli twinarginine translocase with folding variants of a de novodesigned heme protein
 Escape of Hepatitis C Virus from Epitope I Neutralization Increases Sensitivity of Other Neutralization Epitopes
 A Polyamide Inhibits Replication of Vesicular Stomatitis Virus by Targeting RNA in the Nucleocapsid
 CryoEM structure of a herpesvirus capsid at 3.1 A
 Structure of a eukaryotic cytoplasmic pre40S ribosomal subunit
 Shaping Substrate Selectivity in a BroadSpectrum Metallo{beta}Lactamase
 A peptidylic inhibitor for neutralizing expanded CAG RNAinduced nucleolar stress in polyglutamine diseases
 Structural Study of the CTerminal Domain of Nonstructural Protein 1 from Japanese Encephalitis Virus
 Molecular mechanism for inhibition of twinfilin by phosphoinositides
 A crossdomain charge interaction governs the activity of NO synthase
 Interaction between the Caenorhabditis elegans centriolar protein SAS5 and microtubules facilitates organelle assembly
 Structure and Membrane Topography of the VibrioType Secretin Complex from the Type 2 Secretion System of Enteropathogenic Escherichia coli
 An intrinsic lipidbinding interface controls sphingosine kinase 1 function
 CRISPR RNA and antiCRISPR protein binding to the Xanthomonas albilineans Csy1Csy2 heterodimer in the type IF CRISPRCas system
 Thiosulfate sulfurtransferaselike domaincontaining 1 protein interacts with thioredoxin
 PvdO is required for the oxidation of dihydropyoverdine as the last step of fluorophore formation in Pseudomonas fluorescens
 G4941K substitution in the porelining S6 helix of the skeletal muscle ryanodine receptor increases RyR1 sensitivity to cytosolic and luminal Ca2+
 Structural insights into ankyrin repeatmediated recognition of the kinesin motor protein KIF21A by KANK1, a scaffold protein in focal adhesion
 Probing the interaction between the histone methyltransferase/deacetylase subunit RBBP4/7 and the transcription factor BCL11A in epigenetic complexes
 An ultraprocessive, accurate reverse transcriptase encoded by a metazoan group II intron
 Molecular Basis for Necitumumab Inhibition of EGFR Variants Associated with Acquired Cetuximab Resistance
 StructureFunction Dissection of Pseudorabies Virus Glycoprotein B Fusion Loops
 Peptide Uptake Is Essential for Borrelia burgdorferi Viability and Involves Structural and Regulatory Complexity of its Oligopeptide Transporter
 Characterization of SPP inhibitors suppressing propagation of HCV and protozoa
 The periplasmic transaminase PtaA of Pseudomonas fluorescens converts the glutamic acid residue at the pyoverdine fluorophore to {alpha}ketoglutaric acid
 Crystal structures of Mmm1 and Mdm12Mmm1 reveal mechanistic insight into phospholipid trafficking at ERmitochondria contact sites
 The bacterial metacleavage hydrolase LigY belongs to the amidohydrolase superfamily, not to the {alpha}/{beta}hydrolase superfamily
 Structural analyses of the Haemophilus influenzae peptidoglycan synthase activator LpoA suggest multiple conformations in solution
 The intrinsically disordered Nterminal domain of galectin3 dynamically mediates multisite selfassociation of the protein through fuzzy interactions
 Asymmetric Trimeric Ring Structure of the Nucleocapsid Protein of Tospovirus
 TBK1 mutations in Italian patients with amyotrophic lateral sclerosis: genetic and functional characterisation
 Copper homeostasis networks in the bacterium Pseudomonas aeruginosa
 Pleiotropic Effects of ResistanceBreaking Mutations on Particle Stability Provide Insight into Life History Evolution of a Plant RNA Virus
 Structural insights into transcription initiation by yeast RNA polymerase I
 Structure of a symmetric photosynthetic reaction centerphotosystem
 Homeodomainlike DNA binding proteins control the haploidtodiploid transition in Dictyostelium
 Structure of the competence pilus major pilin ComGC in Streptococcus pneumoniae
 Structural and biochemical analyses indicate that a bacterial persulfide dioxygenaserhodanese fusion protein functions in sulfur assimilation
 Interaction between the PH and START domains of ceramide transfer protein competes with phosphatidylinositol 4phosphate binding by the PH domain
 Nearatomic resolution cryoelectron microscopy structure of the 30fold homooligomeric SpoIIIAG channel essential to spore formation in Bacillus subtilis
 An improved Akt reporter reveals intra and intercellular heterogeneity and oscillations in signal transduction
 Structural Insight into Nucleoprotein Conformation Change Chaperoned by VP35 Peptide in Marburg Virus
 Electrostatic control of DNA intersegmental translocation by the ETS transcription factor ETV6
 Mutation of the {alpha}tubulin Tuba1a leads to straighter microtubules and perturbs neuronal migration
 Singlemolecule nanopore enzymology
 Ionpulling simulations provide insights into the mechanisms of channel opening of the skeletal muscle ryanodine receptor
 Symmetryrelated proton transfer pathways in respiratory complex I
 MCTBI: a web server for predicting metal ion effects in RNA structures
 Transmembrane helix hydrophobicity is an energetic barrier during the retrotranslocation of integral membrane ERAD substrates
 Energetics and dynamics of a lightdriven sodiumpumping rhodopsin
 Structural studies of viperin, an antiviral radical SAM enzyme
 The Lon proteaselike domain in the bacterial RecA paralog RadA is required for DNA binding and repair
 Two glycosaminoglycanbinding domains of the mouse cytomegalovirusencoded chemokine MCK2 are critical for oligomerization of the fulllength protein
 Mechanistic insight into the nucleusvacuole junction based on the Vac8pNvj1p crystal structure
 TOGtubulin binding specificity promotes microtubule dynamics and mitotic spindle formation
 NLRP3 inflammasome assembly is regulated by phosphorylation of the pyrin domain
 A biochemical framework for eIF4Edependent mRNA export and nuclear recycling of the export machinery
 The RxLR Motif of the Host Targeting Effector AVR3a of Phytophthora infestans Is Cleaved before Secretion