Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism.

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

Mechanism for the Regulated Control of Bacterial Transcription Termination by a Universal Adaptor Protein.

NusG/Spt5 proteins are the only transcription factors utilized by all cellular organisms. In enterobacteria, NusG antagonizes the transcription termination activity of Rho, a hexameric helicase, during the synthesis of ribosomal and actively translated mRNAs. Paradoxically, NusG helps Rho act on untranslated transcripts, including non-canonical antisense RNAs and those arising from translational stress; how NusG fulfills these disparate functions is unknown. Here, we demonstrate that NusG activates Rho by assisting helicase isomerization from an open-ring, RNA-loading state to a closed-ring, catalytically active translocase. A crystal structure of closed-ring Rho in complex with NusG reveals the physical basis for this activation and further explains how Rho is excluded from translationally competent RNAs. This study demonstrates how a universally conserved transcription factor acts to modulate the activity of a ring-shaped ATPase motor and establishes how the innate sequence bias of a termination factor can be modulated to silence pervasive, aberrant transcription.

Fast pressure-jump all-atom simulations and experiments reveal site-specific protein dehydration-folding dynamics.

As theory and experiment have shown, protein dehydration is a major contributor to protein folding. Dehydration upon folding can be characterized directly by all-atom simulations of fast pressure drops, which create desolvated pockets inside the nascent hydrophobic core. Here, we study pressure-drop refolding of three λ-repressor fragment (λ6-85) mutants computationally and experimentally. The three mutants report on tertiary structure formation via different fluorescent helix-helix contact pairs. All-atom simulations of pressure drops capture refolding and unfolding of all three mutants by a similar mechanism, thus validating the nonperturbative nature of the fluorescent contact probes. Analysis of simulated interprobe distances shows that the α-helix 1-3 pair distance displays a slower characteristic time scale than the 1-2 or 3-2 pair distance. To see whether slow packing of α-helices 1 and 3 is reflected in the rate-limiting folding step, fast pressure-drop relaxation experiments captured refolding on a millisecond time scale. These experiments reveal that refolding monitored by 1-3 contact formation indeed is much slower than when monitored by 1-2 or 3-2 contact formation. Unlike the case of the two-state folder [three-α-helix bundle (α3D)], whose drying and core formation proceed in concert, λ6-85 repeatedly dries and rewets different local tertiary contacts before finally forming a solvent-excluded core, explaining the non-two-state behavior observed during refolding in molecular dynamics simulations. This work demonstrates that proteins can explore desolvated pockets and dry globular states numerous times before reaching the native conformation.

Theoretical Description of the Primary Proton-Coupled Electron Transfer Reaction in the Cytochrome bc1 Complex.

The cytochrome bc1 complex is a transmembrane enzymatic protein complex that plays a central role in cellular energy production and is present in both photosynthetic and respiratory chain organelles. Its reaction mechanism is initiated by the binding of a quinol molecule to an active site, followed by a series of charge transfer reactions between the quinol and protein subunits. Previous work hypothesized that the primary reaction was a concerted proton-coupled electron transfer (PCET) reaction because of the apparent absence of intermediate states associated with single proton or electron transfer reactions. In the present study, the kinetics of the primary bc1 complex PCET reaction is investigated with a vibronically nonadiabatic PCET theory in conjunction with all-atom molecular dynamics simulations and electronic structure calculations. The computed rate constants and relatively high kinetic isotope effects are consistent with experimental measurements on related biomimetic systems. The analysis implicates a concerted PCET mechanism with significant hydrogen tunneling and nonadiabatic effects in the bc1 complex. Moreover, the employed theoretical framework is shown to serve as a general strategy for describing PCET reactions in bioenergetic systems.

NAMD goes quantum: an integrative suite for hybrid simulations.

