Lipid binding in rice nonspecific lipid transfer protein-1 complexes from Oryza sativa

Nonspecific lipid transfer proteins (nsLTPs) facilitate the transfer of phospholipids, glycolipids, fatty acids and steroids between membranes, with wide-ranging binding affinities. Three crystal structures of rice nsLTP1 from Oryza sativa, complexed with myristic (MYR), palmitic (PAL) or stearic acid (STE) were determined. The overall structures of the rice nsLTP1 complexes belong to the four-helix bundle folding with a long C-terminal loop. The nsLTP1-MYR and the nsLTP1-STE complexes bind a single fatty acid while the nsLTP1-PAL complex binds two molecules of fatty acids. The C-terminal loop region is elastic in order to accommodate a diverse range of lipid molecules. The lipid molecules interact with the nsLTP1-binding cavity mainly with hydrophobic interactions. Significant conformational changes were observed in the binding cavity and the C-terminal loop of the rice nsLTP1 upon lipid binding.     

Mutagenesis study of rice nonspecific lipid transfer protein 2 reveals residues that contribute to structure and ligand binding

Plant nonspecific lipid transfer protein 2 (nsLTP2) is a small (7 kDa) protein that binds lipid-like ligands. An inner hydrophobic cavity surrounded by alpha-helices is the defining structural feature of nsLTP2. Although nsLTP2 structures have been reported earlier, the detailed mechanisms of ligand binding and lipid transfer remain unclear. In this study, we used site-directed mutagenesis to determine the role of various hydrophobic residues (L8, I15, F36, F39, Y45, Y48, and V49) in the structure, stability, ligand binding, and lipid transfer activity of rice nsLTP2. Three single mutations (L8A, F36A, and V49A) drastically alter the native tertiary structure and perturb ligand binding and lipid transfer activity. Therefore, these three residues are structurally important. The Y45A mutant, however, retains a native-like structure but has decreased lipid binding affinity and lipid transfer activity, implying that this aromatic residue is critical for these biological functions. The mutants, I15A and Y48A, exhibit quite different ligand binding affinities. Y48 is involved in planar sterol binding but not linear lysophospholipid association. As for I15A, it had the highest dehydroergosterol binding affinity in spite of the lower lipid binding and transfer abilities. Our results suggest that the long alkyl side chain of I15 would restrict the flexibility of loop I (G13-A19) for sterol entry. Finally, F39A can markedly increase the exposed hydrophobic surface to maintain its transfer efficiency despite reduced ligand binding affinity. These findings suggest that the residues forming the hydrophobic cavity play various important roles in the structure and function of rice nsLTP2.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Comparison of solution and crystal structures of maize nonspecific lipid transfer protein: A model for a potential in vivo lipid carrier protein

The three-dimensional solution structure of maize nonspecific lipid transfer protein (nsLTP) obtained by nuclear magnetic resonance (NMR) is compared to the X-ray structure. Although both structures are very similar, some local structural differences are observed in the first and the fourth helices and in several side-chain conformations. These discrepancies arise partly from intermolecular contacts in the crystal lattice. The main characteristic of nsLTP structures is the presence of an internal hydrophobic cavity whose volume was found to vary from 237 to 513 ų without major variations in the 15 solution structures. Comparison of crystal and NMR structures shows the existence of another small hollow at the periphery of the protein containing a water molecule in the X-ray structure, which could play an important structural role. A model of the complexed form of maize nsLTP by α-lysopalmitoylphosphatidylcholine was built by docking the lipid inside the protein cavity of the NMR structure. The main structural feature is a hydrogen bond found also in the X-ray structure of the complex maize nsLTP/palmitate between the hydroxyl of Tyr81 and the carbonyl of the lipid. Comparison of 12 primary sequences of nsLTPs emphasizes that all residues delineating the cavities calculated on solution and X-ray structures are conserved, which suggests that this large cavity is a common feature of all compared plant nsLTPs. Furthermore several conserved basic residues seem to be involved in the stabilization of the protein architecture. Proteins 31:160–171, 1998. © 1998 Wiley-Liss, Inc.     

Crystal structure of nonspecific lipid transfer protein from Solanum melongena

