Characterizing the interactions of two lipid modifications with lipid rafts: farnesyl anchors vs. palmitoyl anchors

Farnesyl (Far) and palmitoyl (Pal) anchors play important roles in the traffic of many lipidated proteins. Herein, we show the distinctive interactions and influences of the two lipid modifications on lipid rafts (LRs) and non-raft-like membranes using molecular dynamics simulations. Palmitoyl anchors behave in a more ordered fashion, pack tighter with the lipids of LRs and diffuse at a slower rate than farnesyl anchors in LRs. When interacting with non-raft-like membranes these two types of anchors become less ordered, pack more loosely with lipids, and diffuse at a higher rate. By calculating both the number of contacts per chain and the number of contact atoms per carbon of the two anchors with the lipid components, we found that the palmitoyl chains preferred to associate with the saturated chains of lipids and cholesterol molecules in LRs, while farnesyl chains favored association with saturated chains and unsaturated chains. For non-raft-like membranes, these two lipid anchors had roughly the same preference for the three types of contact lipid chains. Additionally, palmitoyl anchors caused cholesterol to orient more perpendicular to the membrane surface, surrounding lipids to become more ordered, and lipid lateral fluidity to reduce significantly, compared to farnesyl anchors in LRs. By contrast, the POPE and DSPC became much less ordered, cholesterol became more tilted, and lipids became more fluid, when the two types lipid anchors were inserted in non-raft-like membranes. These findings are useful for understanding the traffic mechanisms of lipidated proteins with farnesyl and palmitoyl modifications in cell membranes.

Probing the disparate effects of arginine and lysine residues on antimicrobial peptide/bilayer association

Antimicrobial peptides (AMPs) are key components of the innate immune response and represent promising templates for the development of broad-spectrum alternatives to conventional antibiotics. Most AMPs are short, cationic peptides that interact more strongly with negatively charged prokaryotic membranes than net neutral eukaryotic ones. Both AMPs and synthetic analogues with arginine-like side chains are more active against bacteria than those with lysine-like amine groups, though the atomistic mechanism for this increase in potency remains unclear. To examine this, we conducted comparative molecular dynamics simulations of a model negatively-charged membrane system interacting with two mutants of the AMP KR-12: one with lysine residues mutated to arginines (R-KR12) and one with arginine residues mutated to lysine (K-KR12). Simulations show that both partition analogously to the bilayer and display similar preferences for hydrogen bonding with the anionic POPGs. However, R-KR12 binds stronger to the bilayer than K-KR12 and forms significantly more hydrogen bonds, leading to considerably longer interaction times. Additional simulations with methylated R-KR12 and charge-modified K-KR12 mutants show that the extensive interaction seen in the R-KR12 system is partly due to arginine's strong atomic charge distribution, rather than being purely an effect of the greater number of hydrogen bond donors. Finally, free energy simulations reveal that both peptides are disordered in solution but form an amphipathic α-helix when inserted into the bilayer headgroup region. Overall, these results highlight the role of charge and hydrogen bond strength in peptide bilayer insertion, and offer potential insights for designing more potent analogues in the future.

Transcriptional and computational study of expansins differentially expressed in response to inclination in radiata pine

Plants have the ability to reorient their vertical growth when exposed to inclination. This response can be as quick as 2 h in inclined young pine (Pinus radiata D. Don) seedlings, with over accumulation of lignin observed after 9 days s. Several studies have identified expansins involved in cell expansion among other developmental processes in plants. Six putative expansin genes were identified in cDNA libraries isolated from inclined pine stems. A differential transcript abundance was observed by qPCR analysis over a time course of inclination. Five genes changed their transcript accumulation in both stem sides in a spatial and temporal manner compared with non-inclined stem. To compare these expansin genes, and to suggest a possible mechanism of action at molecular level, the structures of the predicted proteins were built by comparative modeling methodology. An open groove on the surface of the proteins composed of conserved zresidues was observed. Using a cellulose polymer as ligand the protein-ligand interaction was evaluated, with the results showing differences in the protein-ligand interaction mode. Differences in the binding energy interaction can be explained by changes in some residues that generate differences in electrostatic surface in the open groove region, supporting the participation of six members of multifamily proteins in this specific process. The data suggests participation of different expansin proteins in the dissembling and remodeling of the complex cell wall matrix during the reorientation response to inclination.

Understanding the colloidal dispersion stability of 1D and 2D materials: Perspectives from molecular simulations and theoretical modeling

The colloidal dispersion stability of 1D and 2D materials in the liquid phase is critical for scalable nano-manufacturing, chemical modification, composites production, and deployment as conductive inks or nanofluids. Here, we review recent computational and theoretical studies carried out by our group to model the dispersion stability of 1D and 2D materials, including single-walled carbon nanotubes, graphene, and graphene oxide in aqueous surfactant solutions or organic solvents. All-atomistic (AA) molecular dynamics (MD) simulations can probe the molecular level details of the adsorption morphology of surfactants and solvents around these materials, as well as quantify the interaction energy between the nanomaterials mediated by surfactants or solvents. Utilizing concepts from reaction kinetics and diffusion, one can directly predict the rate constants for the aggregation kinetics and dispersion life times using MD outputs. Furthermore, the use of coarse-grained (CG) MD simulations allows quantitative prediction of surfactant adsorption isotherms. Combined with the Poisson-Boltzmann equation, the Langmuir isotherm, and the DLVO theory, one can directly use CGMD outputs to: (i) predict electrostatic potentials around the nanomaterial, (ii) correlate surfactant surface coverages with surfactant concentrations in the bulk dispersion medium, and (iii) determine energy barriers against coagulation. Finally, we discuss challenges associated with studying emerging 2D materials, such as, hexagonal boron nitride (h-BN), phosphorene, and transition metal dichalcogenides (TMDCs), including molybdenum disulfide (MoS₂). An outlook is provided to address these challenges with plans to develop force-field parameters for MD simulations to enable predictive modeling of emerging 2D materials in the liquid phase.

Spectroscopy analysis and molecular dynamics studies on the binding of penicillin V and sulbactam to beta-lactamase II from Bacillus cereus

The molecular recognition and interaction of beta-lactamase II from Bacillus cereus (Bc II) with penicillin V (PV) and sulbactam (Sul) especially conformational changes of Bc II in the binding process were studied through spectroscopy analysis in combination with molecular dynamics (MD) simulation. The results show that in the binding process, a new coordination bond is observed between the Zn₂ of Bc II and the carboxyl-O of PV or Sul by replacing His204. Electrostatic interaction between Zn₂ and the ligand provide main driving force for the binding affinity. Compared with apo Bc II, there are mainly four loops showing significant conformational changes in ligand-bound Bc II. A weak conformational transformation from β-sheets to random coils is observed in the loop2 of ligand-bound Bc II. The conformational transformation may depend on the functional group and binding pose of the ligand, giving the binding pocket greater flexibility and accordingly allowing for an induced fit of the enzyme-ligand binding site around the newly introduced ligand. The change in the loop2 of ligand-bound Bc II may lead to the opening of the binding pocket of Bc II. Therefore, loop2 can be considered a gate for control of ligand access in Bc II, hence its dynamic response should be considered in new drug design and development.

Simulated pressure denaturation thermodynamics of ubiquitin

