Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Molecular mechanism and energetics of coupling between substrate binding and product release in the F1-ATPase catalytic cycle.

F1-ATPase is a motor protein that couples the rotation of its rotary [Formula: see text] subunit with ATP synthesis or hydrolysis. Single-molecule experiments indicate that nucleotide binding and release events occur almost simultaneously during the synthesis cycle, allowing the energy gain due to spontaneous binding of ADP to one catalytic [Formula: see text] subunit to be directly harnessed for driving the release of ATP from another rather than being dissipated as heat. Here, we examine the unknown mechanism of this coupling that is critical for an exceptionally high mechanochemical efficiency of F1-ATPase by means of all-atom free-energy simulations. We find that nondissipative and kinetically fast progression of the motor in the synthesis direction requires a concerted conformational change involving the closure of the ADP-binding [Formula: see text] subunit followed by the gradual opening of the ATP-releasing [Formula: see text] subunit over the course of the 30 to 40° rotary substep of the [Formula: see text] subunit. This rotary substep, preceding the ATP-dependent metastable state, allows for the recovery of a large portion of the ADP binding energy in the conformation of ATP-bound [Formula: see text] that gradually adopts the low-affinity conformation, captured also by the recent cryo-EM structure of this elusive state. The release of ATP from this nearly open conformation leads to its further opening, which enables the progression of the motor to the next catalytic metastable state. Our simulations explain this energy conversion mechanism in terms of intersubunit and ligand-protein interactions.

Structural Basis for Polyproline-Mediated Ribosome Stalling and Rescue by the Translation Elongation Factor EF-P.

Ribosomes synthesizing proteins containing consecutive proline residues become stalled and require rescue via the action of uniquely modified translation elongation factors, EF-P in bacteria, or archaeal/eukaryotic a/eIF5A. To date, no structures exist of EF-P or eIF5A in complex with translating ribosomes stalled at polyproline stretches, and thus structural insight into how EF-P/eIF5A rescue these arrested ribosomes has been lacking. Here we present cryo-EM structures of ribosomes stalled on proline stretches, without and with modified EF-P. The structures suggest that the favored conformation of the polyproline-containing nascent chain is incompatible with the peptide exit tunnel of the ribosome and leads to destabilization of the peptidyl-tRNA. Binding of EF-P stabilizes the P-site tRNA, particularly via interactions between its modification and the CCA end, thereby enforcing an alternative conformation of the polyproline-containing nascent chain, which allows a favorable substrate geometry for peptide bond formation.

Bending-torsional elasticity and energetics of the plus-end microtubule tip.

SignificanceThe mechanochemical basis of microtubule growth, which is essential for the normal function and division of eukaryotic cells, has remained elusive and controversial, despite extensive work. In particular, recent findings have created the paradox that the microtubule plus-end tips look very similar during both growing and shrinking phases, thereby challenging the traditional textbook picture. Our large-scale atomistic simulations resolve this paradox and explain microtubule growth and shrinkage dynamics as a process governed by energy barriers between protofilament conformations, the heights of which are in turn fine-tuned by different nucleotide states, thus implementing an information-driven Brownian ratchet.

Folding of VemP into translation-arresting secondary structure is driven by the ribosome exit tunnel.

The ribosome is a fundamental biomolecular complex that synthesizes proteins in cells. Nascent proteins emerge from the ribosome through a tunnel, where they may interact with the tunnel walls or small molecules such as antibiotics. These interactions can cause translational arrest with notable physiological consequences. Here, we studied the arrest caused by the regulatory peptide VemP, which is known to form α-helices inside the ribosome tunnel near the peptidyl transferase center under specific conditions. We used all-atom molecular dynamics simulations of the entire ribosome and circular dichroism spectroscopy to study the driving forces of helix formation and how VemP causes the translational arrest. To that aim, we compared VemP dynamics in the ribosome tunnel with its dynamics in solution. We show that the VemP peptide has a low helical propensity in water and that the propensity is higher in mixtures of water and trifluorethanol. We propose that helix formation within the ribosome is driven by the interactions of VemP with the tunnel and that a part of VemP acts as an anchor. This anchor might slow down VemP progression through the tunnel enabling α-helix formation, which causes the elongation arrest.

