Structure, energetics, and dynamics of lipid-protein interactions: A molecular dynamics study of the gramicidin A channel in a DMPC bilayer

Cancer and Autism: How PTEN Mutations Degrade Function at the Membrane and Isoform Expression in the Human Brain

Mutations causing loss of PTEN lipid phosphatase activity can promote cancer, benign tumors (PHTS), and neurodevelopmental disorders (NDDs). Exactly how they preferentially trigger distinct phenotypic outcomes has been puzzling. Here, we demonstrate that PTEN mutations differentially allosterically bias P loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations of the functional wild-type state, while sampled conformations for the strong, cancer-related driver mutation hotspots favor catalysis-primed conformations, suggesting that NDD mutations are likely to be weaker, and our large-scale simulations show why. Prenatal PTEN isoform expression data suggest exons 5 and 7, which harbor NDD mutations, as cancer-risk carriers. Since cancer requires more than a single mutation, our conformational and genomic analysis helps discover how same protein mutations can foster different clinical manifestations, articulates a role for co-occurring background latent driver mutations, and uncovers relationships of splicing isoform expression to life expectancy.

Hyperoside inhibits pancreatic lipase activity in vitro and reduces fat accumulation in vivo

Hyperoside, the main component of many anti-obesity plants, might exhibit a lipase inhibition effect to reduce fat accumulation. The anti-obesity effect of hyperoside was investigated by studying its inhibitory effect and mechanism on pancreatic lipase in vitro and evaluating its ability to reduce lipid accumulation in vivo. Hyperoside is a mixed-type inhibitor of lipase with an IC₅₀ of 0.67 ± 0.02 mmol L^-in vitro. Hyperoside changed the secondary conformation of lipase, increased the α-helix content, and changed the microenvironment of lipase through static quenching. The interaction between hyperoside and lipase results from a strong binding spontaneous exothermic reaction, mainly through hydrogen bonding, van der Waals force and electrostatic force. Hyperoside protected hepatic lipid accumulation and adipose tissue hypertrophy and reduced the expression of inflammatory factors in high-fat diet-induced rats. Moreover, hyperoside had a good inhibitory effect on lipase activity in serum and increased fecal fat excretion, thereby reducing lipid absorption in vivo. This study provides theoretical support for the research and development of hyperoside in fat-reducing functional foods.

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.

Mutations in G6PC2 gene with increased risk for development of type 2 diabetes: Understanding via computational approach

An increase in the fast blood glucose (FBG) levels has been linked to an increased risk of developing a chronic condition, type 2 diabetes (T2D). The mutation in the G6PC2 gene was identified to have a lead role in the modulation of FBG levels. The abnormal regulation of this enzyme influences glucose-stimulated insulin secretion (GSIS), which controls the insulin levels corresponding to the system's glucose level. This study focuses on the mutations at the G6PC2 gene, which cause the variation from normal expression levels and increase the risk of T2D. We examined the non-synonymous single nucleotide polymorphisms (nsSNPs) present in the G6PC2 and subjected them to pathogenicity, stability, residue conservation, and membrane simulation. The individual representation of surrounding amino acids in the mutant (I63T) model showed the loss of hydrophobic interactions compared to the native G6PC2. In addition, the trajectory results from the membrane simulation exhibited reduced stability, and the least compactness was identified for the I63T mutant model. Our study shed light on the structural and conformational changes at the transmembrane region due to the I63T mutation in G6PC2. Additionally, the Gibbs free energy landscape analysis against the two principal components showed structural differences and decreased the conformational stability of the I63T mutant model compared to the native. Like those presented in this study, dynamical simulations may indeed be crucial to comprehending the structural insights of G6PC2 mutations in cardiovascular-associated mortality and T2D.

Structural determinants of peptide-dependent TAP1-TAP2 transit passage targeted by viral proteins and altered by cancer-associated mutations

The TAP1-TAP2 complex transports antigenic peptide substrates into the endoplasmic reticulum (ER). In ER, the peptides are further processed and loaded on the major histocompatibility class (MHC) I molecules by the peptide loading complex (PLC). The TAP transporters are linked with the PLC; a target for cancers and viral immune evasion. But the mechanisms whereby the cancer-derived mutations in TAP1-TAP2 or viral factors targeting the PLC, interfere peptide transport are only emerging. This study describes that transit of peptides through TAP can take place via two different channels (4 or 8 helices) depending on peptide length and sequence. Molecular dynamics and binding affinity predictions of peptide-transporters demonstrated that smaller peptides (8-10 mers; e.g. AAGIGILTV, SIINFEKL) can transport quickly through the transport tunnel compared to longer peptides (15-mer; e.g. ENPVVHFFKNIVTPR). In line with a regulated and selective peptide transport by TAPs, the immunopeptidome upon IFN-γ treatment in melanoma cells induced the shorter length (9-mer) peptide presentation over MHC-I that exhibit a relatively weak binding affinity with TAP. A conserved distance between N and C terminus residues of the studied peptides in the transport tunnel were reported. Furthermore, by adversely interacting with the TAP transport passage or affecting TAPNBD domains tilt movement, the viral proteins and cancer-derived mutations in TAP1-TAP2 may induce allosteric effects in TAP that block conformation of the tunnel (closed towards ER lumen). Interestingly, some cancer-associated mutations (e.g. TAP1R372Q and TAP2R373H) can specifically interfere with selective transport channels (i.e. for longer-peptides). These results provide a model for how viruses and cancer-associated mutations targeting TAP interfaces can affect MHC-I antigen presentation, and how the IFN-γ pathway alters MHC-I antigen presentation via the kinetics of peptide transport.

