Accuracy of Alternate Nonpolarizable Force Fields for the Determination of Protein-Ligand Binding Affinities Dominated by Cation-π Interactions

Synthesis, insertion, and characterization of SARS-CoV-2 membrane protein within lipid bilayers

Throughout history, coronaviruses have posed challenges to both public health and the global economy; nevertheless, methods to combat them remain rudimentary, primarily due to the absence of experiments to understand the function of various viral components. Among these, membrane (M) proteins are one of the most elusive because of their small size and challenges with expression. Here, we report the development of an expression system to produce tens to hundreds of milligrams of M protein per liter of Escherichia coli culture. These large yields render many previously inaccessible structural and biophysical experiments feasible. Using cryo-electron microscopy and atomic force microscopy, we image and characterize individual membrane-incorporated M protein dimers and discover membrane thinning in the vicinity, which we validated with molecular dynamics simulations. Our results suggest that the resulting line tension, along with predicted induction of local membrane curvature, could ultimately drive viral assembly and budding.

Binding Thermodynamics of Fourth-Generation EGFR Inhibitors Revealed by Absolute Binding Free Energy Calculations

The overexpression or mutation of the kinase domain of the epidermal growth factor receptor (EGFR) is strongly associated with non-small-cell lung cancer (NSCLC). EGFR tyrosine kinase inhibitors (TKIs) have proven to be effective in treating NSCLC patients. However, EGFR mutations can result in drug resistance. To elucidate the mechanisms underlying this resistance and inform future drug development, we examined the binding affinities of BLU-945, a recently reported fourth-generation TKI, to wild-type EGFR (EGFRWT) and its double-mutant (L858R/T790M; EGFRDM) and triple-mutant (L858R/T790M/C797S; EGFRTM) forms. We compared the binding affinities of BLU-945, BLU-945 analogues, CH7233163 (another fourth-generation TKI), and erlotinib (a first-generation TKI) using absolute binding free energy calculations. Our findings reveal that BLU-945 and CH7233163 exhibit binding affinities to both EGFRDM and EGFRTM stronger than those of erlotinib, corroborating experimental data. We identified K745 and T854 as the key residues in the binding of fourth-generation EGFR TKIs. Electrostatic forces were the predominant driving force for the binding of fourth-generation TKIs to EGFR mutants. Furthermore, we discovered that the incorporation of piperidinol and sulfone groups in BLU-945 substantially enhanced its binding capacity to EGFR mutants. Our study offers valuable theoretical insights for optimizing fourth-generation EGFR TKIs.

Structural basis for ligand recognition and signaling of hydroxy-carboxylic acid receptor 2

Hydroxycarboxylic acid receptors (HCAR1, HCAR2, and HCAR3) transduce Gi/o signaling upon biding to molecules such as lactic acid, butyric acid and 3-hydroxyoctanoic acid, which are associated with lipolytic and atherogenic activity, and neuroinflammation. Although many reports have elucidated the function of HCAR2 and its potential as a therapeutic target for treating not only dyslipidemia but also neuroimmune disorders such as multiple sclerosis and Parkinson's disease, the structural basis of ligand recognition and ligand-induced Gi-coupling remains unclear. Here we report three cryo-EM structures of the human HCAR2-Gi signaling complex, each bound with different ligands: niacin, acipimox or GSK256073. All three agonists are held in a deep pocket lined by residues that are not conserved in HCAR1 and HCAR3. A distinct hairpin loop at the HCAR2 N-terminus and extra-cellular loop 2 (ECL2) completely enclose the ligand. These structures also reveal the agonist-induced conformational changes propagated to the G-protein-coupling interface during activation. Collectively, the structures presented here are expected to help in the design of ligands specific for HCAR2, leading to new drugs for the treatment of various diseases such as dyslipidemia and inflammation.

Structure and dynamics of differential ligand binding in the human ρ-type GABAA receptor