Conformational Dynamics of the Most Efficient Carboxylase Contributes to Efficient CO2 Fixation.

Crotonyl-CoA carboxylase/reductase (Ccr) is one of the fastest CO2 fixing enzymes and has become part of efficient artificial CO2-fixation pathways in vitro, paving the way for future applications. The underlying mechanism of its efficiency, however, is not yet completely understood. X-ray structures of different intermediates in the catalytic cycle reveal tetramers in a dimer of dimers configuration with two open and two closed active sites. Upon binding a substrate, this active site changes its conformation from the open state to the closed state. It is challenging to predict how these coupled conformational changes will alter the CO2 binding affinity to the reaction's active site. To determine whether the open or closed conformations of Ccr affect binding of CO2 to the active site, we performed all-atom molecular simulations of the various conformations of Ccr. The open conformation without a substrate showed the highest binding affinity. The CO2 binding sites are located near the catalytic relevant Asn81 and His365 residues and in an optimal position for CO2 fixation. Furthermore, they are unaffected by substrate binding, and CO2 molecules stay in these binding sites for a longer time. Longer times at these reactive binding sites facilitate CO2 fixation through the nucleophilic attack of the reactive enolate in the closed conformation. We previously demonstrated that the Asn81Leu variant cannot fix CO2. Simulations of the Asn81Leu variant explain the loss of activity through the removal of the Asn81 and His365 binding sites. Overall, our findings show that the conformational dynamics of the enzyme controls CO2 binding. Conformational changes in Ccr increase the level of CO2 in the open subunit before the substrate is bound, the active site closes, and the reaction starts. The full catalytic Ccr cycle alternates among CO2 addition, conformational change, and chemical reaction in the four subunits of the tetramer coordinated by communication between the two dimers.

Choice of fluorophore affects dynamic DNA nanostructures.

The ability to dynamically remodel DNA origami structures or functional nanodevices is highly desired in the field of DNA nanotechnology. Concomitantly, the use of fluorophores to track and validate the dynamics of such DNA-based architectures is commonplace and often unavoidable. It is therefore crucial to be aware of the side effects of popular fluorophores, which are often exchanged without considering the potential impact on the system. Here, we show that the choice of fluorophore can strongly affect the reconfiguration of DNA nanostructures. To this end, we encapsulate a triple-stranded DNA (tsDNA) into water-in-oil compartments and functionalize their periphery with a single-stranded DNA handle (ssDNA). Thus, the tsDNA can bind and unbind from the periphery by reversible opening of the triplex and subsequent strand displacement. Using a combination of experiments, molecular dynamics (MD) simulations, and reaction-diffusion modelling, we demonstrate for 12 different fluorophore combinations that it is possible to alter or even inhibit the DNA nanostructure formation-without changing the DNA sequence. Besides its immediate importance for the design of pH-responsive switches and fluorophore labelling, our work presents a strategy to precisely tune the energy landscape of dynamic DNA nanodevices.

GROMACS in the Cloud: A Global Supercomputer to Speed Up Alchemical Drug Design.

We assess costs and efficiency of state-of-the-art high-performance cloud computing and compare the results to traditional on-premises compute clusters. Our use case is atomistic simulations carried out with the GROMACS molecular dynamics (MD) toolkit with a particular focus on alchemical protein-ligand binding free energy calculations. We set up a compute cluster in the Amazon Web Services (AWS) cloud that incorporates various different instances with Intel, AMD, and ARM CPUs, some with GPU acceleration. Using representative biomolecular simulation systems, we benchmark how GROMACS performs on individual instances and across multiple instances. Thereby we assess which instances deliver the highest performance and which are the most cost-efficient ones for our use case. We find that, in terms of total costs, including hardware, personnel, room, energy, and cooling, producing MD trajectories in the cloud can be about as cost-efficient as an on-premises cluster given that optimal cloud instances are chosen. Further, we find that high-throughput ligand-screening can be accelerated dramatically by using global cloud resources. For a ligand screening study consisting of 19 872 independent simulations or ∼200 µs of combined simulation trajectory, we made use of diverse hardware available in the cloud at the time of the study. The computations scaled-up to reach peak performance using more than 4 000 instances, 140 000 cores, and 3 000 GPUs simultaneously. Our simulation ensemble finished in about 2 days in the cloud, while weeks would be required to complete the task on a typical on-premises cluster consisting of several hundred nodes.