Simulations of protein thermodynamics are generally difficult to perform and provide limited information. It is desirable to increase the degree of detail provided by simulation and thereby the potential insight into the thermodynamic properties of proteins. In this study, we outline how to analyze simulation trajectories to decompose conformation-specific, parameter free, thermodynamically defined protein volumes into residue-based contributions. The total volumes are obtained using established methods from Fluctuation Solution Theory, while the volume decomposition is new and is performed using a simple proximity method. Native and fully extended ubiquitin are used as the test conformations. Changes in the protein volumes are then followed as a function of pressure, allowing for conformation-specific protein compressibility values to also be obtained. Residue volume and compressibility values indicate significant contributions to protein denaturation thermodynamics from nonpolar and coil residues, together with a general negative compressibility exhibited by acidic residues.

Integrative modelling of TIR domain-containing adaptor molecule inducing interferon-β (TRIF) provides insights into its autoinhibited state

BACKGROUND

TRIF is a key protein in antiviral innate immunity, operating downstream of TLRs. TRIF activation leads to the production of interferon-β and pro-inflammatory cytokines. There is evidence from experiments to suggest that the N-terminal domain of TRIF binds to its TIR domain to avoid constitutive activation. However, no structure of a complex between the N-terminal domain and the TIR domain exists till date. The disordered nature of the region connecting the N-terminal domain and the TIR domain compounds the issue of elucidating the mechanism of autoinhibition of TRIF. In this study, we have employed an integrative approach consisting of mutual information analysis, docking, molecular dynamics simulations and residue network analysis, in combination with existing experimental data to provide a glimpse of TRIF in its autoinhibited state.

RESULTS

Our extensive docking approach reveals that the N-terminal domain binds to the BB loop-B helix region of the TIR domain, consistent with experimental observations. Long length molecular dynamics simulations of 1 microsecond performed on the docked model highlights residues participating in hydrogen bonding and hydrophobic interactions at the interface. A pair of residues present in the vicinity of the interface is also predicted by mutual information analysis, to co-evolve. Residues mediating long-range interactions within the TIR domain of TRIF were identified using residue network analysis.

CONCLUSIONS

Based on the results of the modelling and residue network analysis, we propose that the N-terminal domain binds to the BB loop region of the TIR domain, thereby preventing its homodimersation. The binding of TRIF to TLR3 or TRAM could induce a slight conformational change, causing the interactions between the N-terminal domain and TIR domain to disrupt, thereby exposing the BB loop and rendering it amenable for higher-order oligomerisation.

REVIEWERS

This article was reviewed by Michael Gromiha, Srikrishna Subramaniam and Peter Bond (nominated by Chandra Verma).

The effect of β-glucan and its potential analog on the structure of Dectin-1 receptor

Dectin-1 is a recently discovered pattern-recognition receptor that plays an important role in antifungal innate immunity, which acts a specific receptor for β-glucan (BG). The present study, aimed at clarifying effect of BG and a new analog, maltotriose (MT) on Dectin-1 receptor. We implemented molecular docking of MT on Dectin-1 along with model-independent all-atom-molecular dynamics simulations. Simulations were carried out at three levels of complexity: (1) Apo-Dectin-1; (2) BG:Dectin-1; (3) MT:Dectin-1. All three system complexes were undergone stability check before showing a comparative analysis. A characteristic feature, noted for the MT:Dectin-1, is a shifting of loops (loop1 and loop2) orientation towards atoms of MT, a broad interaction suggested a robust and tight binding on comparison with BG:Dectin-1. Free energy estimation corroborated the observation, which furthermore, made a close agreement by revealing contribution of energy components of interacting residues. In addition, cluster analysis of complexes exhibit a smooth continuous transition to a new confirmation, represented by a series of clusters each having a longer lifetime. Principal component analysis revealed a broken pipe at binding site of BG:Dectin-1 during movement of atoms whereas in MT:Dectin-1 exhibited wide band and high amplitude motion of atoms in trajectory, was due to loop orientation toward MT. Observation was further shown by measuring distances and hydrogen binding calculation. Simulations of the BG:Dectin-1 and MT:Dectin-1 complex revealed first time the influence of BG and MT ligands. This study might extend the knowledge of the BG and MT interaction on Dectin-1 and proposed further potential bioassay of MT.

The pH-dependent assembly of Chaplin E from Streptomyces coelicolor

Chaplin E, is one of five self-assembling peptides secreted by Streptomyces coelicolor that assist aerial growth by lowering the surface tension of water. Although the surface activity of a mixture of chaplin peptides has observed to depend on pH, it is unclear how the solvent environment (i.e. pH) influences the structure, assembly and subsequent functionality of these individual peptides. In this study, the conformation and fibril forming propensity of the Chaplin E peptide was assessed as a function of pH using a combination of experimental measurements and molecular dynamics simulations. At an acidic pH of 3.0, Chaplin E retained a random coil structure, whereas at the isoelectric point of 6.7 or a basic pH of 10.0, Chaplin E rapidly formed amyloid fibrils rich in β-sheet structure with high efficiency (>93%). Molecular dynamics simulations indicate the persistence of greater α-helical content at the N-terminus at high pH; this is likely partly due to the lack of electrostatic repulsion between residues His6 and Lys10. Since fibril formation was observed at high but not at low pH, we propose that the presence of an N-terminal α-helix in the monomeric form of Chaplin E is required for aggregation and conversion to β-amyloid fibrils. The pH sensitivity of Chaplin E peptide structure provides a route to control peptide assembly and may be important for the physiological function of this peptide, as a surface active agent in the transition from vegetative to aerial growth and could assist Streptomyces coelicolor in response to environmental fluctuations in pH.

Influence of membrane composition on the binding and folding of a membrane lytic peptide from the non-enveloped flock house virus

Using a combination of coarse-grained and atomistic molecular dynamics simulations we have investigated the membrane binding and folding properties of the membrane lytic peptide of Flock House virus (FHV). FHV is an animal virus and an excellent model system for studying cell entry mechanisms in non-enveloped viruses. FHV undergoes a maturation event where the 44 C-terminal amino acids are cleaved from the major capsid protein, forming the membrane lytic (γ) peptides. Under acidic conditions, γ is released from the capsid interior allowing the peptides to bind and disrupt membranes. The first 21 N-terminal residues of γ, termed γ₁, have been resolved in the FHV capsid structure and γ₁ has been the subject of in vitro studies. γ₁ is structurally dynamic as it adopts helical secondary structure inside the capsid and on membranes, but it is disordered in solution. In vitro studies have shown the binding free energies to POPC or POPG membranes are nearly equivalent, but binding to POPC is enthalpically driven, while POPG binding is entropically driven. Through coarse-grained and multiple microsecond all-atom simulations the membrane binding and folding properties of γ₁ are investigated against homogeneous and heterogeneous bilayers to elucidate the dependence of the microenvironment on the structural properties of γ₁. Our studies provide a rationale for the thermodynamic data and suggest binding of γ₁ to POPG bilayers occurs in a disordered state, but γ₁ must adopt a helical conformation when binding POPC bilayers.

Molecular dynamics simulation studies of novel β-lactamase inhibitor