Lipids are considered to protect protein allergens from proteolysis and are generally seen to exist in a bound form. One of the well‐known plant protein families with bound lipids is non‐specific lipid transfer proteins (nsLTPs). Structure‐function relationships in the case of the members of non‐specific lipid transfer protein family are not clearly understood. As part of exploring the seed proteome, we have analyzed the proteome of a member of Solanaceae family, Solanum melongena (eggplant) and a non‐specific lipid transfer protein from S. melongena, SM80.2 was purified, crystallized and the structure was determined at 1.87 Å resolution. Overall, the tertiary structure is a cluster of α‐helices forming an internal hydrophobic cavity. Absence of conserved Tyr79, known to govern the plasticity of hydrophobic cavity, and formation of hydrogen bond between Asn79 and Asn36 further reduced the pocket size. Structural analysis of SM80.2 thus gives insight about a new hydrogen bond mediated mechanism followed in closure of the binding pocket. Extra electron densities observed at two different places on the protein surface and not in the cavity could provide interesting physiological relevance. In light of allergenic properties, probably overlapping of epitopic region and ligand binding on surface could be a main reason. This work shows first crystal structure of A‐like nsLTP with a close binding pocket and extra density on the surface suggesting a plausible intermediate state during transfer.     

The mechanism of high affinity pentasaccharide binding to antithrombin,insights from Gaussian accelerated molecular dynamics simulations

Abstract

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from ‘static’ structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a ‘non-activated’ AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. ‘Generalized correlation’ matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.

Communicated by Ramaswamy H. Sarma     

Magnetization transfer from laser-polarized xenon to protons located in the hydrophobic cavity of the wheat nonspecific lipid transfer protein

Nonspecific lipid transfer protein from wheat is studied by liquid-state NMR in the presence of xenon. The gas-protein interaction is indicated by the dependence of the protein proton chemical shifts on the xenon pressure and formally confirmed by the first observation of magnetization transfer from laser-polarized xenon to the protein protons. Twenty-six heteronuclear nOes have allowed the characterization of four interaction sites inside the wheat ns-LTP cavity. Their locations are in agreement with the variations of the chemical shifts under xenon pressure and with solvation simulations. The richness of the information obtained by the noble gas with a nuclear polarization multiplied by approximately 12,000 makes this approach based on dipolar cross-relaxation with laser-polarized xenon promising for probing protein hydrophobic pockets at ambient pressure.     

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein,retinoic acid ligand,and water molecules

Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.     

Conformational states of the glucocorticoid receptor DNA-binding domain from molecular dynamics simulations

Molecular dynamics simulations (MD) have been performed on variant crystal and NMR-derived structures of the glucocorticoid receptor DNA-binding domain (GR DBD). A loop region five residues long, the so-called D-box, exhibits significant flexibility, and transient perturbations of the tetrahedral geometry of two structurally important Cys4 zinc finger are seen, coupled to conformational changes in the D-box. In some cases, one of the Cys ligands to zinc exchanges with water, although no global distortion of the protein structure is observed. Thus, from MD simulation, dynamics of the D-box could partly be explained by solvent effects in conjunction with structural reformation of the zinc finger.     

Understanding peptide competitive inhibition of botulinum neurotoxin a binding to SV2 protein via molecular dynamics simulations

Botulinum neurotoxins (BoNTs) are known as the most toxic natural substances. Synaptic vesicle protein 2 (SV2) has been proposed to be a protein receptor for BoNT/A. Recently, two short peptides (BoNT/A‐A2 and SV2C‐A3) were designed to inhibit complex formation between the BoNT/A receptor‐binding domain (BoNT/A‐RBD) and the synaptic vesicle protein 2C luminal domain (SV2C‐LD). In this article, the two peptide complex systems are studied by molecular dynamics (MD) simulations. The structural stability analysis indicates that BoNT/A‐A2 system is more stable than SV2C‐A3 system. The conformational analysis implies that the β‐sheet in BoNT/A‐A2 system maintains its secondary structure but the two β‐strands in SV2C‐A3 system have remarkable conformational changes. Based on the calculation of hydrogen bonds, hydrophobic interactions and cation‐π interactions, it is found that the internal hydrogen bonds play crucial roles in the structural stability of the peptides. Because of the stable secondary structure, the β‐sheet in BoNT/A‐A2 system establishes effective interactions at the interface and inhibits BoNT/A‐RBD binding to SV2C‐LD. In contrast, without other β‐strands forming internal hydrogen bonds, the two isolated β‐strands in SV2C‐A3 system become the random coil. This conformational change breaks important hydrogen bonds and weakens cation‐π interaction in the interface, so the complex formation is only partially inhibited by the two β‐strands. These results are consistent with experimental studies and may be helpful in understanding the inhibition mechanisms of peptide inhibitors. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 597–608, 2015.     