The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane

Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP₃ produced by PI3Ks. Excess PIP₃ promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP₂ and PIP₃. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP₃. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP₃ allosterically unfolding the N-terminal PIP₂ binding domain, positioning it favorably for the polybasic motif interaction with PIP₂. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer.

•

PTEN localizes on membrane microdomains enriched in phosphoinositides, as PI3K does

•

Full PTEN activation requires both signaling lipids, PIP₂ and PIP₃

•

Strong salt bridge interactions sustain stable PTEN membrane localization

•

Substrate-induced P loop conformational change implicates PTEN catalytic activity

Modifications on the Tetrahydroquinoline Scaffold Targeting a Phenylalanine Cluster on GPER as Antiproliferative Compounds against Renal, Liver and Pancreatic Cancer Cells

The implementation of chemo- and bioinformatics tools is a crucial step in the design of structure-based drugs, enabling the identification of more specific and effective molecules against cancer without side effects. In this study, three new compounds were designed and synthesized with suitable absorption, distribution, metabolism, excretion and toxicity (ADME-tox) properties and high affinity for the G protein-coupled estrogen receptor (GPER) binding site by in silico methods, which correlated with the growth inhibitory activity tested in a cluster of cancer cell lines. Docking and molecular dynamics (MD) simulations accompanied by a molecular mechanics/generalized Born surface area (MMGBSA) approach yielded the binding modes and energetic features of the proposed compounds on GPER. These in silico studies showed that the compounds reached the GPER binding site, establishing interactions with a phenylalanine cluster (F206, F208 and F278) required for GPER molecular recognition of its agonist and antagonist ligands. Finally, a 3-(4,5-dimethylthiazol-2-yl)2,5-diphenyltetrazolium bromide (MTT) assay showed growth inhibitory activity of compounds 4, 5 and 7 in three different cancer cell lines-MIA Paca-2, RCC4-VA and Hep G2-at micromolar concentrations. These new molecules with specific chemical modifications of the GPER pharmacophore open up the possibility of generating new compounds capable of reaching the GPER binding site with potential growth inhibitory activities against nonconventional GPER cell models.

In Silico Prediction of the Binding, Folding, Insertion, and Overall Stability of Membrane-Active Peptides

Membrane-active peptides (MAPs) are short-length peptides used for potential biomedical applications in diagnostic imaging of tissues, targeted drug delivery, gene delivery, and antimicrobials and antibiotics. The broad appeal of MAPs is that they are infinitely variable, relatively low cost, and biocompatible. However, experimentally characterizing the specific properties of a MAP or its many variants is a low-resolution and potentially time-consuming endeavor; molecular dynamics (MD) simulations have emerged as an invaluable tool in identifying the biophysical interactions that are fundamental to the function of MAPs. In this chapter, a step-by-step approach to discreetly model the binding, folding, and insertion of a membrane-active peptide to a model lipid bilayer using MD simulations is described. Detailed discussion is devoted to the critical aspects of running these types of simulations: prior knowledge of the system, understanding the strengths and weaknesses of molecular mechanics force fields, proper construction and equilibration of the system, realistically estimating both experimental and computational timescales, and leveraging analysis to make direct comparisons to experimental results as often as possible.

Assessing the extent of the structural and dynamic modulation of membrane lipids due to pore forming toxins: insights from molecular dynamics simulations

Infections caused by many virulent bacterial strains are triggered by the release of pore forming toxins (PFTs), which form oligomeric transmembrane pore complexes on the target plasma membrane. The spatial extent of the perturbation to the surrounding lipids during pore formation is relatively unexplored. Using all-atom molecular dynamics simulations, we investigate the changes in the structure and dynamics of lipids in a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid bilayer in the presence of contrasting PFTs. Cytolysin A (ClyA), an α toxin with its inserted wedge shaped bundle of inserted α helices, induces significant asymmetry across the membrane leaflets in comparison with α hemolysin (AHL), a β toxin. Despite the differences in hydrophobic mismatch and uniquely different topologies of the two oligomers, perturbations to lipid order as reflected in the tilt angle and order parameters and membrane thinning are short ranged, lying within ∼2.5 nm from the periphery of either pore complex, and commensurate with distances typically associated with van der Waals forces. In contrast, the spatial extent of perturbations to the lipid dynamics extends outward to at least 4 nm for both proteins, and the continuous survival probabilities reveal the presence of a tightly bound shell of lipids in this region. Displacement probability distributions show long tails and the distinctly non-Gaussian features reflect the induced dynamic heterogeneity. A detailed profiling of the protein-lipid contacts with tyrosine, tryptophan, lysine and arginine residues shows increased non-polar contacts in the cytoplasmic leaflet for both PFTs, with a higher number of atomic contacts in the case of AHL in the extracellular leaflet due to the mushroom-like topology of the pore complex. The short ranged nature of the perturbations observed in this simple one component membrane suggests inherent plasticity of membrane lipids enabling the recovery of the structure and membrane fluidity even in the presence of these large oligomeric transmembrane protein assemblies. This observation has implications in membrane repair processes such as budding or vesicle fusion events used to mitigate PFT virulence, where the underlying lipid dynamics and fluidity in the vicinity of the pore complex are expected to play an important role.