The neurotransmitter γ-aminobutyric acid (GABA) drives critical inhibitory processes in and beyond the nervous system, partly via ionotropic type-A receptors (GABAARs). Pharmacological properties of ρ-type GABAARs are particularly distinctive, yet the structural basis for their specialization remains unclear. Here, we present cryo-EM structures of a lipid-embedded human ρ1 GABAAR, including a partial intracellular domain, under apo, inhibited, and desensitized conditions. An apparent resting state, determined first in the absence of modulators, was recapitulated with the specific inhibitor (1,2,5,6-tetrahydropyridin-4-yl)methylphosphinic acid and blocker picrotoxin and provided a rationale for bicuculline insensitivity. Comparative structures, mutant recordings, and molecular simulations with and without GABA further explained the sensitized but slower activation of ρ1 relative to canonical subtypes. Combining GABA with picrotoxin also captured an apparent uncoupled intermediate state. This work reveals structural mechanisms of gating and modulation with applications to ρ-specific pharmaceutical design and to our biophysical understanding of ligand-gated ion channels.

Modern non-polarizable force fields diverge in modeling the enzyme-substrate complex of a canonical serine protease

Classical molecular dynamics simulation is a powerful and established method of modern computational chemistry. Being able to obtain accurate information on molecular behavior is crucial to get valuable insights into structure-function relationships that translate into fundamental findings and practical applications. Active sites of enzymes are known to be particularly intricate, therefore, simpler non-polarizable force fields may provide an inaccurate description. In this work, we addressed this hypothesis in a case of a canonical serine triad protease trypsin in its complex with a substrate-mimicking inhibitor. We tested six modern and popular force fields to find that significantly diverging results may be obtained. Amber FB-15 and OPLS-AA/M turned out to model the active site incorrectly. Amber ff19sb and ff15ipq demonstrated mixed performance. The best performing force fields were CHARMM36m and Amber ff99sb-ildn, therefore, they are recommended for use with this and related systems. We speculate that a similar lack of cross-force field convergence may be characteristic of other enzymatic systems. Therefore, we advocate for careful consideration of different force fields in any study within the field of computational enzymology.

Conformational dynamics and putative substrate extrusion pathways of the N-glycosylated outer membrane factor CmeC from Campylobacter jejuni

The outer membrane factor CmeC of the efflux machinery CmeABC plays an important role in conferring antibiotic and bile resistance to Campylobacter jejuni. Curiously, the protein is N-glycosylated, with the glycans playing a key role in the effective function of this system. In this work we have employed atomistic equilibrium molecular dynamics simulations of CmeC in a representative model of the C. jejuni outer membrane to characterise the dynamics of the protein and its associated glycans. We show that the glycans are more conformationally labile than had previously been thought. The extracellular loops of CmeC visit the open and closed states freely suggesting the absence of a gating mechanism on this side, while the narrow periplasmic entrance remains tightly closed, regulated via coordination to solvated cations. We identify several cation binding sites on the interior surface of the protein. Additionally, we used steered molecular dynamics simulations to elucidate translocation pathways for a bile acid and a macrolide antibiotic. These, and additional equilibrium simulations suggest that the anionic bile acid utilises multivalent cations to climb the ladder of acidic residues that line the interior surface of the protein.

Meta-Analysis Reveals That Absolute Binding Free-Energy Calculations Approach Chemical Accuracy

Systematic and quantitative analysis of the reliability of formally exact methods that in silico calculate absolute protein-ligand binding free energies remains lacking. Here, we provide, for the first time, evidence-based information on the reliability of these methods by statistically studying 853 cases from 34 different research groups through meta-analysis. The results show that formally exact methods approach chemical accuracy (error = 1.58 kcal/mol), even if people are challenging difficult tasks such as blind drug screening in recent years. The geometrical-pathway-based methods prove to possess a better convergence ability than the alchemical ones, while the latter have a larger application range. We also reveal the importance of always using the latest force fields to guarantee reliability and discuss the pros and cons of turning to an implicit solvent model in absolute binding free-energy calculations. Moreover, based on the meta-analysis, an evidence-based guideline for in silico binding free-energy calculations is provided.