New Delhi Metallo-β-Lactamase-1 (NDM-1) has drawn great attention due to its diverse antibiotic resistant activity. It can hydrolyze almost all clinically available β-lactam antibiotics. To inhibit the activity of NDM-1 a new strategy is proposed using computational methods. Molecular dynamics (MD) simulations are used to analyze the molecular interactions between selected inhibitor candidates and NDM-1 structure. The enzyme-ligand complex is subject to binding free energy calculations using MM(PB/GB)SA methods. The role of each residue of the active site contributing in ligand binding affinity is explored using energy decomposition analysis. Furthermore, a hydrogen bonding network between ligand and enzyme active site is observed and key residues are identified ensuring that the ligand stays inside the active site and maintains its movement towards the active site pocket. A production run of 150ns is carried out and results are analyzed using root mean square deviation (RMSD), root mean square fluctuation (RMSF), and radius of gyration (Rg) to explain the stability of enzyme ligand complex. Important active site residue e.g. PHE70, VAL73, TRP93, HIS122, GLN123, ASP124, HIS189, LYS216, CYS208, LYS211, ALA215, HIS250, and SER251 were observed to be involved in ligand attachemet inside the active site pocket, hence depicting its inhibitor potential. Hydrogen bonds involved in structural stability are analyzed through radial distribution function (RDF) and contribution of important residues involved in ligand movement is explained using a novel analytical tool, axial frequency distribution (AFD) to observe the role of important hydrogen bonding partners between ligand atoms and active site residues.

Analysis of Species-Selectivity of Human, Mouse and Rat Cytochrome P450 1A and 2B Subfamily Enzymes using Molecular Modeling, Docking and Dynamics Simulations

Cytochrome P450 (CYP) 1A and 2B subfamily enzymes are important drug metabolizing enzymes, and are highly conserved across species in terms of sequence homology. However, there are major to minor structural and macromolecular differences which provide for species-selectivity and substrate-selectivity. Therefore, species-selectivity of CYP1A and CYP2B subfamily proteins across human, mouse and rat was analyzed using molecular modeling, docking and dynamics simulations when the chiral molecules quinine and quinidine were used as ligands. The three-dimensional structures of 17 proteins belonging to CYP1A and CYP2B subfamilies of mouse and rat were predicted by adopting homology modeling using the available structures of human CYP1A and CYP2B proteins as templates. Molecular docking and dynamics simulations of quinine and quinidine with CYP1A subfamily proteins revealed the existence of species-selectivity across the three species. On the other hand, in the case of CYP2B subfamily proteins, no role for chirality of quinine and quinidine in forming complexes with CYP2B subfamily proteins of the three species was indicated. Our findings reveal the roles of active site amino acid residues of CYP1A and CYP2B subfamily proteins and provide insights into species-selectivity of these enzymes across human, mouse, and rat.

Single-probe multistate detection of DNA via aggregation-induced emission on a graphene oxide platform

Graphene and graphene oxides (GO), or their reduced forms, have been introduced in a variety of biosensing platforms and have exhibited enhanced performance levels in these forms. We herein report a DNA sensing platform consisting of aggregation-induced emission (AIE) molecules and complementary DNA (comDNA) adsorbed on GO. We experimentally turned the AIE molecule on and off by adjusting its distance, which correlates with DNA structures as shown in our computational results, from the GO sheet, which quenches depending on its distance from the graphene plane. The changes in florescence are reproducible, which demonstrates the probe's ability to identify the binding state of the DNA. Our molecular dynamics simulation results reveal strong π-π interactions between single-strand DNA (ssDNA) and GO, which enable the ssDNA molecule to move closer to the graphene oxide. This reduces the center of mass and binding free energies in the simulation. When hybridized with comDNA, the increased distance, evidenced by the reduced interaction, eliminates the quenching effect and turns on the AIE molecule. Our protocol use of the AIE molecule as a probe thus avoids the complicated steps involved in covalent functionalization and allows the rapid and label-free detection of DNA molecules.

STATEMENT OF SIGNIFICANCE

A simple, rapid method of fluorescent measurement of DNA hybridization in the presence of graphene (oxide) is presented. Conventional fluorescent dyes offer high performance in biosensors. However, labeling procedures are synthetically demanding in time and resources making it less cost-effective. Molecules with aggregation-induced-emission (AIE) property have advantages over traditional fluorescent molecules because of their intrinsic preference for detection as a turn-on probe and their single-molecule detection ability. Previous work has shown AIE dyes act as excellent "label-free" bioprobes with high sensitivity but with limited selectivity. Graphene oxide (GO) with its unique optical properties and affinity to different kinds of biomolecules can be used as an auxiliary to enhance selectivity of AIE dyes. In this work, we report a label-free strategy to detect DNA of particular sequence by water-soluble AIE probes with the aid of GO, supported by the computational explanations for this phenomenon.

Insights into proton translocation in cbb3 oxidase from MD simulations

Heme-copper oxidases are membrane protein complexes that catalyse the final step of the aerobic respiration, namely the reduction of oxygen to water. The energy released during catalysis is coupled to the active translocation of protons across the membrane, which contributes to the establishment of an electrochemical gradient that is used for ATP synthesis. The distinctive C-type (or cbb₃) cytochrome c oxidases, which are mostly present in proteobacteria, exhibit a number of unique structural and functional features, including high catalytic activity at low oxygen concentrations. At the moment, the functioning mechanism of C-type oxidases, in particular the proton transfer/pumping mechanism presumably via a single proton channel, is still poorly understood. In this work we used all-atom molecular dynamics simulations and continuum electrostatics calculations to obtain atomic-level insights into the hydration and dynamics of a cbb₃ oxidase. We provide the details of the water dynamics and proton transfer pathways for both the "chemical" and "pumped" protons, and show that formation of protonic connections is strongly affected by the protonation state of key residues, namely H243, E323 and H337.

Performance of a docking/molecular dynamics protocol for virtual screening of nutlin-class inhibitors of Mdmx

A virtual screening protocol involving docking and molecular dynamics has been tested against the results of fluorescence polarization assays testing the potency of a series of compounds of the nutlin class for inhibition of the interaction between p53 and Mdmx, an interaction identified as a driver of certain cancers. The protocol uses a standard docking method (AutoDock) with a cutoff based on the AutoDock score (ADscore), followed by molecular dynamics simulation with a cutoff based on root-mean-square-deviation (RMSD) from the docked pose. An analysis of the experimental and computational results shows modest performance of ADscore alone, but dramatically improved performance when RMSD is also used.

Thermal buckling behavior of defective CNTs under pre-load: A molecular dynamics study

Current study is concentrated on the extraordinary properties of defective carbon nanotubes (CNTs). The role of vacancy defects in thermal buckling response of precompressed CNTs is explored via molecular dynamics (MD) simulations. Defective CNTs are initially compressed at a certain ratio of their critical buckling strain and then undergo a uniform temperature rise. Comprehensive study is implemented on both armchair and zigzag CNTs with different vacancy defects including monovacancy, symmetric bivacancy and asymmetric bivacancy. The results reveal that defects have a pronounced impact on the buckling behavior of CNTs; interestingly, defective CNTs under compressive pre-load show higher resistance to thermal buckling than pristine ones. In the following, the buckling response of defective CNTs is shown to be dependent on the vacancy defects, location of defects and chirality.

Effect of methyl groups on conformational properties of small ionized comb-like polyelectrolytes at the atomic level

