Integrating Hydrogen Deuterium Exchange-Mass Spectrometry with Molecular Simulations Enables Quantification of the Conformational Populations of the Sugar Transporter XylE.

A yet unresolved challenge in structural biology is to quantify the conformational states of proteins underpinning function. This challenge is particularly acute for membrane proteins owing to the difficulties in stabilizing them for in vitro studies. To address this challenge, we present an integrative strategy that combines hydrogen deuterium exchange-mass spectrometry (HDX-MS) with ensemble modeling.

Structural dynamics in the evolution of SARS-CoV-2 spike glycoprotein.

SARS-CoV-2 spike glycoprotein mediates receptor binding and subsequent membrane fusion. It exists in a range of conformations, including a closed state unable to bind the ACE2 receptor, and an open state that does so but displays more exposed antigenic surface. Spikes of variants of concern (VOCs) acquired amino acid changes linked to increased virulence and immune evasion.

Insights into autoregulation of a membrane protein complex by its cytoplasmic domains.

Membrane transporters mediate the passage of molecules across membranes and are essential for cellular function. While the transmembrane region of these proteins is responsible for substrate transport, often the cytoplasmic regions are required for modulating their activity. However, it can be difficult to obtain atomic-resolution descriptions of these autoregulatory domains by classical structural biology techniques, especially if they lack a single, defined structure.

Modeling the native ensemble of PhuS using enhanced sampling MD and HDX-ensemble reweighting.

The cytoplasmic heme binding protein from Pseudomonas aeruginosa, PhuS, plays two essential roles in regulating heme uptake and iron homeostasis. First, PhuS shuttles exogenous heme to heme oxygenase (HemO) for degradation and iron release. Second, PhuS binds DNA and modulates the transcription of the prrF/H small RNAs (sRNAs) involved in the iron-sparing response.

Testing the Limitations of MD-Based Local Electric Fields Using the Vibrational Stark Effect in Solution: Penicillin G as a Test Case.

Noncovalent interactions underlie nearly all molecular processes in the condensed phase from solvation to catalysis. Their quantification within a physically consistent framework remains challenging. Experimental vibrational Stark effect (VSE)-based solvatochromism can be combined with molecular dynamics (MD) simulations to quantify the electrostatic forces in solute-solvent interactions for small rigid molecules and, by extension, when these solutes bind in enzyme active sites.

Chloride-dependent conformational changes in the GlyT1 glycine transporter.

The human GlyT1 glycine transporter requires chloride for its function. However, the mechanism by which Cl exerts its influence is unknown. To examine the role that Cl plays in the transport cycle, we measured the effect of Cl on both glycine binding and conformational changes.

Interpreting Hydrogen-Deuterium Exchange Experiments with Molecular Simulations: Tutorials and Applications of the HDXer Ensemble Reweighting Software [Article v1.0].

Hydrogen-deuterium exchange (HDX) is a comprehensive yet detailed probe of protein structure and dynamics and, coupled to mass spectrometry, has become a powerful tool for investigating an increasingly large array of systems. Computer simulations are often used to help rationalize experimental observations of exchange, but interpretations have frequently been limited to simple, subjective correlations between microscopic dynamical fluctuations and the observed macroscopic exchange behavior. With this in mind, we previously developed the HDX ensemble reweighting approach and associated software, HDXer, to aid the objective interpretation of HDX data using molecular simulations.

Structural predictions of the functions of membrane proteins from HDX-MS.

HDX-MS has emerged as a powerful tool to interrogate the structure and dynamics of proteins and their complexes. Recent advances in the methodology and instrumentation have enabled the application of HDX-MS to membrane proteins. Such targets are challenging to investigate with conventional strategies.

The Role of Electrostatics in Enzymes: Do Biomolecular Force Fields Reflect Protein Electric Fields?

Preorganization of large, directionally oriented, electric fields inside protein active sites has been proposed as a crucial contributor to catalytic mechanism in many enzymes, and it may be efficiently investigated at the atomistic level with molecular dynamics simulations. Here, we evaluate the ability of the AMOEBA polarizable force field, as well as the additive Amber ff14SB and Charmm C36m models, to describe the electric fields present inside the active site of the peptidyl-prolyl isomerase cyclophilin A. We compare the molecular mechanical electric fields to those calculated with a fully first-principles quantum mechanical (QM) representation of the protein, solvent, and ions, and find that AMOEBA consistently shows far greater correlation with the QM electric fields than either of the additive force fields tested.

