Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework.

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Exploring CRD mobility during RAS/RAF engagement at the membrane.

During the activation of mitogen-activated protein kinase (MAPK) signaling, the RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF bind to active RAS at the plasma membrane. The orientation of RAS at the membrane may be critical for formation of the RAS-RBDCRD complex and subsequent signaling. To explore how RAS membrane orientation relates to the protein dynamics within the RAS-RBDCRD complex, we perform multiscale coarse-grained and all-atom molecular dynamics (MD) simulations of KRAS4b bound to the RBD and CRD domains of RAF-1, both in solution and anchored to a model plasma membrane. Solution MD simulations describe dynamic KRAS4b-CRD conformations, suggesting that the CRD has sufficient flexibility in this environment to substantially change its binding interface with KRAS4b. In contrast, when the ternary complex is anchored to the membrane, the mobility of the CRD relative to KRAS4b is restricted, resulting in fewer distinct KRAS4b-CRD conformations. These simulations implicate membrane orientations of the ternary complex that are consistent with NMR measurements. While a crystal structure-like conformation is observed in both solution and membrane simulations, a particular intermolecular rearrangement of the ternary complex is observed only when it is anchored to the membrane. This configuration emerges when the CRD hydrophobic loops are inserted into the membrane and helices α3-5 of KRAS4b are solvent exposed. This membrane-specific configuration is stabilized by KRAS4b-CRD contacts that are not observed in the crystal structure. These results suggest modulatory interplay between the CRD and plasma membrane that correlate with RAS/RAF complex structure and dynamics, and potentially influence subsequent steps in the activation of MAPK signaling.

Full scale structural, mechanical and dynamical properties of HIV-1 liposomes.

Enveloped viruses are enclosed by a lipid membrane inside of which are all of the components necessary for the virus life cycle; viral proteins, the viral genome and metabolites. Viral envelopes are lipid bilayers that adopt morphologies ranging from spheres to tubes. The envelope is derived from the host cell during viral replication. Thus, the composition of the bilayer depends on the complex constitution of lipids from the host-cell's organelle(s) where assembly and/or budding of the viral particle occurs. Here, molecular dynamics (MD) simulations of authentic, asymmetric HIV-1 liposomes are used to derive a unique level of resolution of its full-scale structure, mechanics and dynamics. Analysis of the structural properties reveal the distribution of thicknesses of the bilayers over the entire liposome as well as its global fluctuations. Moreover, full-scale mechanical analyses are employed to derive the global bending rigidity of HIV-1 liposomes. Finally, dynamical properties of the lipid molecules reveal important relationships between their 3D diffusion, the location of lipid-rafts and the asymmetrical composition of the envelope. Overall, our simulations reveal complex relationships between the rich lipid composition of the HIV-1 liposome and its structural, mechanical and dynamical properties with critical consequences to different stages of HIV-1's life cycle.

Array programming with NumPy.

Array programming provides a powerful, compact and expressive syntax for accessing, manipulating and operating on data in vectors, matrices and higher-dimensional arrays. NumPy is the primary array programming library for the Python language. It has an essential role in research analysis pipelines in fields as diverse as physics, chemistry, astronomy, geoscience, biology, psychology, materials science, engineering, finance and economics. For example, in astronomy, NumPy was an important part of the software stack used in the discovery of gravitational waves1 and in the first imaging of a black hole2. Here we review how a few fundamental array concepts lead to a simple and powerful programming paradigm for organizing, exploring and analysing scientific data. NumPy is the foundation upon which the scientific Python ecosystem is constructed. It is so pervasive that several projects, targeting audiences with specialized needs, have developed their own NumPy-like interfaces and array objects. Owing to its central position in the ecosystem, NumPy increasingly acts as an interoperability layer between such array computation libraries and, together with its application programming interface (API), provides a flexible framework to support the next decade of scientific and industrial analysis.

Author Correction: SciPy 1.0: fundamental algorithms for scientific computing in Python.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SciPy 1.0: fundamental algorithms for scientific computing in Python.