The pathway to GTPase activation of elongation factor SelB on the ribosome.

In all domains of life, selenocysteine (Sec) is delivered to the ribosome by selenocysteine-specific tRNA (tRNASec) with the help of a specialized translation factor, SelB in bacteria. Sec-tRNASec recodes a UGA stop codon next to a downstream mRNA stem-loop. Here we present the structures of six intermediates on the pathway of UGA recoding in Escherichia coli by single-particle cryo-electron microscopy. The structures explain the specificity of Sec-tRNASec binding by SelB and show large-scale rearrangements of Sec-tRNASec. Upon initial binding of SelB-Sec-tRNASec to the ribosome and codon reading, the 30S subunit adopts an open conformation with Sec-tRNASec covering the sarcin-ricin loop (SRL) on the 50S subunit. Subsequent codon recognition results in a local closure of the decoding site, which moves Sec-tRNASec away from the SRL and triggers a global closure of the 30S subunit shoulder domain. As a consequence, SelB docks on the SRL, activating the GTPase of SelB. These results reveal how codon recognition triggers GTPase activation in translational GTPases.

Heterogeneous and rate-dependent streptavidin-biotin unbinding revealed by high-speed force spectroscopy and atomistic simulations.

Receptor-ligand interactions are essential for biological function and their binding strength is commonly explained in terms of static lock-and-key models based on molecular complementarity. However, detailed information on the full unbinding pathway is often lacking due, in part, to the static nature of atomic structures and ensemble averaging inherent to bulk biophysics approaches. Here we combine molecular dynamics and high-speed force spectroscopy on the streptavidin-biotin complex to determine the binding strength and unbinding pathways over the widest dynamic range. Experiment and simulation show excellent agreement at overlapping velocities and provided evidence of the unbinding mechanisms. During unbinding, biotin crosses multiple energy barriers and visits various intermediate states far from the binding pocket, while streptavidin undergoes transient induced fits, all varying with loading rate. This multistate process slows down the transition to the unbound state and favors rebinding, thus explaining the long lifetime of the complex. We provide an atomistic, dynamic picture of the unbinding process, replacing a simple two-state picture with one that involves many routes to the lock and rate-dependent induced-fit motions for intermediates, which might be relevant for other receptor-ligand bonds.

Spatially resolved free-energy contributions of native fold and molten-globule-like Crambin.

The folding stability of a protein is governed by the free-energy difference between its folded and unfolded states, which results from a delicate balance of much larger but almost compensating enthalpic and entropic contributions. The balance can therefore easily be shifted by an external disturbance, such as a mutation of a single amino acid or a change of temperature, in which case the protein unfolds. Effects such as cold denaturation, in which a protein unfolds because of cooling, provide evidence that proteins are strongly stabilized by the solvent entropy contribution to the free-energy balance. However, the molecular mechanisms behind this solvent-driven stability, their quantitative contribution in relation to other free-energy contributions, and how the involved solvent thermodynamics is affected by individual amino acids are largely unclear. Therefore, we addressed these questions using atomistic molecular dynamics simulations of the small protein Crambin in its native fold and a molten-globule-like conformation, which here served as a model for the unfolded state. The free-energy difference between these conformations was decomposed into enthalpic and entropic contributions from the protein and spatially resolved solvent contributions using the nonparametric method Per|Mut. From the spatial resolution, we quantified the local effects on the solvent free-energy difference at each amino acid and identified dependencies of the local enthalpy and entropy on the protein curvature. We identified a strong stabilization of the native fold by almost 500 kJ mol-1 due to the solvent entropy, revealing it as an essential contribution to the total free-energy difference of (53 ± 84) kJ mol-1. Remarkably, more than half of the solvent entropy contribution arose from induced water correlations.