The dynamic properties of the Hepatitis C Virus E2 envelope protein unraveled by molecular dynamics

Hepatitis C Virus (HCV) is one of the most persistent human viruses. Although effective therapeutic approaches have been recently discovered, their use is limited by the elevated costs. Therefore, the development of alternative/complementary strategies is an urgent need. The E2 glycoprotein, the most immunogenic HCV protein, and its variants represent natural candidates to achieve this goal. Here we report an extensive molecular dynamics (MD) analysis of the intrinsic properties of E2. Our data provide interesting clues on the global and local intrinsic dynamic features of the protein. Present MD data clearly indicate that E2 combines a flexible structure with a network of covalent bonds. Moreover, the analysis of the two most important antigenic regions of the protein provides some interesting insights into their intrinsic structural and dynamic properties. Our data indicate that a fluctuating β-hairpin represents a populated state by the region E2^412?423. Interestingly, the analysis of the epitope E2^427?446 conformation, that undergoes a remarkable rearrangement in the simulation, has significant similarities with the structure that the E2^430?442 fragment adopts in complex with a neutralizing antibody. Present data also suggest that the strict conservation of Gly436 in E2 protein of different HCV genotypes is likely dictated by structural restraints. Moreover, the analysis of the E2^412?423 flexibility provides insights into the mechanisms that some antibodies adopt to anchor Trp437 that is fully buried in E2. Finally, the present investigation suggests that MD simulations should systematically complement crystallographic studies on flexible proteins that are studied in combination with antibodies.     

Position-resolved free energy of solvation for amino acids in lipid membranes from molecular dynamics simulations

Studies of insertion and interactions of amino acids in lipid membranes are pivotal to our understanding of membrane protein structure and function. Calculating the insertion cost as a function of transmembrane helix sequence is thus an important step towards improved membrane protein prediction and eventually drug design. Here, we present position-dependent free energies of solvation for all amino acid analogs along the membrane normal. The profiles cover the entire region from bulk water to hydrophobic core, and were produced from all-atom molecular dynamics simulations. Experimental differences corresponding to mutations and costs for entire segments match experimental data well, and in addition the profiles provide the spatial resolution currently not available from experiments. Polar side-chains largely maintain their hydration and assume quite ordered conformations, which indicates the solvation cost is mainly entropic. The cost of solvating charged side-chains is not only significantly lower than for implicit solvation models, but also close to experiments, meaning these could well maintain their protonation states inside the membrane. The single notable exception to the experimental agreement is proline, which is quite expensive to introduce in vivo despite its hydrophobicity--a difference possibly explained by kinks making it harder to insert helices in the translocon.     

The role of dynamics in modulating ligand exchange in intracellular lipid binding proteins

Lipids are essential for many biological processes and crucial in the pathogenesis of several diseases. Intracellular lipid-binding proteins (iLBPs) provide mobile hydrophobic binding sites that allow hydrophobic or amphipathic lipid molecules to penetrate into and across aqueous layers. Thus iLBPs mediate the lipid transport within the cell and participate to a spectrum of tissue-specific pathways involved in lipid homeostasis. Structural studies have shown that iLBPs' binding sites are inaccessible from the bulk, implying that substrate binding should involve a conformational change able to produce a ligand entry portal. Many studies have been reported in the last two decades on iLBPs indicating that their dynamics play a pivotal role in regulating ligand binding and targeted release. The ensemble of reported data has not been reviewed until today. This review is thus intended to summarize and possibly generalize the results up to now described, providing a picture which could help to identify the missing notions necessary to improve our understanding of the role of dynamics in iLBPs' molecular recognition. Such notions would clarify the chemistry of lipid binding to iLBPs and set the basis for the development of new drugs.     

