Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism.

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

A molecular switch controls the impact of cholesterol on a Kir channel.

SignificanceCholesterol is one of the main components found in plasma membranes and is involved in lipid-dependent signaling enabled by integral membrane proteins such as inwardly rectifying potassium (Kir) channels. Similar to other ion channels, most of the Kir channels are down-regulated by cholesterol. One of the very few notable exceptions is Kir3.4, which is up-regulated by this important lipid. Here, we discovered and characterized a molecular switch that controls the impact (up-regulation vs. down-regulation) of cholesterol on Kir3.4. Our results provide a detailed molecular mechanism of tunable cholesterol regulation of a potassium channel.

Molecular dynamics simulations of lipid-protein interactions in SLC4 proteins.

The SLC4 family of secondary bicarbonate transporters is responsible for the transport of HCO3-, CO32-, Cl-, Na+, K+, NH3, and H+, which are necessary for regulation of pH and ion homeostasis. They are widely expressed in numerous tissues throughout the body and function in different cell types with different membrane properties. Potential lipid roles in SLC4 function have been reported in experimental studies, focusing mostly on two members of the family: AE1 (Cl-/HCO3- exchanger) and NBCe1 (Na+-CO32-cotransporter). Previous computational studies of the outward-facing state of AE1 with model lipid membranes revealed enhanced protein-lipid interactions between cholesterol (CHOL) and phosphatidylinositol bisphosphate (PIP2). However, the protein-lipid interactions in other members of the family and other conformation states are still poorly understood and this precludes the detailed studies of a potential regulatory role for lipids in the SLC4 family. In this work, we performed coarse-grained and atomistic molecular dynamics simulations on three members of the SLC4 family with different transport modes: AE1, NBCe1, and NDCBE (an Na+-CO32-/Cl- exchanger), in model HEK293 membranes consisting of CHOL, PIP2, phosphatidylcholine, phosphatidylethanolamine, phosphatidylserine, and sphingomyelin. The recently resolved inward-facing state of AE1 was also included in the simulations. Lipid-protein contact analysis of the simulated trajectories was performed with the ProLint server, which provides a multitude of visualization tools for illustration of areas of enhanced lipid-protein contact and identification of putative lipid binding sites within the protein matrix. We observed enrichment of CHOL and PIP2 around all proteins with subtle differences in their distribution depending on the protein type and conformation state. Putative binding sites were identified for CHOL, PIP2, phosphatidylcholine, and sphingomyelin in the three studied proteins, and their potential roles in the SLC4 transport function, conformational transition, and protein dimerization are discussed.

Martini 3: a general purpose force field for coarse-grained molecular dynamics.

The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

An implementation of the Martini coarse-grained force field in OpenMM.

We describe a complete implementation of Martini 2 and Martini 3 in the OpenMM molecular dynamics software package. Martini is a widely used coarse-grained force field with applications in biomolecular simulation, materials, and broader areas of chemistry. It is implemented as a force field but makes extensive use of facilities unique to the GROMACS software, including virtual sites and bonded terms that are not commonly used in standard atomistic force fields. OpenMM is a flexible molecular dynamics package widely used for methods development and is competitive in speed on GPUs with other commonly used packages. OpenMM has facilities to easily implement new force field terms, external forces and fields, and other nonstandard features, which we use to implement all force field terms used in Martini 2 and Martini 3. This allows Martini simulations, starting with GROMACS topology files that are processed by custom scripts, with all the added flexibility of OpenMM. We provide a GitHub repository with test cases, compare accuracy and performance between GROMACS and OpenMM, and discuss the limitations of our implementation in terms of direct comparison with GROMACS. We describe a use case that implements the Modeling Employing Limited Data method to apply experimental constraints in a Martini simulation to efficiently determine the structure of a protein complex. We also discuss issues and a potential solution with the Martini 2 topology for cholesterol.

Can two wrongs make a right? F508del-CFTR ion channel rescue by second-site mutations in its transmembrane domains.