SciPy is an open-source scientific computing library for the Python programming language. Since its initial release in 2001, SciPy has become a de facto standard for leveraging scientific algorithms in Python, with over 600 unique code contributors, thousands of dependent packages, over 100,000 dependent repositories and millions of downloads per year. In this work, we provide an overview of the capabilities and development practices of SciPy 1.0 and highlight some recent technical developments.

Scalable Analysis of Authentic Viral Envelopes on FRONTERA.

Enveloped viruses, such as SARS-CoV-2, infect cells via fusion of their envelope with the host membrane. By employing molecular simulations to characterize viral envelopes, researchers can gain insights into key determinants of infection. Here, the Frontera supercomputer is leveraged for large-scale modeling and analysis of authentic viral envelopes, whose lipid compositions are complex and realistic. Visual Molecular Dynamics (VMD) with support for MPI is employed, overcoming previous computational limitations and enabling investigation into virus biology at an unprecedented scale. The techniques applied here to an authentic HIV-1 envelope at two levels of spatial resolution (29 million particles and 280 million atoms) are broadly applicable to the study of other viruses. The authors are actively employing these techniques to develop and characterize an authentic SARS-CoV-2 envelope. A general framework for carrying out scalable analysis of simulation trajectories on Frontera is presented, expanding the utility of the machine in humanity's ongoing fight against infectious diseases.

How nanoscale protein interactions determine the mesoscale dynamic organisation of bacterial outer membrane proteins.

The spatiotemporal organisation of membranes is often characterised by the formation of large protein clusters. In Escherichia coli, outer membrane protein (OMP) clustering leads to OMP islands, the formation of which underpins OMP turnover and drives organisation across the cell envelope. Modelling how OMP islands form in order to understand their origin and outer membrane behaviour has been confounded by the inherent difficulties of simulating large numbers of OMPs over meaningful timescales. Here, we overcome these problems by training a mesoscale model incorporating thousands of OMPs on coarse-grained molecular dynamics simulations. We achieve simulations over timescales that allow direct comparison to experimental data of OMP behaviour. We show that specific interaction surfaces between OMPs are key to the formation of OMP clusters, that OMP clusters present a mesh of moving barriers that confine newly inserted proteins within islands, and that mesoscale simulations recapitulate the restricted diffusion characteristics of OMPs.

Protein crowding and lipid complexity influence the nanoscale dynamic organization of ion channels in cell membranes.

Cell membranes are crowded and complex environments. To investigate the effect of protein-lipid interactions on dynamic organization in mammalian cell membranes, we have performed coarse-grained molecular dynamics simulations containing >100 copies of an inwardly rectifying potassium (Kir) channel which forms specific interactions with the regulatory lipid phosphatidylinositol 4,5-bisphosphate (PIP2). The tendency of protein molecules to cluster has the effect of organizing the membrane into dynamic compartments. At the same time, the diversity of lipids present has a marked effect on the clustering behavior of ion channels. Sub-diffusion of proteins and lipids is observed. Protein crowding alters the sub-diffusive behavior of proteins and lipids such as PIP2 which interact tightly with Kir channels. Protein crowding also affects bilayer properties, such as membrane undulations and bending rigidity, in a PIP2-dependent manner. This interplay between the diffusion and the dynamic organization of Kir channels may have important implications for channel function.

Membrane Compartmentalization Reducing the Mobility of Lipids and Proteins within a Model Plasma Membrane.

The cytoskeleton underlying cell membranes may influence the dynamic organization of proteins and lipids within the bilayer by immobilizing certain transmembrane (TM) proteins and forming corrals within the membrane. Here, we present coarse-grained resolution simulations of a biologically realistic membrane model of asymmetrically organized lipids and TM proteins. We determine the effects of a model of cytoskeletal immobilization of selected membrane proteins using long time scale coarse-grained molecular dynamics simulations. By introducing compartments with varying degrees of restraints within the membrane models, we are able to reveal how compartmentalization caused by cytoskeletal immobilization leads to reduced and anomalous diffusional mobility of both proteins and lipids. This in turn results in a reduced rate of protein dimerization within the membrane and of hopping of membrane proteins between compartments. These simulations provide a molecular realization of hierarchical models often invoked to explain single-molecule imaging studies of membrane proteins.