Comb-like polycarboxylate ether (PCE) molecules with different content of methyl groups substituted on backbone and different location of methyl groups substituted on the side chains, respectively, were designed and were studied in explicit salt solutions by all-atom molecular dynamics simulations. Methyl groups substituted on the backbone of PCE have a great effect on the conformation of PCE. Stiffness of charged backbone was not only affected by the rotational freedom but also the electrostatic repulsion between the charged COO^- groups. The interaction of counterions (Na⁺) with COO^- groups for PCE3 (with part of AA substituted by MAA on the backbone) was stronger and the screen effect was great, which decided the smaller size of PCE3. The interaction between water and COO^- groups was strong regardless of the content of AA substituted by MAA on the backbone. The effect of methyl groups substituted on the different location of side chains on the conformation of PCE was less than that of methyl groups substituted on the backbone. The equilibrium sizes of the four PCE molecules with methyl groups substituted on the side chains were similar. Graphical Abstract Effect of methyl groups on conformational properties of small ionized comb-like polyelectrolytes at the atomic level.

Crystal structure and molecular dynamics studies of L-amino acid oxidase from Bothrops atrox

L-amino acid oxidases (LAAOs) are dimeric flavoproteins that catalyze the deamination of L-amino acid to α-keto acid, producing ammonia and hydrogen peroxide. In this study, we report the crystal structure and molecular dynamics simulations of LAAO from the venom of Bothrops atrox (BatroxLAAO). BatroxLAAO presents several biological and pharmacological properties with promising biomedical applications. BatroxLAAO structure contains the highly conserved structural pattern of LAAOs comprising a FAD-binding domain, substrate-binding domain and helical domain, and a dimeric arrangement that can be stabilized by zinc. Also, molecular dynamics results show an asymmetric behavior, and a direct communication between FAD- and substrate-binding domains of counterpart subunits. These findings shed light on the structural role of dimerization to catalytic mechanism of SV-LAAOs.

Joint refinement of FRET measurements using spectroscopic and computational tools

The variability of the orientation factor is a long-standing challenge in converting FRET efficiency measurements into donor-acceptor distances. We propose the use of molecular dynamics (MD) simulations to characterize orientation distributions and thus improve the accuracy of distance measurements. Here, we test this approach by comparing experimental and simulated FRET efficiencies for a model donor-acceptor pair of enhanced cyan and enhanced yellow FPs connected by a flexible linker. Several spectroscopic techniques were used to characterize FRET in solution. In addition, a series of atomistic MD simulations of a total length of 1.5 μs were carried out to calculate the distances and the orientation factor in the FRET-pair. The resulting MD-based and experimentally measured FRET efficiency histograms coincided with each other, allowing for direct comparison of distance distributions. Despite the fact that the calculated average orientation factor was close to 2/3, the application of the average κ² to the entire histogram of FRET efficiencies resulted in a substantial artificial broadening of the calculated distribution of apparent donor-acceptor distances. By combining single pair-FRET measurements with computational tools, we demonstrate that accounting for the donor and acceptor orientation heterogeneity is critical for accurate representation of the donor-acceptor distance distribution from FRET measurements.

Lid domain plasticity and lipid flexibility modulate enzyme specificity in human monoacylglycerol lipase

Human monoacylglycerol lipase (MAGL) is a membrane-interacting enzyme that generates pro-inflammatory signaling molecules. For this reason, MAGL inhibition is a promising strategy to treat pain, cancer, and neuroinflammatory diseases. MAGL can hydrolyze monoacylglycerols bearing an acyl chain of different lengths and degrees of unsaturation, cleaving primarily the endocannabinoid 2-arachidonoylglycerol. Importantly, the enzymatic binding site of MAGL is confined by a 75-amino-acid-long, flexible cap domain, named 'lid domain', which is structurally similar to that found in several other lipases. However, it is unclear how lid domain plasticity affects catalysis in MAGL. By integrating extensive molecular dynamics simulations and free-energy calculations with mutagenesis and kinetic experiments, we here define a lid-domain-mediated mechanism for substrate selection and binding in MAGL catalysis. In particular, we clarify the key role of Phe159 and Ile179, two conserved residues within the lid domain, in regulating substrate specificity in MAGL. We conclude by proposing that other structurally related lipases may share this lid-domain-mediated mechanism for substrate specificity.

Molecular mechanisms of bio-catalysis of heme extraction from hemoglobin

Red blood cell hemolysis in sickle cell disease (SCD) releases free hemoglobin. Extracellular hemoglobin and its degradation products, free heme and iron, are highly toxic due to oxidative stress induction and decrease in nitric oxide availability. We propose an approach that helps to eliminate extracellular hemoglobin toxicity in SCD by employing a bacterial protein system that evolved to extract heme from extracellular hemoglobin. NEAr heme Transporter (NEAT) domains from iron-regulated surface determinant proteins from Staphylococcus aureus specifically bind free heme as well as facilitate its extraction from hemoglobin. We demonstrate that a purified NEAT domain fused with human haptoglobin β-chain is able to remove heme from hemoglobin and reduce heme content and peroxidase activity of hemoglobin. We further use molecular dynamics (MD) simulations to resolve molecular pathway of heme transfer from hemoglobin to NEAT, and to elucidate molecular mechanism of such heme transferring process. Our study is the first of its kind, in which simulations are employed to characterize the process of heme leaving hemoglobin and subsequent rebinding with a NEAT domain. Our MD results highlight important amino acid residues that facilitate heme transfer and will guide further studies for the selection of best NEAT candidate to attenuate free hemoglobin toxicity.

A combined EPR and MD simulation study of a nitroxyl spin label with restricted internal mobility sensitive to protein dynamics

EPR studies combined with fully atomistic Molecular Dynamics (MD) simulations and an MD-EPR simulation method provide evidence for intrinsic low rotameric mobility of a nitroxyl spin label, Rn, compared to the more widely employed label MTSL (R1). Both experimental and modelling results using two structurally different sites of attachment to Myoglobin show that the EPR spectra of Rn are more sensitive to the local protein environment than that of MTSL. This study reveals the potential of using the Rn spin label as a reporter of protein motions.

Dynamic interactions between a membrane binding protein and lipids induce fluctuating diffusivity

Pleckstrin homology (PH) domains are membrane-binding lipid recognition proteins that interact with phosphatidylinositol phosphate (PIP) molecules in eukaryotic cell membranes. Diffusion of PH domains plays a critical role in biological reactions on membrane surfaces. Although diffusivity can be estimated by long-time measurements, it lacks information on the short-time diffusive nature. We reveal two diffusive properties of a PH domain bound to the surface of a PIP-containing membrane using molecular dynamics simulations. One is fractional Brownian motion, attributed to the motion of the lipids with which the PH domain interacts. The other is temporally fluctuating diffusivity; that is, the short-time diffusivity of the bound protein changes substantially with time. Moreover, the diffusivity for short-time measurements is intrinsically different from that for long-time measurements. This fluctuating diffusivity results from dynamic changes in interactions between the PH domain and PIP molecules. Our results provide evidence that the complexity of protein-lipid interactions plays a crucial role in the diffusion of proteins on biological membrane surfaces. Changes in the diffusivity of PH domains and related membrane-bound proteins may in turn contribute to the formation/dissolution of protein complexes in membranes.

The plastoquinol-plastoquinone exchange mechanism in photosystem II: insight from molecular dynamics simulations

