Viral DNA polymerase structures reveal mechanisms of antiviral drug resistance.

DNA polymerases are important drug targets, and many structural studies have captured them in distinct conformations. However, a detailed understanding of the impact of polymerase conformational dynamics on drug resistance is lacking. We determined cryoelectron microscopy (cryo-EM) structures of DNA-bound herpes simplex virus polymerase holoenzyme in multiple conformations and interacting with antivirals in clinical use. These structures reveal how the catalytic subunit Pol and the processivity factor UL42 bind DNA to promote processive DNA synthesis. Unexpectedly, in the absence of an incoming nucleotide, we observed Pol in multiple conformations with the closed state sampled by the fingers domain. Drug-bound structures reveal how antivirals may selectively bind enzymes that more readily adopt the closed conformation. Molecular dynamics simulations and the cryo-EM structure of a drug-resistant mutant indicate that some resistance mutations modulate conformational dynamics rather than directly impacting drug binding, thus clarifying mechanisms that drive drug selectivity.

High-resolution structures with bound Mn2+ and Cd2+ map the metal import pathway in an Nramp transporter.

Transporters of the Nramp (Natural resistance-associated macrophage protein) family import divalent transition metal ions into cells of most organisms. By supporting metal homeostasis, Nramps prevent diseases and disorders related to metal insufficiency or overload. Previous studies revealed that Nramps take on a LeuT fold and identified the metal-binding site. We present high-resolution structures of Deinococcus radiodurans (Dra)Nramp in three stable conformations of the transport cycle revealing that global conformational changes are supported by distinct coordination geometries of its physiological substrate, Mn2+, across conformations, and by conserved networks of polar residues lining the inner and outer gates. In addition, a high-resolution Cd2+-bound structure highlights differences in how Cd2+ and Mn2+ are coordinated by DraNramp. Complementary metal binding studies using isothermal titration calorimetry with a series of mutated DraNramp proteins indicate that the thermodynamic landscape for binding and transporting physiological metals like Mn2+ is different and more robust to perturbation than for transporting the toxic Cd2+ metal. Overall, the affinity measurements and high-resolution structural information on metal substrate binding provide a foundation for understanding the substrate selectivity of essential metal ion transporters like Nramps.

Target-driven machine learning-enabled virtual screening (TAME-VS) platform for early-stage hit identification.

High-throughput screening (HTS) methods enable the empirical evaluation of a large scale of compounds and can be augmented by virtual screening (VS) techniques to save time and money by using potential active compounds for experimental testing. Structure-based and ligand-based virtual screening approaches have been extensively studied and applied in drug discovery practice with proven outcomes in advancing candidate molecules. However, the experimental data required for VS are expensive, and hit identification in an effective and efficient manner is particularly challenging during early-stage drug discovery for novel protein targets. Herein, we present our TArget-driven Machine learning-Enabled VS (TAME-VS) platform, which leverages existing chemical databases of bioactive molecules to modularly facilitate hit finding. Our methodology enables bespoke hit identification campaigns through a user-defined protein target. The input target ID is used to perform a homology-based target expansion, followed by compound retrieval from a large compilation of molecules with experimentally validated activity. Compounds are subsequently vectorized and adopted for machine learning (ML) model training. These machine learning models are deployed to perform model-based inferential virtual screening, and compounds are nominated based on predicted activity. Our platform was retrospectively validated across ten diverse protein targets and demonstrated clear predictive power. The implemented methodology provides a flexible and efficient approach that is accessible to a wide range of users. The TAME-VS platform is publicly available at https://github.com/bymgood/Target-driven-ML-enabled-VS to facilitate early-stage hit identification.

Velcrin-induced selective cleavage of tRNALeu(TAA) by SLFN12 causes cancer cell death.

Velcrin compounds kill cancer cells expressing high levels of phosphodiesterase 3A (PDE3A) and Schlafen family member 12 (SLFN12) by inducing complex formation between these two proteins, but the mechanism of cancer cell killing by the PDE3A-SLFN12 complex is not fully understood. Here, we report that the physiological substrate of SLFN12 RNase is tRNALeu(TAA). SLFN12 selectively digests tRNALeu(TAA), and velcrin treatment promotes the cleavage of tRNALeu(TAA) by inducing PDE3A-SLFN12 complex formation in vitro. We found that distinct sequences in the variable loop and acceptor stem of tRNALeu(TAA) are required for substrate digestion. Velcrin treatment of sensitive cells results in downregulation of tRNALeu(TAA), ribosome pausing at Leu-TTA codons and global inhibition of protein synthesis. Velcrin-induced cleavage of tRNALeu(TAA) by SLFN12 and the concomitant global inhibition of protein synthesis thus define a new mechanism of apoptosis initiation.

