Trimeric HIV-1-Env Structures Define Glycan Shields from Clades A, B, and G.

The HIV-1-envelope (Env) trimer is covered by a glycan shield of ∼90 N-linked oligosaccharides, which comprises roughly half its mass and is a key component of HIV evasion from humoral immunity. To understand how antibodies can overcome the barriers imposed by the glycan shield, we crystallized fully glycosylated Env trimers from clades A, B, and G, visualizing the shield at 3.4-3.7 Å resolution. These structures reveal the HIV-1-glycan shield to comprise a network of interlocking oligosaccharides, substantially ordered by glycan crowding, that encase the protein component of Env and enable HIV-1 to avoid most antibody-mediated neutralization. The revealed features delineate a taxonomy of N-linked glycan-glycan interactions. Crowded and dispersed glycans are differently ordered, conserved, processed, and recognized by antibody. The structures, along with glycan-array binding and molecular dynamics, reveal a diversity in oligosaccharide affinity and a requirement for accommodating glycans among known broadly neutralizing antibodies that target the glycan-shielded trimer.

The inhibitory mechanism of a small protein reveals its role in antimicrobial peptide sensing.

A large number of small membrane proteins have been uncovered in bacteria, but their mechanism of action has remained mostly elusive. Here, we investigate the mechanism of a physiologically important small protein, MgrB, which represses the activity of the sensor kinase PhoQ and is widely distributed among enterobacteria. The PhoQ/PhoP two-component system is a master regulator of the bacterial virulence program and interacts with MgrB to modulate bacterial virulence, fitness, and drug resistance. A combination of cross-linking approaches with functional assays and protein dynamic simulations revealed structural rearrangements due to interactions between MgrB and PhoQ near the membrane/periplasm interface and along the transmembrane helices. These interactions induce the movement of the PhoQ catalytic domain and the repression of its activity. Without MgrB, PhoQ appears to be much less sensitive to antimicrobial peptides, including the commonly used C18G. In the presence of MgrB, C18G promotes MgrB to dissociate from PhoQ, thus activating PhoQ via derepression. Our findings reveal the inhibitory mechanism of the small protein MgrB and uncover its importance in antimicrobial peptide sensing.

Structure and supramolecular organization of the canine distemper virus attachment glycoprotein.

Canine distemper virus (CDV) is an enveloped RNA morbillivirus that triggers respiratory, enteric, and high incidence of severe neurological disorders. CDV induces devastating outbreaks in wild and endangered animals as well as in domestic dogs in countries associated with suboptimal vaccination programs. The receptor-binding tetrameric attachment (H)-protein is part of the morbilliviral cell entry machinery. Here, we present the cryo-electron microscopy (cryo-EM) structure and supramolecular organization of the tetrameric CDV H-protein ectodomain. The structure reveals that the morbilliviral H-protein is composed of three main domains: stalk, neck, and heads. The most unexpected feature was the inherent asymmetric architecture of the CDV H-tetramer being shaped by the neck, which folds into an almost 90° bent conformation with respect to the stalk. Consequently, two non-contacting receptor-binding H-head dimers, which are also tilted toward each other, are located on one side of an intertwined four helical bundle stalk domain. Positioning of the four protomer polypeptide chains within the neck domain is guided by a glycine residue (G158), which forms a hinge point exclusively in two protomer polypeptide chains. Molecular dynamics simulations validated the stability of the asymmetric structure under near physiological conditions and molecular docking showed that two receptor-binding sites are fully accessible. Thus, this spatial organization of the CDV H-tetramer would allow for concomitant protein interactions with the stalk and head domains without steric clashes. In summary, the structure of the CDV H-protein ectodomain provides new insights into the morbilliviral cell entry system and offers a blueprint for next-generation structure-based antiviral drug discovery.

High-resolution structure of the amino acid transporter AdiC reveals insights into the role of water molecules and networks in oligomerization and substrate binding.