Exploring the biotransformation of N-(2-hydroxyphenyl)-2-propylpentanamide (an aryl valproic acid derivative) by CYP2C11, using in silico predictions and in vitro studies

OBJECTIVES

N-(2-hydroxyphenyl)-2-propylpentanamide (HO-AAVPA), a derivative of valproic acid (VPA), has been proposed as a potential anticancer agent due to its improved antiproliferative effects in some cancer cell lines. Although there is evidence that VPA is metabolized by cytochrome P450 2C11 rat isoform, HO-AAVPA CYP-mediated metabolism has not yet been fully explored. Therefore, in this work, the biotransformation of HO-AAVPA by CYP2C11 was investigated.

METHODS

Kinetic parameters and spectral interaction between HO-AAVPA and CYP were evaluated using rat liver microsomes. The participation of CYP2C11 in metabolism of HO-AAVPA was confirmed by cimetidine (CIM) inhibition assay. Docking and molecular dynamics simulations coupled to MMGBSA methods were used in theoretical study.

KEY FINDINGS

HO-AAVPA is metabolized by CYP enzymes (K_M  = 38.94 µm), yielding a hydroxylated metabolite according to its HPLC retention time (5.4 min) and MS analysis (252.2 m/z). In addition, CIM inhibition in rat liver microsomes (K_i  = 59.23 µm) confirmed that CYP2C11 is mainly involved in HO-AAVPA metabolism. Furthermore, HO-AAVPA interacts with CYP2C11 as a type I ligand. HO-AAVPA is stabilized at the CYP2C11 ligand recognition site through a map of interactions similar to other typical CYP2C11 substrates.

CONCLUSION

Therefore, rat liver CYP2C11 isoform is able to metabolize HO-AAVPA.

The quaternary assembly of KRas4B with Raf-1 at the membrane

Proximally located in the membrane, oncogenic Ras dimers (or nanoclusters) can recruit and promote Raf dimerization and MAPK (Raf/MEK/ERK) signaling. Among Ras isoforms, KRas4B is the most frequently mutated. Recent data on the binary KRas4B–Raf-1 complex suggested that Raf-1 CRD not only executes membrane anchorage, but also supports the high-affinity interaction of Raf-1 RBD with KRas4B catalytic domain. For a detailed mechanistic picture of Raf activation at the membrane, we employ explicit MD simulations of the quaternary KRas4B–Raf-1 complex. The complex contains two active GTP-bound KRas4B proteins forming a dimer through the allosteric lobe interface and two tandem RBD-CRD segments of Raf-1 interacting with the effector lobes at both ends of the KRas4B dimer. We show that Raf-1 RBD-CRD supports stable KRas4B dimer at preferred interface and orientation at the membrane, thereby cooperatively enhancing the affinity of the KRas4B–Raf-1 interaction. We propose that a Ras dimer at the membrane can increase the population of proximal Raf kinase domains, promoting kinase domain dimerization in the cytoplasm. Collectively, the dynamic Ras–Raf assembly promotes Raf activation not by allostery; instead, Ras activates Raf by shifting its ensemble toward kinase domain-accessible states through enhanced affinity at the membrane.

Quantitative Characterization of Protein-Lipid Interactions by Free Energy Simulation between Binary Bilayers

Using a recently developed binary bilayer system (BBS) consisting of two patches of laterally contacting bilayers, umbrella sampling molecular dynamics (MD) simulations were performed for quantitative characterization of protein-lipid interactions. The BBS is composed of 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) with an embedded model membrane protein, a gramicidin A (gA) channel. The calculated free energy difference for the transfer of a gA channel from DLPC (hydrophobic thickness ≈ 21.5 Å) to DMPC (hydrophobic thickness ≈ 25.5 Å) bilayers, ΔG(DLPC → DMPC), is -2.2 ± 0.7 kcal/mol. This value appears at odds with the traditional view that the hydrophobic length of the gA channel is ∼22 Å. To understand this discrepancy, we first note that recent MD simulations by different groups have shown that lipid bilayer thickness profiles in the vicinity of a gA channel differ qualitatively from the deformation profile predicted from continuum elastic bilayer models. Our MD simulations at low and high gA:lipid molar ratios and different membrane compositions indicate that the gA channel's effective hydrophobic length is ∼26 Å. Using this effective hydrophobic length, ΔG(DLPC → DMPC) determined here is in excellent agreement with predictions based on continuum elastic models (-3.0 to -2.2 kcal/mol) where the bilayer deformation energy is approximated as a harmonic function of the mismatch between the channel's effective hydrophobic length and the hydrophobic thickness of the bilayer. The free energy profile for gA in the BBS includes a barrier at the interface between the two bilayers which can be attributed to the line tension at the interface between two bilayers with different hydrophobic thicknesses. This observation implies that translation of a peptide between two different regions of a cell membrane (such as between the liquid ordered and disordered phases) may include effects of a barrier at the interface in addition to the relative free energies of the species far from the interface. The BBS allows for direct transfer free energy calculations between bilayers without a need of a reference medium, such as bulk water, and thus provides an efficient simulation protocol for the quantitative characterization of protein-lipid interactions at all-atom resolution.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

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.