Correction to "An Iterative Approach Guides Discovery of the FabI Inhibitor Fabimycin, a Late-Stage Antibiotic Candidate with In Vivo Efficacy against Drug-Resistant Gram-Negative Infections".

[This corrects the article DOI: 10.1021/acscentsci.2c00598.].

Revealing a Hidden Intermediate of Rotatory Catalysis with X-ray Crystallography and Molecular Simulations.

The mechanism of rotatory catalysis in ATP-hydrolyzing molecular motors remains an unresolved puzzle in biological energy transfer. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, knowledge on how the coupling between the chemical and mechanical steps within motors enforces directional rotatory movements remains fragmentary. Even more contentious is to pinpoint the rate-limiting step of a multistep rotation process. Here, using vacuolar or V1-type hexameric ATPase as an exemplary rotational motor, we present a model of the complete 4-step conformational cycle involved in rotatory catalysis. First, using X-ray crystallography, a new intermediate or "dwell" is identified, which enables the release of an inorganic phosphate (or Pi) after ATP hydrolysis. Using molecular dynamics simulations, this new dwell is placed in a sequence with three other crystal structures to derive a putative cyclic rotation path. Free-energy simulations are employed to estimate the rate of the hexameric protein transformations and delineate allosteric effects that allow new reactant ATP entry only after hydrolysis product exit. An analysis of transfer entropy brings to light how the side-chain-level interactions transcend into larger-scale reorganizations, highlighting the role of the ubiquitous arginine-finger residues in coupling chemical and mechanical information. An inspection of all known rates encompassing the 4-step rotation mechanism implicates the overcoming of the ADP interactions with V1-ATPase to be the rate-limiting step of motor action.

An Iterative Approach Guides Discovery of the FabI Inhibitor Fabimycin, a Late-Stage Antibiotic Candidate with In Vivo Efficacy against Drug-Resistant Gram-Negative Infections.

Genomic studies and experiments with permeability-deficient strains have revealed a variety of biological targets that can be engaged to kill Gram-negative bacteria. However, the formidable outer membrane and promiscuous efflux pumps of these pathogens prevent many candidate antibiotics from reaching these targets. One such promising target is the enzyme FabI, which catalyzes the rate-determining step in bacterial fatty acid biosynthesis. Notably, FabI inhibitors have advanced to clinical trials for Staphylococcus aureus infections but not for infections caused by Gram-negative bacteria. Here, we synthesize a suite of FabI inhibitors whose structures fit permeation rules for Gram-negative bacteria and leverage activity against a challenging panel of Gram-negative clinical isolates as a filter for advancement. The compound to emerge, called fabimycin, has impressive activity against >200 clinical isolates of Escherichia coli, Klebsiella pneumoniae, and Acinetobacter baumannii, and does not kill commensal bacteria. X-ray structures of fabimycin in complex with FabI provide molecular insights into the inhibition. Fabimycin demonstrates activity in multiple mouse models of infection caused by Gram-negative bacteria, including a challenging urinary tract infection model. Fabimycin has translational promise, and its discovery provides additional evidence that antibiotics can be systematically modified to accumulate in Gram-negative bacteria and kill these problematic pathogens.

Protein Flexibility and Dissociation Pathway Differentiation Can Explain Onset of Resistance Mutations in Kinases.

Understanding how mutations render a drug ineffective is a problem of immense relevance. Often the mechanism through which mutations cause drug resistance can be explained purely through thermodynamics. However, the more perplexing situation is when two proteins have the same drug binding affinities but different residence times. In this work, we demonstrate how all-atom molecular dynamics simulations using recent developments grounded in statistical mechanics can provide a detailed mechanistic rationale for such variances. We discover dissociation mechanisms for the anti-cancer drug Imatinib (Gleevec) against wild-type and the N368S mutant of Abl kinase. We show how this point mutation triggers far-reaching changes in the protein's flexibility and leads to a different, much faster, drug dissociation pathway. We believe that this work marks an efficient and scalable approach to obtain mechanistic insight into resistance mutations in biomolecular receptors that are hard to explain using a structural perspective.