Multiple binding modes for palmitate to barley lipid transfer protein facilitated by the presence of proline 12

Molecular dynamics simulations have been used to characterise the binding of the fatty acid ligand palmitate in the barley lipid transfer protein 1 (LTP) internal cavity. Two different palmitate binding modes (1 and 2), with similar protein–ligand interaction energies, have been identified using a variety of simulation strategies. These strategies include applying experimental protein–ligand atom–atom distance restraints during the simulation, or protonating the palmitate ligand, or using the vacuum GROMOS 54B7 force‐field parameter set for the ligand during the initial stages of the simulations. In both the binding modes identified the palmitate carboxylate head group hydrogen bonds with main chain amide groups in helix A, residues 4 to 19, of the protein. In binding mode 1 the hydrogen bonds are to Lys 11, Cys 13, and Leu 14 and in binding mode 2 to Thr 15, Tyr 16, Val 17, Ser 24 and also to the OH of Thr 15. In both cases palmitate binding exploits irregularity of the intrahelical hydrogen‐bonding pattern in helix A of barley LTP due to the presence of Pro 12. Simulations of two variants of barley LTP, namely the single mutant Pro12Val and the double mutant Pro12Val Pro70Val, show that Pro 12 is required for persistent palmitate binding in the LTP cavity. Overall, the work identifies key MD simulation approaches for characterizing the details of protein–ligand interactions in complexes where NMR data provide insufficient restraints.     

A comparative study of the backbone dynamics of two closely related lipid binding proteins: Bovine heart fatty acid binding protein and porcine ileal lipid binding protein

The backbone dynamics of bovine heart fatty acid binding protein (H-FABP) and porcine ileal lipid binding protein (ILBP) were studied by 15N NMR relaxation (T1 and T2) and steady state heteronuclear 15N{1H} NOE measurements. The microdynamic parameters characterizing the backbone mobility were determined using the model-free approach. For H-FABP, the non-terminal backbone amide groups display a rather compact protein structure of low flexibility. In contrast, for ILBP an increased number of backbone amide groups display unusually high internal mobility. Furthermore, the data indicate a higher degree of conformational exchange processes in the sec-msec time range for ILBP compared to H-FABP. These results suggest significant differences in the conformational stability for these two structurally highly homologous members of the fatty acid binding protein family.     

CGDB: A database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

Structural and dynamic roles of permanent water molecules in ligand molecular recognition by chicken liver bile acid binding protein

Chicken liver bile acid binding protein (cL-BABP) crystallizes with water molecules in its binding site. To obtain insights on the role of internal water, we performed two 100 ns molecular dynamics (MD) simulations in explicit solvent for cL-BABP, as apo form and as a complex with two molecules of cholic acid, and analyzed in detail the dynamics properties of all water molecules. The diffusion coefficients of the more persistent internal water molecules are significantly different from the bulk, but similar between the two protein forms. A different number of molecules and a different organization are observed for apo- and holo-cL-BABP. Most water molecules identified in the binding site of the apo-crystal diffuse to the bulk during the simulation. In contrast, almost all the internal waters of the holo-crystal maintain the same interactions with internal sidechains and ligands, which suggests they have a relevant role in protein-ligand molecular recognition. Only in the presence of these water molecules we were able to reproduce, by a classical molecular docking approach, the structure of the complex cL-BABP::cholic acid with a low ligand root mean square deviation (RMSD) with respect to its reference positioning. Literature data reported a conserved pattern of hydrogen bonds between a single water molecule and three amino acid residues of the binding site in a series of crystallized FABP. In cL-BABP, the interactions between this conserved water molecule and the three residues are present in the crystal of both apo- and holo-cL-BABP but are lost immediately after the start of molecular dynamics. Copyright (c) 2008 John Wiley & Sons, Ltd.