Challenges facing nanotoxicology and nanomedicine due to cellular diversity

This review examines the interaction of nanomaterials (NMs) with cells from the perspective of major cellular differentiations. The structure and composition of cells reflect their role and function in a particular organ or environment. The normal differentiated-state and diseased cells may respond to NMs very differently. This review progresses with due care on nanotoxicology while emphasizing the potential of NMs in treating stress-associated disorders, including cancer and degeneration. The striking potential of NMs in inducing ROS, scavenging ROS, depleting cellular antioxidants, replenishing antioxidants, mimicking antioxidant enzyme activity, and modulating the immune system all show their considerable potential in treating cancer and other aging-associated disorders. It is now clear that NMs become more active and versatile when they come into contact with biological machinery, surprisingly in some cases, in a manner dependent on cell type. The mechanisms leading to the contrasting bioresponse of NMs ranging from toxicity to anticancer and from cell survival to carcinogenicity followed by their immuno-modulating potential show NMs to be a highly promising agent in biomedical therapy. This first-of-its-kind article seeks the challenges to be addressed that could provide a solid rationale in translating the promises of nanomedicine. A thorough understanding of normal and cancer biology could help to minimize the gap between basic and translational research in nanotechnology-based therapy.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

The syndromic deafness mutation G12R impairs fast and slow gating in Cx26 hemichannels

Mutations in connexin 26 hemichannels that cause syndromic deafness have a gain-of-function phenotype that is poorly understood. García et al. show that one such mutation impairs fast and slow gating in these hemichannels because of an interaction between the N terminus and intracellular loop.

Mutations in connexin 26 (Cx26) hemichannels can lead to syndromic deafness that affects the cochlea and skin. These mutations lead to gain-of-function hemichannel phenotypes by unknown molecular mechanisms. In this study, we investigate the biophysical properties of the syndromic mutant Cx26G12R (G12R). Unlike wild-type Cx26, G12R macroscopic hemichannel currents do not saturate upon depolarization, and deactivation is faster during hyperpolarization, suggesting that these channels have impaired fast and slow gating. Single G12R hemichannels show a large increase in open probability, and transitions to the subconductance state are rare and short-lived, demonstrating an inoperative fast gating mechanism. Molecular dynamics simulations indicate that G12R causes a displacement of the N terminus toward the cytoplasm, favoring an interaction between R12 in the N terminus and R99 in the intracellular loop. Disruption of this interaction recovers the fast and slow voltage-dependent gating mechanisms. These results suggest that the mechanisms of fast and slow gating in connexin hemichannels are coupled and provide a molecular mechanism for the gain-of-function phenotype displayed by the syndromic G12R mutation.

Lipids Alter Rhodopsin Function via Ligand-like and Solvent-like Interactions

Rhodopsin, a prototypical G protein-coupled receptor, is a membrane protein that can sense dim light. This highly effective photoreceptor is known to be sensitive to the composition of its lipidic environment, but the molecular mechanisms underlying this fine-tuned modulation of the receptor's function and structural stability are not fully understood. There are two competing hypotheses to explain how this occurs: 1) lipid modulation occurs via solvent-like interactions, where lipid composition controls membrane properties like hydrophobic thickness, which in turn modulate the protein's conformational equilibrium; or 2) protein-lipid interactions are ligand-like, with specific hot spots and long-lived binding events. By analyzing an ensemble of all-atom molecular dynamics simulations of five different states of rhodopsin, we show that a local ordering effect takes place in the membrane upon receptor activation. Likewise, docosahexaenoic acid acyl tails and phosphatidylethanolamine headgroups behave like weak ligands, preferentially binding to the receptor in inactive-like conformations and inducing subtle but significant structural changes.

Metabolite amyloid-like fibrils interact with model membranes

Metabolite assemblies interaction with membranes further extend the “amyloid hypothesis” to include small metabolites which serve as amyloidogenic building blocks.

Molecular modeling study on resistance of WT/D473H SMO to antagonists LDE-225 and LEQ-506

The smoothened (SMO) receptor, an essential signal transducer in the Hedgehog pathway, was targeted with antagonists to suppress the tumor. It is interesting that SMO D473H mutation confers resistance on inhibitor LDE-225 rather than LEQ-506. In this paper, the binding modes of them against the wild type and mutant SMO receptors were identified to gain insights into the resistant and non-resistant factors, based on a comprehensive protocol involving molecular docking, molecular dynamic simulations, free energy calculation and decomposition. A comparison of resistant LDE-225 and non-resistant LEQ-506 indicates that the volume of the binding cavity decreases seriously in the mutant complex with resistant LDE-225. In addition, the D473H mutation disrupts the hydrogen bond network with residues R400 and Q477, which results in the TM6 conformation inward. Owing to the absence of the hydrogen bond, residues R400 and Q477 make weak contributions to LDE-225. However, the D473H mutation along with TM6 conformational change has no effect on non-resistant LEQ-506. Finally, the resistance ascribes to adverse interaction between the greater polarity of mutant residue H473 and the nonpolar phenmethyl of LDE-225. The elaborate insights into structural and energetic mechanism of drug resistance provide an effective strategy to design rationally non-resistant antagonists.

Molecular Dynamics Methodologies for Probing Cannabinoid Ligand/Receptor Interaction