The Role of the Membrane in the Structure and Biophysical Robustness of the Dengue Virion Envelope.

The dengue virion is surrounded by an envelope of membrane proteins surrounding a lipid bilayer. We have combined the cryoelectron microscopy structures of the membrane proteins (PDB: 3J27) with a lipid bilayer whose composition is based on lipidomics data for insect cell membranes, to obtain a near-atomic resolution computational model of the envelope of the dengue virion. A coarse-grained molecular dynamics simulation on the microsecond timescale enables analysis of key biophysical properties of the dengue outer envelope. Properties analyzed include area per lipid values (for a spherical virion with a mixed lipid composition), bilayer thickness, and lipid diffusion coefficients. Despite the absence of cholesterol from the lipid bilayer, the virion exhibits biophysical robustness (slow lipid diffusion alongside stable bilayer thickness, virion diameter, and shape) that matches the cholesterol-rich membrane of influenza A, with similarly anomalous diffusion of lipids. Biophysical robustness of the envelope may confer resilience to environmental perturbations.

Computational virology: From the inside out.

Viruses typically pack their genetic material within a protein capsid. Enveloped viruses also have an outer membrane made up of a lipid bilayer and membrane-spanning glycoproteins. X-ray diffraction and cryoelectron microscopy provide high resolution static views of viral structure. Molecular dynamics (MD) simulations may be used to provide dynamic insights into the structures of viruses and their components. There have been a number of simulations of viral capsids and (in some cases) of the inner core of RNA or DNA packaged within them. These simulations have generally focussed on the structural integrity and stability of the capsid and/or on the influence of the nucleic acid core on capsid stability. More recently there have been a number of simulation studies of enveloped viruses, including HIV-1, influenza A, and dengue virus. These have addressed the dynamic behaviour of the capsid, the matrix, and/or of the outer envelope. Analysis of the dynamics of the lipid bilayer components of the envelopes of influenza A and of dengue virus reveals a degree of biophysical robustness, which may contribute to the stability of virus particles in different environments. Significant computational challenges need to be addressed to aid simulation of complex viruses and their membranes, including the need to integrate structural data from a range of sources to enable us to move towards simulations of intact virions. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Nothing to sneeze at: a dynamic and integrative computational model of an influenza A virion.

The influenza virus is surrounded by an envelope composed of a lipid bilayer and integral membrane proteins. Understanding the structural dynamics of the membrane envelope provides biophysical insights into aspects of viral function, such as the wide-ranging survival times of the virion in different environments. We have combined experimental data from X-ray crystallography, nuclear magnetic resonance spectroscopy, cryo-electron microscopy, and lipidomics to build a model of the intact influenza A virion. This is the basis of microsecond-scale coarse-grained molecular dynamics simulations of the virion, providing simulations at different temperatures and with varying lipid compositions. The presence of the Forssman glycolipid alters a number of biophysical properties of the virion, resulting in reduced mobility of bilayer lipid and protein species. Reduced mobility in the virion membrane may confer physical robustness to changes in environmental conditions. Our simulations indicate that viral spike proteins do not aggregate and thus are competent for multivalent immunoglobulin G interactions.

Methodologies for the analysis of instantaneous lipid diffusion in MD simulations of large membrane systems.

Interactions between lipids and membrane proteins play a key role in determining the nanoscale dynamic and structural properties of biological membranes. Molecular dynamics (MD) simulations provide a valuable tool for studying membrane models, complementing experimental approaches. It is now possible to simulate large membrane systems, such as simplified models of bacterial and viral envelope membranes. Consequently, there is a pressing need to develop tools to visualize and quantify the dynamics of these immense systems, which typically comprise millions of particles. To tackle this issue, we have developed visual and quantitative analyses of molecular positions and their velocity field using path line, vector field and streamline techniques. This allows us to highlight large, transient flow-like movements of lipids and to better understand crowding within the lipid bilayer. The current study focuses on visualization and analysis of lipid dynamics. However, the methods are flexible and can be readily applied to e.g. proteins and nanoparticles within large complex membranes. The protocols developed here are readily accessible both as a plugin for the molecular visualization program VMD and as a module for the MDAnalysis library.