Structure of SARS-CoV-2 membrane protein essential for virus assembly

The coronavirus membrane protein (M) is the most abundant viral structural protein and plays a central role in virus assembly and morphogenesis. However, the process of M protein-driven virus assembly are largely unknown. Here, we report the cryo-electron microscopy structure of the SARS-CoV-2 M protein in two different conformations. M protein forms a mushroom-shaped dimer, composed of two transmembrane domain-swapped three-helix bundles and two intravirion domains. M protein further assembles into higher-order oligomers. A highly conserved hinge region is key for conformational changes. The M protein dimer is unexpectedly similar to SARS-CoV-2 ORF3a, a viral ion channel. Moreover, the interaction analyses of M protein with nucleocapsid protein (N) and RNA suggest that the M protein mediates the concerted recruitment of these components through the positively charged intravirion domain. Our data shed light on the M protein-driven virus assembly mechanism and provide a structural basis for therapeutic intervention targeting M protein.

Accurate Description of Solvent-Exposed Salt Bridges with a Non-polarizable Force Field Incorporating Solvent Effects

The strength of salt bridges resulting from the interaction of cations and anions is modulated by their environment. However, polarization of the solvent molecules by the charged moieties makes the accurate description of cation-anion interactions in an aqueous solution by means of a pairwise additive potential energy function and classical combination rules particularly challenging. In this contribution, aiming at improving the representation of solvent-exposed salt-bridge interactions with an all-atom non-polarizable force field, we put forth here a parametrization strategy. First, the interaction of a cation and an anion is characterized by hybrid quantum mechanical/molecular mechanics (QM/MM) potential of mean force (PMF) calculations, whereby constantly exchanging solvent molecules around the ions are treated at the quantum mechanical level. The Lennard-Jones (LJ) parameters describing the salt-bridge ion pairs are then optimized to match the reference QM/MM PMFs through the so-called nonbonded FIX, or NBFIX, feature of the CHARMM force field. We apply the new set of parameters, coined CHARMM36m-SBFIX, to the calculation of association constants for the ammonium-acetate and guanidinium-acetate complexes, the osmotic pressures for glycine zwitterions, guanidinium, and acetate ions, and to the simulation of both folded and intrinsically disordered proteins. Our findings indicate that CHARMM36m-SBFIX improves the description of solvent-exposed salt-bridge interactions, both structurally and thermodynamically. However, application of this force field to the standard binding free-energy calculation of a protein-ligand complex featuring solvent-excluded salt-bridge interactions leads to a poor reproduction of the experimental value, suggesting that the parameters optimized in an aqueous solution cannot be readily transferred to describe solvent-excluded salt-bridge interactions. Put together, owing to their sensitivity to the environment, modeling salt-bridge interactions by means of a single, universal set of LJ parameters remains a daunting theoretical challenge.

Accurate determination of protein:ligand standard binding free energies from molecular dynamics simulations

Designing a reliable computational methodology to calculate protein:ligand standard binding free energies is extremely challenging. The large change in configurational enthalpy and entropy that accompanies the association of ligand and protein is notoriously difficult to capture in naive brute-force simulations. Addressing this issue, the present protocol rests upon a rigorous statistical mechanical framework for the determination of protein:ligand binding affinities together with the comprehensive Binding Free-Energy Estimator 2 (BFEE2) application software. With the knowledge of the bound state, available from experiments or docking, application of the BFEE2 protocol with a reliable force field supplies in a matter of days standard binding free energies within chemical accuracy, for a broad range of protein:ligand complexes. Limiting undesirable human intervention, BFEE2 assists the end user in preparing all the necessary input files and performing the post-treatment of the simulations towards the final estimate of the binding affinity.