The cannabinoid type 1 and 2 G-protein-coupled receptors are currently important pharmacological targets with significant drug discovery potential. These receptors have been shown to display functional selectivity or biased agonism, a property currently thought to have substantial therapeutic potential. Although recent advances in crystallization techniques have provided a wealth of structural information about this important class of membrane-embedded proteins, these structures lack dynamical information. In order to fully understand the interplay of structure and function for this important class of proteins, complementary techniques that address the dynamical aspects of their function are required such as NMR as well as a variety of other spectroscopies. Complimentary to these experimental approaches is molecular dynamics, which has been effectively used to help unravel, at the atomic level, the dynamics of ligand binding and activation of these membrane-bound receptors. Here, we discuss and present several representative examples of the application of molecular dynamics simulations to the understanding of the signatures of ligand-binding and -biased signaling at the cannabinoid type 1 and 2 receptors.

CHARMM-GUI 10 years for biomolecular modeling and simulation

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.

Membrane insertion of fusion peptides from Ebola and Marburg viruses studied by replica-exchange molecular dynamics simulations

This work presents replica-exchange molecular dynamics simulations of inserting a 16-residue Ebola virus fusion peptide into a membrane bilayer. A computational approach is applied for modeling the peptide at the explicit all-atom level and the membrane-aqueous bilayer by a generalized Born continuum model with a smoothed switching function (GBSW). We provide an assessment of the model calculations in terms of three metrics: (1) the ability to reproduce the NMR structure of the peptide determined in the presence of SDS micelles and comparable structural data on other fusion peptides; (2) determination of the effects of the mutation Trp-8 to Ala and sequence discrimination of the homologous Marburg virus; and (3) calculation of potentials of mean force for estimating the partitioning free energy and their comparison to predictions from the Wimley-White interfacial hydrophobicity scale. We found the GBSW implicit membrane model to produce results of limited accuracy in conformational properties of the peptide when compared to the NMR structure, yet the model resolution is sufficient to determine the effect of sequence differentiation on peptide-membrane integration. © 2016 Wiley Periodicals, Inc.

Membrane-associated Ras dimers are isoform-specific: K-Ras dimers differ from H-Ras dimers

Are the dimer structures of active Ras isoforms similar? This question is significant since Ras can activate its effectors as a monomer; however, as a dimer, it promotes Raf's activation and MAPK (mitogen-activated protein kinase) cell signalling. In the present study, we model possible catalytic domain dimer interfaces of membrane-anchored GTP-bound K-Ras4B and H-Ras, and compare their conformations. The active helical dimers formed by the allosteric lobe are isoform-specific: K-Ras4B-GTP favours the α3 and α4 interface; H-Ras-GTP favours α4 and α5. Both isoforms also populate a stable β-sheet dimer interface formed by the effector lobe; a less stable β-sandwich interface is sustained by salt bridges of the β-sheet side chains. Raf's high-affinity β-sheet interaction is promoted by the active helical interface. Collectively, Ras isoforms' dimer conformations are not uniform; instead, the isoform-specific dimers reflect the favoured interactions of the HVRs (hypervariable regions) with cell membrane microdomains, biasing the effector-binding site orientations, thus isoform binding selectivity.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

Understanding the molecular basis of agonist/antagonist mechanism of GPER1/GPR30 through structural and energetic analyses

The G-protein coupled receptors (GPCRs) represent the largest superfamily of membrane proteins in charge to pass the cell signaling after binding with their cognate ligands to the cell interior. In breast cancer, a GPCR named GPER1 plays a key role in the process of growth and the proliferation of cancer cells. In a previous study, theoretical methods were applied to construct a model of GPER1, which later was submitted to molecular dynamics (MD) simulations to perform a docking calculation. Based on this preceding work, it is known that GPER1 is sensitive to structural differences in its binding site. However, due to the nature of that past study, conformational changes linked to the ligand binding were not observed. Therefore, in this study, in order to explore the conformational changes coupled to the agonist/antagonist binding, MD simulations of about 0.25μs were performed for the free and bound states, summarizing 0.75μs of MD simulation in total. For the bound states, one agonist (G-1) and antagonist (G-15) were chosen since is widely known that these two molecules cause an impact on GPER1 mobility. Based on the conformational ensemble generated through MD simulations, we found that despite G-1 and G-15 being stabilized by similar map of residues, the structural differences between both ligands impact the hydrogen bond pattern not only at the GPER1 binding site but also along the seven-helix bundle, causing significant differences in the conformational mobility along the extracellular and cytoplasmic domain, and to a lesser degree in the curvatures of helix 2, helix 3 and helix 7 between the free and bound states, which is in agreement with reported literature, and might be linked to microscopic characteristics of the activated-inactivated transition. Furthermore, binding free energy calculations using the MM/GBSA method for the bound states, followed by an alanine scanning analysis allowed us to identify some important residues for the complex stabilization.

The higher level of complexity of K-Ras4B activation at the membrane