Primary and secondary dimer interfaces of the fibroblast growth factor receptor 3 transmembrane domain: characterization via multiscale molecular dynamics simulations.

Receptor tyrosine kinases are single-pass membrane proteins that form dimers within the membrane. The interactions of their transmembrane domains (TMDs) play a key role in dimerization and signaling. Fibroblast growth factor receptor 3 (FGFR3) is of interest as a G380R mutation in its TMD is the underlying cause of ~99% of the cases of achondroplasia, the most common form of human dwarfism. The structural consequences of this mutation remain uncertain: the mutation shifts the position of the TMD relative to the lipid bilayer but does not alter the association free energy. We have combined coarse-grained and all-atom molecular dynamics simulations to study the dimerization of wild-type, heterodimer, and mutant FGFR3 TMDs. The simulations reveal that the helices pack together in the dimer to form a flexible interface. The primary packing mode is mediated by a Gx3G motif. There is also a secondary dimer interface that is more highly populated in heterodimer and mutant configurations that may feature in the molecular mechanism of pathology. Both coarse-grained and atomistic simulations reveal a significant shift of the G380R mutant dimer TMD relative to the bilayer to allow interactions of the arginine side chain with lipid headgroup phosphates.

Structural features of the apelin receptor N-terminal tail and first transmembrane segment implicated in ligand binding and receptor trafficking.

G-protein coupled receptors (GPCRs) comprise a large family of membrane proteins with rich functional diversity. Signaling through the apelin receptor (AR or APJ) influences the cardiovascular system, central nervous system and glucose regulation. Pathophysiological involvement of apelin has been shown in atherosclerosis, cancer, human immunodeficiency virus-1 (HIV-1) infection and obesity. Here, we present the high-resolution nuclear magnetic resonance (NMR) spectroscopy-based structure of the N-terminus and first transmembrane (TM) segment of AR (residues 1-55, AR55) in dodecylphosphocholine micelles. AR55 consists of two disrupted helices, spanning residues D14-K25 and A29-R55(1.59). Molecular dynamics (MD) simulations of AR built from a hybrid of experimental NMR and homology model-based restraints allowed validation of the AR55 structure in the context of the full-length receptor in a hydrated bilayer. AR55 structural features were functionally probed using mutagenesis in full-length AR through monitoring of apelin-induced extracellular signal-regulated kinase (ERK) phosphorylation in transiently transfected human embryonic kidney (HEK) 293A cells. Residues E20 and D23 form an extracellular anionic face and interact with lipid headgroups during MD simulations in the absence of ligand, producing an ideal binding site for a cationic apelin ligand proximal to the membrane-water interface, lending credence to membrane-catalyzed apelin-AR binding. In the TM region of AR55, N46(1.50) is central to a disruption in helical character. G42(1.46), G45(1.49) and N46(1.50), which are all involved in the TM helical disruption, are essential for proper trafficking of AR. In summary, we introduce a new correlative NMR spectroscopy and computational biochemistry methodology and demonstrate its utility in providing some of the first high-resolution structural information for a peptide-activated GPCR TM domain.

Multifaceted substrate capture scheme of a rhomboid protease.

Rhomboid proteases are integral membrane serine proteases that catalyze peptide bond hydrolysis in biological membranes. Little is currently known about the interaction of enzyme and substrate. Coarse-grained molecular dynamics simulations in hydrated lipid bilayers are employed herein to study the interaction of the E. coli rhomboid protease GlpG (ecGlpG) with the transmembrane domain (TMD) of the substrate Spitz. Spitz does not associate with ecGlpG exclusively at the putative substrate gate near TMD 5. Instead, there are six prominent and stable interaction sites, including one between TMDs 1 and 3, with the closest enzyme-substrate proximity occurring at the ends of helical TMDs or in loops. Bilayer thinning is observed proximal to ecGlpG, but there is no evidence of additional thinning of the bilayer upon interaction with substrate. We suggest that the initial interaction between enzyme and substrate, or substrate capture event, is not limited to a single site on the enzyme, and may be driven by juxtamembrane electrostatic interactions. The findings are of additional interest because catalytically inactive rhomboids (iRhoms) are now known to interact with the substrates of their catalytically active counterparts and to antagonize the enzyme-driven pathways.