Hybrid methods that combine quantum mechanics (QM) and molecular mechanics (MM) can be applied to studies of reaction mechanisms in locations ranging from active sites of small enzymes to multiple sites in large bioenergetic complexes. By combining the widely used molecular dynamics and visualization programs NAMD and VMD with the quantum chemistry packages ORCA and MOPAC, we created an integrated, comprehensive, customizable, and easy-to-use suite (http://www.ks.uiuc.edu/Research/qmmm). Through the QwikMD interface, setup, execution, visualization, and analysis are streamlined for all levels of expertise.

Lateral Segregation of Photosystem I in Cyanobacterial Thylakoids.

Photosystem I (PSI) is the dominant photosystem in cyanobacteria and it plays a pivotal role in cyanobacterial metabolism. Despite its biological importance, the native organization of PSI in cyanobacterial thylakoid membranes is poorly understood. Here, we use atomic force microscopy (AFM) to show that ordered, extensive macromolecular arrays of PSI complexes are present in thylakoids from Thermosynechococcus elongatus, Synechococcus sp PCC 7002, and Synechocystis sp PCC 6803. Hyperspectral confocal fluorescence microscopy and three-dimensional structured illumination microscopy of Synechocystis sp PCC 6803 cells visualize PSI domains within the context of the complete thylakoid system. Crystallographic and AFM data were used to build a structural model of a membrane landscape comprising 96 PSI trimers and 27,648 chlorophyll a molecules. Rather than facilitating intertrimer energy transfer, the close associations between PSI primarily maximize packing efficiency; short-range interactions with Complex I and cytochrome b6f are excluded from these regions of the membrane, so PSI turnover is sustained by long-distance diffusion of the electron donors at the membrane surface. Elsewhere, PSI-photosystem II contact zones provide sites for docking phycobilisomes and the formation of megacomplexes. PSI-enriched domains in cyanobacteria might foreshadow the partitioning of PSI into stromal lamellae in plants, similarly sustained by long-distance diffusion of electron carriers.

Scalable molecular dynamics on CPU and GPU architectures with NAMD.

NAMDis a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

Structural insights into the functional cycle of the ATPase module of the 26S proteasome.

In eukaryotic cells, the ubiquitin-proteasome system (UPS) is responsible for the regulated degradation of intracellular proteins. The 26S holocomplex comprises the core particle (CP), where proteolysis takes place, and one or two regulatory particles (RPs). The base of the RP is formed by a heterohexameric AAA+ ATPase module, which unfolds and translocates substrates into the CP. Applying single-particle cryo-electron microscopy (cryo-EM) and image classification to samples in the presence of different nucleotides and nucleotide analogs, we were able to observe four distinct conformational states (s1 to s4). The resolution of the four conformers allowed for the construction of atomic models of the AAA+ ATPase module as it progresses through the functional cycle. In a hitherto unobserved state (s4), the gate controlling access to the CP is open. The structures described in this study allow us to put forward a model for the 26S functional cycle driven by ATP hydrolysis.

Angular measurements of the dynein ring reveal a stepping mechanism dependent on a flexible stalk.

The force-generating mechanism of dynein differs from the force-generating mechanisms of other cytoskeletal motors. To examine the structural dynamics of dynein's stepping mechanism in real time, we used polarized total internal reflection fluorescence microscopy with nanometer accuracy localization to track the orientation and position of single motors. By measuring the polarized emission of individual quantum nanorods coupled to the dynein ring, we determined the angular position of the ring and found that it rotates relative to the microtubule (MT) while walking. Surprisingly, the observed rotations were small, averaging only 8.3°, and were only weakly correlated with steps. Measurements at two independent labeling positions on opposite sides of the ring showed similar small rotations. Our results are inconsistent with a classic power-stroke mechanism, and instead support a flexible stalk model in which interhead strain rotates the rings through bending and hinging of the stalk. Mechanical compliances of the stalk and hinge determined based on a 3.3-µs molecular dynamics simulation account for the degree of ring rotation observed experimentally. Together, these observations demonstrate that the stepping mechanism of dynein is fundamentally different from the stepping mechanisms of other well-studied MT motors, because it is characterized by constant small-scale fluctuations of a large but flexible structure fully consistent with the variable stepping pattern observed as dynein moves along the MT.

Mature HIV-1 capsid structure by cryo-electron microscopy and all-atom molecular dynamics.

Retroviral capsid proteins are conserved structurally but assemble into different morphologies. The mature human immunodeficiency virus-1 (HIV-1) capsid is best described by a 'fullerene cone' model, in which hexamers of the capsid protein are linked to form a hexagonal surface lattice that is closed by incorporating 12 capsid-protein pentamers. HIV-1 capsid protein contains an amino-terminal domain (NTD) comprising seven α-helices and a ß-hairpin, a carboxy-terminal domain (CTD) comprising four α-helices, and a flexible linker with a 310-helix connecting the two structural domains. Structures of the capsid-protein assembly units have been determined by X-ray crystallography; however, structural information regarding the assembled capsid and the contacts between the assembly units is incomplete. Here we report the cryo-electron microscopy structure of a tubular HIV-1 capsid-protein assembly at 8 Å resolution and the three-dimensional structure of a native HIV-1 core by cryo-electron tomography. The structure of the tubular assembly shows, at the three-fold interface, a three-helix bundle with critical hydrophobic interactions. Mutagenesis studies confirm that hydrophobic residues in the centre of the three-helix bundle are crucial for capsid assembly and stability, and for viral infectivity. The cryo-electron-microscopy structures enable modelling by large-scale molecular dynamics simulation, resulting in all-atom models for the hexamer-of-hexamer and pentamer-of-hexamer elements as well as for the entire capsid. Incorporation of pentamers results in closer trimer contacts and induces acute surface curvature. The complete atomic HIV-1 capsid model provides a platform for further studies of capsid function and for targeted pharmacological intervention.

Conserved methionine dictates substrate preference in Nramp-family divalent metal transporters.

Natural resistance-associated macrophage protein (Nramp) family transporters catalyze uptake of essential divalent transition metals like iron and manganese. To discriminate against abundant competitors, the Nramp metal-binding site should favor softer transition metals, which interact either covalently or ionically with coordinating molecules, over hard calcium and magnesium, which interact mainly ionically. The metal-binding site contains an unusual, but conserved, methionine, and its sulfur coordinates transition metal substrates, suggesting a vital role in their transport. Using a bacterial Nramp model system, we show that, surprisingly, this conserved methionine is dispensable for transport of the physiological manganese substrate and similar divalents iron and cobalt, with several small amino acid replacements still enabling robust uptake. Moreover, the methionine sulfur's presence makes the toxic metal cadmium a preferred substrate. However, a methionine-to-alanine substitution enables transport of calcium and magnesium. Thus, the putative evolutionary pressure to maintain the Nramp metal-binding methionine likely exists because it-more effectively than any other amino acid-increases selectivity for low-abundance transition metal transport in the presence of high-abundance divalents like calcium and magnesium.

Structure of the human 26S proteasome at a resolution of 3.9 Å.

Protein degradation in eukaryotic cells is performed by the Ubiquitin-Proteasome System (UPS). The 26S proteasome holocomplex consists of a core particle (CP) that proteolytically degrades polyubiquitylated proteins, and a regulatory particle (RP) containing the AAA-ATPase module. This module controls access to the proteolytic chamber inside the CP and is surrounded by non-ATPase subunits (Rpns) that recognize substrates and deubiquitylate them before unfolding and degradation. The architecture of the 26S holocomplex is highly conserved between yeast and humans. The structure of the human 26S holocomplex described here reveals previously unidentified features of the AAA-ATPase heterohexamer. One subunit, Rpt6, has ADP bound, whereas the other five have ATP in their binding pockets. Rpt6 is structurally distinct from the other five Rpt subunits, most notably in its pore loop region. For Rpns, the map reveals two main, previously undetected, features: the C terminus of Rpn3 protrudes into the mouth of the ATPase ring; and Rpn1 and Rpn2, the largest proteasome subunits, are linked by an extended connection. The structural features of the 26S proteasome observed in this study are likely to be important for coordinating the proteasomal subunits during substrate processing.

PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations.

Molecular dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth has led to the production of a large amount of simulation data that need to be filtered for relevant information to address specific biomedical and biochemical questions. One of the most relevant molecular properties that can be investigated by all-atom MD simulations is the time-dependent evolution of the complex noncovalent interaction networks governing such fundamental aspects as molecular recognition, binding strength, and mechanical and structural stability. Extracting, evaluating, and visualizing noncovalent interactions is a key task in the daily work of structural biologists. We have developed PyContact, an easy-to-use, highly flexible, and intuitive graphical user interface-based application, designed to provide a toolkit to investigate biomolecular interactions in MD trajectories. PyContact is designed to facilitate this task by enabling identification of relevant noncovalent interactions in a comprehensible manner. The implementation of PyContact as a standalone application enables rapid analysis and data visualization without any additional programming requirements, and also preserves full in-program customization and extension capabilities for advanced users. The statistical analysis representation is interactively combined with full mapping of the results on the molecular system through the synergistic connection between PyContact and VMD. We showcase the capabilities and scientific significance of PyContact by analyzing and visualizing in great detail the noncovalent interactions underlying the ion permeation pathway of the human P2X3 receptor. As a second application, we examine the protein-protein interaction network of the mechanically ultrastable cohesin-dockering complex.

Nascent peptide assists the ribosome in recognizing chemically distinct small molecules.

Regulation of gene expression in response to the changing environment is critical for cell survival. For instance, binding of macrolide antibiotics to the ribosome promotes translation arrest at the leader open reading frames ermCL and ermBL, which is necessary for inducing the antibiotic resistance genes ermC and ermB. Cladinose-containing macrolides such as erythromycin (ERY), but not ketolides such as telithromycin (TEL), arrest translation of ermCL, whereas either ERY or TEL stall ermBL translation. How the ribosome distinguishes between chemically similar small molecules is unknown. We show that single amino acid changes in the leader peptide switch the specificity of recognition of distinct molecules, triggering gene activation in response to ERY alone, to TEL alone or to both antibiotics or preventing stalling altogether. Thus, the ribosomal response to chemical signals can be modulated by minute changes in the nascent peptide, suggesting that protein sequences could have been optimized for rendering translation sensitive to environmental cues.

The ribosome can discriminate the chirality of amino acids within its peptidyl-transferase center.

The cellular translational machinery (TM) synthesizes proteins using exclusively L- or achiral aminoacyl-tRNAs (aa-tRNAs), despite the presence of D-amino acids in nature and their ability to be aminoacylated onto tRNAs by aa-tRNA synthetases. The ubiquity of L-amino acids in proteins has led to the hypothesis that D-amino acids are not substrates for the TM. Supporting this view, protein engineering efforts to incorporate D-amino acids into proteins using the TM have thus far been unsuccessful. Nonetheless, a mechanistic understanding of why D-aa-tRNAs are poor substrates for the TM is lacking. To address this deficiency, we have systematically tested the translation activity of D-aa-tRNAs using a series of biochemical assays. We find that the TM can effectively, albeit slowly, accept D-aa-tRNAs into the ribosomal aa-tRNA binding (A) site, use the A-site D-aa-tRNA as a peptidyl-transfer acceptor, and translocate the resulting peptidyl-D-aa-tRNA into the ribosomal peptidyl-tRNA binding (P) site. During the next round of continuous translation, however, we find that ribosomes carrying a P-site peptidyl-D-aa-tRNA partition into subpopulations that are either translationally arrested or that can continue translating. Consistent with its ability to arrest translation, chemical protection experiments and molecular dynamics simulations show that P site-bound peptidyl-D-aa-tRNA can trap the ribosomal peptidyl-transferase center in a conformation in which peptidyl transfer is impaired. Our results reveal a novel mechanism through which D-aa-tRNAs interfere with translation, provide insight into how the TM might be engineered to use D-aa-tRNAs, and increase our understanding of the physiological role of a widely distributed enzyme that clears D-aa-tRNAs from cells.

Dynamic allostery governs cyclophilin A-HIV capsid interplay.

Host factor protein Cyclophilin A (CypA) regulates HIV-1 viral infectivity through direct interactions with the viral capsid, by an unknown mechanism. CypA can either promote or inhibit viral infection, depending on host cell type and HIV-1 capsid (CA) protein sequence. We have examined the role of conformational dynamics on the nanosecond to millisecond timescale in HIV-1 CA assemblies in the escape from CypA dependence, by magic-angle spinning (MAS) NMR and molecular dynamics (MD). Through the analysis of backbone (1)H-(15)N and (1)H-(13)C dipolar tensors and peak intensities from 3D MAS NMR spectra of wild-type and the A92E and G94D CypA escape mutants, we demonstrate that assembled CA is dynamic, particularly in loop regions. The CypA loop in assembled wild-type CA from two strains exhibits unprecedented mobility on the nanosecond to microsecond timescales, and the experimental NMR dipolar order parameters are in quantitative agreement with those calculated from MD trajectories. Remarkably, the CypA loop dynamics of wild-type CA HXB2 assembly is significantly attenuated upon CypA binding, and the dynamics profiles of the A92E and G94D CypA escape mutants closely resemble that of wild-type CA assembly in complex with CypA. These results suggest that CypA loop dynamics is a determining factor in HIV-1's escape from CypA dependence.

Chemomechanical Coupling in Hexameric Protein-Protein Interfaces Harnesses Energy within V-Type ATPases.

ATP synthase is the most prominent bioenergetic macromolecular motor in all life forms, utilizing the proton gradient across the cell membrane to fuel the synthesis of ATP. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, the precise molecular mechanism whereby vacuolar (V-type) ATP synthase fulfills its biological function remains largely fragmentary. Recently, crystallographers provided the first high-resolution view of ATP activity in Enterococcus hirae V1-ATPase. Employing a combination of transition-path sampling and high-performance free-energy methods, the sequence of conformational transitions involved in a functional cycle accompanying ATP hydrolysis has been investigated in unprecedented detail over an aggregate simulation time of 65 µs. Our simulated pathways reveal that the chemical energy produced by ATP hydrolysis is harnessed via the concerted motion of the protein-protein interfaces in the V1-ring, and is nearly entirely consumed in the rotation of the central stalk. Surprisingly, in an ATPase devoid of a central stalk, the interfaces of this ring are perfectly designed for inducing ATP hydrolysis. However, in a complete V1-ATPase, the mechanical property of the central stalk is a key determinant of the rate of ATP turnover. The simulations further unveil a sequence of events, whereby unbinding of the hydrolysis product (ADP + Pi) is followed by ATP uptake, which, in turn, leads to the torque generation step and rotation of the center stalk. Molecular trajectories also bring to light multiple intermediates, two of which have been isolated in independent crystallography experiments.

Contributions of Charged Residues in Structurally Dynamic Capsid Surface Loops to Rous Sarcoma Virus Assembly.

UNLABELLED: Extensive studies of orthoretroviral capsids have shown that many regions of the CA protein play unique roles at different points in the virus life cycle. The N-terminal domain (NTD) flexible-loop (FL) region is one such example: exposed on the outer capsid surface, it has been implicated in Gag-mediated particle assembly, capsid maturation, and early replication events. We have now defined the contributions of charged residues in the FL region of the Rous sarcoma virus (RSV) CA to particle assembly. Effects of mutations on assembly were assessed in vivo and in vitro and analyzed in light of new RSV Gag lattice models. Virus replication was strongly dependent on the preservation of charge at a few critical positions in Gag-Gag interfaces. In particular, a cluster of charges at the beginning of FL contributes to an extensive electrostatic network that is important for robust Gag assembly and subsequent capsid maturation. Second-site suppressor analysis suggests that one of these charged residues, D87, has distal influence on interhexamer interactions involving helix α7. Overall, the tolerance of FL to most mutations is consistent with current models of Gag lattice structures. However, the results support the interpretation that virus evolution has achieved a charge distribution across the capsid surface that (i) permits the packing of NTD domains in the outer layer of the Gag shell, (ii) directs the maturational rearrangements of the NTDs that yield a functional core structure, and (iii) supports capsid function during the early stages of virus infection. IMPORTANCE: The production of infectious retrovirus particles is a complex process, a choreography of protein and nucleic acid that occurs in two distinct stages: formation and release from the cell of an immature particle followed by an extracellular maturation phase during which the virion proteins and nucleic acids undergo major rearrangements that activate the infectious potential of the virion. This study examines the contributions of charged amino acids on the surface of the Rous sarcoma virus capsid protein in the assembly of appropriately formed immature particles and the maturational transitions that create a functional virion. The results provide important biological evidence in support of recent structural models of the RSV immature virions and further suggest that immature particle assembly and virion maturation are controlled by an extensive network of electrostatic interactions and long-range communication across the capsid surface.

Computational de novo design of antibodies binding to a peptide with high affinity.

Antibody drugs play a critical role in infectious diseases, cancer, autoimmune diseases, and inflammation. However, experimental methods for the generation of therapeutic antibodies such as using immunized mice or directed evolution remain time consuming and cannot target a specific antigen epitope. Here, we describe the application of a computational framework called OptMAVEn combined with molecular dynamics to de novo design antibodies. Our reference system is antibody 2D10, a single-chain antibody (scFv) that recognizes the dodecapeptide DVFYPYPYASGS, a peptide mimic of mannose-containing carbohydrates. Five de novo designed scFvs sharing less than 75% sequence similarity to all existing natural antibody sequences were generated using OptMAVEn and their binding to the dodecapeptide was experimentally characterized by biolayer interferometry and isothermal titration calorimetry. Among them, three scFvs show binding affinity to the dodecapeptide at the nM level. Critically, these de novo designed scFvs exhibit considerably diverse modeled binding modes with the dodecapeptide. The results demonstrate the potential of OptMAVEn for the de novo design of thermally and conformationally stable antibodies with high binding affinity to antigens and encourage the targeting of other antigen targets in the future. Biotechnol. Bioeng. 2017;114: 1331-1342. © 2017 Wiley Periodicals, Inc.

Advances in the molecular dynamics flexible fitting method for cryo-EM modeling.

Molecular Dynamics Flexible Fitting (MDFF) is an established technique for fitting all-atom structures of molecules into corresponding cryo-electron microscopy (cryo-EM) densities. The practical application of MDFF is simple but requires a user to be aware of and take measures against a variety of possible challenges presented by each individual case. Some of these challenges arise from the complexity of a molecular structure or the limited quality of available structural models and densities to be interpreted, while others stem from the intricacies of MDFF itself. The current article serves as an overview of the strategies that have been developed since MDFF's inception to overcome common challenges and successfully perform MDFF simulations.