In the photosystem II (PSII) of oxygenic photosynthetic organisms, the reaction center (RC) core mediates the light-induced electron transfer leading to water splitting and production of reduced plastoquinone molecules. The reduction of plastoquinone to plastoquinol lowers PSII affinity for the latter and leads to its release. However, little is known about the role of protein dynamics in this process. Here, molecular dynamics simulations of the complete PSII complex embedded in a lipid bilayer have been used to investigate the plastoquinol release mechanism. A distinct dynamic behavior of PSII in the presence of plastoquinol is observed which, coupled to changes in charge distribution and electrostatic interactions, causes disruption of the interactions seen in the PSII-plastoquinone complex and leads to the "squeezing out" of plastoquinol from the binding pocket. Displacement of plastoquinol closes the second water channel, recently described in a 2.9 Å resolution PSII structure (Guskov et al. in Nat Struct Mol Biol 16:334-342, 2009), allowing to rule out the proposed "alternating" mechanism of plastoquinol-plastoquinone exchange, while giving support to the "single-channel" one. The performed simulations indicated a pivotal role of D₁-Ser264 in modulating the dynamics of the plastoquinone binding pocket and plastoquinol-plastoquinone exchange via its interaction with D₁-His252 residue. The effects of the disruption of this hydrogen bond network on the PSII redox reactions were experimentally assessed in the D₁ site-directed mutant Ser264Lys.

Information-Driven, Ensemble Flexible Peptide Docking Using HADDOCK

Modeling protein-peptide interactions remains a significant challenge for docking programs due to the inherent highly flexible nature of peptides, which often adopt different conformations whether in their free or bound forms. We present here a protocol consisting of a hybrid approach, combining the most frequently found peptide conformations in complexes with representative conformations taken from molecular dynamics simulations of the free peptide. This approach intends to broaden the range of conformations sampled during docking. The resulting ensemble of conformations is used as a starting point for information-driven flexible docking with HADDOCK. We demonstrate the performance of this protocol on six cases of increasing difficulty, taken from a protein-peptide benchmark set. In each case, we use knowledge of the binding site on the receptor to drive the docking process. In the majority of cases where MD conformations are added to the starting ensemble for docking, we observe an improvement in the quality of the resulting models.

Protocols for Molecular Dynamics Simulations of RNA Nanostructures

Molecular dynamics (MD) simulations have been used as one of the main research tools to study a wide range of biological systems and bridge the gap between X-ray crystallography or NMR structures and biological mechanism. In the field of RNA nanostructures, MD simulations have been used to fix steric clashes in computationally designed RNA nanostructures, characterize the dynamics, and investigate the interaction between RNA and other biomolecules such as delivery agents and membranes.In this chapter we present examples of computational protocols for molecular dynamics simulations in explicit and implicit solvent using the Amber Molecular Dynamics Package. We also show examples of post-simulation analysis steps and briefly mention selected tools beyond the Amber package. Limitations of the methods, tools, and protocols are also discussed. Most of the examples are illustrated for a small RNA duplex (helix), but the protocols are applicable to any nucleic acid structure, subject only to the computational speed and memory limitations of the hardware available to the user.

Hybrid In Silico/In Vitro Approaches for the Identification of Functional Cholesterol-Binding Domains in Membrane Proteins

In eukaryotic cells, cholesterol is an important regulator of a broad range of membrane proteins, including receptors, transporters, and ion channels. Understanding how cholesterol interacts with membrane proteins is a difficult task because structural data of these proteins complexed with cholesterol are scarce. Here, we describe a dual approach based on in silico studies of protein-cholesterol interactions, combined with physico-chemical measurements of protein insertion into cholesterol-containing monolayers. Our algorithm is validated through careful analysis of the effect of key mutations within and outside the predicted cholesterol-binding site. Our method is illustrated by a complete analysis of cholesterol-binding to Alzheimer's β-amyloid peptide, a protein that penetrates the plasma membrane of brain cells through a cholesterol-dependent process.

Absolute binding free energies for octa-acids and guests in SAMPL5 : Evaluating binding free energies for octa-acid and guest complexes in the SAMPL5 blind challenge

As part of the SAMPL5 blind prediction challenge, we calculate the absolute binding free energies of six guest molecules to an octa-acid (OAH) and to a methylated octa-acid (OAMe). We use the double decoupling method via thermodynamic integration (TI) or Hamiltonian replica exchange in connection with the Bennett acceptance ratio (HREM-BAR). We produce the binding poses either through manual docking or by using GalaxyDock-HG, a docking software developed specifically for this study. The root mean square deviations for our most accurate predictions are 1.4 kcal mol-1 for OAH with TI and 1.9 kcal mol-1 for OAMe with HREM-BAR. Our best results for OAMe were obtained for systems with ionic concentrations corresponding to the ionic strength of the experimental solution. The most problematic system contains a halogenated guest. Our attempt to model the σ-hole of the bromine using a constrained off-site point charge, does not improve results. We use results from molecular dynamics simulations to argue that the distinct binding affinities of this guest to OAH and OAMe are due to a difference in the flexibility of the host. We believe that the results of this extensive analysis of host-guest complexes will help improve the protocol used in predicting binding affinities for larger systems, such as protein-substrate compounds.

Interaction between Pin1 and its natural product inhibitor epigallocatechin-3-gallate by spectroscopy and molecular dynamics simulations

The binding of epigallocatechin-3-gallate (EGCG) to wild type Pin1 in solution was studied by spectroscopic methods and molecular dynamics simulations in this research to explore the binding mode and inhibition mechanism. The binding constants and number of binding sites per Pin1 for EGCG were calculated through the Stern-Volmer equation. The values of binding free energy and thermodynamic parameters were calculated and indicated that hydrogen bonds, electrostatic interaction and Van der Waals interaction played the major role in the binding process. The alterations of Pin1 secondary structure in the presence of EGCG were confirmed by far-UV circular dichroism spectra. The binding model at atomic-level revealed that EGCG was bound to the Glu12, Lys13, Arg14, Met15 and Arg17 in WW domain. Furthermore, EGCG could also interact with Arg69, Asp112, Cys113 and Ser114 in PPIase domain.

Buckling analysis of defective cross-linked functionalized single- and double-walled carbon nanotubes with polyethylene chains using molecular dynamics simulations

Functionalized carbon nanotubes (CNTs) can be used for improving the mechanical properties and load transfer in nanocomposites. In this research, the buckling behavior of perfect and defective cross-linked functionalized CNTs with polyethylene (PE) chains is studied employing molecular dynamics (MD) simulations. Two different configurations with the consideration of vacancy defects, namely mapped and wrapped, are selected. According to the results, critical buckling force of cross-linked functionalized CNTs with PE chains increases as compared to pure CNTs, especially in the case of double-walled carbon nanotubes (DWCNTs). By contrast, it is demonstrated that critical strain of cross-linked functionalized CNTs decreases as compared to that of pristine CNTs. Also, it is observed that increasing the weight percentage leads to the higher increase and the decrease in critical buckling force and strain of cross-linked functionalized CNTs, respectively. Moreover, the presence of defect considerably reduces both critical buckling force and strain of cross-linked functionalized CNTs. Finally, it is shown that the critical buckling strain is more sensitive to the presence of defects as compared to critical buckling force.

Understanding protein domain-swapping using structure-based models of protein folding

In domain-swapping, two or more identical protein monomers exchange structural elements and fold into dimers or multimers whose units are structurally similar to the original monomer. Domain-swapping is of biotechnological interest because inhibiting domain-swapping can reduce disease-causing fibrillar protein aggregation. To achieve such inhibition, it is important to understand both the energetics that stabilize the domain-swapped structure and the protein dynamics that enable the swapping. Structure-based models (SBMs) encode the folded structure of the protein in their potential energy functions. SBMs have been successfully used to understand diverse aspects of monomer folding. Symmetrized SBMs model interactions between two identical protein chains using only intra-monomer interactions. Molecular dynamics simulations of such symmetrized SBMs have been used to correctly predict the domain-swapped structure and to understand the mechanism of domain-swapping. Here, we review such models and illustrate that monomer topology determines key aspects of domain-swapping. However, in some proteins, specifics of local energetic interactions modulate domain-swapping and these need to be added to the symmetrized SBMs. We then summarize some general principles of the mechanism of domain-swapping that emerge from the symmetrized SBM simulations. Finally, using our own results, we explore how symmetrized SBMs could be used to design domain-swapping in proteins.