Poor Person's pH Simulation of Membrane Proteins.

pH conditions are central to the functioning of all biomolecules. However, implications of pH changes are nontrivial on a molecular scale. Though a rigorous microscopic definition of pH exists, its implementation in classical molecular dynamics (MD) simulations is cumbersome, and more so in large integral membrane systems. In this chapter, an integrative pipeline is described that combines Multi-Conformation Continuum Electrostatics (MCCE) computations with MD simulations to capture the effect of transient protonation states on the coupled conformational changes in transmembrane proteins. The core methodologies are explained, and all the software required to set up this pipeline are outlined with their key parameters. All associated analyses of structure and function are provided using two case studies, namely those of bioenergetic complexes: NADH dehydrogenase (complex I) and Vo domain of V-type ATPase. The hybrid MCCE-MD pipeline has allowed the discovery of hydrogen bond networks, ligand binding pathways, and disease-causing mutations.

Cryo-EM model validation recommendations based on outcomes of the 2019 EMDataResource challenge.

This paper describes outcomes of the 2019 Cryo-EM Model Challenge. The goals were to (1) assess the quality of models that can be produced from cryogenic electron microscopy (cryo-EM) maps using current modeling software, (2) evaluate reproducibility of modeling results from different software developers and users and (3) compare performance of current metrics used for model evaluation, particularly Fit-to-Map metrics, with focus on near-atomic resolution. Our findings demonstrate the relatively high accuracy and reproducibility of cryo-EM models derived by 13 participating teams from four benchmark maps, including three forming a resolution series (1.8 to 3.1 Å). The results permit specific recommendations to be made about validating near-atomic cryo-EM structures both in the context of individual experiments and structure data archives such as the Protein Data Bank. We recommend the adoption of multiple scoring parameters to provide full and objective annotation and assessment of the model, reflective of the observed cryo-EM map density.

CryoFold: determining protein structures and data-guided ensembles from cryo-EM density maps.

Cryo-electron microscopy (EM) requires molecular modeling to refine structural details from data. Ensemble models arrive at low free-energy molecular structures, but are computationally expensive and limited to resolving only small proteins that cannot be resolved by cryo-EM. Here, we introduce CryoFold - a pipeline of molecular dynamics simulations that determines ensembles of protein structures directly from sequence by integrating density data of varying sparsity at 3-5 Å resolution with coarse-grained topological knowledge of the protein folds. We present six examples showing its broad applicability for folding proteins between 72 to 2000 residues, including large membrane and multi-domain systems, and results from two EMDB competitions. Driven by data from a single state, CryoFold discovers ensembles of common low-energy models together with rare low-probability structures that capture the equilibrium distribution of proteins constrained by the density maps. Many of these conformations, unseen by traditional methods, are experimentally validated and functionally relevant. We arrive at a set of best practices for data-guided protein folding that are controlled using a Python GUI.

Hydrogen-deuterium exchange mass spectrometry captures distinct dynamics upon substrate and inhibitor binding to a transporter.

Proton-coupled transporters use transmembrane proton gradients to power active transport of nutrients inside the cell. High-resolution structures often fail to capture the coupling between proton and ligand binding, and conformational changes associated with transport. We combine HDX-MS with mutagenesis and MD simulations to dissect the molecular mechanism of the prototypical transporter XylE. We show that protonation of a conserved aspartate triggers conformational transition from outward-facing to inward-facing state. This transition only occurs in the presence of substrate xylose, while the inhibitor glucose locks the transporter in the outward-facing state. MD simulations corroborate the experiments by showing that only the combination of protonation and xylose binding, and not glucose, sets up the transporter for conformational switch. Overall, we demonstrate the unique ability of HDX-MS to distinguish between the conformational dynamics of inhibitor and substrate binding, and show that a specific allosteric coupling between substrate binding and protonation is a key step to initiate transport.

Cryo-EM and MD infer water-mediated proton transport and autoinhibition mechanisms of Vo complex.