Is nucleotide exchange sufficient to activate K-Ras4B? To signal, oncogenic rat sarcoma (Ras) anchors in the membrane and recruits effectors by exposing its effector lobe. With the use of NMR and molecular dynamics (MD) simulations, we observed that in solution, farnesylated guanosine 5'-diphosphate (GDP)-bound K-Ras4B is predominantly autoinhibited by its hypervariable region (HVR), whereas the GTP-bound state favors an activated, HVR-released state. On the anionic membrane, the catalytic domain adopts multiple orientations, including parallel (∼180°) and perpendicular (∼90°) alignments of the allosteric helices, with respect to the membrane surface direction. In the autoinhibited state, the HVR is sandwiched between the effector lobe and the membrane; in the active state, with membrane-anchored farnesyl and unrestrained HVR, the catalytic domain fluctuates reinlessly, exposing its effector-binding site. Dimerization and clustering can reduce the fluctuations. This achieves preorganized, productive conformations. Notably, we also observe HVR-autoinhibited K-Ras4B-GTP states, with GDP-bound-like orientations of the helices. Thus, we propose that the GDP/GTP exchange may not be sufficient for activation; instead, our results suggest that the GDP/GTP exchange, HVR sequestration, farnesyl insertion, and orientation/localization of the catalytic domain at the membrane conjointly determine the active or inactive state of K-Ras4B. Importantly, K-Ras4B-GTP can exist in active and inactive states; on its own, GTP binding may not compel K-Ras4B activation.-Jang, H., Banerjee, A., Chavan, T. S, Lu, S., Zhang, J., Gaponenko, V., Nussinov, R. The higher level of complexity of K-Ras4B activation at the membrane.

Predicting peptide vaccine candidates against H1N1 influenza virus through theoretical approaches

Identification of potential epitopes that might activate the immune system has been facilitated by the employment of algorithms that use experimental data as templates. However, in order to prove the affinity and the map of interactions between the receptor (major histocompatibility complex, MHC, or T-cell receptor) and the potential epitope, further computational studies are required. Docking and molecular dynamics (MDs) simulations have been an effective source of generating structural information at molecular level in immunology. Herein, in order to provide a detailed understanding of the origins of epitope recognition and to select the best peptide candidate to develop an epitope-based vaccine, docking and MDs simulations in combination with MMGBSA free energy calculations and per-residue free energy decomposition were performed, taking as starting complexes those formed between four designed epitopes (P1-P4) from hemagglutinin (HA) of the H1N1 influenza virus and MHC-II anchored in POPC membrane. Our results revealed that the energetic contributions of individual amino acids within the pMHC-II complexes are mainly dictated by van der Waals interactions and the nonpolar part of solvation energy, whereas the electrostatic interactions corresponding to hydrogen bonds and salt bridges determine the binding specificity, being the most favorable interactions formed between p4 and MHC-II. Then, P1-P4 epitopes were synthesized and tested experimentally to compare theoretical and experimental results. Experimental results show that P4 elicited the highest strong humoral immune response to HA of the H1N1 and may induce antibodies that are cross-reactive to other influenza subtypes, suggesting that it could be a good candidate for the development of a peptide-based vaccine.

Computational Methods for Structural and Functional Studies of Alzheimer's Amyloid Ion Channels

Aggregation can be studied by a range of methods, experimental and computational. Aggregates form in solution, across solid surfaces, and on and in the membrane, where they may assemble into unregulated leaking ion channels. Experimental probes of ion channel conformations and dynamics are challenging. Atomistic molecular dynamics (MD) simulations are capable of providing insight into structural details of amyloid ion channels in the membrane at a resolution not achievable experimentally. Since data suggest that late stage Alzheimer's disease involves formation of toxic ion channels, MD simulations have been used aiming to gain insight into the channel shapes, morphologies, pore dimensions, conformational heterogeneity, and activity. These can be exploited for drug discovery. Here we describe computational methods to model amyloid ion channels containing the β-sheet motif at atomic scale and to calculate toxic pore activity in the membrane.

Characterization of a novel water pocket inside the human Cx26 hemichannel structure

Connexins (Cxs) are a family of vertebrate proteins constituents of gap junction channels (GJCs) that connect the cytoplasm of adjacent cells by the end-to-end docking of two Cx hemichannels. The intercellular transfer through GJCs occurs by passive diffusion allowing the exchange of water, ions, and small molecules. Despite the broad interest to understand, at the molecular level, the functional state of Cx-based channels, there are still many unanswered questions regarding structure-function relationships, perm-selectivity, and gating mechanisms. In particular, the ordering, structure, and dynamics of water inside Cx GJCs and hemichannels remains largely unexplored. In this work, we describe the identification and characterization of a believed novel water pocket-termed the IC pocket-located in-between the four transmembrane helices of each human Cx26 (hCx26) monomer at the intracellular (IC) side. Using molecular dynamics (MD) simulations to characterize hCx26 internal water structure and dynamics, six IC pockets were identified per hemichannel. A detailed characterization of the dynamics and ordering of water including conformational variability of residues forming the IC pockets, together with multiple sequence alignments, allowed us to propose a functional role for this cavity. An in vitro assessment of tracer uptake suggests that the IC pocket residue Arg-143 plays an essential role on the modulation of the hCx26 hemichannel permeability.

The Melibiose Transporter of Escherichia coli: CRITICAL CONTRIBUTION OF LYS-377 TO THE STRUCTURAL ORGANIZATION OF THE INTERACTING SUBSTRATE BINDING SITES