Energetics of Glutathione Binding to Human Eukaryotic Elongation Factor 1 Gamma: Isothermal Titration Calorimetry and Molecular Dynamics Studies

The energetics of ligand binding to human eukaryotic elongation factor 1 gamma (heEF1γ) was investigated using reduced glutathione (GSH), oxidised glutathione (GSSG), glutathione sulfonate and S-hexylglutathione as ligands. The experiments were conducted using isothermal titration calorimetry, and the findings were supported using computational studies. The data show that the binding of these ligands to heEF1γ is enthalpically favourable and entropically driven (except for the binding of GSSG). The full length heEF1γ binds GSSG with lower affinity (K _d = 115 μM), with more hydrogen-bond contacts (ΔH = -73.8 kJ/mol) and unfavourable entropy (-TΔS = 51.7 kJ/mol) compared to the glutathione transferase-like N-terminus domain of heEF1γ, which did not show preference to any specific ligand. Computational free binding energy calculations from the 10 ligand poses show that GSSG and GSH consistently bind heEF1γ, and that both ligands bind at the same site with a folded bioactive conformation. This study reveals the possibility that heEF1γ is a glutathione-binding protein.

Molecular mechanism of R-bicalutamide switching from androgen receptor antagonist to agonist induced by amino acid mutations using molecular dynamics simulations and free energy calculation

R-bicalutamide, a first generation antiandrogen, was used to treat prostate cancer for decades. Although it is very effective at the beginning, resistance appears after 2-3 years of treatment. Mutation of androgen receptor (AR) is considered a main reason for drug resistance. It is reported that AR W741C, W741L, W741C_T877A, T877A, F876L, F876L_T877A and L701H mutations can convert R-bicalutamide from AR antagonist to agonist, but the switching mechanisms are not clear. In this study, molecular dynamics simulations and molecular mechanics generalized Born surface area (MM-GBSA) calculations were performed to analyze the interaction mechanisms between R-bicalutamide and wild type/mutant ARs. The results indicate that helix H12, which lies on the top of AR LBD like a cover, plays a vital role in R-bicalutamide binding. When interacting with AR, the B-ring of R-bicalutamide pushes H12 aside, distorting the coactivator binding site (AF2) resulting in the inactivation of transcription. Several residue mutations appear to enlarge the distance between the B-ring of R-bicalutamide and H12, reducing steric clash, which is conducive to a closed H12 conformation, leading to the formation of the coactivator binding site AF2 and increased transcription. Hydrogen bond and per-residue free energy decomposition analyses are also investigated to explore the interacting mechanisms, and M895 is found to be a key residue in the antagonist mechanism. The obtained molecular mechanisms will aid rational screening and design of novel AR antagonists, even to mutant AR.

Elucidating a chemical defense mechanism of Antarctic sponges: A computational study

In 2000, a novel secondary metabolite (erebusinone, Ereb) was isolated from the Antarctic sea sponge, Isodictya erinacea. The bioactivity of Ereb was investigated, and it was found to inhibit molting when fed to the arthropod species Orchomene plebs. Xanthurenic acid (XA) is a known endogenous molt regulator present in arthropods. Experimental studies have confirmed that XA inhibits molting by binding to either (or both) of two P450 enzymes (CYP315a1 or CYP314a1) that are responsible for the final two hydroxylations in the production of the molt-inducing hormone, 20-hydroxyecdysone (20E). The lack of crystal structures and biochemical assays for CYP315a1 or CYP314a1, has prevented further experimental exploration of XA and Ereb's molt inhibition mechanisms. Herein, a wide array of computational techniques - homology modeling, molecular dynamics simulations, binding site bioinformatics, flexible receptor-flexible ligand docking, and molecular mechanics-generalized Born surface area calculations - have been employed to elucidate the structure-function relationships between the aforementioned P450s and the two described small molecule inhibitors (Ereb and XA). Results indicate that Ereb likely targets CYP315a1 by interacting with a network of aromatic residues in the binding site, while XA may inhibit both CYP315a1 and CYP314a1 because of its aromatic, as well as charged nature.

Molecular docking and molecular dynamics simulation analyses of urea with ammoniated and ammoxidized lignin

Ammoniated lignin, prepared through the Mannich reaction of lignin, has more advantages as a slow-release carrier of urea molecules than ammoxidized lignin and lignin. The advantages of the ammoniated lignin include its amine groups added and its high molecular mass kept as similar as that of lignin. Three organic molecules including guaiacyl, 2-hydroxybenzylamine and 5-carbamoylpentanoic acid are monomers respectively in lignin, ammoniated lignin and ammoxidized lignin. We studied the difference between the interactions of lignin, ammoniated lignin and ammoxidized lignin with respect to urea, based on radial distribution functions (RDFs) results from molecular dynamics (MD) simulations. Glass transition temperature (T_g) and solubility parameter (δ) of ammoniated and ammoxidized lignin have been calculated by MD simulations in the constant-temperature and constant-pressure ensemble (NPT). Molecular docking results showed the interaction sites of the urea onto the ammoniated and ammoxidized lignin and three different interaction modes were identified. Root mean square deviation (RMSD) values could indicate the mobilities of the urea molecule affected by the three different interaction modes. A series of MD simulations in the constant-temperature and constant-volume ensemble (NVT) helped us to calculate the diffusivity of urea which was affected by the content of urea in ammoniated and ammoxidized lignin.

Behavior of sphingomyelin and ceramide in a tear film lipid layer model

Tear film lipid layer is a complex lipid mixture forming the outermost interface between eye and environment. Its key characteristics, such as surface tension and structural stability, are governed by the presence of polar lipids. The origin of these lipids and exact composition of the mixture are still elusive. We focus on two minor polar lipid components of the tear film lipid later: sphingomyelin and ceramide. By employing coarse grain molecular dynamics in silico simulations accompanied by Langmuir balance experiments we provide molecular-level insight into behavior of these two lipids in a tear film lipid layer model. Sphingomyelin headgroups are significantly exposed at the water-lipids boundary while ceramide molecules are incorporated between other lipids frequently interacting with nonpolar lipids. Even though these two lipids increase surface tension of the film, their molecular-level behavior suggests that they have a stabilizing effect on the tear film lipid layer.

Structural Basis for Differences in Dynamics Induced by Leu Versus Ile Residues in the CD Loop of Kir Channels

The effect of the conserved Leu/Ile site in the CD loop on the gating dynamics of Kir channels and corresponding micro-structural mechanism remains unclear. Molecular dynamics simulations were performed to investigate the structural mechanism of chicken Kir2.2. Compared to WT, the I223L mutant channel bound to PIP₂ more strongly, was activated more rapidly, and maintained the activation state more stably after PIP₂ dissociation. Cellular electrophysiology assays of mouse Kir2.1 and human Kir2.2 indicated that, consistent with simulations, the Leu residue increased the channel responses to PIP₂ through increased binding affinity and faster activation kinetics, and the deactivation kinetics decreased upon PIP₂ inhibition. The Ile residue induced the opposite responses. This difference was attributed to the distinct hydrophobic side chain symmetries of Leu and Ile; switching between these residues caused the interaction network to redistribute and offered effective conformation transduction in the Leu systems, which had more rigid and independent subunits.