Interpretation of biomolecular NMR spin relaxation parameters.

Biomolecular nuclear magnetic resonance (NMR) spin relaxation experiments provide exquisite information on the picosecond to nanosecond timescale motions of bond vectors. Spin-lattice (T1) and spin-spin (T2) relaxation times and the steady-state nuclear Overhauser effect (NOE) are the first set of parameters extracted from typical 15N or 13C NMR relaxation experiments. Therefore, verifying that T1, T2, and NOE are consistent with theoretical predictions is an important step before carrying out the more detailed model-free and reduced spectral density mapping analyses commonly employed. In this mini-review, we discuss the essential motional parameters used to describe biomolecular dynamics in the context of a variety of examples of folded and intrinsically disordered proteins and peptides in aqueous and membrane mimetic environments. Estimates of these parameters can be used as input for an online interface, introduced herein, allowing plotting of trends of T1, T2, and NOE with magnetic field strength. The plots may serve as a first-check to the spectroscopist preparing to embark on a detailed NMR relaxation analysis.

Correlating structure, dynamics, and function in transmembrane segment VII of the Na+/H+ exchanger isoform 1.

We place (15)N nuclear magnetic resonance relaxation analysis and functional mutagenesis studies in the context of our previous structural and mutagenesis work to correlate structure, dynamics and function for the seventh transmembrane segment of the human Na(+)/H(+) exchanger isoform 1. Although G261-S263 was previously identified as an interruption point in the helical structure of this isolated transmembrane peptide in dodecylphosphocholine micelles, and rapid conformational exchange was implicated in the NOE measurements, the six (15)N labelled residues examined in this study all have similar dynamics on the ps-ns time scale. A mathematical model incorporating chemical exchange is the best fit for residues G261, L264, and A268. This implies that a segment of residues from G261 to A268 samples different conformations on the mus-ms time scale. Chemical exchange on an intermediate time scale is consistent with an alternating-access cycle where E262 is bent away from the cytosol during proton translocation by the exchanger. The functional importance of chemical exchange at G261-A268 is corroborated by the abrogated activity of the full-length exchanger with the bulky and restricting Ile substitutions F260I, G261I, E262I, S263I, and A268I.

Structural insight into G-protein coupled receptor binding by apelin.

Apelin peptides are the cognate ligands for the G-protein coupled receptor APJ, with functions in the cardiovascular and central nervous systems, in glucose metabolism and as a human immunodeficiency virus (HIV-1) coreceptor. Apelin is found in 13-36 residue forms in vivo. The structures of five isoforms of apelin at physiological versus low (5-6 degrees C) temperature are compared here using circular dichroism (CD) and nuclear magnetic resonance (NMR) spectroscopy, demonstrating increased structure at low temperature. Far-ultraviolet (UV) CD spectra are predominantly random coil for apelin isoforms, but are convoluted by unusual bands from the C-terminal phenylalanine side chain. These bands, assigned using F13A-apelin-13, are accentuated at 5 degrees C and imply conformational restriction. At 35 degrees C, the R6-L9 region of apelin-17 is well structured, consistent with previous mutagenesis results showing necessity of this segment for apelin-APJ binding and activation. At 5 degrees C, R6-L9 retains its structuring while the functionally critical C-terminal G13-F17 region also becomes highly structured. Type IV beta-turns and some polyproline-II structure alongside F17 side chain motional reduction correlate well with CD spectral properties. Cis-trans peptide bond isomerization at P14 and P16 produces two sequentially assignable conformers (both trans:both cis approximately 4:1) alongside less populated conformers. Chemical shift assignment of apelin-12, -13 and pyroglutamate-apelin-13 implies highly similar structuring and the same isomerization at the C-terminus. Based on the apelin-17 structure, a two-step binding and activation mechanism is hypothesized.