Determination of Hydrophobic Lengths of Membrane Proteins with the HDGB Implicit Membrane Model

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Common synaptic phenotypes arising from diverse mutations in the human NMDA receptor subunit GluN2A

Dominant mutations in the human gene GRIN2A, encoding NMDA receptor (NMDAR) subunit GluN2A, make a significant and growing contribution to the catalogue of published single-gene epilepsies. Understanding the disease mechanism in these epilepsy patients is complicated by the surprising diversity of effects that the mutations have on NMDARs. Here we have examined the cell-autonomous effect of five GluN2A mutations, 3 loss-of-function and 2 gain-of-function, on evoked NMDAR-mediated synaptic currents (NMDA-EPSCs) in CA1 pyramidal neurons in cultured hippocampal slices. Despite the mutants differing in their functional incorporation at synapses, prolonged NMDA-EPSC current decays (with only marginal changes in charge transfer) were a common effect for both gain- and loss-of-function mutants. Modelling NMDA-EPSCs with mutant properties in a CA1 neuron revealed that the effect of GRIN2A mutations can lead to abnormal temporal integration and spine calcium dynamics during trains of concerted synaptic activity. Investigations beyond establishing the molecular defects of GluN2A mutants are much needed to understand their impact on synaptic transmission.

The cell-autonomous effect of five severe loss- or gain-of-function GluN2A (NMDA receptor) mutations is assessed on evoked NMDAR-mediated synaptic currents in CA1 pyramidal neurons in cultured mouse hippocampal slices. Data and modelling suggest that mutant-like NMDA-EPSCs can lead to abnormal temporal summation and spine calcium dynamics.

Spatial arrangement of proteins in planar and curved membranes by PPM 3.0

Cellular protrusions, invaginations, and many intracellular organelles have strongly curved membrane regions. Transmembrane and peripheral membrane proteins that induce, sense, or stabilize such regions cannot be properly fitted into a single flat bilayer. To treat such proteins, we developed a new method and a web tool, PPM 3.0, for positioning proteins in curved or planar, single or multiple membranes. This method determines the energetically optimal spatial position, the hydrophobic thickness, and the radius of intrinsic curvature of a membrane-deforming protein structure by arranging it in a single or several sphere-shaped or planar membrane sections. In addition, it can define the lipid-embedded regions of a protein that simultaneously spans several membranes or determine the optimal position of a peptide in a spherical micelle. The PPM 3.0 web server operates with 17 types of biological membranes and 4 types of artificial bilayers. It is publicly available at https://opm.phar.umich.edu/ppm_server3. PPM 3.0 was applied to identify and characterize arrangements in membranes of 128 proteins with a significant intrinsic curvature, such as BAR domains, annexins, Piezo, and MscS mechanosensitive channels, cation-chloride cotransporters, as well as mitochondrial ATP synthases, calcium uniporters, and TOM complexes. These proteins form large complexes that are mainly localized in mitochondria, plasma membranes, and endosomes. Structures of bacterial drug efflux pumps, AcrAB-TolC, MexAB-OrpM, and MacAB-TolC, were positioned in both membranes of the bacterial cell envelop, while structures of multimeric gap-junction channels were arranged in two opposed cellular membranes.

Diverse Scientific Benchmarks for Implicit Membrane Energy Functions

Energy functions are fundamental to biomolecular modeling. Their success depends on robust physical formalisms, efficient optimization, and high-resolution data for training and validation. Over the past 20 years, progress in each area has advanced soluble protein energy functions. Yet, energy functions for membrane proteins lag behind due to sparse and low-quality data, leading to overfit tools. To overcome this challenge, we assembled a suite of 12 tests on independent data sets varying in size, diversity, and resolution. The tests probe an energy function's ability to capture membrane protein orientation, stability, sequence, and structure. Here, we present the tests and use the franklin2019 energy function to demonstrate them. We then identify areas for energy function improvement and discuss potential future integration with machine-learning-based optimization methods. The tests are available through the Rosetta Benchmark Server (https://benchmark.graylab.jhu.edu/) and GitHub (https://github.com/rfalford12/Implicit-Membrane-Energy-Function-Benchmark).

Protein Structure Prediction and Design in a Biologically Realistic Implicit Membrane

Protein design is a powerful tool for elucidating mechanisms of function and engineering new therapeutics and nanotechnologies. Although soluble protein design has advanced, membrane protein design remains challenging because of difficulties in modeling the lipid bilayer. In this work, we developed an implicit approach that captures the anisotropic structure, shape of water-filled pores, and nanoscale dimensions of membranes with different lipid compositions. The model improves performance in computational benchmarks against experimental targets, including prediction of protein orientations in the bilayer, ΔΔG calculations, native structure discrimination, and native sequence recovery. When applied to de novo protein design, this approach designs sequences with an amino acid distribution near the native amino acid distribution in membrane proteins, overcoming a critical flaw in previous membrane models that were prone to generating leucine-rich designs. Furthermore, the proteins designed in the new membrane model exhibit native-like features including interfacial aromatic side chains, hydrophobic lengths compatible with bilayer thickness, and polar pores. Our method advances high-resolution membrane protein structure prediction and design toward tackling key biological questions and engineering challenges.

Generalized Born Implicit Solvent Models for Biomolecules

It would often be useful in computer simulations to use an implicit description of solvation effects, instead of explicitly representing the individual solvent molecules. Continuum dielectric models often work well in describing the thermodynamic aspects of aqueous solvation and can be very efficient compared to the explicit treatment of the solvent. Here, we review a particular class of so-called fast implicit solvent models, generalized Born (GB) models, which are widely used for molecular dynamics (MD) simulations of proteins and nucleic acids. These approaches model hydration effects and provide solvent-dependent forces with efficiencies comparable to molecular-mechanics calculations on the solute alone; as such, they can be incorporated into MD or other conformational searching strategies in a straightforward manner. The foundations of the GB model are reviewed, followed by examples of newer, emerging models and examples of important applications. We discuss their strengths and weaknesses, both for fidelity to the underlying continuum model and for the ability to replace explicit consideration of solvent molecules in macromolecular simulations.

Structure refinement of membrane proteins via molecular dynamics simulations

A refinement protocol based on physics-based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge-based or implicit membrane-based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid-facing residues. Scoring with knowledge-based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane-based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.