Interpretation of HDX Data by Maximum-Entropy Reweighting of Simulated Structural Ensembles.

Hydrogen-deuterium exchange combined with mass spectrometry (HDX-MS) is a widely applied biophysical technique that probes the structure and dynamics of biomolecules without the need for site-directed modifications or bio-orthogonal labels. The mechanistic interpretation of HDX data, however, is often qualitative and subjective, owing to a lack of quantitative methods to rigorously translate observed deuteration levels into atomistic structural information. To help address this problem, we have developed a methodology to generate structural ensembles that faithfully reproduce HDX-MS measurements.

Evaluating Anti-CD32b F(ab) Conformation Using Molecular Dynamics and Small-Angle X-Ray Scattering.

Complementary strategies of small-angle x-ray scattering (SAXS) and crystallographic analysis are often used to determine atomistic three-dimensional models of macromolecules and their variability in solution. This combination of techniques is particularly valuable when applied to macromolecular complexes to detect changes within the individual binding partners. Here, we determine the x-ray crystallographic structure of a F(ab) fragment in complex with CD32b, the only inhibitory Fc-γ receptor in humans, and compare the structure of the F(ab) from the crystal complex to SAXS data for the F(ab) alone in solution.

Evaluation of solvation free energies for small molecules with the AMOEBA polarizable force field.

The effects of electronic polarization in biomolecular interactions will differ depending on the local dielectric constant of the environment, such as in solvent, DNA, proteins, and membranes. Here the performance of the AMOEBA polarizable force field is evaluated under nonaqueous conditions by calculating the solvation free energies of small molecules in four common organic solvents. Results are compared with experimental data and equivalent simulations performed with the GAFF pairwise-additive force field.

Advanced Potential Energy Surfaces for Molecular Simulation.

Advanced potential energy surfaces are defined as theoretical models that explicitly include many-body effects that transcend the standard fixed-charge, pairwise-additive paradigm typically used in molecular simulation. However, several factors relating to their software implementation have precluded their widespread use in condensed-phase simulations: the computational cost of the theoretical models, a paucity of approximate models and algorithmic improvements that can ameliorate their cost, underdeveloped interfaces and limited dissemination in computational code bases that are widely used in the computational chemistry community, and software implementations that have not kept pace with modern high-performance computing (HPC) architectures, such as multicore CPUs and modern graphics processing units (GPUs). In this Feature Article we review recent progress made in these areas, including well-defined polarization approximations and new multipole electrostatic formulations, novel methods for solving the mutual polarization equations and increasing the MD time step, combining linear-scaling electronic structure methods with new QM/MM methods that account for mutual polarization between the two regions, and the greatly improved software deployment of these models and methods onto GPU and CPU hardware platforms.

Evaluating Parametrization Protocols for Hydration Free Energy Calculations with the AMOEBA Polarizable Force Field.

Hydration free energy (HFE) calculations are often used to assess the performance of biomolecular force fields and the quality of assigned parameters. The AMOEBA polarizable force field moves beyond traditional pairwise additive models of electrostatics and may be expected to improve upon predictions of thermodynamic quantities such as HFEs over and above fixed-point-charge models. The recent SAMPL4 challenge evaluated the AMOEBA polarizable force field in this regard but showed substantially worse results than those using the fixed-point-charge GAFF model.

Comparing experimental and computational alanine scanning techniques for probing a prototypical protein-protein interaction.

The central role of protein-protein interactions in a wide range of cellular processes makes them a target for research and drug discovery. A variety of methods, both experimental and theoretical, exist for probing protein interfaces for residues that affect activity and binding affinity. Using as an example a protein-protein complex between trypsin and a nine-residue synthetic peptide, we experimentally assay-binding affinities for a variety of mutants and determine their relative free energy of binding, ΔΔG, to rank the importance of interface residues to binding.