We examine the role of Lys-377, the only charged residue in helix XI, on the functional mechanism of the Na(+)-sugar melibiose symporter from Escherichia coli. Intrinsic fluorescence, FRET, and Fourier transform infrared difference spectroscopy reveal that replacement of Lys-377 with either Cys, Val, Arg, or Asp disables both Na(+) and melibiose binding. On the other hand, molecular dynamics simulations extending up to 200-330 ns reveal that Lys-377 (helix XI) interacts with the anionic side chains of two of the three putative ligands for cation binding (Asp-55 and Asp-59 in helix II). When Asp-59 is protonated during the simulations, Lys-377 preferentially interacts with Asp-55. Interestingly, when a Na(+) ion is positioned in the Asp-55-Asp-59 environment, Asp-124 in helix IV (a residue essential for melibiose binding) reorients and approximates the Asp-55-Asp-59 pair, and all three acidic side chains act as Na(+) ligands. Under these conditions, the side chain of Lys-377 interacts with the carboxylic moiety of these three Asp residues. These data highlight the crucial role of the Lys-377 residue in the spatial organization of the Na(+) binding site. Finally, the analysis of the second-site revertants of K377C reveals that mutation of Ile-22 (in helix I) preserves Na(+) binding, whereas that of melibiose is largely abolished according to spectroscopic measurements. This amino acid is located in the border of the sugar-binding site and might participate in sugar binding through apolar interactions.

Alchembed: A Computational Method for Incorporating Multiple Proteins into Complex Lipid Geometries

A necessary step prior to starting any membrane protein computer simulation is the creation of a well-packed configuration of protein(s) and lipids. Here, we demonstrate a method, alchembed, that can simultaneously and rapidly embed multiple proteins into arrangements of lipids described using either atomistic or coarse-grained force fields. During a short simulation, the interactions between the protein(s) and lipids are gradually switched on using a soft-core van der Waals potential. We validate the method on a range of membrane proteins and determine the optimal soft-core parameters required to insert membrane proteins. Since all of the major biomolecular codes include soft-core van der Waals potentials, no additional code is required to apply this method. A tutorial is included in the Supporting Information.

Mechanisms of membrane binding of small GTPase K-Ras4B farnesylated hypervariable region

K-Ras4B belongs to a family of small GTPases that regulates cell growth, differentiation and survival. K-ras is frequently mutated in cancer. K-Ras4B association with the plasma membrane through its farnesylated and positively charged C-terminal hypervariable region (HVR) is critical to its oncogenic function. However, the structural mechanisms of membrane association are not fully understood. Here, using confocal microscopy, surface plasmon resonance, and molecular dynamics simulations, we observed that K-Ras4B can be distributed in rigid and loosely packed membrane domains. Its membrane binding domain interaction with phospholipids is driven by membrane fluidity. The farnesyl group spontaneously inserts into the disordered lipid microdomains, whereas the rigid microdomains restrict the farnesyl group penetration. We speculate that the resulting farnesyl protrusion toward the cell interior allows oligomerization of the K-Ras4B membrane binding domain in rigid microdomains. Unlike other Ras isoforms, K-Ras4B HVR contains a single farnesyl modification and positively charged polylysine sequence. The high positive charge not only modulates specific HVR binding to anionic phospholipids but farnesyl membrane orientation. Phosphorylation of Ser-181 prohibits spontaneous farnesyl membrane insertion. The mechanism illuminates the roles of HVR modifications in K-Ras4B targeting microdomains of the plasma membrane and suggests an additional function for HVR in regulation of Ras signaling.

Molecular dynamics simulations of membrane proteins

Membrane protein structures are underrepresented in the Protein Data Bank (PDB) due to difficulties associated with expression and crystallization. As such, it is one area where computational studies, particularly Molecular Dynamics (MD) simulations, can provide useful additional information. Recently, there has been substantial progress in the simulation of lipid bilayers and membrane proteins embedded within them. Initial efforts at simulating membrane proteins embedded within a lipid bilayer were relatively slow and interactive processes, but recent advances now mean that the setup and running of membrane protein simulations is somewhat more straightforward, though not without its problems. In this chapter, we outline practical methods for setting up and running MD simulations of a membrane protein embedded within a lipid bilayer and discuss methodologies that are likely to contribute future improvements.

The mechanosensitive channel of small conductance (MscS) functions as a Jack-in-the box

Phenotypical analysis of the lipid interacting residues in the closed state of the mechanosensitive channel of small conductance (MscS) from Escherichia coli (E. coli) has previously shown that these residues are critical for channel function. In the closed state, mutation of individual hydrophobic lipid lining residues to alanine, thus reducing the hydrophobicity, resulted in phenotypic changes that were observable using in vivo assays. Here, in an analogous set of experiments, we identify eleven residues in the first transmembrane domain of the open state of MscS that interact with the lipid bilayer. Each of these residues was mutated to alanine and leucine to modulate their hydrophobic interaction with the lipid tail-groups in the open state. The effects of these changes on channel function were analyzed using in vivo bacterial assays and patch clamp electrophysiology. Mutant channels were found to be functionally indistinguishable from wildtype MscS. Thus, mutation of open-state lipid interacting residues does not differentially stabilize or destabilize the open, closed, intermediate, or transition states of MscS. Based on these results and other data from the literature, we propose a new gating paradigm for MscS where MscS acts as a "Jack-In-The-Box" with the intrinsic bilayer lateral pressure holding the channel in the closed state. In this model, upon application of extrinsic tension the channel springs into the open state due to relief of the intrinsic lipid bilayer pressure.