Rotary vacuolar adenosine triphosphatases (V-ATPases) drive transmembrane proton transport through a Vo proton channel subcomplex. Despite recent high-resolution structures of several rotary ATPases, the dynamic mechanism of proton pumping remains elusive. Here, we determined a 2.7-Å cryo-electron microscopy (cryo-EM) structure of yeast Vo proton channel in nanodisc that reveals the location of ordered water molecules along the proton path, details of specific protein-lipid interactions, and the architecture of the membrane scaffold protein. Moreover, we uncover a state of Vo that shows the c-ring rotated by ~14°. Molecular dynamics simulations demonstrate that the two rotary states are in thermal equilibrium and depict how the protonation state of essential glutamic acid residues couples water-mediated proton transfer with c-ring rotation. Our cryo-EM models and simulations also rationalize a mechanism for inhibition of passive proton transport as observed for free Vo that is generated as a result of V-ATPase regulation by reversible disassembly in vivo.

Retrieving functional pathways of biomolecules from single-particle snapshots.

A primary reason for the intense interest in structural biology is the fact that knowledge of structure can elucidate macromolecular functions in living organisms. Sustained effort has resulted in an impressive arsenal of tools for determining the static structures. But under physiological conditions, macromolecules undergo continuous conformational changes, a subset of which are functionally important. Techniques for capturing the continuous conformational changes underlying function are essential for further progress. Here, we present chemically-detailed conformational movies of biological function, extracted data-analytically from experimental single-particle cryo-electron microscopy (cryo-EM) snapshots of ryanodine receptor type 1 (RyR1), a calcium-activated calcium channel engaged in the binding of ligands. The functional motions differ substantially from those inferred from static structures in the nature of conformationally active structural domains, the sequence and extent of conformational motions, and the way allosteric signals are transduced within and between domains. Our approach highlights the importance of combining experiment, advanced data analysis, and molecular simulations.

Molecular Dynamics Flexible Fitting: All You Want to Know About Resolution Exchange.