Can We Rely on Computational Predictions To Correctly Identify Ligand Binding Sites on Novel Protein Drug Targets? Assessment of Binding Site Prediction Methods and a Protocol for Validation of Predicted Binding Sites

In the field of medicinal chemistry there is increasing focus on identifying key proteins whose biochemical functions can firmly be linked to serious diseases. Such proteins become targets for drug or inhibitor molecules that could treat or halt the disease through therapeutic action or by blocking the protein function respectively. The protein must be targeted at the relevant biologically active site for drug or inhibitor binding to be effective. As insufficient experimental data is available to confirm the biologically active binding site for novel protein targets, researchers often rely on computational prediction methods to identify binding sites. Presented herein is a short review on structure-based computational methods that (i) predict putative binding sites and (ii) assess the druggability of predicted binding sites on protein targets. This review briefly covers the principles upon which these methods are based, where they can be accessed and their reliability in identifying the correct binding site on a protein target. Based on this review, we believe that these methods are useful in predicting putative binding sites, but as they do not account for the dynamic nature of protein-ligand binding interactions, they cannot definitively identify the correct site from a ranked list of putative sites. To overcome this shortcoming, we strongly recommend using molecular docking to predict the most likely protein-ligand binding site(s) and mode(s), followed by molecular dynamics simulations and binding thermodynamics calculations to validate the docking results. This protocol provides a valuable platform for experimental and computational efforts to design novel drugs and inhibitors that target disease-related proteins.

Polyunsaturated chains in asymmetric lipids disorder raft mixtures and preferentially associate with α-Synuclein

Using molecular dynamics simulations, we have explored the effect of asymmetric lipids-specifically those that contain one polyunsaturated (PUFA) and one saturated fatty acid chain-on phase separation in heterogeneous membranes. These lipids are prevalent in neuronal membranes, particularly in synaptic membranes, where the Parkinson's Disease protein α-Synuclein (αS) is found. We have therefore explored the relationship between asymmetric, PUFA-containing lipids, and αS. The simulations show that asymmetric lipids partition to the liquid disordered (Ld) phase of canonical raft mixtures because of the highly disordered PUFA chain. In the case of a membrane built to mimic the lipid composition of a synaptic vesicle, the PUFA-containing asymmetric lipids completely disrupt phase separation. Because αS is positively charged, we show that it partitions with negatively charged lipids, regardless of the saturation state of the chains. Additionally, αS preferentially associates with the polyunsaturated fatty acid tails of both charged and neutral lipids. This is a consequence of those chains' ability to accommodate the void beneath the amphipathic helix. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.

The Integrin Receptor in Biologically Relevant Bilayers: Insights from Molecular Dynamics Simulations

Integrins are heterodimeric (αβ) cell surface receptors that are potential therapeutic targets for a number of diseases. Despite the existence of structural data for all parts of integrins, the structure of the complete integrin receptor is still not available. We have used available structural data to construct a model of the complete integrin receptor in complex with talin F2-F3 domain. It has been shown that the interactions of integrins with their lipid environment are crucial for their function but details of the integrin/lipid interactions remain elusive. In this study an integrin/talin complex was inserted in biologically relevant bilayers that resemble the cell plasma membrane containing zwitterionic and charged phospholipids, cholesterol and sphingolipids to study the dynamics of the integrin receptor and its effect on bilayer structure and dynamics. The results of this study demonstrate the dynamic nature of the integrin receptor and suggest that the presence of the integrin receptor alters the lipid organization between the two leaflets of the bilayer. In particular, our results suggest elevated density of cholesterol and of phosphatidylserine lipids around the integrin/talin complex and a slowing down of lipids in an annulus of ~30 Å around the protein due to interactions between the lipids and the integrin/talin F2-F3 complex. This may in part regulate the interactions of integrins with other related proteins or integrin clustering thus facilitating signal transduction across cell membranes.

Energetic and conformational features linked to the monomeric and dimeric states of bovine BLG

Bovine β-lactoglobulin (BLG) belong to the lipocalin family. This is a group of proteins involved in the binding and transporting of hydrophobic molecules. Experimental and theoretical reports have stated its complex structural behavior in solution, with coupled effects between homodimerization and ligand recognition. Nonetheless, structural evidence at the atomic level about the cause of this coupled effect has not been reported to date. To address this issue microsecond molecular dynamics (MD) simulations were combined with the molecular mechanics generalized Born surface area (MM/GBSA) approach, clustering analysis and principal component analysis (PCA), to explore the conformational complexity of BLG protein-protein self-association and palmitic acid (PLM) or dodecyl sulfate (SDS) ligand recognition in the monomeric and dimeric state. MD simulations, coupled to the MM/GBSA method, revealed that dimerization exerts contrasting effects on the ligand-binding capacity of BLG. Protein dimerization decreases PLM affinity, promoting dimer association. For SDS the dimeric state increases affinity, enhancing dimer dissociation. MD simulations based on PCA revealed that while few differences in the conformational subspace are observed between the free and bound monomer and dimer coupling for PLM, substantial changes are observed between the free and bound monomer and dimer coupling for SDS.

In-silico structural analysis of E509K mutation in LARGE and T192M mutation in Alpha Dystroglycan in the inhibition of glycosylation of Alpha Dystroglycan by LARGE

Impaired glycosylation of cellular receptor Alpha Dystroglycan (α-DG) leads to dystroglycanopathy. Glycoprotein α-DG is the receptor protein in the Dystrophin Associated Protein Complex (DAPC), a macromolecular gathering on muscle cell membrane to form a bridge between extracellular matrix (ECM) and cellular actin cytoskeleton. Proper glycosylation of α-DG is mediated by the glycosylating enzyme LARGE. Mutations either in α-DG or in LARGE lead to improper glycosylations of α-DG thereby hampering the formation of final Laminin binding form α-DG resulting in dystroglycanopathy. In our current work, we explored the structural changes associated with the presence of mutations in α-DG as well as in the enzyme LARGE. We further extended our research to understand the effect of the mutations onto protein-enzyme interactions. Moreover, since LARGE transfers the sugar moiety (glucuronic acid; GlcA) onto α-DG, we tried to analyze what effect the mutation in LARGE confers on this enzyme ligand interaction. This work for the first time addressed the molecular changes occurring in the structures α-DG, LARGE and their interactions and shed lights on the as yet poorly understood mechanism behind the dystroglycanopathy onset.

Probing the binding of Cu(2+) ions to a fragment of the Aβ(1-42) polypeptide using fluorescence spectroscopy, isothermal titration calorimetry and molecular dynamics simulations

