Assembly of HIV-1 particles is initiated by the trafficking of viral Gag polyproteins from the cytoplasm to the plasma membrane (PM), where they co-localize and bud to form immature particles. Membrane targeting is mediated by the N-terminally myristoylated matrix (MA) domain of Gag and is dependent on the PM marker phosphatidylinositol-4,5-bisphosphate [PI(4,5)P2]. Recent studies revealed that PI(4,5)P2 molecules containing truncated acyl chains [tr-PI(4,5)P2] are capable of binding MA in an “extended lipid” conformation and promoting myristoyl exposure. Here we report that tr-PI(4,5)P2 molecules also readily bind to non-membrane proteins, including HIV-1 capsid (CA), which prompted us to re-examine MA-PI(4,5)P2 interactions using native lipids and membrane mimetic liposomes and bicelles. Liposome binding trends observed using a recently developed NMR approach paralleled results of flotation assays, although the affinities measured under the equilibrium conditions of NMR experiments were significantly higher. Native PI(4,5)P2 enhanced MA binding to liposomes designed to mimic non-raft-like regions of the membrane, supporting proposals that Gag binding may nucleate raft formation. Studies with bicelles revealed a subset of surface and myr-associated MA residues that are sensitive to native PI(4,5)P2, but cleft residues that interact with the 2′-acyl chains of tr-PI(4,5)P2 molecules in aqueous solution were insensitive to native PI(4,5)P2 in bicelles. Our findings call to question extended-lipid MA:membrane binding models, and instead support a model put forward from coarse-grained simulations indicating that binding is mediated predominantly by dynamic, electrostatic interactions between conserved basic residues of MA and multiple PI(4,5)P2 and phosphatidylserine molecules.

In the accepted model for human immunodeficiency virus preassembly in infected host cells, the anchoring to the intracellular leaflet of the membrane of the matrix domain (MA) that lies at the N-terminus of the viral Gag protein precursor appears to be one of the crucial steps for particle assembly. In this study, we simulated the membrane anchoring of human immunodeficiency virus-1 myristoylated MA protein using a coarse-grained representation of both the protein and the membrane. Our calculations first suggest that the myristoyl group could spontaneously release from its initial hydrophobic pocket before MA protein interacts with the lipid membrane. All-atom simulations confirmed this possibility with a related energy cost estimated to be ~5 kcal.mol(-1). The phosphatidylinositol (4,5) bisphosphate (PI(4,5)P2) head binds preferentially to the MA highly basic region as described in available NMR data, but interestingly without flipping of its 2' acyl chain into the MA protein. Moreover, MA was able to confine PI(4,5)P2 lipids all around its molecular surface after having found a stable orientation at the membrane surface. Our results suggest that this orientation is dependent on Myr anchoring and that this confinement induces a lateral segregation of PI(4,5)P2 in domains. This is consistent with a PI(4,5)P2 enrichment of the virus envelope as compared to the host cell membrane.

The molecular features that dominate the binding mode of agonists by a broadly tuned olfactory receptor are analyzed through a joint approach combining cell biology, calcium imaging, and molecular modeling. The odorant/receptor affinities, estimated through statistics accrued during molecular dynamics simulations, are in accordance with the experimental ranking. Although in many systems receptors recognize their target through a network of oriented interactions, such as H-bonding, the binding by broadly tuned olfactory receptors is dominated by non-polar terms. We show how such a feature allows chemicals belonging to different chemical families to similarly activate the receptors through compensations of interactions within the binding site.

A modeling of lignified secondary plant cell walls adapted to grass has been achieved, using molecular dynamics for time up to 180 ns, applied to systems composed of cellulose, xylan, water, and lignin. The overall model, which was 70 nm thick for a volume of 74.4 nm 3 , consisted of two crystalline cellulose layers, each being two molecules deep, separated by an interlayer space where the three other components were located. Whereas the cellulose and xylan chemistry was fixed, 18 lignin systems were considered that varied not only in guaiacyl, syringyl, and p-hydroxyphenyl composition, but also in chain length, linkage types, and the presence or absence of coumaryl units. The stabilized models showed a well-defined interface between xylan and cellulose, but some interpenetration of xylan into the lignin part of the models. A survey of the 18 models showed that their lignin component was amorphous and that their density profile was very variable and essentially model dependent. This variability was also reflected in the co-orientation of the lignin phenyl rings with respect to the cellulose surfaces, some systems showing some orientation whereas others did not. The pattern of void distribution accessible to water varied from one system to the next, but the overall void volume was systematically established at around 3%, accepting around 200 water molecules. The estimated mobility of the water molecules interacting with lignin was 1.5 times greater than that interacting with carbohydrates.

In order to model the adsorption of xylan on cellulose, we have simulated, at the atomic level, the gas phase adsorption of small xylan fragments having 5 skeletal b (1 ? 4) xylosyl residues (X 5 ), using molecular dynamics simulations. A first regime was considered, corresponding to a low surface coverage, with the adsorption of isolated X 5 in various initial orientations. In this regime, the simulation indicated that X 5 moved toward extended conformations, some of them being helical, with the possibility of either 2 1 or left-handed 3 1 helices. During the simulation, the X 5 fragments became preferentially oriented, parallel or anti parallel with respect to the cellulose chain axis. Substitution of the X 5 backbone by either GlcA and/or Araf side chains had no major influence on either the conformation or the efficiency of the interaction. However, the presence of side chains favored orientations of the X 5 backbone inclined with respect to the cellulose chain axis. In a second regime corresponding to monolayer coverage, the geometrical features of the adsorption of the xylan fragments on cellulose was roughly the same as that in the individual coverage situation. In this case, the monolayer became equilibrated at 0.14 g of xylan fragments for each g of cellulose, a figure that compared favourably with the values obtained in experimental adsorption of xylan on bacterial cellulose.

Providing a rationale that associates a chemical structure of an odorant to its induced perception has been sought for a long time. To achieve this, a detailed atomic structure of both the odorant and the olfactory receptor must be known. State-of-the-art techniques to model the 3D structure of an olfactory receptor in complex with various odorants are presented here. These range from sequence alignment with known structures to molecular dynamics simulations in a realistic environment.

We present a comparison of various computational approaches aiming at predicting the binding free energy in ligand-protein systems where the ligand is located within a highly hydrophobic cavity. The relative binding free energy between similar ligands is obtained by means of the thermodynamic integration (TI) method and compared to experimental data obtained through isothermal titration calorimetry measurements. The absolute free energy of binding prediction was obtained on a similar system (a pyrazine derivative bound to a lipocalin) by TI, potential of mean force (PMF) and also by means of the MMPBSA protocols. Although the TI protocol performs poorly either with an explicit or an implicit solvation scheme, the PMF calculation using an implicit solvation scheme leads to encouraging results, with a prediction of the binding affinity being 2 kcal mol(-1) lower than the experimental value. The use of an implicit solvation scheme appears to be well suited for the study of such hydrophobic systems, due to the lack of water molecules within the binding site.

A joint approach combining free-energy calculations and calcium-imaging assays on the broadly tuned human 1G1 olfactory receptor is reported. The free energy of binding of ten odorants was computed by means of molecular-dynamics simulations. This state function allows separating the experimentally determined eight agonists from the two non-agonists. This study constitutes a proof-of-principle for the computational deorphanization of olfactory receptors.