Deletion of phenylalanine 508 (F508del) in the cystic fibrosis transmembrane conductance regulator (CFTR) anion channel is the most common cause of cystic fibrosis. The F508 residue is located on nucleotide-binding domain 1 (NBD1) in contact with the cytosolic extensions of the transmembrane helices, in particular intracellular loop 4 (ICL4). To investigate how absence of F508 at this interface impacts the CFTR protein, we carried out a mutagenesis scan of ICL4 by introducing second-site mutations at 11 positions in cis with F508del. Using an image-based fluorescence assay, we measured how each mutation affected membrane proximity and ion-channel function. The scan strongly validated the effectiveness of R1070W at rescuing F508del defects. Molecular dynamics simulations highlighted two features characterizing the ICL4/NBD1 interface of F508del/R1070W-CFTR: flexibility, with frequent transient formation of interdomain hydrogen bonds, and loosely stacked aromatic sidechains (F1068, R1070W, and F1074, mimicking F1068, F508, and F1074 in WT CFTR). F508del-CFTR displayed a distorted aromatic stack, with F1068 displaced toward the space vacated by F508, while in F508del/R1070F-CFTR, which largely retained F508del defects, R1070F could not form hydrogen bonds and the interface was less flexible. Other ICL4 second-site mutations which partially rescued F508del-CFTR included F1068M and F1074M. Methionine side chains allow hydrophobic interactions without the steric rigidity of aromatic rings, possibly conferring flexibility to accommodate the absence of F508 and retain a dynamic interface. These studies highlight how both hydrophobic interactions and conformational flexibility might be important at the ICL4/NBD1 interface, suggesting possible structural underpinnings of F508del-induced dysfunction.

ProLint: a web-based framework for the automated data analysis and visualization of lipid-protein interactions.

The functional activity of membrane proteins is carried out in a complex lipid environment. Increasingly, it is becoming clear that lipids are an important player in regulating or generally modulating their activity. A routinely used method to gain insight into this interplay between lipids and proteins are Molecular Dynamics (MD) simulations, since they allow us to study interactions at atomic or near-atomic detail as a function of time. A major bottleneck, however, is analyzing and visualizing lipid-protein interactions, which, in practice, is a time-demanding task. Here, we present ProLint (www.prolint.ca), a webserver that completely automates analysis of MD generated files and visualization of lipid-protein interactions. Analysis is modular allowing users to select their preferred method, and visualization is entirely interactive through custom built applications that enable a detailed qualitative and quantitative exploration of lipid-protein interactions. ProLint also includes a database of published MD results that have been processed through the ProLint workflow and can be visualized by anyone regardless of their level of experience with MD. The automated analysis, feature-rich visualization, database integration, and open-source distribution with an easy to install process, will allow ProLint to become a routine workflow in lipid-protein interaction studies.

Curvature-based sorting of eight lipid types in asymmetric buckled plasma membrane models.

Curvature is a fundamental property of biological membranes and has essential roles in cellular function. Bending of membranes can be induced by their lipid and protein compositions, as well as peripheral proteins, such as those that make up the cytoskeleton. An important aspect of membrane function is the grouping of lipid species into microdomains, or rafts, which serve as platforms for specific biochemical processes. The fluid mosaic model of membranes has evolved to recognize the importance of curvature and leaflet asymmetry, and there are efforts toward evaluating their functional roles. This work investigates the effect of curvature on the sorting of lipids in buckled asymmetric bilayers containing eight lipid types, approximating an average mammalian plasma membrane, through coarse-grained (CG) molecular dynamics (MD) simulations with the Martini force field. The simulations reveal that 1) leaflet compositional asymmetry can induce curvature asymmetry, 2) lipids are sorted by curvature to different extents, and 3) curvature-based partitioning trends show moderate to strong correlations with lipid molecular volumes and head to tail bead ratios, respectively. The findings provide unique insights into the role of curvature in membrane organization, and the curvature-based sorting trends should be useful references for later investigations and potentially interpreting the functional roles of specific lipids.

Effects of cholesterol and PIP2 on interactions between glycophorin A and Band 3 in lipid bilayers.

In the erythrocyte membrane, the interactions between glycophorin A (GPA) and Band 3 are associated strongly with the biological function of the membrane and several blood disorders. In this work, using coarse-grained molecular-dynamics simulations, we systematically investigate the effects of cholesterol and phosphatidylinositol-4,5-bisphosphate (PIP2) on the interactions of GPA with Band 3 in the model erythrocyte membranes. We examine the dynamics of the interactions of GPA with Band 3 in different lipid bilayers on the microsecond time scale and calculate the binding free energy between GPA and Band 3. The results indicate that cholesterols thermodynamically favor the binding of GPA to Band 3 by increasing the thickness of the lipid bilayer and by producing an effective attraction between the proteins due to the depletion effect. Cholesterols also slow the kinetics of the binding of GPA to Band 3 by reducing the lateral mobility of the lipids and proteins and may influence the binding sites between the proteins. The anionic PIP2 lipids prefer binding to the surface of the proteins through electrostatic attraction between the PIP2 headgroup and the positively charged residues on the protein surface. Ions in the solvent facilitate PIP2 aggregation, which promotes the binding of GPA to Band 3.

