Solution structure of two xenoantigens: alpha Gal-LacNAc and alpha Gal-Lewis X

Organ hyperacute rejection, a phenomenon occurring during discordant xenotransplantation, is due to the recognition of an oligosaccharide epitope by human xenoreactive natural antibodies. In addition to the alpha Gal(1-3)beta Gal(1-4)GlcNAc trisaccharide, a fucosylated structure, alpha Gal-Lewis X, has been shown to be recognized by the antibodies. Both the trisaccharide and the tetrasaccharide have been synthesized by chemical methods. A complete nuclear magnetic resonance characterization of the two compounds has been performed, including the measurements of two-dimensional nuclear Overhauser effect spectroscopy data. Molecular dynamics simulations were run for several ns in the presence of explicit water molecules. The combination of experimental and theoretical approaches revealed the effect of an additional fucose residue on the conformational behavior of the xenoantigen. This branched fucose strongly rigidifies the N-acetyllactosamine. The effect on the alpha Gal(1-3)Gal fragment is less marked. In the presence of fucose, the terminal alpha Gal residue can still adopt two different conformations, but the equilibrium populations are modified.     

HMGN dynamics and chromatin function.

Recent studies indicate that most nuclear proteins, including histone H1 and HMG are highly mobile and their interaction with chromatin is transient. These findings suggest that the structure of chromatin is dynamic and the protein composition at any particular chromatin site is not fixed. Here we discuss how the dynamic behavior of the nucleosome binding HMGN proteins affects the structure and function of chromatin. The high intranuclear mobility of HMGN insures adequate supply of protein throughout the nucleus and serves to target these proteins to their binding sites. Transient interactions of the proteins with nucleosomes destabilize the higher order chromatin, enhance the access to nucleosomal DNA, and impart flexibility to the chromatin fiber. While roaming the nucleus, the HMGN proteins encounter binding partners and form metastable multiprotein complexes, which modulate their chromatin interactions. Studies with HMGN proteins underscore the important role of protein dynamics in chromatin function.     

Rapid generation of a representative ensemble of N-glycan conformations