All-atom molecular dynamics simulation of a photosystem i/detergent complex

All-atom molecular dynamics (MD) simulation was used to investigate the solution structure and dynamics of the photosynthetic pigment-protein complex photosystem I (PSI) from Thermosynechococcus elongatus embedded in a toroidal belt of n-dodecyl-β-d-maltoside (DDM) detergent. Evaluation of root-mean-square deviations (RMSDs) relative to the known crystal structure show that the protein complex surrounded by DDM molecules is stable during the 200 ns simulation time, and root-mean-square fluctuation (RMSF) analysis indicates that regions of high local mobility correspond to solvent-exposed regions such as turns in the transmembrane α-helices and flexible loops on the stromal and lumenal faces. Comparing the protein-detergent complex to a pure detergent micelle, the detergent surrounding the PSI trimer is found to be less densely packed but with more ordered detergent tails, contrary to what is seen in most lipid bilayer models. We also investigated any functional implications for the observed conformational dynamics and protein-detergent interactions, discovering interesting structural changes in the psaL subunits associated with maintaining the trimeric structure of the protein. Importantly, we find that the docking of soluble electron mediators such as cytochrome c6 and ferredoxin to PSI is not significantly impacted by the solubilization of PSI in detergent.

Smoothened (SMO) receptor mutations dictate resistance to vismodegib in basal cell carcinoma

Basal cell carcinomas (BCCs) and a subset of medulloblastomas are characterized by loss-of-function mutations in the tumor suppressor gene, PTCH1. PTCH1 normally functions by repressing the activity of the Smoothened (SMO) receptor. Inactivating PTCH1 mutations result in constitutive Hedgehog pathway activity through uncontrolled SMO signaling. Targeting this pathway with vismodegib, a novel SMO inhibitor, results in impressive tumor regression in patients harboring genetic defects in this pathway. However, a secondary mutation in SMO has been reported in medulloblastoma patients following relapse on vismodegib to date. This mutation preserves pathway activity, but appears to confer resistance by interfering with drug binding. Here we report for the first time on the molecular mechanisms of resistance to vismodegib in two BCC cases. The first case, showing progression after 2 months of continuous vismodegib (primary resistance), exhibited the new SMO G497W mutation. The second case, showing a complete clinical response after 5 months of treatment and a subsequent progression after 11 months on vismodegib (secondary resistance), exhibited a PTCH1 nonsense mutation in both the pre- and the post-treatment specimens, and the SMO D473Y mutation in the post-treatment specimens only. In silico analysis demonstrated that SMO(G497W) undergoes a conformational rearrangement resulting in a partial obstruction of the protein drug entry site, whereas the SMO D473Y mutation induces a direct effect on the binding site geometry leading to a total disruption of a stabilizing hydrogen bond network. Thus, the G497W and D473Y SMO mutations may represent two different mechanisms leading to primary and secondary resistance to vismodegib, respectively.

Computational investigation of cholesterol binding sites on mitochondrial VDAC

The mitochondrial voltage-dependent anion channel (VDAC) allows passage of ions and metabolites across the mitochondrial outer membrane. Cholesterol binds mammalian VDAC, and we investigated the effects of binding to human VDAC1 with atomistic molecular dynamics simulations that totaled 1.4 μs. We docked cholesterol to specific sites on VDAC that were previously identified with NMR, and we tested the reliability of multiple docking results in each site with simulations. The most favorable binding modes were used to build a VDAC model with cholesterol occupying five unique sites, and during multiple 100 ns simulations, cholesterol stably and reproducibly remained bound to the protein. For comparison, VDAC was simulated in systems with identical components but with cholesterol initially unbound. The dynamics of loops that connect adjacent β-strands were most affected by bound cholesterol, with the averaged root-mean-square fluctuation (RMSF) of multiple residues altered by 20–30%. Cholesterol binding also stabilized charged residues inside the channel and localized the surrounding electrostatic potentials. Despite this, ion diffusion through the channel was not significantly affected by bound cholesterol, as evidenced by multi-ion potential of mean force measurements. Although we observed modest effects of cholesterol on the open channel, our model will be particularly useful in experiments that investigate how cholesterol affects VDAC function under applied electrochemical forces and also how other ligands and proteins interact with the channel.

Structure-based engineering of a minimal porin reveals loop-independent channel closure

Porins, like outer membrane protein G (OmpG) of Escherichia coli, are ideal templates among ion channels for protein and chemical engineering because of their robustness and simple architecture. OmpG shows fast transitions between open and closed states, which were attributed to loop 6 (L6). As flickering limits single-channel-based applications, we pruned L6 by either 8 or 12 amino acids. While the open probabilities of both L6 variants resemble that of native OmpG, their gating frequencies were reduced by 63 and 81%, respectively. Using the 3.2 Å structure of the shorter L6 variant in the open state, we engineered a minimal porin (220 amino acids), where all remaining extramembranous loops were truncated. Unexpectedly, this minimized porin still exhibited gating, but it was 5-fold less frequent than in OmpG. The residual gating of the minimal pore is hence independent of L6 rearrangements and involves narrowing of the ion conductance pathway most probably driven by global stretching-flexing deformations of the membrane-embedded β-barrel.