Spatiotemporal Resolution of Conformational Changes in Biomolecules by Combining Pulsed Electron-Electron Double Resonance Spectroscopy with Microsecond Freeze-Hyperquenching.

The function of proteins is linked to their conformations that can be resolved with several high-resolution methods. However, only a few methods can provide the temporal order of intermediates and conformational changes, with each having its limitations. Here, we combine pulsed electron-electron double resonance spectroscopy with a microsecond freeze-hyperquenching setup to achieve spatiotemporal resolution in the angstrom range and lower microsecond time scale. We show that the conformational change of the Cα-helix in the cyclic nucleotide-binding domain of the Mesorhizobium loti potassium channel occurs within about 150 µs and can be resolved with angstrom precision. Thus, this approach holds great promise for obtaining 4D landscapes of conformational changes in biomolecules.

Microtubule instability driven by longitudinal and lateral strain propagation.

Tubulin dimers associate longitudinally and laterally to form metastable microtubules (MTs). MT disassembly is preceded by subtle structural changes in tubulin fueled by GTP hydrolysis. These changes render the MT lattice unstable, but it is unclear exactly how they affect lattice energetics and strain. We performed long-time atomistic simulations to interrogate the impacts of GTP hydrolysis on tubulin lattice conformation, lateral inter-dimer interactions, and (non-)local lateral coordination of dimer motions. The simulations suggest that most of the hydrolysis energy is stored in the lattice in the form of longitudinal strain. While not significantly affecting lateral bond stability, the stored elastic energy results in more strongly confined and correlated dynamics of GDP-tubulins, thereby entropically destabilizing the MT lattice.

How proteins open fusion pores: insights from molecular simulations.

Fusion proteins can play a versatile and involved role during all stages of the fusion reaction. Their roles go far beyond forcing the opposing membranes into close proximity to drive stalk formation and fusion. Molecular simulations have played a central role in providing a molecular understanding of how fusion proteins actively overcome the free energy barriers of the fusion reaction up to the expansion of the fusion pore. Unexpectedly, molecular simulations have revealed a preference of the biological fusion reaction to proceed through asymmetric pathways resulting in the formation of, e.g., a stalk-hole complex, rim-pore, or vertex pore. Force-field based molecular simulations are now able to directly resolve the minimum free-energy path in protein-mediated fusion as well as quantifying the free energies of formed reaction intermediates. Ongoing developments in Graphics Processing Units (GPUs), free energy calculations, and coarse-grained force-fields will soon gain additional insights into the diverse roles of fusion proteins.

Sequential Water and Headgroup Merger: Membrane Poration Paths and Energetics from MD Simulations.

Membrane topology changes such as poration, stalk formation, and hemifusion rupture are essential to cellular function, but their molecular details, energetics, and kinetics are still not fully understood. Here, we present a unified energetic and mechanistic picture of metastable pore defects in tensionless lipid membranes. We used an exhaustive committor analysis to test and select optimal reaction coordinates and also to determine the nucleation mechanism. These reaction coordinates were used to calculate free-energy landscapes that capture the full process and end states. The identified barriers agree with the committor analysis. To enable sufficient sampling of the complete transition path for our molecular dynamics simulations, we developed a "gizmo" potential biasing scheme. The simulations suggest that the essential step in the nucleation is the initial merger of lipid headgroups at the nascent pore center. To facilitate this event, an indentation pathway is energetically preferred to a hydrophobic defect. Continuous water columns that span the indentation were determined to be on-path transients that precede the nucleation barrier. This study gives a quantitative description of the nucleation mechanism and energetics of small metastable pores and illustrates a systematic approach to uncover the mechanisms of diverse cellular membrane remodeling processes.

tRNA Dissociation from EF-Tu after GTP Hydrolysis: Primary Steps and Antibiotic Inhibition.