The voltage-sensing domain of a hERG1 mutant is a cation-selective channel.

A cationic leak current known as an "omega current" may arise from mutations of the first charged residue in the S4 of the voltage sensor domains of sodium and potassium voltage-gated channels. The voltage-sensing domains (VSDs) in these mutated channels act as pores allowing nonspecific passage of cations, such as Li+, K+, Cs+, and guanidinium. Interestingly, no omega currents have been previously detected in the nonswapped voltage-gated potassium channels such as the human-ether-a-go-go-related (hERG1), hyperpolarization-activated cyclic nucleotide-gated, and ether-a-go-go channels. In this work, we discovered a novel omega current by mutating the first charged residue of the S4 of the hERG1, K525 to serine. To characterize this omega current, we used various probes, including the hERG1 pore domain blocker, dofetilide, to show that the omega current does not require cation flux via the canonical pore domain. In addition, the omega flux does not cross the conventional selectivity filter. We also show that the mutated channel (K525S hERG1) conducts guanidinium. These data are indicative of the formation of an omega current channel within the VSD. Using molecular dynamics simulations with replica-exchange umbrella sampling simulations of the wild-type hERG1 and the K525S hERG1, we explored the molecular underpinnings governing the cation flow in the VSD of the mutant. We also show that the wild-type hERG1 may form water crevices supported by the biophysical surface accessibility data. Overall, our multidisciplinary study demonstrates that the VSD of hERG1 may act as a cation-selective channel wherein a mutation of the first charged residue in the S4 generates an omega current. Our simulation uncovers the atomistic underpinning of this mechanism.

Computational Insights into the Role of Cholesterol in Inverted Hexagonal Phase Stabilization and Endosomal Drug Release.

Cholesterol is a major component of many lipid-based drug delivery systems, including cationic lipid nanoparticles. Despite its critical role in the drug release stage, the underlying molecular mechanism by which cholesterol assists in endosomal escape remains unclear. An efficient drug release from the endosome requires endosomal disruption. This disruption is believed to involve a lamellar-to-inverted hexagonal (Lα-HII) phase transition upon fusion of the lipid nanoparticle with the endosomal membrane. We used molecular dynamics simulations to study the structural properties of HII systems composed of an anionic lipid distearoyl phosphatidylserine (DSPS), an ionizable cationic lipid (KC2H), and cholesterol for several hydration levels and molar ratios. This system corresponds to the lipid mixtures in the hypothesized HII structure formed upon fusion and is of interest for the rational design of ionizable cationic lipids, including KC2, for an optimal drug release. Simulations suggest a geometry- and symmetry-driven lipid sorting and cholesterol-DSPS co-location around the water cores. Cholesterol preferentially co-locates with negatively charged saturated DSPS lipids at interstitial angles. The observed cholesterol-DSPS co-location results in an overall increase in the DSPS acyl chains' order parameters, which we propose to assist in stabilizing the HII phase by stretching the DSPS acyl chains for filling the voids formed by three adjacent lipid tubules. Furthermore, a systematic increase in the cholesterol concentration increased the lattice plane spacing and the water core radius but decreased the undulations along the lipid tubule axis. We propose that cholesterol and the degree of saturation/polyunsaturation of the lipid acyl chains, and not the lipid charge, are the main contributors in facilitating the Lα-HII phase transition and stabilizing/destabilizing the formed HII phase, whereas the positive charge of the ionizable cationic lipid promotes the LNP-endosomal membrane adhesion and assists in initiating the fusion process at the local contact area. We also propose that the effect of cholesterol on the HII structure and curvature is the main underlying reason for the well-documented HII stabilization and destabilization at low and high molar concentrations of cholesterol, respectively.

Supramolecular Organization of SARS-CoV and SARS-CoV-2 Virions Revealed by Coarse-Grained Models of Intact Virus Envelopes.