BACKGROUND: The L-arginine/agmatine transporter AdiC is part of the arginine-dependent extreme acid resistance system of the bacterium Escherichia coli and its pathogenic varieties such as strain E. coli O157:H7. At the present time, there is a lack of knowledge concerning the role of water molecules and networks for the structure and function of AdiC, and solute transporters in general. RESULTS: The structure of the L-arginine/agmatine transporter AdiC was determined at 1.7 Å resolution by X-ray crystallography. This high resolution allowed for the identification of numerous water molecules buried in the structure. In combination with molecular dynamics (MD) simulations, we demonstrate that water molecules play an important role for stabilizing the protein and key residues, and act as placeholders for atoms of the AdiC substrates L-arginine and agmatine. MD simulations unveiled flexibility and restrained mobility of gating residues W202 and W293, respectively. Furthermore, a water-filled cavity was identified at the dimer interface of AdiC. The two monomers formed bridging interactions through water-mediated hydrogen bonds. The accessibility and presence of water molecules in this cavity was confirmed with MD simulations. Point mutations disrupting the interfacial water network validated the importance of water molecules for dimer stabilization. CONCLUSIONS: This work gives new insights into the role and importance of water molecules in the L-arginine/agmatine transporter AdiC for protein stabilization and substrate-binding site shaping and as placeholders of substrate atoms. Furthermore, and based on the observed flexibility and restrained mobility of gating residues, a mechanistic role of the gate flexibility in the transport cycle was proposed. Finally, we identified a water-filled cavity at the dimeric interface that contributes to the stability of the amino acid transporter oligomer.

Structural heterogeneity and intersubject variability of Aß in familial and sporadic Alzheimer's disease.

Point mutations in the amyloid-ß (Aß) coding region produce a combination of mutant and WT Aß isoforms that yield unique clinicopathologies in familial Alzheimer's disease (fAD) and cerebral amyloid angiopathy (fCAA) patients. Here, we report a method to investigate the structural variability of amyloid deposits found in fAD, fCAA, and sporadic AD (sAD). Using this approach, we demonstrate that mutant Aß determines WT Aß conformation through prion template-directed misfolding. Using principal component analysis of multiple structure-sensitive fluorescent amyloid-binding dyes, we assessed the conformational variability of Aß deposits in fAD, fCAA, and sAD patients. Comparing many deposits from a given patient with the overall population, we found that intrapatient variability is much lower than interpatient variability for both disease types. In a given brain, we observed one or two structurally distinct forms. When two forms coexist, they segregate between the parenchyma and cerebrovasculature, particularly in fAD patients. Compared with sAD samples, deposits from fAD patients show less intersubject variability, and little overlap exists between fAD and sAD deposits. Finally, we examined whether E22G (Arctic) or E22Q (Dutch) mutants direct the misfolding of WT Aß, leading to fAD-like plaques in vivo. Intracerebrally injecting mutant Aß40 fibrils into transgenic mice expressing only WT Aß induced the deposition of plaques with many biochemical hallmarks of fAD. Thus, mutant Aß40 prions induce a conformation of WT Aß similar to that found in fAD deposits. These findings indicate that diverse AD phenotypes likely arise from one or more initial Aß prion conformations, which kinetically dominate the spread of prions in the brain.

Antibodyomics: bioinformatics technologies for understanding B-cell immunity to HIV-1.