Glycosylated proteins are ubiquitous components of extracellular matrices and cellular surfaces where their oligosaccharide moieties are implicated in a wide range of cell-cell and cell-matrix recognition events. Glycans constitute highly flexible molecules. Only a small number of glycan X-ray structures is available for which sufficient electron density for an entire oligosaccharide chain has been observed. An unambiguous structure determination based on NMR-derived geometric constraints alone is often not possible. Time consuming computational approaches such as Monte Carlo calculations and molecular dynamics simulations have been widely used to explore the conformational space accessible to complex carbohydrates. The generation of a comprehensive data base for N-glycan fragments based on long time molecular dynamics simulations is presented. The fragments are chosen in such a way that the effects of branched N-glycan structures are taken into account. The prediction database constitutes the basis of a procedure to generate a complete set of all possible conformations for a given N-glycan. The constructed conformations are ranked according to their energy content. The resulting conformations are in reasonable agreement with experimental data. A web interface has been established (http://www.dkfz.de/spec/glydict/), which enables to input any N-glycan of interest and to receive an ensemble of generated conformations within a few minutes.     

MoDEL (Molecular Dynamics Extended Library): a database of atomistic molecular dynamics trajectories

More than 1700 trajectories of proteins representative of monomeric soluble structures in the protein data bank (PDB) have been obtained by means of state-of-the-art atomistic molecular dynamics simulations in near-physiological conditions. The trajectories and analyses are stored in a large data warehouse, which can be queried for dynamic information on proteins, including interactions. Here, we describe the project and the structure and contents of our database, and provide examples of how it can be used to describe the global flexibility properties of proteins. Basic analyses and trajectories stripped of solvent molecules at a reduced resolution level are available from our web server.     

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .     

NMR: prediction of molecular alignment from structure using the PALES software

Orientational restraints such as residual dipolar couplings promise to overcome many of the problems that traditionally limited liquid-state nuclear magnetic resonance spectroscopy. Recently, we developed methods to predict a molecular alignment tensor and thus residual dipolar couplings for a given molecular structure. This provides many new opportunities for the study of the structure and dynamics of proteins, nucleic acids, oligosaccharides and small molecules. This protocol details the use of the software PALES (Prediction of AlignmEnt from Structure) for prediction of an alignment tensor from a known three-dimensional (3D) coordinate file of a solute. The method is applicable to alignment of molecules in many neutral and charged orienting media and takes into account the molecular shape and 3D charge distribution of the molecule.     

Hydroxymethyluracil modifications enhance the flexibility and hydrophilicity of double-stranded DNA

Oxidation of a DNA thymine to 5-hydroxymethyluracil is one of several recently discovered epigenetic modifications. Here, we report the results of nanopore translocation experiments and molecular dynamics simulations that provide insight into the impact of this modification on the structure and dynamics of DNA. When transported through ultrathin solid-state nanopores, short DNA fragments containing thymine modifications were found to exhibit distinct, reproducible features in their transport characteristics that differentiate them from unmodified molecules. Molecular dynamics simulations suggest that 5-hydroxymethyluracil alters the flexibility and hydrophilicity of the DNA molecules, which may account for the differences observed in our nanopore translocation experiments. The altered physico-chemical properties of DNA produced by the thymine modifications may have implications for recognition and processing of such modifications by regulatory DNA-binding proteins.     

Binding of the estrogen receptor to DNA. The role of waters.

Molecular dynamics simulations are carried out to investigate the binding of the estrogen receptor, a member of the nuclear hormone receptor family, to specific and non-specific DNA. Two systems have been simulated, each based on the crystallographic structure of a complex of a dimer of the estrogen receptor DNA binding domain with DNA. One structure includes the dimer and a consensus segment of DNA, ds(CCAGGTCACAGTGACCTGG); the other structure includes the dimer and a nonconsensus segment of DNA, ds(CCAGAACACAGTGACCTGG). The simulations involve an atomic model of the protein-DNA complex, counterions, and a sphere of explicit water with a radius of 45 A. The molecular dynamics package NAMD was used to obtain 100 ps of dynamics for each system with complete long-range electrostatic interactions. Analysis of the simulations revealed differences in the protein-DNA interactions for consensus and nonconsensus sequences, a bending and unwinding of the DNA, a slight rearrangement of several amino acid side chains, and inclusion of water molecules at the protein-DNA interface region. Our results indicate that binding specificity and stability is conferred by a network of direct and water mediated protein-DNA hydrogen bonds. For the consensus sequence, the network involves three water molecules, residues Glu-25, Lys-28, Lys-32, Arg-33, and bases of the DNA. The binding differs for the nonconsensus DNA sequence in which case the fluctuating network of hydrogen bonds allows water molecules to enter the protein-DNA interface. We conclude that water plays a role in furnishing DNA binding specificity to nuclear hormone receptors.     

Rotational decoupling between the hydrophilic and hydrophobic regions in lipid membranes

Cells use homeostatic mechanisms to ensure an optimal composition of distinct types of lipids in cellular membranes. The hydrophilic region of biological lipid membranes is mainly composed of several types of phospholipid headgroups that interact with incoming molecules, nanoparticles, and viruses, whereas the hydrophobic region consists of a distribution of acyl chains and sterols affecting membrane fluidity/rigidity related properties and forming an environment for membrane-bound molecules such as transmembrane proteins. A fundamental open question is to what extent the motions of these regions are coupled and, consequently, how strongly the interactions of phospholipid headgroups with other molecules depend on the properties and composition of the membrane hydrophobic core. We combine advanced solid-state nuclear magnetic resonance spectroscopy with high-fidelity molecular dynamics simulations to demonstrate how the rotational dynamics of choline headgroups remain nearly unchanged (slightly faster) with incorporation of cholesterol into a phospholipid membrane, contrasting the well-known extreme slowdown of the other phospholipid segments. Notably, our results suggest a new paradigm in which phospholipid dipole headgroups interact as quasi-freely rotating flexible dipoles at the interface, independent of the properties in the hydrophobic region.     

Molecular dynamics simulations of glycoclusters and glycodendrimers

Protein-carbohydrate recognition plays a crucial role in a wide range of biological processes, required both for normal physiological functions and the onset of disease. Nature uses multivalency in carbohydrate-protein interactions as a strategy to overcome the low affinity found for singular binding of an individual saccharide epitope to a single carbohydrate recognition domain of a lectin. To mimic the complex multi-branched oligosaccharides found in glycoconjugates, which form the structural basis of multivalent carbohydrate-protein interactions, so-called glycoclusters and glycodendrimers have been designed to serve as high-affinity ligands of the respective receptor proteins. To allow a rational design of glycodendrimer-type molecules with regard to the receptor structures involved in carbohydrate recognition, a deeper knowledge of the dynamics of such molecules is desirable. Most glycodendrimers have to be considered highly flexible molecules with their conformational preferences most difficult to elucidate by experimental methods. Longtime molecular dynamics (MD) simulations with inclusion of explicit solvent molecules are suited to explore the conformational space accessible to glycodendrimers. Here, a detailed geometric and conformational analysis of 15 glycodendrimers and glycoclusters has been accomplished, which differ with regard to their core moieties, spacer characteristics and the type of terminal carbohydrate units. It is shown that the accessible conformational space depends strongly on the structural features of the core and spacer moieties and even on the type of terminating sugars. The obtained knowledge about possible spatial distributions of the sugar epitopes exposed on the investigated hyperbranched neoglycoconjugates is detailed for all examples and forms important information for the interpretation and prediction of affinity data, which can be deduced from biological testing of these multivalent neoglycoconjugates.     

Spectrins and the Golgi

Several isoforms of spectrin membrane skeleton proteins have been localized to the Golgi complex. Golgi-specific membrane skeleton proteins associate with the Golgi as a detergent-resistant cytoskeletal structure that likely undergoes a dynamic assembly process that accommodates Golgi membrane dynamics. This review discusses the potential roles for this molecule in Golgi functions. In particular, it will focus on a recently identified distant cousin to conventional erythroid spectrin variously named Syne-1, Nesprin, myne, Enaptin, MSP-300, and Ank-1. Syne-1 has the novel ability to bind to both the Golgi and the nuclear envelope, a property that raises several intriguing and novel insights into Golgi structure and function. These include (1) the facilitation of interactions between Golgi and transitional ER sites on the nuclear envelope of muscle cells, and (2) an ability to impart localized specificity to the secretory pathway within large multinucleate syncytia such as skeletal muscle fibers.     

Structure of a protein–detergent complex: the balance between detergent cohesion and binding

Despite the major interest in membrane proteins at functional, genomic, and therapeutic levels, their biochemical and structural study remains challenging, as they require, among other things, solubilization in detergent micelles. The complexity of this task derives from the dependence of membrane protein structure on their anisotropic environment, influenced by a delicate balance between many different physicochemical properties. To study such properties in a small protein–detergent complex, we used fluorescence measurements and molecular dynamics (MD) simulations on the transmembrane part of glycophorin A (GpAtm) solubilized in micelles of dihexanoylphosphatidylcholine (DHPC) detergent. Fluorescence measurements show that DHPC has limited ability to solubilize the peptide, while MD provides a possible molecular explanation for this. We observe that the detergent molecules are balanced between two different types of interactions: cohesive interactions between detergent molecules that hold the micelle together, and adhesive interactions with the peptide. While the cohesive interactions are detergent mediated, the adhesion to the peptide depends on the specific interactions between the hydrophobic parts of the detergent and the topography of the peptide dictated by the amino acids. The balance between these two parameters results in a certain frustration of the system and rather slow equilibration. These observations suggest how molecular properties of detergents could influence membrane protein stabilization and solubilization.     

Molecular dynamics simulation of the RNA complex of a double-stranded RNA-binding domain reveals dynamic features of the intermolecular interface and its hydration

The interaction between double-stranded RNA (dsRNA) and the third double-stranded domain (dsRBD) from Drosophila Staufen protein represents a paradigm to understand how the dsRBD protein family, one of the most common RNA-binding protein units, binds dsRNA. The nuclear magnetic resonance (NMR) structure of this complex and the x-ray structure of another family member revealed the stereochemical basis for recognition, but also raised new questions. Although the crystallographic studies revealed a highly ordered interface containing numerous water-mediated contacts, NMR suggested extensive residual motion at the interface. To address how interfacial motion contributes to molecular recognition in the dsRBD-dsRNA system, we conducted a 2-ns molecular dynamics simulation of the complex derived from Staufen protein and of the separate protein and RNA components. The results support the observation that a high degree of conformational flexibility is retained upon complex formation and that this involves interfacial residues that are critical for dsRBD-dsRNA binding. The structural origin of this residual flexibility is revealed by the analysis of the trajectory of motion. Individual basic side chains switch continuously from one RNA polar group to another with a residence time seldom exceeding 100 ps, while retaining favorable interaction with RNA throughout much of the simulation. Short-lived water molecules mediate some of these interactions for a large fraction of the trajectory studied here. This result indicates that water molecules are not statically associated with the interface, but continuously exchange with the bulk solvent on a 1-10-ps time scale. This work provides new insight into dsRBD-dsRNA recognition and builds upon a growing body of evidence, suggesting that short-lived dynamic interactions play important roles in protein-nucleic acid interactions.     

Molecular dynamics simulations on pars intercerebralis major peptide-C (PMP-C) reveal the role of glycosylation and disulfide bonds in its enhanced structural stability and function

Fucosylation of Thr 9 in pars intercerebralis major peptide-C (PMP-C) enhances its structural stability and functional ability as a serine protease inhibitor. In order to understand the role of disulfide bonds and glycosylation on the structure and function of PMP-C, we have carried out multiple explicit solvent molecular dynamics (MD) simulations on fucosylated and non-fucosylated forms of PMP-C, both in the presence and absence of the disulfide bonds. Our simulations revealed that there were no significant structural changes in the native disulfide bonded forms of PMP-C due to fucosylation. On the other hand, the non-fucosylated form of PMP-C without disulfide bonds showed larger deviations from the starting structure than the fucosylated form. However, the structural deviations were restricted to the terminal regions while core β-sheet retained its hydrogen bonded structure even in absence of disulfide bonds as well as fucosylation. Interestingly, fucosylation of disulfide bonded native PMP-C led to a decreased thermal flexibility in the residue stretch 29-32 which is known to interact with the active site of the target proteases. Our analysis revealed that disulfide bonds covalently connect the residue stretch 29-32 to the central β-sheet of PMP-C and using a novel network of side chain interactions and disulfide bonds fucosylation at Thr 9 is altering the flexibility of the stretch 29-32 located at a distal site. Thus, our simulations explain for the first time, how presence of disulfide bonds between conserved cysteines and fucosylation enhance the function of PMP-C as a protease inhibitor.     

Role of flexibility and polarity as determinants of the hydration of internal cavities and pockets in proteins

Molecular dynamics simulations of Staphylococcal nuclease and of 10 variants with internal polar or ionizable groups were performed to investigate systematically the molecular determinants of hydration of internal cavities and pockets in proteins. In contrast to apolar cavities in rigid carbon structures, such as nanotubes or buckeyballs, internal cavities in proteins that are large enough to house a few water molecules will most likely be dehydrated unless they contain a source of polarity. The water content in the protein interior can be modulated by the flexibility of protein elements that interact with water, which can impart positional disorder to water molecules, or bias the pattern of internal hydration that is stabilized. This might explain differences in the patterns of hydration observed in crystal structures obtained at cryogenic and room temperature conditions. The ability of molecular dynamics simulations to determine the most likely sites of water binding in internal pockets and cavities depends on its efficiency in sampling the hydration of internal sites and alternative protein and water conformations. This can be enhanced significantly by performing multiple molecular dynamics simulations as well as simulations started from different initial hydration states.     

Molecular dynamics simulations of the minor ampullate spidroin modular amino acid sequence from Parawixia bistriatra: insights into silk tertiary structure and fibre formation

Spider fibres are primarily composed of proteins that are secreted by specialised glands found in different groups of arthropods. Because of their unique mechanical characteristics, it is of great interest to understand how the influence of repetitive modules within the fibres affects the final protein structure. Because each fibre is composed of a diverse set of repeated modular sequences, the differences between fibres allow for their structural comparison and, thereby, their functional comparison. Herein, we present molecular dynamics simulations of partial sequences from minor ampullate Spidroin (MiSp) silk of the Brazilian species Parawixia bistriata. Our data show that the formation of β-sheet structures is directly related to the N-termini alignment of the modules. The N-terminal alignment gives rise to a high number of hydrogen bonds whose formation is driven by repeated alanine (Ala) sequences, which, in turn, lead to increased fibre strength. This increased fibre strength contributes significantly to the final tertiary structure of the silk.     

Nucleoporins’ exclusive amino acid sequence features regulate their transient interaction with and selectivity of cargo complexes in the nuclear pore

Nucleocytoplasmic traffic of nucleic acids and proteins across the nuclear envelop via the nuclear pore complexes (NPCs) is vital for eukaryotic cells. NPCs screen transported macromolecules based on their morphology and surface chemistry. This selective nature of the NPC-mediated traffic is essential for regulating the fundamental functions of the nucleus, such as gene regulation, protein synthesis, and mechanotransduction. Despite the fundamental role of the NPC in cell and nuclear biology, the detailed mechanisms underlying how the NPC works have remained largely unknown. The critical components of NPCs enabling their selective barrier function are the natively unfolded phenylalanine- and glycine-rich proteins called “FG-nucleoporins” (FG Nups). These intrinsically disordered proteins are tethered to the inner wall of the NPC, and together form a highly dynamic polymeric meshwork whose physicochemical conformation has been the subject of intense debate. We observed that specific sequence features (called largest positive like-charge regions, or lpLCRs), characterized by extended subsequences that only possess positively charged amino acids, significantly affect the conformation of FG Nups inside the NPC. Here we investigate how the presence of lpLCRs affects the interactions between FG Nups and their interactions with the cargo complex. We combine coarse-grained molecular dynamics simulations with time-resolved force distribution analysis to disordered proteins to explore the behavior of the system. Our results suggest that the number of charged residues in the lpLCR domain directly governs the average distance between Phe residues and the intensity of interaction between them. As a result, the number of charged residues within lpLCR determines the balance between the hydrophobic interaction and the electrostatic repulsion and governs how dense and disordered the hydrophobic network formed by FG Nups is. Moreover, changing the number of charged residues in an lpLCR domain can interfere with ultrafast and transient interactions between FG Nups and the cargo complex.     

Carbohydrates as recognition molecules in infection and immunity

Consideration of host-parasite interactions encompasses a wide range of phenomena from adhesion to epithelial surfaces to interactions with cells of the immune system. This review focuses on the role of carbohydrates as recognition molecules in these complex interactions. The abundant glycoproteins and glycolipids of cell surfaces of both prokaryotic and eukaryotic cells have the ability to exist in a variety of spatial configurations through alpha- and beta-linkages and the formation of branched structures. This ability carries with it the opportunity of acting as informational molecules greater than that possible for proteins or nucleic acids. The blood group substances are probably the best characterized of these carbohydrate containing molecules. Whilst at present a detailed understanding of the importance of these molecules in host-parasite interaction is lacking, the material covered in this discussion emphasizes the way in which carbohydrate based recognition has been shown to be involved and may provide the basis for further understanding.     

The carbohydrate-binding site in galectin-3 is preorganized to recognize a sugarlike framework of oxygens: ultra-high-resolution structures and water dynamics

The recognition of carbohydrates by proteins is a fundamental aspect of communication within and between living cells. Understanding the molecular basis of carbohydrate-protein interactions is a prerequisite for the rational design of synthetic ligands. Here we report the high- to ultra-high-resolution crystal structures of the carbohydrate recognition domain of galectin-3 (Gal3C) in the ligand-free state (1.08 ? at 100 K, 1.25 ? at 298 K) and in complex with lactose (0.86 ?) or glycerol (0.9 ?). These structures reveal striking similarities in the positions of water and carbohydrate oxygen atoms in all three states, indicating that the binding site of Gal3C is preorganized to coordinate oxygen atoms in an arrangement that is nearly optimal for the recognition of β-galactosides. Deuterium nuclear magnetic resonance (NMR) relaxation dispersion experiments and molecular dynamics simulations demonstrate that all water molecules in the lactose-binding site exchange with bulk water on a time scale of nanoseconds or shorter. Nevertheless, molecular dynamics simulations identify transient water binding at sites that agree well with those observed by crystallography, indicating that the energy landscape of the binding site is maintained in solution. All heavy atoms of glycerol are positioned like the corresponding atoms of lactose in the Gal3C complexes. However, binding of glycerol to Gal3C is insignificant in solution at room temperature, as monitored by NMR spectroscopy or isothermal titration calorimetry under conditions where lactose binding is readily detected. These observations make a case for protein cryo-crystallography as a valuable screening method in fragment-based drug discovery and further suggest that identification of water sites might inform inhibitor design.