The coronavirus disease 19 (COVID-19) pandemic is causing a global health crisis and has already caused a devastating societal and economic burden. The pathogen, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has a high sequence and architecture identity with SARS-CoV, but far more people have been infected by SARS-CoV-2. Here, combining the structural data from cryo-electron microscopy and structure prediction, we constructed bottom-up Martini coarse-grained models of intact SARS-CoV and SARS-CoV-2 envelopes. Microsecond molecular dynamics simulations were performed, allowing us to explore their dynamics and supramolecular organization. Both SARS-CoV and SARS-CoV-2 envelopes present a spherical morphology, with structural proteins forming multiple string-like islands in the membrane and clusters between the heads of spike proteins. Critical differences between the SARS-CoV and SARS-CoV-2 envelopes are the interaction pattern between the spike proteins and the flexibility of the spike proteins. Our models provide structural and dynamic insights into the SARS virus envelopes and could be used for further investigation, such as drug design and membrane fusion and fission processes.

In vitro analyses of suspected arrhythmogenic thin filament variants as a cause of sudden cardiac death in infants.

Sudden unexpected death of an infant (SUDI) is a devastating occurrence for families. To investigate the genetic pathogenesis of SUDI, we sequenced >70 genes from 191 autopsy-negative SUDI victims. Ten infants sharing a previously unknown variant in troponin I (TnI) were identified. The mutation (TNNI1 R37C+/-) is in the fetal/neonatal paralog of TnI, a gene thought to be expressed in the heart up to the first 24 months of life. Using phylogenetic analysis and molecular dynamics simulations, it was determined that arginine at residue 37 in TNNI1 may play a critical functional role, suggesting that the variant may be pathogenic. We investigated the biophysical properties of the TNNI1 R37C mutation in human reconstituted thin filaments (RTFs) using fluorometry. RTFs reconstituted with the mutant R37C TnI exhibited reduced Ca2+-binding sensitivity due to an increased Ca2+ off-rate constant. Furthermore, we generated TNNI1 R37C+/- mutants in human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) using CRISPR-Cas9. In monolayers of hiPSC-CMs, we simultaneously monitored voltage and Ca2+ transients through optical mapping and compared them to their isogenic controls. We observed normal intrinsic beating patterns under control conditions in TNNI1 R37C+/- at stimulation frequencies of 55 beats/min (bpm), but these cells showed no restitution with increased stimulation frequency to 65 bpm and exhibited alternans at >75 bpm. The WT hiPSC-CMs did not exhibit any sign of arrhythmogenicity even at stimulation frequencies of 120 bpm. The approach used in this study provides critical physiological and mechanistic bases to investigate sarcomeric mutations in the pathogenesis of SUDI.

Refinement of a cryo-EM structure of hERG: Bridging structure and function.

The human-ether-a-go-go-related gene (hERG) encodes the voltage-gated potassium channel (KCNH2 or Kv11.1, commonly known as hERG). This channel plays a pivotal role in the stability of phase 3 repolarization of the cardiac action potential. Although a high-resolution cryo-EM structure is available for its depolarized (open) state, the structure surprisingly did not feature many functionally important interactions established by previous biochemical and electrophysiology experiments. Using molecular dynamics flexible fitting (MDFF), we refined the structure and recovered the missing functionally relevant salt bridges in hERG in its depolarized state. We also performed electrophysiology experiments to confirm the functional relevance of a novel salt bridge predicted by our refinement protocol. Our work shows how refinement of a high-resolution cryo-EM structure helps to bridge the existing gap between the structure and function in the voltage-sensing domain (VSD) of hERG.

Structural basis for antibacterial peptide self-immunity by the bacterial ABC transporter McjD.

Certain pathogenic bacteria produce and release toxic peptides to ensure either nutrient availability or evasion from the immune system. These peptides are also toxic to the producing bacteria that utilize dedicated ABC transporters to provide self-immunity. The ABC transporter McjD exports the antibacterial peptide MccJ25 in Escherichia coli Our previously determined McjD structure provided some mechanistic insights into antibacterial peptide efflux. In this study, we have determined its structure in a novel conformation, apo inward-occluded and a new nucleotide-bound state, high-energy outward-occluded intermediate state, with a defined ligand binding cavity. Predictive cysteine cross-linking in E. coli membranes and PELDOR measurements along the transport cycle indicate that McjD does not undergo major conformational changes as previously proposed for multi-drug ABC exporters. Combined with transport assays and molecular dynamics simulations, we propose a novel mechanism for toxic peptide ABC exporters that only requires the transient opening of the cavity for release of the peptide. We propose that shielding of the cavity ensures that the transporter is available to export the newly synthesized peptides, preventing toxic-level build-up.