Numerous antibodies have been identified from HIV-1-infected donors that neutralize diverse strains of HIV-1. These antibodies may provide the basis for a B cell-mediated HIV-1 vaccine. However, it has been unclear how to elicit similar antibodies by vaccination. To address this issue, we have undertaken an informatics-based approach to understand the genetic and immunologic processes controlling the development of HIV-1-neutralizing antibodies. As DNA sequencing comprises the fastest growing database of biological information, we focused on incorporating next-generation sequencing of B-cell transcripts to determine the origin, maturation pathway, and prevalence of broadly neutralizing antibody lineages (Antibodyomics1, 2, 4, and 6). We also incorporated large-scale robotic analyses of serum neutralization to identify and quantify neutralizing antibodies in donor cohorts (Antibodyomics3). Statistical analyses furnish another layer of insight (Antibodyomics5), with physical characteristics of antibodies and their targets through molecular dynamics simulations (Antibodyomics7) and free energy perturbation analyses (Antibodyomics8) providing information-rich output. Functional interrogation of individual antibodies (Antibodyomics9) and synthetic antibody libraries (Antibodyomics10) also yields multi-dimensional data by which to understand and improve antibodies. Antibodyomics, described here, thus comprise resolution-enhancing tools, which collectively embody an information-driven discovery engine aimed toward the development of effective B cell-based vaccines.

RosENet: Improving Binding Affinity Prediction by Leveraging Molecular Mechanics Energies with an Ensemble of 3D Convolutional Neural Networks.

The worldwide increase and proliferation of drug resistant microbes, coupled with the lag in new drug development, represents a major threat to human health. In order to reduce the time and cost for exploring the chemical search space, drug discovery increasingly relies on computational biology approaches. One key step in these approaches is the need for the rapid and accurate prediction of the binding affinity for potential leads. Here, we present RosENet (Rosetta Energy Neural Networks), an ensemble of three-dimensional (3D) Convolutional Neural Networks (CNNs), which combines voxelized molecular mechanics energies and molecular descriptors for predicting the absolute binding affinity of protein-ligand complexes. By leveraging the physicochemical properties captured by the molecular force field, our ensemble model achieved a Root Mean Square Error (RMSE) of 1.24 on the PDBBind v2016 core set. We also explored some limitations and the robustness of the PDBBind data set and our approach on nearly 500 structures, including structures determined by Nuclear Magnetic Resonance and virtual screening experiments. Our study demonstrated that molecular mechanics energies can be voxelized and used to help improve the predictive power of the CNNs. In the future, our framework can be extended to features extracted from other biophysical and biochemical models, such as molecular dynamics simulations.

De novo design of covalently constrained mesosize protein scaffolds with unique tertiary structures.

The folding of natural proteins typically relies on hydrophobic packing, metal binding, or disulfide bond formation in the protein core. Alternatively, a 3D structure can be defined by incorporating a multivalent cross-linking agent, and this approach has been successfully developed for the selection of bicyclic peptides from large random-sequence libraries. By contrast, there is no general method for the de novo computational design of multicross-linked proteins with predictable and well-defined folds, including ones not found in nature. Here we use Rosetta and Tertiary Motifs (TERMs) to design small proteins that fold around multivalent cross-linkers. The hydrophobic cross-linkers stabilize the fold by macrocyclic restraints, and they also form an integral part of a small apolar core. The designed CovCore proteins were prepared by chemical synthesis, and their structures were determined by solution NMR or X-ray crystallography. These mesosized proteins, lying between conventional proteins and small peptides, are easily accessible either through biosynthetic precursors or chemical synthesis. The unique tertiary structures and ease of synthesis of CovCore proteins indicate that they should provide versatile templates for developing inhibitors of protein-protein interactions.

XFEL structures of the influenza M2 proton channel: Room temperature water networks and insights into proton conduction.

The M2 proton channel of influenza A is a drug target that is essential for the reproduction of the flu virus. It is also a model system for the study of selective, unidirectional proton transport across a membrane. Ordered water molecules arranged in "wires" inside the channel pore have been proposed to play a role in both the conduction of protons to the four gating His37 residues and the stabilization of multiple positive charges within the channel. To visualize the solvent in the pore of the channel at room temperature while minimizing the effects of radiation damage, data were collected to a resolution of 1.4 Å using an X-ray free-electron laser (XFEL) at three different pH conditions: pH 5.5, pH 6.5, and pH 8.0. Data were collected on the Inwardopen state, which is an intermediate that accumulates at high protonation of the His37 tetrad. At pH 5.5, a continuous hydrogen-bonded network of water molecules spans the vertical length of the channel, consistent with a Grotthuss mechanism model for proton transport to the His37 tetrad. This ordered solvent at pH 5.5 could act to stabilize the positive charges that build up on the gating His37 tetrad during the proton conduction cycle. The number of ordered pore waters decreases at pH 6.5 and 8.0, where the Inwardopen state is less stable. These studies provide a graphical view of the response of water to a change in charge within a restricted channel environment.