In recent years, owing to the advances in instrumentation, cryo-EM has emerged as the go-to tool for obtaining high-resolution structures of biomolecular systems. However, building three-dimensional atomic structures of biomolecules from these high-resolution maps remains a concern for the traditional map-guided structure-determination schemes. Recently, we developed a computational tool, Resolution Exchange Molecular Dynamics Flexible Fitting (ReMDFF) to address this problem by re-refining a search model against a series of maps of progressively higher resolutions, which ends with the original experimental resolution (Wang et al., J Struct Biol 204(2):319-328, 2018). In this chapter, we present a step-by-step outline for preparing, executing, and analyzing ReMDFF refinements of simple proteins and multimeric complexes. The structure determination of carbon monoxide dehydrogenase and Mg2+-channel CorA are employed as case studies. All scripts are provided via GitHub (Vant, Resolution exchange molecular dynamics flexible fitting (ReMDFF) all you want to know about flexible fitting, 2019, https://github.com/jvant/ReMDFF_Singharoy_Group.git ).

Charge Transfer and Chemo-Mechanical Coupling in Respiratory Complex I.

The mitochondrial respiratory chain, formed by five protein complexes, utilizes energy from catabolic processes to synthesize ATP. Complex I, the first and the largest protein complex of the chain, harvests electrons from NADH to reduce quinone, while pumping protons across the mitochondrial membrane. Detailed knowledge of the working principle of such coupled charge-transfer processes remains, however, fragmentary due to bottlenecks in understanding redox-driven conformational transitions and their interplay with the hydrated proton pathways. Complex I from Thermus thermophilus encases 16 subunits with nine iron-sulfur clusters, reduced by electrons from NADH. Here, employing the latest crystal structure of T. thermophilus complex I, we have used microsecond-scale molecular dynamics simulations to study the chemo-mechanical coupling between redox changes of the iron-sulfur clusters and conformational transitions across complex I. First, we identify the redox switches within complex I, which allosterically couple the dynamics of the quinone binding pocket to the site of NADH reduction. Second, our free-energy calculations reveal that the affinity of the quinone, specifically menaquinone, for the binding-site is higher than that of its reduced, menaquinol form-a design essential for menaquinol release. Remarkably, the barriers to diffusive menaquinone dynamics are lesser than that of the more ubiquitous ubiquinone, and the naphthoquinone headgroup of the former furnishes stronger binding interactions with the pocket, favoring menaquinone for charge transport in T. thermophilus. Our computations are consistent with experimentally validated mutations and hierarchize the key residues into three functional classes, identifying new mutation targets. Third, long-range hydrogen-bond networks connecting the quinone-binding site to the transmembrane subunits are found to be responsible for proton pumping. Put together, the simulations reveal the molecular design principles linking redox reactions to quinone turnover to proton translocation in complex I.

Flexible Fitting of Small Molecules into Electron Microscopy Maps Using Molecular Dynamics Simulations with Neural Network Potentials.

Despite significant advances in resolution, the potential for cryo-electron microscopy (EM) to be used in determining the structures of protein-drug complexes remains unrealized. Determination of accurate structures and coordination of bound ligands necessitates simultaneous fitting of the models into the density envelopes, exhaustive sampling of the ligand geometries, and, most importantly, concomitant rearrangements in the side chains to optimize the binding energy changes. In this article, we present a flexible-fitting pipeline where molecular dynamics flexible fitting (MDFF) is used to refine structures of protein-ligand complexes from 3 to 5 Å electron density data. Enhanced sampling is employed to explore the binding pocket rearrangements. To provide a model that can accurately describe the conformational dynamics of the chemically diverse set of small-molecule drugs inside MDFF, we use QM/MM and neural-network potential (NNP)/MM models of protein-ligand complexes, where the ligand is represented using the QM or NNP model, and the protein is represented using established molecular mechanical force fields (e.g., CHARMM). This pipeline offers structures commensurate to or better than recently submitted high-resolution cryo-EM or X-ray models, even when given medium to low-resolution data as input. The use of the NNPs makes the algorithm more robust to the choice of search models, offering a radius of convergence of 6.5 Å for ligand structure determination. The quality of the predicted structures was also judged by density functional theory calculations of ligand strain energy. This strain potential energy is found to systematically decrease with better fitting to density and improved ligand coordination, indicating correct binding interactions. A computationally inexpensive protocol for computing strain energy is reported as part of the model analysis protocol that monitors both the ligand fit as well as model quality.

XFEL and NMR Structures of Francisella Lipoprotein Reveal Conformational Space of Drug Target against Tularemia.

Francisella tularensis is the causative agent for the potentially fatal disease tularemia. The lipoprotein Flpp3 has been identified as a virulence determinant of tularemia with no sequence homology outside the Francisella genus. We report a room temperature structure of Flpp3 determined by serial femtosecond crystallography that exists in a significantly different conformation than previously described by the NMR-determined structure. Furthermore, we investigated the conformational space and energy barriers between these two structures by molecular dynamics umbrella sampling and identified three low-energy intermediate states, transitions between which readily occur at room temperature. We have also begun to investigate organic compounds in silico that may act as inhibitors to Flpp3. This work paves the road to developing targeted therapeutics against tularemia and aides in our understanding of the disease mechanisms of tularemia.

Computational Dissection of Membrane Transport at a Microscopic Level.

Membrane transporters are key gatekeeper proteins at cellular membranes that closely control the traffic of materials. Their function relies on structural rearrangements of varying degrees that facilitate substrate translocation across the membrane. Characterizing these functionally important molecular events at a microscopic level is key to our understanding of membrane transport, yet challenging to achieve experimentally. Recent advances in simulation technology and computing power have rendered molecular dynamics (MD) simulation a powerful biophysical tool to investigate a wide range of dynamical events spanning multiple spatial and temporal scales. Here, we review recent studies of diverse membrane transporters using computational methods, with an emphasis on highlighting the technical challenges, key lessons learned, and new opportunities to illuminate transporter structure and function.

Integrating hydrogen-deuterium exchange mass spectrometry with molecular dynamics simulations to probe lipid-modulated conformational changes in membrane proteins.

Biological membranes define the boundaries of cells and are composed primarily of phospholipids and membrane proteins. It has become increasingly evident that direct interactions of membrane proteins with their surrounding lipids play key roles in regulating both protein conformations and function. However, the exact nature and structural consequences of these interactions remain difficult to track at the molecular level. Here, we present a protocol that specifically addresses this challenge. First, hydrogen-deuterium exchange mass spectrometry (HDX-MS) of membrane proteins incorporated into nanodiscs of controlled lipid composition is used to obtain information on the lipid species that are involved in modulating the conformational changes in the membrane protein. Then molecular dynamics (MD) simulations in lipid bilayers are used to pinpoint likely lipid-protein interactions, which can be tested experimentally using HDX-MS. By bringing together the MD predictions with the conformational readouts from HDX-MS, we have uncovered key lipid-protein interactions implicated in stabilizing important functional conformations. This protocol can be applied to virtually any integral membrane protein amenable to classic biophysical studies and for which a near-atomic-resolution structure or homology model is available. This protocol takes ~4 d to complete, excluding the time for data analysis and MD simulations, which depends on the size of the protein under investigation.