CHARMM-GUI Nanodisc Builder for modeling and simulation of various nanodisc systems

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.

CHARMM-GUI PDB Reader and Manipulator: Covalent Ligand Modeling and Simulation

Molecular modeling and simulation serve an important role in exploring biological functions of proteins at the molecular level, which is complementary to experiments. CHARMM-GUI (https://www.charmm-gui.org) is a web-based graphical user interface that generates complex molecular simulation systems and input files, and we have been continuously developing and expanding its functionalities to facilitate various complex molecular modeling and make molecular dynamics simulations more accessible to the scientific community. Currently, covalent drug discovery emerges as a popular and important field. Covalent drug forms a chemical bond with specific residues on the target protein, and it has advantages in potency for its prolonged inhibition effects. Even though there are higher demands in modeling PDB protein structures with various covalent ligand types, proper modeling of covalent ligands remains challenging. This work presents a new functionality in CHARMM-GUI PDB Reader & Manipulator that can handle a diversity of ligand-amino acid linkage types, which is validated by a careful benchmark study using over 1,000 covalent ligand structures in RCSB PDB. We hope that this new functionality can boost the modeling and simulation study of covalent ligands.

The structure, self-assembly and dynamics of lipid nanodiscs revealed by computational approaches

Nanodisc technology is increasingly being used in structural, biochemical and biophysical studies of membrane proteins. The computational approaches have revealed many important features of nanodisc assembly, structures and dynamics. Therefore, we reviewed the application of computational approaches, especially molecular modeling and molecular dyncamics (MD) simulations, to characterize nanodiscs, including the structural models, assembly and disassembly, protocols for modeling, structural properties and dynamics, and protein-lipid interactions in nanodiscs. More amazing computational studies about nanodiscs are looked forward to in the future.

Electron videography of a lipid-protein tango

Biological phenomena, from enzymatic catalysis to synaptic transmission, originate in the structural transformations of biomolecules and biomolecular assemblies in liquid water. However, directly imaging these nanoscopic dynamics without probes or labels has been a fundamental methodological challenge. Here, we developed an approach for "electron videography"-combining liquid phase electron microscopy with molecular modeling-with which we filmed the nanoscale structural fluctuations of individual, suspended, and unlabeled membrane protein nanodiscs in liquid. Systematic comparisons with biochemical data and simulation indicate the graphene encapsulation involved can afford sufficiently reduced effects of the illuminating electron beam for these observations to yield quantitative fingerprints of nanoscale lipid-protein interactions. Our results suggest that lipid-protein interactions delineate dynamically modified membrane domains across unexpectedly long ranges. Moreover, they contribute to the molecular mechanics of the nanodisc as a whole in a manner specific to the protein within. Overall, this work illustrates an experimental approach to film, quantify, and understand biomolecular dynamics at the nanometer scale.

Lipid nanodisc scaffold and size alter the structure of a pentameric ligand-gated ion channel

Lipid nanodiscs have become a standard tool for studying membrane proteins, including using single particle cryo-electron microscopy (cryo-EM). We find that reconstituting the pentameric ligand-gated ion channel (pLGIC), Erwinia ligand-gated ion channel (ELIC), in different nanodiscs produces distinct structures by cryo-EM. The effect of the nanodisc on ELIC structure extends to the extracellular domain and agonist binding site. Additionally, molecular dynamic simulations indicate that nanodiscs of different size impact ELIC structure and that the nanodisc scaffold directly interacts with ELIC. These findings suggest that the nanodisc plays a crucial role in determining the structure of pLGICs, and that reconstitution of ion channels in larger nanodiscs may better approximate a lipid membrane environment.

Alchemical Free Energy Workflows for the Computation of Protein-Ligand Binding Affinities

Alchemical free energy methods can be used for the efficient computation of relative binding free energies during preclinical drug discovery stages. In recent years, this has been facilitated further by the implementation of workflows that enable non-experts to quickly and consistently set up the required simulations. Given the correct input structures, workflows handle the difficult aspects of setting up perturbations, including consistently defining the perturbable molecule, its atom mapping and topology generation, perturbation network generation, running of the simulations via different sampling methods, and analysis of the results. Different academic and commercial workflows are discussed, including FEW, FESetup, FEPrepare, CHARMM-GUI, Transformato, PMX, QLigFEP, TIES, ProFESSA, PyAutoFEP, BioSimSpace, FEP+, Flare, and Orion. These workflows differ in various aspects, such as mapping algorithms or enhanced sampling methods. Some workflows can accommodate more than one molecular dynamics (MD) engine and use external libraries for tasks. Differences between workflows can present advantages for different use cases, however a lack of interoperability of the workflows' components hinders systematic comparisons.

Anti-nanodisc antibodies specifically capture nanodiscs and facilitate molecular interaction kinetics studies for membrane protein

Nanodisc technology has dramatically advanced the analysis of molecular interactions for membrane proteins. A nanodisc is designed as a vehicle for membrane proteins that provide a native-like phospholipid environment and better thermostability in a detergent-free buffer. This enables the determination of the thermodynamic and kinetic parameters of small molecule binding by surface plasmon resonance. In this study, we generated a nanodisc specific anti-MSP (membrane scaffold protein) monoclonal antibody biND5 for molecular interaction analysis of nanodiscs. The antibody, biND5 bound to various types of nanodiscs with sub-nanomolar to nanomolar affinity. Epitope mapping analysis revealed specific recognition of 8 amino acid residues in the exposed helix-4 structure of MSP. Further, we performed kinetics binding analysis between adenosine A_2a receptor reconstituted nanodiscs and small molecule antagonist ZM241385 using biND5 immobilized sensor chips. These results show that biND5 facilitates the molecular interaction kinetics analysis of membrane proteins substituted in nanodiscs.

pyCHARMM: Embedding CHARMM Functionality in a Python Framework

CHARMM is rich in methodology and functionality as one of the first programs addressing problems of molecular dynamics and modeling of biological macromolecules and their partners, e.g., small molecule ligands. When combined with the highly developed CHARMM parameters for proteins, nucleic acids, small molecules, lipids, sugars, and other biologically relevant building blocks, and the versatile CHARMM scripting language, CHARMM has been a trendsetting platform for modeling studies of biological macromolecules. To further enhance the utility of accessing and using CHARMM functionality in increasingly complex workflows associated with modeling biological systems, we introduce pyCHARMM, Python bindings, functions, and modules to complement and extend the extensive set of modeling tools and methods already available in CHARMM. These include access to CHARMM function-generated variables associated with the system (psf), coordinates, velocities and forces, atom selection variables, and force field related parameters. The ability to augment CHARMM forces and energies with energy terms or methods derived from machine learning or other sources, written in Python, CUDA, or OpenCL and expressed as Python callable routines is introduced together with analogous functions callable during dynamics calculations. Integration of Python-based graphical engines for visualization of simulation models and results is also accessible. Loosely coupled parallelism is available for workflows such as free energy calculations, using MBAR/TI approaches or high-throughput multisite λ-dynamics (MSλD) free energy methods, string path optimization calculations, replica exchange, and molecular docking with a new Python-based CDOCKER module. CHARMM accelerated platform kernels through the CHARMM/OpenMM API, CHARMM/DOMDEC, and CHARMM/BLaDE API are also readily integrated into this Python framework. We anticipate that pyCHARMM will be a robust platform for the development of comprehensive and complex workflows utilizing Python and its extensive functionality as well as an optimal platform for users to learn molecular modeling methods and practices within a Python-friendly environment such as Jupyter Notebooks.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Modeling of flexible membrane-bound biomolecular complexes for solution small-angle scattering

Recent advances in protein expression protocols, sample handling, and experimental set up of small-angle scattering experiments have allowed users of the technique to structurally investigate biomolecules of growing complexity and structural disorder. Notable examples include intrinsically disordered proteins, multi-domain proteins and membrane proteins in suitable carrier systems. Here, we outline a modeling scheme for calculating the scattering profiles from such complex samples. This kind of modeling is necessary for structural information to be refined from the corresponding data. The scheme bases itself on a hybrid of classical form factor based modeling and the well-known spherical harmonics-based formulation of small-angle scattering amplitudes. Our framework can account for flexible domains alongside other structurally elaborate components of the molecular system in question. We demonstrate the utility of this modeling scheme through a recent example of a structural model of the growth hormone receptor membrane protein in a phospholipid bilayer nanodisc which is refined against experimental SAXS data. Additionally we investigate how the scattering profiles from the complex would appear under different scattering contrasts. For each contrast situation we discuss what structural information is contained and the related consequences for modeling of the data.

Capturing Lipid Nanodisc Shape and Properties Using a Continuum Elastic Theory

Lipid nanodiscs are nanometric bilayer patches enveloped by confining structures, commonly composed of membrane scaffolding proteins (MSPs). To resolve the interplay between MSP geometry, lipid confinement, and membrane material properties on the nanodisc shape, we apply a continuum elastic theory accounting for lipid bending, tilting, and area deformations. The equilibrium nanodisc shape is then determined by minimizing the elastic free energy functional. Analytic expressions derived under simplifying assumptions demonstrate that the nanodisc shape is sensitive to its size, lipid density, and the lipid tilt and thickness imposed at the contact with the MSP. Under matching physical parameters, these expressions quantitatively reproduce the shape of nanodiscs seen in molecular dynamics simulations, but only if lipid tilt is explicitly considered. We further demonstrate how the bending rigidity can be extracted from the membrane shape profile by fitting the numerically minimized full elastic functional to the membrane shape found in simulations. This fitting procedure faithfully informs on the bending rigidity of nanodiscs larger than ca. 5 nm in radius. The fitted profiles accurately reproduce the increase in bending modulus found using real-space fluctuation analysis of simulated nanodiscs and, for large nanodiscs, also accurately resolve its spatial variations. Our study shows how deformations in lipid patches confined in nanodiscs can be well described by a continuum elastic theory and how this fit can be used to determine local material properties from shape analysis of nanodiscs in simulations. This methodology could potentially allow direct determination of lipid properties from experiments, for example cryo-electron microscopy images of lipid nanodiscs, thereby allowing to guide the development of future nanodisc formulations with desired properties.

Reconfigurable Peptide Analogs of Apolipoprotein A-I Reveal Tunable Features of Nanodisc Assembly

Nanodisc (ND)-forming membrane scaffold proteins or peptides developed from apolipoprotein A-I (apoA-I) have led to considerable promise in structural biology and therapeutic applications. However, the rationale and regularity characteristics in peptide sequence design remain inconclusive. Here, we proposed a consensus-based normalization approach through the reversed engineering of apoA-IΔ1-45 to design reconfigurable apoA-I peptide analogs (APAs) for tunable ND assembly. We present extensive morphological validations and computational simulation analyses on divergent APA-NDs that are generated by our method. Fifteen divergent APAs were generated accordingly to study the assembly machinery of NDs. We show that APA designs exhibit multifactorial influence in terms of varying APA tandem repeats, sequence composition, and lipid-to-APA ratio to form tunable diameters of NDs. There is a strong positive correlation between DMPC-to-APA ratios and ND diameters. Longer APA with more tandem repeats tends to yield higher particle size homogeneity. Our results also suggest proline is a dispensable residue for the APA-ND formation. Interestingly, proline-rich substitution not only provides an inward-bending effect in forming smaller NDs but also induces the cumulative chain flexibility that enables larger ND formation at higher lipid ratios. Additionally, proline-tryptophan residues in APAs play a dominant role in forming larger NDs. Molecular simulation shows that enriched basic and acidic residues in APAs evoke abundant hydrogen-bond and salt bridge networks to reinforce the structural stability of APA-NDs. Together, our findings provide a rational basis for understanding APA design. The proposed model could be extended to other apolipoproteins for desired ND engineering.

Pharmacokinetic Study of Triptolide Nanocarrier in Transdermal Drug Delivery System-Combination of Experiment and Mathematical Modeling

Compared with traditional oral and injection administration, the transdermal administration of traditional Chinese medicine has distinctive characteristics and advantages, which can avoid the "first pass effect" of the liver and the destruction of the gastrointestinal tract, maintain a stable blood concentration, and prolong drug action time. However, the basic theory and technology research in transdermal drug delivery are relatively limited at present, especially regarding research on new carriers of transdermal drug delivery and pharmacokinetic studies of the skin, which has become a bottleneck of transdermal drug delivery development. Triptolide is one of the main active components of Tripterygium wilfordii, which displays activities against mouse models of polycystic kidney disease and pancreatic cancer but its physical properties and severe toxicity limit its therapeutic potential. Due to the previously mentioned advantages of transdermal administration, in this study, we performed a detail analysis of the pharmacokinetics of a new transdermal triptolide delivery system. Triptolide nanoemulsion gels were prepared and served as new delivery systems, and the ex vivo characteristics were described. The metabolic characteristics of the different triptolide transdermal drug delivery formulations were investigated via skin-blood synchronous microdialysis combined with LC/MS. A multiscale modeling framework, molecular dynamics and finite element modeling were adopted to simulate the transport process of triptolide in the skin and to explore the pharmacokinetics and mathematical patterns. This study shows that the three-layer model can be used for transdermal drug delivery system drug diffusion research. Therefore, it is profitable for transdermal drug delivery system design and the optimization of the dosage form. Based on the drug concentration of the in vivo microdialysis measurement technology, the diffusion coefficient of drugs in the skin can be more accurately measured, and the numerical results can be verified. Therefore, the microdialysis technique combined with mathematical modeling provides a very good platform for the further study of transdermal delivery systems. This research will provide a new technology and method for the study of the pharmacokinetics of traditional Chinese medicine transdermal drug delivery. It has important theoretical and practical significance in clarifying the metabolic transformation of percutaneous drug absorption and screening for appropriate drugs and dosage forms of transdermal drug delivery.

Nanodisc Lipids Exhibit Singular Behaviors Implying Critical Phenomena

Nanodiscs are broadly used for characterization of membrane proteins as they are generally assumed to provide a near-native environment. In fact, it is an open question whether the physical properties of lipids in nanodiscs and membrane vesicles of the same lipid composition are identical. Here, we investigate the properties of lipids (1,2-dipalmitoyl-sn-glycero-3-phosphocholine, 1,2-dilauroyl-sn-glycero-3-phosphocholine, and their mixtures) in two different sample types, nanodiscs and multilamellar vesicles, by means of spin-label electron spin resonance techniques. Our results provide a quantitative description of lipid dynamics and ordering, elucidating the molecular details of how lipids in the two sample types behave differently in response to temperature and lipid composition. We show that the properties of lipids are altered in nanodiscs such that the dissimilarity of the fluid and gel lipid phases is reduced, and the first-order phase transitions are largely abolished in nanodiscs. We unveil that the ensemble of lipids in the middle of a nanodisc bilayer, as probed by the end-chain spin-label 16-PC, is promoted to a state close to a miscibility critical point, thereby rendering the phase transitions continuous. Critical phenomena have recently been proposed to explain features of the heterogeneity in native cell membranes. Our results lay the groundwork for how to establish a near-native environment in nanodiscs with simple organization of lipid components.

Structures and Dynamics of Anionic Lipoprotein Nanodiscs

Nanolipoprotein particles known as nanodiscs (NDs) have emerged as versatile and powerful tools for the stabilization of membrane proteins permitting a plethora of structural and biophysical studies. Part of their allure is their flexibility to accommodate many types of lipids and precise control of the composition. However, little is known about how variations in lipid composition impact their structures and dynamics. Herein, we investigate how the introduction of the anionic lipid POPG into POPC NDs impacts these features. Small-angle X-ray and neutron scattering (SAXS and SANS) of variable-composition NDs are complemented with molecular dynamics simulations to interrogate how increasing the concern of POPG impacts the ND shape, structure of the lipid core, and the dynamics of the popular membrane scaffold protein, MSP1D1(-). A convenient benefit of including POPG is that it eliminates D₂O-induced aggregation observed in pure POPC NDs, permitting studies by SANS at multiple contrasts. SAXS and SANS data could be globally fit to a stacked elliptical cylinder model as well as an extension of the model that accounts for membrane curvature. Fitting to both models supports that the introduction of POPG results in strongly elliptical NDs; however, MD simulations predict the curvature of the membrane, thereby supporting the use of the latter model. Trends in the model-independent parameters suggest that increases in POPG reduce the conformational heterogeneity of the MSP1D1(-), which is in agreement with MD simulations that show that the incorporation of sufficient POPG suppresses disengagement of the N-terminal helix from the lipid core. These studies highlight novel structural changes in NDs in response to an anionic lipid and will inform the interpretation of future structural studies of membrane proteins embedded in NDs of mixed lipid composition.

Global fitting of multiple data frames from SEC-SAXS to investigate the structure of next-generation nanodiscs

The combination of online size-exclusion chromatography and small-angle X-ray scattering (SEC-SAXS) is rapidly becoming a key technique for structural investigations of elaborate biophysical samples in solution. Here, a novel model-refinement strategy centred around the technique is outlined and its utility is demonstrated by analysing data series from several SEC-SAXS experiments on phospholipid bilayer nanodiscs. Using this method, a single model was globally refined against many frames from the same data series, thereby capturing the frame-to-frame tendencies of the irradiated sample. These are compared with models refined in the traditional manner, in which refinement is based on the average profile of a set of consecutive frames from the same data series without an in-depth comparison of individual frames. This is considered to be an attractive model-refinement scheme as it considerably lowers the total number of parameters refined from the data series, produces tendencies that are automatically consistent between frames, and utilizes a considerably larger portion of the recorded data than is often performed in such experiments. Additionally, a method is outlined for correcting a measured UV absorption signal by accounting for potential peak broadening by the experimental setup.

Membrane Mixer: A Toolkit for Efficient Shuffling of Lipids in Heterogeneous Biological Membranes

Molecular dynamics (MD) simulations of biological membranes have achieved such levels of sophistication that are commonly used to predict unresolved structures and various properties of lipids and to substantiate experimental data. While achieving sufficient sampling of lipid dynamics remains a major challenge, a commonly used method to improve lipid sampling, e.g., in terms of specific interactions with membrane-associated proteins, is to randomize the initial arrangement of lipid constituents in multiple replicas of simulations, without changing the overall lipid composition of the membrane of interest. Here, we introduce a method that can rapidly generate multiple replicas of lipid bilayers with different spatial and conformational configurations for any given lipid composition. The underlying algorithm, which allows one to shuffle lipids at any desired level, relies on the application of an external potential, here referred to as the "carving potential", that removes clashes/entanglements before lipid positions are exchanged (shuffled), thereby minimizing the energy penalty due to abrupt lipid repositioning. The method is implemented as "Membrane Mixer Plugin (MMP) 1.0" in VMD, with a convenient graphical user interface that guides the user in setting various options and parameters. The plugin is fully automated and generates new membrane replicas more rapidly and conveniently than other analogous tools. The plugin and its capabilities introduced here can be extended to include additional features in future versions.

An in Silico Approach to Reveal the Nanodisc Formulation of Doxorubicin

Molecular dynamic behaviors of nanodisc (ND) formulations of free doxorubicin (DOX) and DOX conjugated lipid prodrug molecules were investigated by molecular dynamics (MD) simulations. We have unveiled how formulation design affects the drug release profile and conformational stability of ND assemblies. Our simulation results indicate that free DOX molecules loaded in the ND system experienced rapid dissociation due to the unfavorable orientation of DOX attached to the lipid surface. It is found that DOX tends to form aggregates with higher drug quantities. In contrast, lipidated DOX-prodrugs incorporated in ND formulations exhibited sufficient ND conformational stability. The drug loading capacity is dependent on the type of lipid molecules grafted on the DOX-prodrug, and the drug loading quantities in a fixed area of NDs follow the order: DOX-BMPH-MP > DOX-BMPH-TC > DOX-BMPH-PTE. To gain further insight into the dynamic characteristics of ND formulations governed by different kinds of lipidation, we investigated the conformational variation of ND components, intermolecular interactions, the solvent accessible surface area, and individual MSP1 residue flexibility. We found that the global conformational stability of DOX-prodrug-loaded ND assemblies is influenced by the molecular flexibility and lipidated forms of DOX-prodrug. We also found that the spontaneous self-aggregation of DOX-prodrugs with increasing quantities on ND could reduce the membrane fluidity and enhance the conformational stability of ND formulations.

CHARMM-GUI Ligand Designer for Template-Based Virtual Ligand Design in a Binding Site

Rational drug design involves a task of finding ligands that would bind to a specific target protein. This work presents CHARMM-GUI Ligand Designer that is an intuitive and interactive web-based tool to design virtual ligands that match the shape and chemical features of a given protein binding site. Ligand Designer provides ligand modification capabilities with 3D visualization that allow researchers to modify and redesign virtual ligands while viewing how the protein-ligand interactions are affected. Virtual ligands can also be parameterized for further molecular dynamics (MD) simulations and free energy calculations. Using 8 targets from 8 different protein classes in the directory of useful decoys, enhanced (DUD-E) data set, we show that Ligand Designer can produce similar ligands to the known active ligands in the crystal structures. Ligand Designer also produces stable protein-ligand complex structures when tested using short MD simulations. We expect that Ligand Designer can be a useful and user-friendly tool to design small molecules in any given potential ligand binding site on a protein of interest.

CHARMM-GUI Free Energy Calculator for Practical Ligand Binding Free Energy Simulations with AMBER

Alchemical free energy methods, such as free energy perturbation (FEP) and thermodynamic integration (TI), become increasingly popular and crucial for drug design and discovery. However, the system preparation of alchemical free energy simulation is an error-prone, time-consuming, and tedious process for a large number of ligands. To address this issue, we have recently presented CHARMM-GUI Free Energy Calculator that can provide input and postprocessing scripts for NAMD and GENESIS FEP molecular dynamics systems. In this work, we extended three submodules of Free Energy Calculator to work with the full suite of GPU-accelerated alchemical free energy methods and tools in AMBER, including input and postprocessing scripts. The BACE1 (β-secretase 1) benchmark set was used to validate the AMBER-TI simulation systems and scripts generated by Free Energy Calculator. The overall results of relatively large and diverse systems are almost equivalent with different protocols (unified and split) and with different timesteps (1, 2, and 4 fs), with R² > 0.9. More importantly, the average free energy differences between two protocols are small and reliable with four independent runs, with a mean unsigned error (MUE) below 0.4 kcal/mol. Running at least four independent runs for each pair with AMBER20 (and FF19SB/GAFF2.1/OPC force fields), we obtained a MUE of 0.99 kcal/mol and root-mean-square error of 1.31 kcal/mol for 58 alchemical transformations in comparison with experimental data. In addition, a set of ligands for T4-lysozyme was used to further validate our free energy calculation protocol whose results are close to experimental data (within 1 kcal/mol). In summary, Free Energy Calculator provides a user-friendly web-based tool to generate the AMBER-TI system and input files for high-throughput binding free energy calculations with access to the full set of GPU-accelerated alchemical free energy, enhanced sampling, and analysis methods in AMBER.

CHARMM-GUI LBS Finder & Refiner for Ligand Binding Site Prediction and Refinement

A protein performs its task by binding a variety of ligands in its local region that is also known as the ligand-binding-site (LBS). Therefore, accurate prediction, characterization, and refinement of LBS can facilitate protein functional annotations and structure-based drug design. In this work, we present CHARMM-GUI LBS Finder & Refiner (https://www.charmm-gui.org/input/lbsfinder) that predicts potential LBS, offers interactive features for local LBS structure analysis, and prepares various molecular dynamics (MD) systems and inputs by setting up distance restraint potentials for LBS structure refinement. LBS Finder & Refiner supports 5 different commonly used simulation programs, such as NAMD, AMBER, GROMACS, GENESIS, and OpenMM, for LBS structure refinement together with hydrogen mass repartitioning. The capability of LBS Finder & Refiner is illustrated through LBS structure predictions and refinements of 48 modeled and 20 apo benchmark target proteins. Overall, successful LBS structure predictions and refinements are seen in our benchmark tests. We hope that LBS Finder & Refiner is useful to predict, characterize, and refine potential LBS on any given protein of interest.

Membrane perturbation by lipidated Atg8 underlies autophagosome biogenesis

Autophagosome biogenesis is an essential feature of autophagy. Lipidation of Atg8 plays a critical role in this process. Previous in vitro studies identified membrane tethering and hemi-fusion/fusion activities of Atg8, yet definitive roles in autophagosome biogenesis remained controversial. Here, we studied the effect of Atg8 lipidation on membrane structure. Lipidation of Saccharomyces cerevisiae Atg8 on nonspherical giant vesicles induced dramatic vesicle deformation into a sphere with an out-bud. Solution NMR spectroscopy of Atg8 lipidated on nanodiscs identified two aromatic membrane-facing residues that mediate membrane-area expansion and fragmentation of giant vesicles in vitro. These residues also contribute to the in vivo maintenance of fragmented vacuolar morphology under stress in fission yeast, a moonlighting function of Atg8. Furthermore, these aromatic residues are crucial for the formation of a sufficient number of autophagosomes and regulate autophagosome size. Together, these data demonstrate that Atg8 can cause membrane perturbations that underlie efficient autophagosome biogenesis.

General Protocol for Constructing Molecular Models of Nanodiscs

Nanodisc technology is increasingly being applied for structural and biophysical studies of membrane proteins. In this work, we present a general protocol for constructing molecular models of nanodiscs for molecular dynamics simulations. The protocol is written in python and based on geometric equations, making it fast and easy to modify, enabling automation and customization of nanodiscs in silico. The novelty being the ability to construct any membrane scaffold protein (MSP) variant fast and easy given only an input sequence. We validated and tested the protocol by simulating seven different nanodiscs of various sizes and with different membrane scaffold proteins, both circularized and noncircularized. The structural and biophysical properties were analyzed and shown to be in good agreement with previously reported experimental data and simulation studies.

Can di-4-ANEPPDHQ reveal the structural differences between nanodiscs and liposomes?

The potential-sensitive di-4-ANEPPDHQ dye is presently gaining popularity in structural studies of the lipid bilayer. Within the bilayer, dye environmental sensitivity originates from the excitation induced charge redistribution and is usually attributed to solvent relaxation. Here, di-4-ANEPPDHQ is utilized to compare the structure of neutral and negatively charged lipid bilayers between two model systems: the nanodiscs and the liposomes. Using the well-established approach of measuring solvatochromic shifts of the steady-state spectra to study the bilayer structural changes has proved insufficient in this case. By applying an in-depth analysis of time-resolved fluorescence decays and emission spectra, we distinguished and characterized two and three distinct emissive di-4-ANEPPDHQ species in the liposomes and the nanodiscs, respectively. These emissive species were ascribed to the dual emission of the dye rather than to solvent relaxation. An additional, long-lived component present in the nanodiscs was associated with a unique domain of high order, postulated recently. Our results reveal that the di-4-ANEPPDHQ steady-state fluorescence should be interpreted with caution. With the experimental approach presented here, the di-4-ANEPPDHQ sensitivity was improved. We confirmed that the bilayer structure is, indeed, altered in the nanodiscs. Moreover, molecular dynamic simulations showed a distribution of the probe in the nanodiscs plane, which is sensitive to lipid composition. In POPC nanodiscs, probe frequently interacts with MSP, while in POPC-POPG nanodiscs, such interactions are rare. We did not observe, however, any impact of those interactions on the probe fluorescence.

DNA-Mediated Stack Formation of Nanodiscs

Membrane-scaffolding proteins (MSPs) derived from apolipoprotein A-1 have become a versatile tool in generating nano-sized discoidal membrane mimetics (nanodiscs) for membrane protein research. Recent efforts have aimed at exploiting their controlled lipid protein ratio and size distribution to arrange membrane proteins in regular supramolecular structures for diffraction studies. Thereby, direct membrane protein crystallization, which has remained the limiting factor in structure determination of membrane proteins, would be circumvented. We describe here the formation of multimers of membrane-scaffolding protein MSP1D1-bounded nanodiscs using the thiol reactivity of engineered cysteines. The mutated positions N42 and K163 in MSP1D1 were chosen to support chemical modification as evidenced by fluorescent labeling with pyrene. Minimal interference with the nanodisc formation and structure was demonstrated by circular dichroism spectroscopy, differential light scattering and size exclusion chromatography. The direct disulphide bond formation of nanodiscs formed by the MSP1D1_N42C variant led to dimers and trimers with low yield. In contrast, transmission electron microscopy revealed that the attachment of oligonucleotides to the engineered cysteines of MSP1D1 allowed the growth of submicron-sized tracts of stacked nanodiscs through the hybridization of nanodisc populations carrying complementary strands and a flexible spacer.

MPBuilder: A PyMOL Plugin for Building and Refinement of Solubilized Membrane Proteins Against Small Angle X-ray Scattering Data

•

SAXS is a unique tool to study MP systems in the native-like environment.

•

The lack of user-friendly software hinders SAXS data interpretation of MPs.

•

We present a novel plugin for the widely used molecular visualization system PyMol.

•

MPBuilder allows to build and refine models against SAXS data.

Membrane proteins (MPs) are the target of numerous structural and functional studies in biological and medical/pharmaceutical sciences. Strategies for the high-throughput structural analysis of MPs and of their perturbations driven by ligands having potential therapeutic applications are uncommon, often requiring scaled up crystallization, electron microscopy, and nuclear magnetic resonance (NMR) efforts. Small-angle X-ray scattering (SAXS) provides a rapid means to study low resolution structures and conformational changes of native MPs in solution without cumbersome sample preparations/treatment. The method requires the MPs solubilized in an appropriate medium (eg. detergents, mixed micelles and nanodiscs) and reliable and robust models are needed to describe the relevant complexes. Here we present MPBuilder, a simple and versatile tool for the generation and refinement of all-atom MP systems in the popular software PyMOL, an environment familiar to most biologists. MPBuilder provides building capability for protein-detergent, bicelle, and lipid-scaffold (saposin nanoparticles, nanodiscs) complexes and links this to the ATSAS software package modules for model refinement and validation against the SAXS data.

Properties of DNA- and Protein-Scaffolded Lipid Nanodiscs

The properties of natural lipid bilayers are vital to the regulation of many membrane proteins. Scaffolded nanodiscs provide an in vitro lipid bilayer platform to host membrane proteins in an environment that approximates native lipid bilayers. However, the properties of scaffold-enclosed bilayers may depart significantly from those of bulk cellular membranes. Therefore, to improve the usefulness of nanodiscs it is essential to understand the properties of lipids restricted by scaffolds. We used computational molecular dynamics and modeling approaches to understand the effects of nanodisc size, scaffold type (DNA or protein), and hydrophobic modification of DNA scaffolds on bilayer stability and degree to which the properties of enclosed bilayers approximate bulk bilayers. With respect to achieving bulk bilayer behavior, we found that charge neutralization of DNA scaffolds was more important than the total hydrophobic content of their modifications: bilayer properties were better for scaffolds having a large number of short alkyl chains than those having fewer long alkyl chains. Further, complete charge neutralization of DNA scaffolds enabled better lipid binding, and more stable bilayers, as shown by steered molecular dynamics simulations that measured the force required to dislodge scaffolds from lipid bilayer patches. Considered together, our simulations provide a guide to the design of DNA-scaffolded nanodiscs suitable for studying membrane proteins.

Distortion of the bilayer and dynamics of the BAM complex in lipid nanodiscs

The β-barrel assembly machinery (BAM) catalyses the folding and insertion of β-barrel outer membrane proteins (OMPs) into the outer membranes of Gram-negative bacteria by mechanisms that remain unclear. Here, we present an ensemble of cryoEM structures of the E. coli BamABCDE (BAM) complex in lipid nanodiscs, determined using multi-body refinement techniques. These structures, supported by single-molecule FRET measurements, describe a range of motions in the BAM complex, mostly localised within the periplasmic region of the major subunit BamA. The β-barrel domain of BamA is in a 'lateral open' conformation in all of the determined structures, suggesting that this is the most energetically favourable species in this bilayer. Strikingly, the BAM-containing lipid nanodisc is deformed, especially around BAM's lateral gate. This distortion is also captured in molecular dynamics simulations, and provides direct structural evidence for the lipid 'disruptase' activity of BAM, suggested to be an important part of its functional mechanism.

Nanodiscs: A toolkit for membrane protein science

Membrane proteins are involved in numerous vital biological processes, including transport, signal transduction and the enzymes in a variety of metabolic pathways. Integral membrane proteins account for up to 30% of the human proteome and they make up more than half of all currently marketed therapeutic targets. Unfortunately, membrane proteins are inherently recalcitrant to study using the normal toolkit available to scientists, and one is most often left with the challenge of finding inhibitors, activators and specific antibodies using a denatured or detergent solubilized aggregate. The Nanodisc platform circumvents these challenges by providing a self-assembled system that renders typically insoluble, yet biologically and pharmacologically significant, targets such as receptors, transporters, enzymes, and viral antigens soluble in aqueous media in a native-like bilayer environment that maintain a target's functional activity. By providing a bilayer surface of defined composition and structure, Nanodiscs have found great utility in the study of cellular signaling complexes that assemble on a membrane surface. Nanodiscs provide a nanometer scale vehicle for the in vivo delivery of amphipathic drugs, therapeutic lipids, tethered nucleic acids, imaging agents and active protein complexes. This means for generating nanoscale lipid bilayers has spawned the successful use of numerous other polymer and peptide amphipathic systems. This review, in celebration of the Anfinsen Award, summarizes some recent results and provides an inroad into the current and historical literature.

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

Alchemical free energy simulations have long been utilized to predict free energy changes for binding affinity and solubility of small molecules. However, while the theoretical foundation of these methods is well established, seamlessly handling many of the practical aspects regarding the preparation of the different thermodynamic end states of complex molecular systems and the numerous processing scripts often remains a burden for successful applications. In this work, we present CHARMM-GUI Free Energy Calculator (http://www.charmm-gui.org/input/fec) that provides various alchemical free energy perturbation molecular dynamics (FEP/MD) systems with input and post-processing scripts for NAMD and GENESIS. Four submodules are available: Absolute Ligand Binder (for absolute ligand binding FEP/MD), Relative Ligand Binder (for relative ligand binding FEP/MD), Absolute Ligand Solvator (for absolute ligand solvation FEP/MD), and Relative Ligand Solvator (for relative ligand solvation FEP/MD). Each module is designed to build multiple systems of a set of selected ligands at once for high-throughput FEP/MD simulations. The capability of Free Energy Calculator is illustrated by absolute and relative solvation FEP/MD of a set of ligands and absolute and relative binding FEP/MD of a set of ligands for T4-lysozyme in solution and the adenosine A_2A receptor in a membrane. The calculated free energy values are overall consistent with the experimental and published free energy results (within ∼1 kcal/mol). We hope that Free Energy Calculator is useful to carry out high-throughput FEP/MD simulations in the field of biomolecular sciences and drug discovery.

Thermally induced structural organization of nanodiscs by coarse-grained molecular dynamics simulations

Membrane scaffold proteins (MSP) nanodiscs have been extensively used in structural study of membrane proteins. In cryo-EM, an incorporation of target proteins into nanodiscs is conducted under a rapid change from cryogenic to ambient temperatures. We present a coarse-grained molecular dynamics (CGMD) study for investigating an effect of temperature on the structural organization of DPPC-nanodisc and POPC-nanodisc. A non-monotonic response of physical quantities (i.e. the lipid order parameter, nanodisc flatness, structural change, solvation property, radius of gyration) with increase in temperature (T = 200-350 K) is found to be associated with the gel-ripple-liquid crystalline phase change within nanodiscs. The reorganization of lipids upon temperature variation induced conformational changes of MSP to minimize hydrophobic exposure of the lipid membrane to an aqueous environment. Structural response to temperature is different to a certain extent between the saturated DPPC and unsaturated POPC.

Confinement in Nanodiscs Anisotropically Modifies Lipid Bilayer Elastic Properties

Lipid nanodiscs are small synthetic lipid bilayer structures that are stabilized in solution by special circumscribing (or scaffolding) proteins or polymers. Because they create native-like environments for transmembrane proteins, lipid nanodiscs have become a powerful tool for structural determination of this class of systems when combined with cryo-electron microscopy or nuclear magnetic resonance. The elastic properties of lipid bilayers determine how the lipid environment responds to membrane protein perturbations, and how the lipid in turn modifies the conformational state of the embedded protein. However, despite the abundant use of nanodiscs in determining membrane protein structure, the elastic material properties of even pure lipid nanodiscs (i.e., without embedded proteins) have not yet been quantitatively investigated. A major hurdle is due to the inherently nonlocal treatment of the elastic properties of lipid systems implemented by most existing methods, both experimental and computational. In addition, these methods are best suited for very large “infinite” size lipidic assemblies, or ones that contain periodicity, in the case of simulations. We have previously described a computational analysis of molecular dynamics simulations designed to overcome these limitations, so it allows quantification of the bending rigidity (K_C) and tilt modulus (κ_t) on a local scale even for finite, nonperiodic systems, such as lipid nanodiscs. Here we use this computational approach to extract values of K_C and κ_t for a set of lipid nanodisc systems that vary in size and lipid composition. We find that the material properties of lipid nanodiscs are different from those of infinite bilayers of corresponding lipid composition, highlighting the effect of nanodisc confinement. Nanodiscs tend to show higher stiffness than their corresponding macroscopic bilayers, and moreover, their material properties vary spatially within them. For small-size MSP1 nanodiscs, the stiffness decreases radially, from a value that is larger in their center than the moduli of the corresponding bilayers by a factor of ∼2–3. The larger nanodiscs (MSP1E3D1 and MSP2N2) show milder spatial changes of moduli that are composition dependent and can be maximal in the center or at some distance from it. These trends in moduli correlate with spatially varying structural properties, including the area per lipid and the nanodisc thickness. Finally, as has previously been reported, nanodiscs tend to show deformations from perfectly flat circular geometries to varying degrees, depending on size and lipid composition. The modulations of lipid elastic properties that we find should be carefully considered when making structural and functional inferences concerning embedded proteins.

Membrane lipids are both the substrates and a mechanistically responsive environment of TMEM16 scramblase proteins

Recent discoveries about functional mechanisms of proteins in the TMEM16 family of phospholipid scramblases have illuminated the dual role of the membrane as both the substrate and a mechanistically responsive environment in the wide range of physiological processes and genetic disorders in which they are implicated. This is highlighted in the review of recent findings from our collaborative investigations of molecular mechanisms of TMEM16 scramblases that emerged from iterative functional, structural, and computational experimentation. In the context of this review, we present new MD simulations and trajectory analyses motivated by the fact that new structural information about the TMEM16 scramblases is emerging from cryo-EM determinations in lipid nanodiscs. Because the functional environment of these proteins in in vivo and in in vitro is closer to flat membranes, we studied comparatively the responses of the membrane to the TMEM16 proteins in flat membranes and nanodiscs. We find that bilayer shapes in the nanodiscs are very different from those observed in the flat membrane systems, but the function-related slanting of the membrane observed at the nhTMEM16 boundary with the protein is similar in the nanodiscs and in the flat bilayers. This changes, however, in the bilayer composed of longer-tail lipids, which is thicker near the phospholipid translocation pathway, which may reflect an enhanced tendency of the long tails to penetrate the pathway and create, as shown previously, a nonconductive environment. These findings support the correspondence between the mechanistic involvement of the lipid environment in the flat membranes, and the nanodiscs. © 2019 Wiley Periodicals, Inc.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

CHARMM-GUI DEER facilitator for spin-pair distance distribution calculations and preparation of restrained-ensemble molecular dynamics simulations

The double electron-electron resonance (DEER) is a powerful structural biology technique to obtain distance information in the range of 18 to 80 å by measuring the dipolar coupling between two unpaired electron spins. The distance distributions obtained from the experiment provide valuable structural information about the protein in its native environment that can be exploited using restrained ensemble molecular dynamics (reMD) simulations. We present a new tool DEER Facilitator in CHARMM-GUI that consists of two modules Spin-Pair Distributor and reMD Prepper to setup simulations that utilize information from DEER experiments. Spin-Pair Distributor provides a web-based interface to calculate the spin-pair distance distribution of labeled sites in a protein using MD simulations. The calculated distribution can be used to guide the selection of the labeling sites in experiments as well as validate different protein structure models. reMD Prepper facilities the setup of reMD simulations using different types of spin labels in four different environments including vacuum, solution, micelle, and bilayer. The applications of these two modules are demonstrated with several test cases. Spin-Pair Distributor and reMD Prepper are available at http://www.charmm-gui.org/input/deer and http://www.charmm-gui.org/input/deerre. DEER Facilitator is expected to facilitate advanced biomolecular modeling and simulation, thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems based on experimental DEER data. © 2019 Wiley Periodicals, Inc.

Frequency selective coherence transfer NMR spectroscopy to study the structural dynamics of high molecular weight proteins

Multidimensional nuclear magnetic resonance (NMR) spectroscopy has enabled detailed characterizations of protein structures and dynamics that are closely linked to functions. However, it leads to a large sensitivity loss in applications to high molecular weight proteins, which is caused by spin relaxation during the frequency discrimination period in the indirect dimension. Here, we describe a selective coherence transfer scheme, which enables us to selectively observe ¹H nuclei bonded to ¹⁵N or ¹³C nuclei with specified resonance frequencies. By utilizing this scheme, we achieved a 2.5- to 6-fold increase in signal height per unit of time with this scheme by avoiding the relaxation loss in the indirect dimension, as compared to the conventional two-dimensional heteronuclear correlation spectroscopy. We also demonstrated the effectiveness of this approach with applications to the membrane protein KirBac1.1, and characterized the functionally relevant conformational exchange process in both detergent micelles and a reconstituted membrane environment, corresponding to the apparent molecular masses of 220 kDa and 300 kDa, respectively.

CHARMM-GUI Glycan Modeler for modeling and simulation of carbohydrates and glycoconjugates

Characterizing glycans and glycoconjugates in the context of three-dimensional structures is important in understanding their biological roles and developing efficient therapeutic agents. Computational modeling and molecular simulation have become an essential tool complementary to experimental methods. Here, we present a computational tool, Glycan Modeler for in silico N-/O-glycosylation of the target protein and generation of carbohydrate-only systems. In our previous study, we developed Glycan Reader, a web-based tool for detecting carbohydrate molecules from a PDB structure and generation of simulation system and input files. As integrated into Glycan Reader in CHARMM-GUI, Glycan Modeler (Glycan Reader & Modeler) enables to generate the structures of glycans and glycoconjugates for given glycan sequences and glycosylation sites using PDB glycan template structures from Glycan Fragment Database (http://glycanstructure.org/fragment-db). Our benchmark tests demonstrate the universal applicability of Glycan Reader & Modeler to various glycan sequences and target proteins. We also investigated the structural properties of modeled glycan structures by running 2-μs molecular dynamics simulations of HIV envelope protein. The simulations show that the modeled glycan structures built by Glycan Reader & Modeler have the similar structural features compared to the ones solved by X-ray crystallography. We also describe the representative examples of glycoconjugate modeling with video demos to illustrate the practical applications of Glycan Reader & Modeler. Glycan Reader & Modeler is freely available at http://charmm-gui.org/input/glycan.