Lipid distributions and transleaflet cholesterol migration near heterogeneous surfaces in asymmetric bilayers.

Specific and nonspecific protein-lipid interactions in cell membranes have important roles in an abundance of biological functions. We have used coarse-grained (CG) molecular dynamics (MD) simulations to assess lipid distributions and cholesterol flipping dynamics around surfaces in a model asymmetric plasma membrane containing one of six structurally distinct entities: aquaporin-1 (AQP1), the bacterial ß-barrel outer membrane proteins OmpF and OmpX, the KcsA potassium channel, the WALP23 peptide and a carbon nanotube (CNT). Our findings revealed varied lipid partitioning and cholesterol flipping times around the different solutes and putative cholesterol binding sites in AQP1 and KcsA. The results suggest that protein-lipid interactions can be highly variable, and that surface-dependent lipid profiles are effectively manifested in CG simulations with the Martini force field.

Emerging Diversity in Lipid-Protein Interactions.

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Computational Modeling of Realistic Cell Membranes.

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Mutagenic Analysis of the Putative ABCC6 Substrate-Binding Cavity Using a New Homology Model.

Inactivating mutations in ABCC6 underlie the rare hereditary mineralization disorder pseudoxanthoma elasticum. ABCC6 is an ATP-binding cassette (ABC) integral membrane protein that mediates the release of ATP from hepatocytes into the bloodstream. The released ATP is extracellularly converted into pyrophosphate, a key mineralization inhibitor. Although ABCC6 is firmly linked to cellular ATP release, the molecular details of ABCC6-mediated ATP release remain elusive. Most of the currently available data support the hypothesis that ABCC6 is an ATP-dependent ATP efflux pump, an un-precedented function for an ABC transporter. This hypothesis implies the presence of an ATP-binding site in the substrate-binding cavity of ABCC6. We performed an extensive mutagenesis study using a new homology model based on recently published structures of its close homolog, bovine Abcc1, to characterize the substrate-binding cavity of ABCC6. Leukotriene C4 (LTC4), is a high-affinity substrate of ABCC1. We mutagenized fourteen amino acid residues in the rat ortholog of ABCC6, rAbcc6, that corresponded to the residues in ABCC1 found in the LTC4 binding cavity. Our functional characterization revealed that most of the amino acids in rAbcc6 corresponding to those found in the LTC4 binding pocket in bovine Abcc1 are not critical for ATP efflux. We conclude that the putative ATP binding site in the substrate-binding cavity of ABCC6/rAbcc6 is distinct from the bovine Abcc1 LTC4-binding site.

Lipid-Protein Interactions Are a Unique Property and Defining Feature of G Protein-Coupled Receptors.

G protein-coupled receptors (GPCRs) are membrane-bound proteins that depend on their lipid environment to carry out their physiological function. Combined efforts from many theoretical and experimental studies on the lipid-protein interaction profile of several GPCRs hint at an intricate relationship of these receptors with their surrounding membrane environment, with several lipids emerging as particularly important. Using coarse-grained molecular dynamics simulations, we explore the lipid-protein interaction profiles of 28 different GPCRs, spanning different levels of classification and conformational states and totaling to 1 ms of simulation time. We find a close relationship with lipids for all GPCRs simulated, in particular, cholesterol and phosphatidylinositol phosphate (PIP) lipids, but the number, location, and estimated strength of these interactions is dependent on the specific GPCR as well as its conformational state. Although both cholesterol and PIP lipids bind specifically to GPCRs, they utilize distinct mechanisms. Interactions with PIP lipids are mediated by charge-charge interactions with intracellular loop residues and stabilized by one or both of the transmembrane helices linked by the loop. Interactions with cholesterol, on the other hand, are mediated by a hydrophobic environment, usually made up of residues from more than one helix, capable of accommodating its ring structure and stabilized by interactions with aromatic and charged/polar residues. Cholesterol binding to GPCRs occurs in a small number of sites, some of which (like the binding site on the extracellular side of transmembrane 6/7) are shared among many class A GPCRs. Combined with a thorough investigation of the local membrane structure, our results provide a detailed picture of GPCR-lipid interactions. Additionally, we provide an accompanying website to interactively explore the lipid-protein interaction profile of all GPCRs simulated to facilitate analysis and comparison of our data.