Membrane Curvature Sensing by Amphipathic Helices: Insights from Implicit Membrane Modeling

Sensing and generation of lipid membrane curvature, mediated by the binding of specific proteins onto the membrane surface, play crucial roles in cell biology. A number of mechanisms have been proposed, but the molecular understanding of these processes is incomplete. All-atom molecular dynamics simulations have offered valuable insights but are extremely demanding computationally. Implicit membrane simulations could provide a viable alternative, but current models apply only to planar membranes. In this work, the implicit membrane model 1 is extended to spherical and tubular membranes. The geometric change from planar to curved shapes is straightforward but insufficient for capturing the full curvature effect, which includes changes in lipid packing. Here, these packing effects are taken into account via the lateral pressure profile. The extended implicit membrane model 1 is tested on the wild-types and mutants of the antimicrobial peptide magainin, the ALPS motif of arfgap1, α-synuclein, and an ENTH domain. In these systems, the model is in qualitative agreement with experiments. We confirm that favorable electrostatic interactions tend to weaken curvature sensitivity in the presence of strong hydrophobic interactions but may actually have a positive effect when those are weak. We also find that binding to vesicles is more favorable than binding to tubes of the same diameter and that the long helix of α-synuclein tends to orient along the axis of tubes, whereas shorter helices tend to orient perpendicular to it. Adoption of a specific orientation could provide a mechanism for coupling protein oligomerization to tubule formation.     

The Structural Basis of IKs Ion-Channel Activation: Mechanistic Insights from Molecular Simulations

Relating ion channel (iCh) structural dynamics to physiological function remains a challenge. Current experimental and computational techniques have limited ability to explore this relationship in atomistic detail over physiological timescales. A framework associating iCh structure to function is necessary for elucidating normal and disease mechanisms. We formulated a modeling schema that overcomes the limitations of current methods through applications of artificial intelligence machine learning. Using this approach, we studied molecular processes that underlie human IKs voltage-mediated gating. IKs malfunction underlies many debilitating and life-threatening diseases. Molecular components of IKs that underlie its electrophysiological function include KCNQ1 (a pore-forming tetramer) and KCNE1 (an auxiliary subunit). Simulations, using the IKs structure-function model, reproduced experimentally recorded saturation of gating-charge displacement at positive membrane voltages, two-step voltage sensor (VS) movement shown by fluorescence, iCh gating statistics, and current-voltage relationship. Mechanistic insights include the following: 1) pore energy profile determines iCh subconductance; 2) the entire protein structure, not limited to the pore, contributes to pore energy and channel subconductance; 3) interactions with KCNE1 result in two distinct VS movements, causing gating-charge saturation at positive membrane voltages and current activation delay; and 4) flexible coupling between VS and pore permits pore opening at lower VS positions, resulting in sequential gating. The new modeling approach is applicable to atomistic scale studies of other proteins on timescales of physiological function.     

Real-Time Insights into Biological Events: In-Cell Processes and Protein-Ligand Interactions

FlowNMR has the aim of continuously monitoring processes that occur in conditions that are not compatible with being carried out within a closed tube. However, it is sample intensive and not suitable for samples, such as proteins or living cells, that are often available in limited volumes and possibly low concentrations. We here propose a dialysis-based modification of a commercial flowNMR setup that allows for recycling the medium while confining the sample (proteins and cells) within the active volume of the tube. This approach is demonstrated in the specific cases of in-cell NMR and protein-based ligand studies.     

Integrating NMR,SAXS, and Atomistic Simulations: Structure and Dynamics of a Two-Domain Protein

Multidomain proteins with two or more independently folded functional domains are prevalent in nature. Whereas most multidomain proteins are linked linearly in sequence, roughly one-tenth possess domain insertions where a guest domain is implanted into a loop of a host domain, such that the two domains are connected by a pair of interdomain linkers. Here, we characterized the influence of the interdomain linkers on the structure and dynamics of a domain-insertion protein in which the guest LysM domain is inserted into a central loop of the host CVNH domain. Expanding upon our previous crystallographic and NMR studies, we applied SAXS in combination with NMR paramagnetic relaxation enhancement to construct a structural model of the overall two-domain system. Although the two domains have no fixed relative orientation, certain orientations were found to be preferred over others. We also assessed the accuracies of molecular mechanics force fields in modeling the structure and dynamics of tethered multidomain proteins by integrating our experimental results with microsecond-scale atomistic molecular dynamics simulations. In particular, our evaluation of two different combinations of the latest force fields and water models revealed that both combinations accurately reproduce certain structural and dynamical properties, but are inaccurate for others. Overall, our study illustrates the value of integrating experimental NMR and SAXS studies with long timescale atomistic simulations for characterizing structural ensembles of flexibly linked multidomain systems.     

Molecular Insights into Variable Electron Transfer in Amphibian Cryptochrome