Glycosylator: a Python framework for the rapid modeling of glycans.

BACKGROUND: Carbohydrates are a class of large and diverse biomolecules, ranging from a simple monosaccharide to large multi-branching glycan structures. The covalent linkage of a carbohydrate to the nitrogen atom of an asparagine, a process referred to as N-linked glycosylation, plays an important role in the physiology of many living organisms. Most software for glycan modeling on a personal desktop computer requires knowledge of molecular dynamics to interface with specialized programs such as CHARMM or AMBER. There are a number of popular web-based tools that are available for modeling glycans (e.g., GLYCAM-WEB (http:// https://dev.glycam.org/gp/ ) or Glycosciences.db ( http://www.glycosciences.de/ )). However, these web-based tools are generally limited to a few canonical glycan conformations and do not allow the user to incorporate glycan modeling into their protein structure modeling workflow. RESULTS: Here, we present Glycosylator, a Python framework for the identification, modeling and modification of glycans in protein structure that can be used directly in a Python script through its application programming interface (API) or through its graphical user interface (GUI). The GUI provides a straightforward two-dimensional (2D) rendering of a glycoprotein that allows for a quick visual inspection of the glycosylation state of all the sequons on a protein structure. Modeled glycans can be further refined by a genetic algorithm for removing clashes and sampling alternative conformations. Glycosylator can also identify specific three-dimensional (3D) glycans on a protein structure using a library of predefined templates. CONCLUSIONS: Glycosylator was used to generate models of glycosylated protein without steric clashes. Since the molecular topology is based on the CHARMM force field, new complex sugar moieties can be generated without modifying the internals of the code. Glycosylator provides more functionality for analyzing and modeling glycans than any other available software or webserver at present. Glycosylator will be a valuable tool for the glycoinformatics and biomolecular modeling communities.

Structure and Function of the Transmembrane Domain of NsaS, an Antibiotic Sensing Histidine Kinase in Staphylococcus aureus.

NsaS is one of four intramembrane histidine kinases (HKs) in Staphylococcus aureus that mediate the pathogen's response to membrane active antimicrobials and human innate immunity. We describe the first integrative structural study of NsaS using a combination of solution state NMR spectroscopy, chemical-cross-linking, molecular modeling and dynamics. Three key structural features emerge: First, NsaS has a short N-terminal amphiphilic helix that anchors its transmembrane (TM) bundle into the inner leaflet of the membrane such that it might sense neighboring proteins or membrane deformations. Second, the transmembrane domain of NsaS is a 4-helix bundle with significant dynamics and structural deformations at the membrane interface. Third, the intracellular linker connecting the TM domain to the cytoplasmic catalytic domains of NsaS is a marginally stable helical dimer, with one state likely to be a coiled-coil. Data from chemical shifts, heteronuclear NOE, H/D exchange measurements and molecular modeling suggest that this linker might adopt different conformations during antibiotic induced signaling.

Insights from NMR Spectroscopy into the Conformational Properties of Man-9 and Its Recognition by Two HIV Binding Proteins.

Man9 GlcNAc2 (Man-9) present at the surface of HIV makes up the binding sites of several HIV-neutralizing agents and the mammalian lectin DC-SIGN, which is involved in cellular immunity and trans-infections. We describe the conformational properties of Man-9 in its free state and when bound by the HIV entry-inhibitor protein microvirin (MVN), and define the minimum epitopes of both MVN and DC-SIGN by using NMR spectroscopy. To facilitate the implementation of 3D 13 C-edited spectra to deconvolute spectral overlap and to determine the solution structure of Man-9, we developed a robust expression system for the production of 13 C,15 N-labeled glycans in mammalian cells. The studies reveal that Man-9 interacts with HIV-binding proteins through distinct epitopes and adopts diverse conformations in the bound state. In combination with molecular dynamics simulations we observed receptor-bound conformations to be sampled by Man-9 in the free state, thus suggesting a conformational selection mechanism for diverse recognition.

Critical roles of isoleucine-364 and adjacent residues in a hydrophobic gate control of phospholipid transport by the mammalian P4-ATPase ATP8A2.

P4-ATPases (flippases) translocate specific phospholipids such as phosphatidylserine from the exoplasmic leaflet of the cell membrane to the cytosolic leaflet, upholding an essential membrane asymmetry. The mechanism of flipping this giant substrate has remained an enigma. We have investigated the importance of amino acid residues in transmembrane segment M4 of mammalian P4-ATPase ATP8A2 by mutagenesis. In the related ion pumps Na(+),K(+)-ATPase and Ca(2+)-ATPase, M4 moves during the enzyme cycle, carrying along the ion bound to a glutamate. In ATP8A2, the corresponding residue is an isoleucine, which recently was found mutated in patients with cerebellar ataxia, mental retardation, and dysequilibrium syndrome. Our analyses of the lipid substrate concentration dependence of the overall and partial reactions of the enzyme cycle in mutants indicate that, during the transport across the membrane, the phosphatidylserine head group passes near isoleucine-364 (I364) and that I364 is critical to the release of the transported lipid into the cytosolic leaflet. Another M4 residue, N359, is involved in recognition of the lipid substrate on the exoplasmic side. Our functional studies are supported by structural homology modeling and molecular dynamics simulations, suggesting that I364 and adjacent hydrophobic residues function as a hydrophobic gate that separates the entry and exit sites of the lipid and directs sequential formation and annihilation of water-filled cavities, thereby enabling transport of the hydrophilic phospholipid head group in a groove outlined by the transmembrane segments M1, M2, M4, and M6, with the hydrocarbon chains following passively, still in the membrane lipid phase.

Structural Polymorphism of Alzheimer's ß-Amyloid Fibrils as Controlled by an E22 Switch: A Solid-State NMR Study.

The amyloid-ß (Aß) peptide of Alzheimer's disease (AD) forms polymorphic fibrils on the micrometer and molecular scales. Various fibril growth conditions have been identified to cause polymorphism, but the intrinsic amino acid sequence basis for this polymorphism has been unclear. Several single-site mutations in the center of the Aß sequence cause different disease phenotypes and fibrillization properties. The E22G (Arctic) mutant is found in familial AD and forms protofibrils more rapidly than wild-type Aß. Here, we use solid-state NMR spectroscopy to investigate the structure, dynamics, hydration and morphology of Arctic E22G Aß40 fibrils. (13)C, (15)N-labeled synthetic E22G Aß40 peptides are studied and compared with wild-type and Osaka E22Δ Aß40 fibrils. Under the same fibrillization conditions, Arctic Aß40 exhibits a high degree of polymorphism, showing at least four sets of NMR chemical shifts for various residues, while the Osaka and wild-type Aß40 fibrils show a single or a predominant set of chemical shifts. Thus, structural polymorphism is intrinsic to the Arctic E22G Aß40 sequence. Chemical shifts and inter-residue contacts obtained from 2D correlation spectra indicate that one of the major Arctic conformers has surprisingly high structural similarity with wild-type Aß42. (13)C-(1)H dipolar order parameters, (1)H rotating-frame spin-lattice relaxation times and water-to-protein spin diffusion experiments reveal substantial differences in the dynamics and hydration of Arctic, Osaka and wild-type Aß40 fibrils. Together, these results strongly suggest that electrostatic interactions in the center of the Aß peptide sequence play a crucial role in the three-dimensional fold of the fibrils, and by inference, fibril-induced neuronal toxicity and AD pathogenesis.

A coiled coil switch mediates cold sensing by the thermosensory protein DesK.

The thermosensor histidine kinase DesK from Bacillus subtilis senses changes in membrane fluidity initiating an adaptive response. Structural changes in DesK have been implicated in transmembrane signaling, but direct evidence is still lacking. On the basis of structure-guided mutagenesis, we now propose a mechanism of DesK-mediated signal sensing and transduction. The data indicate that stabilization/destabilization of a 2-helix coiled coil, which connects the transmembrane sensory domain of DesK to its cytosolic catalytic region, is crucial to control its signaling state. Computational modeling and simulations reveal couplings between protein, water and membrane mechanics. We propose that membrane thickening is the main driving force for signal sensing and that it acts by inducing helix stretching and rotation prompting an asymmetric kinase-competent state. Overall, the known structural changes of the sensor kinase, as well as further dynamic rearrangements that we now predict, consistently link structure determinants to activity modulation.

Palmitoylated calnexin is a key component of the ribosome-translocon complex.

A third of the human genome encodes N-glycosylated proteins. These are co-translationally translocated into the lumen/membrane of the endoplasmic reticulum (ER) where they fold and assemble before they are transported to their final destination. Here, we show that calnexin, a major ER chaperone involved in glycoprotein folding is palmitoylated and that this modification is mediated by the ER palmitoyltransferase DHHC6. This modification leads to the preferential localization of calnexin to the perinuclear rough ER, at the expense of ER tubules. Moreover, palmitoylation mediates the association of calnexin with the ribosome-translocon complex (RTC) leading to the formation of a supercomplex that recruits the actin cytoskeleton, leading to further stabilization of the assembly. When formation of the calnexin-RTC supercomplex was affected by DHHC6 silencing, mutation of calnexin palmitoylation sites or actin depolymerization, folding of glycoproteins was impaired. Our findings thus show that calnexin is a stable component of the RTC in a manner that is exquisitely dependent on its palmitoylation status. This association is essential for the chaperone to capture its client proteins as they emerge from the translocon, acquire their N-linked glycans and initiate folding.

Infrared and fluorescence assessment of the hydration status of the tryptophan gate in the influenza A M2 proton channel.

The M2 proton channel of the influenza A virus has been the subject of extensive studies because of its critical role in viral replication. As such, we now know a great deal about its mechanism of action, especially how it selects and conducts protons in an asymmetric fashion. The conductance of this channel is tuned to conduct protons at a relatively low biologically useful rate, which allows acidification of the viral interior of a virus entrapped within an endosome, but not so great as to cause toxicity to the infected host cell prior to packaging of the virus. The dynamic, structural and chemical features that give rise to this tuning are not fully understood. Herein, we use a tryptophan (Trp) analog, 5-cyanotryptophan, and various methods, including linear and nonlinear infrared spectroscopies, static and time-resolved fluorescence techniques, and molecular dynamics simulations, to site-specifically interrogate the structure and hydration dynamics of the Trp41 gate in the transmembrane domain of the M2 proton channel. Our results suggest that the Trp41 sidechain adopts the t90 rotamer, the χ2 dihedral angle of which undergoes an increase of approximately 35° upon changing the pH from 7.4 to 5.0. Furthermore, we find that Trp41 is situated in an environment lacking bulk-like water, and somewhat surprisingly, the water density and dynamics do not show a measurable difference between the high (7.4) and low (5.0) pH states. Since previous studies have shown that upon channel opening water flows into the cavity above the histidine tetrad (His37), the present finding thus provides evidence indicating that the lack of sufficient water molecules near Trp41 needed to establish a continuous hydrogen bonding network poses an additional energetic bottleneck for proton conduction.

Perturbations of the straight transmembrane α-helical structure of the amyloid precursor protein affect its processing by γ-secretase.

The amyloid precursor protein (APP) is a widely expressed type I transmembrane (TM) glycoprotein present at the neuronal synapse. The proteolytic cleavage by γ-secretase of its C-terminal fragment produces amyloid-ß (Aß) peptides of different lengths, the deposition of which is an early indicator of Alzheimer disease. At present, there is no consensus on the conformation of the APP-TM domain at the biological membrane. Although structures have been determined by NMR in detergent micelles, their conformation is markedly different. Here we show by using molecular simulations that the APP-TM region systematically prefers a straight α-helical conformation once embedded in a membrane bilayer. However, APP-TM is highly flexible, and its secondary structure is strongly influenced by the surrounding lipid environment, as when enclosed in detergent micelles. This behavior is confirmed when analyzing in silico the atomistic APP-TM population observed by residual dipolar couplings and double electron-electron resonance spectroscopy. These structural and dynamic features are critical in the proteolytic processing of APP by the γ-secretase enzyme, as suggested by a series of Gly(700) mutants. Affecting the hydration and flexibility of APP-TM, these mutants invariantly show an increase in the production of Aß38 compared with Aß40 peptides, which is reminiscent of the effect of γ-secretase modulators inhibitors.

LipidBuilder: A Framework To Build Realistic Models for Biological Membranes.

The physical and chemical characterization of biological membranes is of fundamental importance for understanding the functional role of lipid bilayers in shaping cells and organelles, steering vesicle trafficking and promoting membrane-protein signaling. Molecular dynamics simulations stand as a powerful tool to probe the properties of membranes at atomistic level. However, the biological membrane is highly complex, and closely mimicking its physiological constitution in silico is not a straightforward task. Here, we present LipidBuilder, a framework for creating and storing models of biologically relevant phospholipid species with acyl tails of heterogeneous composition. LipidBuilder also enables the assembly of these database-stored lipids into realistic bilayers featuring asymmetric distribution on layer leaflets and concentration of given membrane constituents as defined, for example, by lipidomics experiments. The ability of LipidBuilder to assemble robust membrane models was validated by simulating membranes of homogeneous lipid composition for which experimental data are available. Furthermore, taking advantage of the extensive lipid headgroup repertoire, we assembled models of membranes of heterogeneous nature as naturally found in viral (phage PRD1), bacterial (Salmonella enterica, Laurinavicius , S. ; Kakela , R. ; Somerharju , P. ; Bamford , D. H. ; Virology 2004 , 322 , 328 - 336 ) and plant (Chlorella kessleri, Rezanka , T. ; Podojil , M. ; J. Chromatogr. 1989 , 463 , 397 - 408 ) organisms. These realistic membrane models were built using a near-exact lipid composition revealed from analytical chemistry experiments. We suggest LipidBuilder as a useful tool to model biological membranes of near-biological complexity, and as a robust complement to the current efforts to characterize the biophysical properties of biological membrane using molecular simulation.

Assembly of the transmembrane domain of E. coli PhoQ histidine kinase: implications for signal transduction from molecular simulations.

The PhoQP two-component system is a signaling complex essential for bacterial virulence and cationic antimicrobial peptide resistance. PhoQ is the histidine kinase chemoreceptor of this tandem machine and assembles in a homodimer conformation spanning the bacterial inner membrane. Currently, a full understanding of the PhoQ signal transduction is hindered by the lack of a complete atomistic structure. In this study, an atomistic model of the key transmembrane (TM) domain is assembled by using molecular simulations, guided by experimental cross-linking data. The formation of a polar pocket involving Asn202 in the lumen of the tetrameric TM bundle is crucial for the assembly and solvation of the domain. Moreover, a concerted displacement of the TM helices at the periplasmic side is found to modulate a rotation at the cytoplasmic end, supporting the transduction of the chemical signal through a combination of scissoring and rotational movement of the TM helices.