Steady-state and time-resolved fluorescence quenching measurements supported by isothermal titration calorimetry (ITC) and molecular dynamics simulations (MD), with the NMR-derived restraints, were used to investigate the interactions of Cu(2+) ions with a fragment of the Aβ(1-42) polypeptide, Aβ(5-16) with the following sequence: Ac-Arg-His-Asp-Ser-Gly-Tyr-Glu-Val-His-His-Gln-Lys-NH2, denoted as HZ1. The studies presented in this paper, when compared with our previous results (Makowska et al., Spectrochim. Acta A 153: 451-456), show that the affinity of the peptide to metal ions is conformation-dependent. All the measurements were carried out in 20mM 2-(N-morpholino)ethanesulfonic acid (MES) buffer solution, pH6.0. The Stern-Volmer equations, along with spectroscopic observations, were used to determine the quenching and binding parameters. The obtained results unequivocally suggest that Cu(2+) ions quench the fluorescence of HZ1 only through a static quenching mechanism, in contrast to the fragment from the N-terminal part of the FPB28 protein, with sequence Ac-Tyr-Lys-Thr-Ala-Asp-Gly-Lys-Thr-Tyr- NH2 (D9) and its derivative with a single point mutation: Ac-Tyr-Lys-Thr-Ala-Asn-Gly-Lys-Thr-Tyr- NH2 (D9_M), where dynamic quenching occurred. The thermodynamic parameters (ΔITCH, ΔITCS) for the interactions between Cu(2+) ions and the HZ1 peptide were determined from the calorimetric data. The conditional thermodynamic parameters suggest that, under the experimental conditions, the formation of the Cu(2+)-HZ1 complex is both an enthalpy and entropy driven process.

Distribution and dynamics of quinones in the lipid bilayer mimicking the inner membrane of mitochondria

Quinone and its analogues (Q) constitute an important class of compounds that perform key electron transfer reactions in oxidative- and photo-phosphorylation. In the inner membrane of mitochondria, ubiquinone molecules undergo continuous redox transitions enabling electron transfer between the respiratory complexes. In such a dynamic system undergoing continuous turnover for ATP synthesis, an uninterrupted supply of substrate molecules is absolutely necessary. In the current work, we have performed atomistic molecular dynamics simulations and free energy calculations to assess the structure, dynamics, and localization of quinone and its analogues in a lipid bilayer, whose composition mimics the one in the inner mitochondrial membrane. The results show that there is a strong tendency of both quinone and quinol molecules to localize in the vicinity of the lipids' acyl groups, right under the lipid head group region. Additionally, we observe a second location in the middle of the bilayer where quinone molecules tend to stabilize. Translocation of quinone through a lipid bilayer is very fast and occurs in 10-100ns time scale, whereas the translocation of quinol is at least an order of magnitude slower. We suggest that this has important mechanistic implications given that the localization of Q ensures maximal occupancy of the Q-binding sites or Q-entry points in electron transport chain complexes, thereby maintaining an optimal turnover rate for ATP synthesis.

How to Run FAST Simulations

Molecular dynamics (MD) simulations are a powerful tool for understanding enzymes' structures and functions with full atomistic detail. These physics-based simulations model the dynamics of a protein in solution and store snapshots of its atomic coordinates at discrete time intervals. Analysis of the snapshots from these trajectories provides thermodynamic and kinetic properties such as conformational free energies, binding free energies, and transition times. Unfortunately, simulating biologically relevant timescales with brute force MD simulations requires enormous computing resources. In this chapter we detail a goal-oriented sampling algorithm, called fluctuation amplification of specific traits, that quickly generates pertinent thermodynamic and kinetic information by using an iterative series of short MD simulations to explore the vast depths of conformational space.

On the dynamics of molecular self-assembly and the structural analysis of bilayer membranes using coarse-grained molecular dynamics simulations

We present a molecular dynamics simulation study of the self-assembly of coarse-grained lipid molecules from unbiased random initial configurations. Our lipid model is based on a well-tried CG polymer model with an additional potential that mimics the hydrophobic properties of lipid tails. We find that several stages of self-organization of lipid clusters are involved in the dynamics of bilayer formation and that the resulting equilibrium structures sensitively depend on the strength of hydrophobic interactions hc of the lipid tails and on temperature T. The obtained stable lipid membranes are quantitatively analyzed with respect to their local structure and their degree of order. At equilibrium, we obtain self-stabilizing bilayer membrane structures that exhibit a bending stiffness κB and compression modulus KC comparable to experimental measurements under physiological conditions. We present a phase diagram of our lipid model which covers a sol-gel transition, a liquid (or gel-like) phase including stable bilayer structures and vesicle formation, as well as a quasi-crystalline phase. We also determine the exact conditions for temperature T and degree of hydrophobicity hc for stable bilayer formation including closed vesicles.

Exploring the membrane fusion mechanism through force-induced disassembly of HIV-1 six-helix bundle

Enveloped virus, such as HIV-1, employs membrane fusion mechanism to invade into host cell. HIV-1 gp41 ectodomain uses six-helix bundle configuration to accomplish this process. Using molecular dynamic simulations, we confirmed the stability of this six-helix bundle by showing high occupancy of hydrogen bonds and hydrophobic interactions. Key residues and interactions important for the bundle integration were characterized by force-induced unfolding simulations of six-helix bundle, exhibiting the collapse order of these groups of interactions. Moreover, our results in some way concerted with a previous theory that the formation of coiled-coil choose a route which involved cooperative interactions between the N-terminal and C-terminal helix.

Protein stability, conformational change and binding mechanism of human serum albumin upon binding of embelin and its role in disease control

Here, we present the inclusive binding mode of phytochemical embelin, an anticancer drug with human serum albumin (HSA) established under physiological condition. Also, to understand the pharmacological role of embelin molecule, here, we have studied the anti-cancer activity of embelin on human cervical cancer cell line (HeLa cell line), which revealed that embelin showed dose dependent inhibition in the growth of cancer cells and also induces 26.3% of apoptosis at an IC50 value of 29μM. Further, embelin was titrated with HSA and the fluorescence emission quenching of HSA due to the formation of the HSA-embelin complex was observed. The binding constant of this complex is 5.9±.01×10(4)M(-1) and the number of bound embelin molecules is approximately 1.0. Consequently, molecular displacement and computational docking experiments show that the embelin is binding to subdomain IB to HSA. Further evidence from microTOF-Q mass spectrometry showed an increase in mass from 66,563Da to 66,857Da observed for free HSA and HSA+embelin complex, signifying that there is robust binding of embelin with HSA. In addition, the variations of HSA secondary structural elements in presence of embelin were confirmed by circular dichroism which indicates partial unfolding of protein. Furthermore, the transmission electron micrographs established that complex formation leads to aggregation of HSA plus embelin. Molecular dynamics simulations revealed that the stability of the HSA-embelin complexes and results suggests that at around 3500ps the complex reaches equilibration state which clearly contributes to the understanding of the stability of the HSA-embelin complexes.

Asymmetric dynamics of ion channel forming proteins - Hepatitis C virus (HCV) p7 bundles

Protein p7 of hepatitis C virus (HCV) is a short 63 amino acid membrane protein which homo-oligomerises in the lipid membrane to form ion and proton conducting bundles. Two different genotypes (GTs) of p7, 1a and 5a, are used to simulate hexameric bundles of the protein embedded in a fully hydrated lipid bilayer during 400 ns molecular dynamics (MD) simulations. Whilst the bundle of GT 1a is based on a fully computational derived structure, the bundle of GT 5a is based on NMR spectroscopic data. Results of a full correlation analysis (FCA) reveal that albeit structural differences both bundles screen local minima during the simulation. The collective motion of the protein domains is asymmetric. No 'breathing-mode'-like dynamics is observed. The presence of divalent ions, such as Ca-ions affects the dynamics of especially solvent exposed parts of the protein, but leaves the asymmetric domain motion unaffected.

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.