Cryptochrome proteins are activated by the absorption of blue light, leading to the formation of radical pairs through electron transfer in the active site. Recent experimental studies have shown that once some of the amino acid residues in the active site of Xenopus laevis cryptochrome DASH are mutated, radical-pair formation is still observed. In this study, we computationally investigate electron-transfer pathways in the X. laevis cryptochrome DASH by extensively equilibrating a previously established homology model using molecular dynamics simulations and then mutating key amino acids involved in the electron transfer. The electron-transfer pathways are then probed by using tight-binding density-functional theory. We report the alternative electron-transfer pathways resolved at the molecular level and, through comparison of amino acid sequences for cryptochromes from different species, we demonstrate that one of these alternative electron-transfer pathways could be general for all cryptochrome DASH proteins.     

Molecular Insights into the Mechanisms of SUN1 Oligomerization in the Nuclear Envelope

The LINC complex is found in a wide variety of organisms and is formed by the transluminal interaction between outer- and inner-nuclear-membrane KASH and SUN proteins, respectively. Most extensively studied are SUN1 and SUN2 proteins, which are widely expressed in mammals. Although SUN1 and SUN2 play functionally redundant roles in several cellular processes, more recent studies have revealed diverse and distinct functions for SUN1. While several recent in vitro structural studies have revealed the molecular details of various fragments of SUN2, no such structural information is available for SUN1. Herein, we conduct a systematic analysis of the molecular relationships between SUN1 and SUN2, highlighting key similarities and differences that could lead to clues into their distinct functions. We use a wide range of computational tools, including multiple sequence alignments, homology modeling, molecular docking, and molecular dynamic simulations, to predict structural differences between SUN1 and SUN2, with the goal of understanding the molecular mechanisms underlying SUN1 oligomerization in the nuclear envelope. Our simulations suggest that the structural model of SUN1 is stable in a trimeric state and that SUN1 trimers can associate through their SUN domains to form lateral complexes. We also ask whether SUN1 could adopt an inactive monomeric conformation as seen in SUN2. Our results imply that the KASH binding domain of SUN1 is also inhibited in monomeric SUN1 but through weaker interactions than in monomeric SUN2.     

Molecular Insights on the Release of Avibactam from the Acyl-Enzyme Complex

Avibactam is a non-β-lactam β-lactamase inhibitor for treating complicated urinary tract and respiratory infections caused by multidrug-resistant bacterial pathogens, a serious public health threat. Despite its importance, the release mechanism of avibactam from the enzyme-inhibitor complex has been scarcely studied from first principles, considering the total protein environment. This information at the molecular level is essential for the rational design of new antibiotics and inhibitors. In this article, we addressed the release of avibactam from the complex CTX-M-15 by means of molecular dynamics simulations and quantum mechanics/molecular mechanics calculations. This study provides molecular information not available earlier, including exploration of the potential energy surfaces, characterization of the observed intermediate, and their critical points, as well. Our results show that unlike that observed in the acylation reaction, the residues Glu166 and Lys73 would be in their neutral forms. Release of avibactam follows a stepwise mechanism in which the first stage corresponds to the formation of a tetrahedral intermediate, whereas the second stage corresponds to the cleavage of the Ser70-C7 bond, mediated by Lys73, either directly or through Ser130.     

Dynamics of Membrane-Bound G12V-KRAS from Simulations and Single-Molecule FRET in Native Nanodiscs

Recent studies have shown that the small GTPase KRAS adopts multiple orientations with respect to the plane of anionic model membranes, whereby either the three C-terminal helices or the three N-terminal β-strands of the catalytic domain face the membrane. This has functional implications because, in the latter, the membrane occludes the effector-interacting surface. However, it remained unclear how membrane reorientation occurs and, critically, whether it occurs in the cell in which KRAS operates as a molecular switch in signaling pathways. Herein, using data from a 20 μs-long atomistic molecular dynamics simulation of the oncogenic G12V-KRAS mutant in a phosphatidylcholine/phosphatidylserine bilayer, we first show that internal conformational fluctuations of flexible regions in KRAS result in three distinct membrane orientations. We then show, using single-molecule fluorescence resonance energy transfer measurements in native lipid nanodiscs derived from baby hamster kidney cells, that G12V-KRAS samples three conformational states that correspond to the predicted orientations. The combined results suggest that relatively small energy barriers separate orientation states and that signaling-competent conformations dominate the overall population.     

Rotamer Dynamics: Analysis of Rotamers in Molecular Dynamics Simulations of Proteins

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.     

Partitioning and Enhanced Self-Assembly of Actin in Polypeptide Coacervates

Biomolecules exist and function in cellular microenvironments that control their spatial organization, local concentration, and biochemical reactivity. Due to the complexity of native cytoplasm, the development of artificial bioreactors and cellular mimics to compartmentalize, concentrate, and control the local physico-chemical properties is of great interest. Here, we employ self-assembling polypeptide coacervates to explore the partitioning of the ubiquitous cytoskeletal protein actin into liquid polymer-rich droplets. We find that actin spontaneously partitions into coacervate droplets and is enriched by up to ～30-fold. Actin polymerizes into micrometer-long filaments and, in contrast to the globular protein BSA, these filaments localize predominately to the droplet periphery. We observe up to a 50-fold enhancement in the actin filament assembly rate inside coacervate droplets, consistent with the enrichment of actin within the coacervate phase. Together these results suggest that coacervates can serve as a versatile platform in which to localize and enrich biomolecules to study their reactivity in physiological environments.