In each round of ribosomal translation, the translational GTPase elongation factor Tu (EF-Tu) delivers a transfer RNA (tRNA) to the ribosome. After successful decoding, EF-Tu hydrolyzes GTP, which triggers a conformational change that ultimately results in the release of the tRNA from EF-Tu. To identify the primary steps of these conformational changes and how they are prevented by the antibiotic kirromycin, we employed all-atom explicit-solvent molecular dynamics simulations of the full ribosome-EF-Tu complex. Our results suggest that after GTP hydrolysis and Pi release, the loss of interactions between the nucleotide and the switch 1 loop of EF-Tu allows domain D1 of EF-Tu to rotate relative to domains D2 and D3 and leads to an increased flexibility of the switch 1 loop. This rotation induces a closing of the D1-D3 interface and an opening of the D1-D2 interface. We propose that the opening of the D1-D2 interface, which binds the CCA tail of the tRNA, weakens the crucial EF-Tu-tRNA interactions, which lowers tRNA binding affinity, representing the first step of tRNA release. Kirromycin binds within the D1-D3 interface, sterically blocking its closure, but does not prevent hydrolysis. The resulting increased flexibility of switch 1 explains why it is not resolved in kirromycin-bound structures.

CHARMM36m: an improved force field for folded and intrinsically disordered proteins.

The all-atom additive CHARMM36 protein force field is widely used in molecular modeling and simulations. We present its refinement, CHARMM36m (http://mackerell.umaryland.edu/charmm_ff.shtml), with improved accuracy in generating polypeptide backbone conformational ensembles for intrinsically disordered peptides and proteins.

How accurate is circular dichroism-based model validation?

Circular dichroism (CD) spectroscopy is highly sensitive to the secondary structure (SS) composition of proteins. Several methods exist to either estimate the SS composition of a protein or to validate existing structural models using its CD spectrum. The accuracy and precision of these methods depend on the quality of both the measured CD spectrum and the used reference structure. Using a large reference protein set with high-quality CD spectra and synthetic data derived from this set, we quantified deviations from both ideal spectra and reference structures due to experimental limitations. We also determined the impact of these deviations on SS estimation, CD prediction, and SS validation methods of the SESCA analysis package. With regard to the CD spectra, our results suggest intensity scaling errors and non-SS contributions as the main causes of inaccuracies. These factors also can lead to overestimated model errors during validation. The errors of the used reference structures combine non-additively with errors caused by the CD spectrum, which increases the uncertainty of model validation. We have further shown that the effects of scaling errors in the CD spectrum can be nearly eliminated by appropriate re-scaling, and that the accuracy of model validation methods can be improved by accounting for typical non-SS contributions. These improvements have now been implemented within the SESCA package and are available at: https://www.mpibpc.mpg.de/sesca .

More bang for your buck: Improved use of GPU nodes for GROMACS 2018.

We identify hardware that is optimal to produce molecular dynamics (MD) trajectories on Linux compute clusters with the GROMACS 2018 simulation package. Therefore, we benchmark the GROMACS performance on a diverse set of compute nodes and relate it to the costs of the nodes, which may include their lifetime costs for energy and cooling. In agreement with our earlier investigation using GROMACS 4.6 on hardware of 2014, the performance to price ratio of consumer GPU nodes is considerably higher than that of CPU nodes. However, with GROMACS 2018, the optimal CPU to GPU processing power balance has shifted even more toward the GPU. Hence, nodes optimized for GROMACS 2018 and later versions enable a significantly higher performance to price ratio than nodes optimized for older GROMACS versions. Moreover, the shift toward GPU processing allows to cheaply upgrade old nodes with recent GPUs, yielding essentially the same performance as comparable brand-new hardware. © 2019 Wiley Periodicals, Inc.

Engineering fatty acid synthases for directed polyketide production.

In this study, we engineered fatty acid synthases (FAS) for the biosynthesis of short-chain fatty acids and polyketides, guided by a combined in vitro and in silico approach. Along with exploring the synthetic capability of FAS, we aim to build a foundation for efficient protein engineering, with the specific goal of harnessing evolutionarily related megadalton-scale polyketide synthases (PKS) for the tailored production of bioactive natural compounds.