Ordered Water Molecules Research Articles

Monte Carlo Simulation of the Local Ordering of Water Molecules. II. Spatial Correlations and Hydrogen Bonds

The Monte Carlo method is used to calculate spatial distribution functions of oxygen and hydrogen atoms within a large-size water model (33666 SPC/E water molecules) under atmospheric pressure at room temperature. The work focuses on structural interpretation of local densities of water at the distances of about 3–5 A from its molecules. The distribution of the distances between water molecules connected by chains of two or more hydrogen bonds indicates that the molecules between the first and second peaks of the radial distribution function (RDF) are mainly second and, to a lesser extent, third neighbors along the chain of bonds.

The effect of different alcohols on the Asp23-Lys28 and Asp23-Ala42 salt bridges of the most effective peptide in Alzheimer's disease: Molecular dynamics viewpoints

The beta-amyloid peptide Aβ1-42 is the most effective peptide in the process of forming plaque and creating Alzheimer's. After the separation of Aβ1-42 from APP membrane protein, the membrane surface is transmitted to the extracellular environment, which is a crowded environment. On the other hand, stability of salt bridges Asp23-Lys28 and Lys28-Ala42 is important for monomer toxicity and fibrillation formation. In this work, the effects of ethanol, propanol, butanol, pentanol, hexanol, heptanol and octanol on the Asp23-Lys28 and Lys28-Ala42 salt bridges of the Aβ1-42 have been investigated by molecular dynamics simulation. The radial distribution function of the oxygen atoms of the water around the atoms Cγ-Asp23, Nξ-Lys28 and O-Ala42 was calculated in the presence of the alcohols. The results show that the peak height of the radial distribution function around the Cγ-Asp23 atom is larger than the other two atoms. Also, the numbering of water molecules in the interval corresponding to the first peak in the radial distribution function for all atoms involved in the two salt bridges Asp23-Lys28 and Lys28-Ala42 was calculated. The results show that the obtained coordinate numbers are within the range of experimental numbers reported for water. The results also show that the order of water molecules around the O-Ala42 is lower. The results of solvent accessible surface area of Aβ1-42 show that the Lys28-Ala42 Salt bridge stability is more important for toxicity of monomer.

Structural Insight Into Conformational Changes Induced by ATP Binding in a Type III Secretion-Associated ATPase From Shigella flexneri.

Gram-negative bacteria utilize the type III secretion system (T3SS) to inject effector proteins into the host cell cytoplasm, where they subvert cellular functions and assist pathogen invasion. The conserved type III-associated ATPase is critical for the separation of chaperones from effector proteins, the unfolding of effector proteins and translocating them through the narrow channel of the secretion apparatus. However, how ATP hydrolysis is coupled to the mechanical work of the enzyme remains elusive. Herein, we present a complete description of nucleoside triphosphate binding by surface presentation antigens 47 (Spa47) from Shigella flexneri, based on crystal structures containing ATPγS, a catalytic magnesium ion and an ordered water molecule. Combining the crystal structures of Spa47-ATPγS and unliganded Spa47, we propose conformational changes in Spa47 associated with ATP binding, the binding of ATP induces a conformational change of a highly conserved luminal loop, facilitating ATP hydrolysis by the Spa47 ATPase. Additionally, we identified a specific hydrogen bond critical for ATP recognition and demonstrated that, while ATPγS is an ideal analog for probing ATP binding, AMPPNP is a poor ATP mimic. Our findings provide structural insight pertinent for inhibitor design.

Femtosecond Hydration Map of Intrinsically Disordered α-Synuclein

Protein hydration water plays a fundamentally important role in protein folding, binding, assembly, and function. Little is known about the hydration water in intrinsically disordered proteins that challenge the conventional sequence-structure-function paradigm. Here, by combining experiments and simulations, we show the existence of dynamical heterogeneity of hydration water in an intrinsically disordered presynaptic protein, namely α-synuclein, implicated in Parkinson’s disease. We took advantage of nonoccurrence of cysteine in the sequence and incorporated a number of cysteine residues at the N-terminal segment, the central amyloidogenic nonamyloid-β component (NAC) domain, and the C-terminal end of α-synuclein. We then labeled these cysteine variants using environment-sensitive thiol-active fluorophore and monitored the solvation dynamics using femtosecond time-resolved fluorescence. The site-specific femtosecond time-resolved experiments allowed us to construct the hydration map of α-synuclein. Our results show the presence of three dynamically distinct types of water: bulk, hydration, and confined water. The amyloidogenic NAC domain contains dynamically restrained water molecules that are strikingly different from the water molecules present in the other two domains. Atomistic molecular dynamics simulations revealed longer residence times for water molecules near the NAC domain and supported our experimental observations. Additionally, our simulations allowed us to decipher the molecular origin of the dynamical heterogeneity of water in α-synuclein. These simulations captured the quasi-bound water molecules within the NAC domain originating from a complex interplay between the local chain compaction and the sequence composition. Our findings from this synergistic experimental simulation approach suggest longer trapping of interfacial water molecules near the amyloidogenic hotspot that triggers the pathological conversion into amyloids via chain sequestration, chain desolvation, and entropic liberation of ordered water molecules.

Atomic Structure of Water Monolayer on Anatase TiO2(101) Surface

The molecular structure and interactions of interfacial water, i.e., the first monolayers of water adsorbed on surfaces, are critical to determine the microscopic properties of solid–liquid interfaces. In particular, interfacial water on TiO2 anatase (101) is of specific relevance to industrial, environmental, and medical applications. However, owing to the complexity and challenges in sample preparation, direct observation has been limited to individual or short-range ordered water molecules adsorbed on anatase. In this article, we report the direct observation of an unprecedented long-range ordered water monolayer on TiO2 anatase (101). The adsorption structure and inter/intramolecular interactions of the water molecules have been investigated using scanning tunneling microscopy (STM) at 5 K in conjunction with density functional theory calculations. Additionally, the effect of STM tip-induced electric field was calculated to elucidate the ground-state structure of the water monolayer. The water molecules commensurately adsorb on the anatase (101) substrate; each oxygen atom of water forms a dative bond with one of its lone pair electrons to a surface Ti5c atom. The remaining lone pair electron and OH bonds are oriented in the surface plane resulting in hydrophobic properties. The herein presented results are important to understand the complex wetting properties of the TiO2 anatase surface.

The Crystal Structure of a hCA VII Variant Provides Insights into the Molecular Determinants Responsible for Its Catalytic Behavior.

Although important progress has been achieved in understanding the catalytic mechanism of Carbonic Anhydrases, a detailed picture of all factors influencing the catalytic efficiency of the various human isoforms is still missing. In this paper we report a detailed structural study and theoretical pKa calculations on a hCA VII variant. The obtained data were compared with those already known for another thoroughly investigated cytosolic isoform, hCA II. Our structural studies show that in hCA VII the network of ordered water molecules, which connects the zinc bound solvent molecule to the proton shuttle His64, is altered compared to hCA II, causing a reduction of the catalytic efficiency. Theoretical calculations suggest that changes in solvent network are related to the difference in pKa of the proton shuttle in the two enzymes. The residue that plays a major role in determining the diverse pKa values of the proton shuttle is the one in position four, namely His for hCA II and Gly for hCA VII. This residue is located on the protein surface, outside of the active site cavity. These findings are in agreement with our previous studies that highlighted the importance of histidines on the protein surface of hCA II (among which His4) as crucial residues for the high catalytic efficiency of this isoform.

Water-Binding-Mediated Gelation/Crystallization and Thermosensitive Superchirality.

Determination of molecular structural parameters of hydrophobic cholesterol-naphthalimide conjugates for water binding capabilities as well as their moisture-sensitive supramolecular self-assembly were revealed. Water binding was a key factor in leading trace water-induced crystallization against gelation in apolar solvent. Ordered water molecules entrapped in self-assembly arrays revealed by crystal structures behave as hydrogen-bonding linkers to facilitate three-dimensional growth into crystals rather than one-dimensional gel nanofibers. Water binding was also reflected on the supramolecular chirality inversion of vesicle self-assembly in aqueous media via heating-induced dehydration. Structural parameters that favor water binding were evaluated in detail, which could help rationally design organic building units for advancing soft materials, crystal engineering, and chiral recognition.

Water‐Binding‐Mediated Gelation/Crystallization and Thermosensitive Superchirality

AbstractDetermination of molecular structural parameters of hydrophobic cholesterol–naphthalimide conjugates for water binding capabilities as well as their moisture‐sensitive supramolecular self‐assembly were revealed. Water binding was a key factor in leading trace water‐induced crystallization against gelation in apolar solvent. Ordered water molecules entrapped in self‐assembly arrays revealed by crystal structures behave as hydrogen‐bonding linkers to facilitate three‐dimensional growth into crystals rather than one‐dimensional gel nanofibers. Water binding was also reflected on the supramolecular chirality inversion of vesicle self‐assembly in aqueous media via heating‐induced dehydration. Structural parameters that favor water binding were evaluated in detail, which could help rationally design organic building units for advancing soft materials, crystal engineering, and chiral recognition.

From hit to lead: Structure-based discovery of naphthalene-1-sulfonamide derivatives as potent and selective inhibitors of fatty acid binding protein 4

Fatty acid binding protein 4 (FABP4) plays a critical role in metabolism and inflammatory processes and therefore is a potential therapeutic target for immunometabolic diseases such as diabetes and atherosclerosis. Herein, we reported the identification of naphthalene-1-sulfonamide derivatives as novel, potent and selective FABP4 inhibitors by applying a structure-based design strategy. The binding affinities of compounds 16dk, 16do and 16du to FABP4, at the molecular level, are equivalent to or even better than that of BMS309403. The X-ray crystallography complemented by the isothermal titration calorimetry studies revealed the binding mode of this series of inhibitors and the pivotal network of ordered water molecules in the binding pocket of FABP4. Moreover, compounds 16dk and 16do showed good metabolic stabilities in liver microsomes. Further extensive in vivo study demonstrated that 16dk and 16do exhibited a dramatic improvement in glucose and lipid metabolism, by decreasing fasting blood glucose and serum lipid levels, enhancing insulin sensitivity, and ameliorating hepatic steatosis in obese diabetic (db/db) mice.

X-Ray Crystal Structures of the Influenza a M2 Proton Channel Bound to Amantadine, Rimantadine, and Inhibiting Compounds

Inhibition of the M2 proton channel in the influenza A virus prevents viral replication from occurring. Two of the four FDA-approved drugs for the treatment of influenza infections, amantadine and rimantadine, target the M2 channel. However, because M2 is a membrane protein, structural studies of drug binding to the channel using x-ray crystallography have been limited because of the challenging nature of the target, with the one previously solved crystal structure limited to 3.5 Å resolution. Here, we have obtained multiple crystal structures of M2 in the presence of drugs and inhibiting compounds using lipidic cubic phase (LCP) crystallization techniques. We present the first crystal structures of rimantadine bound to M2 in both the Inwardclosed and Inwardopen conformations of the channel (2.0 Å, 2.5 Å), as well as amantadine bound to the Inwardclosed conformation (2.0 Å). At this resolution range the orientation of the bound drug is unambiguous, and the ammonium group of the adamantane drugs can be seen interacting with ordered water molecules present within the channel. Additionally, we report the binding of a dual-inhibiting compound to both the wild type channel (2.6 Å) and the drug-resistant V27A mutant (2.5 Å). The position of the bound inhibitor within the channel shifts in the presence of the V27A mutation. These structures further our understanding of drug binding and inhibition within the M2 proton channel and will help guide the design of compounds for the inhibition of drug-resistant mutants of the influenza M2 channel.

Pushing Size and Resolution Limits of Single Particle Cryo-EM at 200 keV

Technical and methodological advances in single-particle cryo-electron microscopy (cryo-EM) have expanded the technique into a resolution regime that was previously only attainable by X-ray crystallography. Although single-particle cryo-EM has proven to be a useful technique for determining the structures of biomedically relevant molecules at near-atomic resolution, nearly 98% of the structures resolved to better than 4 Å resolution have been determined using 300 keV transmission electron microscopes (TEMs). We demonstrate that it is possible to obtain cryo-EM reconstructions of sub-100 kDa macromolecular complexes to better than 3 Å resolution using a 200 keV TEM. These structures are of sufficient quality to unambiguously assign amino acid rotameric conformations and identify ordered water molecules and bound co-factors, features previously thought only to be resolvable using TEMs operating at 300 keV.

Prediction of Ordered Water Molecules in Protein Binding Sites from Molecular Dynamics Simulations: The Impact of Ligand Binding on Hydration Networks.

Water plays a major role in ligand binding and is attracting increasing attention in structure-based drug design. Water molecules can make large contributions to binding affinity by bridging protein–ligand interactions or by being displaced upon complex formation, but these phenomena are challenging to model at the molecular level. Herein, networks of ordered water molecules in protein binding sites were analyzed by clustering of molecular dynamics (MD) simulation trajectories. Locations of ordered waters (hydration sites) were first identified from simulations of high resolution crystal structures of 13 protein–ligand complexes. The MD-derived hydration sites reproduced 73% of the binding site water molecules observed in the crystal structures. If the simulations were repeated without the cocrystallized ligands, a majority (58%) of the crystal waters in the binding sites were still predicted. In addition, comparison of the hydration sites obtained from simulations carried out in the absence of ligands to those identified for the complexes revealed that the networks of ordered water molecules were preserved to a large extent, suggesting that the locations of waters in a protein–ligand interface are mainly dictated by the protein. Analysis of >1000 crystal structures showed that hydration sites bridged protein–ligand interactions in complexes with different ligands, and those with high MD-derived occupancies were more likely to correspond to experimentally observed ordered water molecules. The results demonstrate that ordered water molecules relevant for modeling of protein–ligand complexes can be identified from MD simulations. Our findings could contribute to development of improved methods for structure-based virtual screening and lead optimization.

Metallopeptidases of Toxoplasma gondii: in silico identification and gene expression.

Metallopeptidases are a family of proteins with domains that remain highly conserved throughout evolution. These hydrolases require divalent metal cation(s) to activate the water molecule in order to carry out their catalytic action on peptide bonds by nucleophilic attack. Metallopeptidases from parasitic protozoa, including Toxoplasma, are investigated because of their crucial role in parasite biology. In the present study, we screened the T. gondii database using PFAM motifs specific for metallopeptidases in association with the MEROPS peptidase Database (release 10.0). In all, 49 genes encoding proteins with metallopeptidase signatures were identified in the Toxoplasma genome. An Interpro Search enabled us to uncover their domain/motif organization, and orthologs with the highest similarity by BLAST were used for annotation. These 49 Toxoplasma metallopeptidases clustered into 15 families described in the MEROPS database. Experimental expression analysis of their genes in the tachyzoite stage revealed transcription for all genes studied. Further research on the role of these peptidases should increase our knowledge of basic Toxoplasma biology and provide opportunities to identify novel therapeutic targets. This type of study would also open a path towards the comparative biology of apicomplexans.

Изучение закономерностей ближнего упорядочения водных молекул методом Монте-Карло. II. Пространственные корреляции и водородные связи

Изучение закономерностей ближнего упорядочения водных молекул методом Монте-Карло. II. Пространственные корреляции и водородные связи

Proton movement and coupling in the POT family of peptide transporters

POT transporters represent an evolutionarily well-conserved family of proton-coupled transport systems in biology. An unusual feature of the family is their ability to couple the transport of chemically diverse ligands to an inwardly directed proton electrochemical gradient. For example, in mammals, fungi, and bacteria they are predominantly peptide transporters, whereas in plants the family has diverged to recognize nitrate, plant defense compounds, and hormones. Although recent structural and biochemical studies have identified conserved sites of proton binding, the mechanism through which transport is coupled to proton movement remains enigmatic. Here we show that different POT transporters operate through distinct proton-coupled mechanisms through changes in the extracellular gate. A high-resolution crystal structure reveals the presence of ordered water molecules within the peptide binding site. Multiscale molecular dynamics simulations confirm proton transport occurs through these waters via Grotthuss shuttling and reveal that proton binding to the extracellular side of the transporter facilitates a reorientation from an inward- to outward-facing state. Together these results demonstrate that within the POT family multiple mechanisms of proton coupling have likely evolved in conjunction with variation of the extracellular gate.

The accretion of the new ice layer on the surface of hexagonal ice crystal and the influence of the local electric field on this process.

The process of creation of a new layer of ice on the basal plane and on the prism plane of a hexagonal ice crystal is analyzed. It is demonstrated that the ordering of water molecules in the already existing crystal affects the freezing. On the basal plane, when the orientations of water molecules in the ice block are random, the arrangement of the new layer in a cubic manner is observed more frequently-approximately 1.7 times more often than in a hexagonal manner. When the water molecules in the ice block are more ordered, it results in the predominance of the oxygen atoms or the hydrogen atoms on the most outer part of the surface of the ice block. In this case, the hexagonal structure is formed more frequently when the supercooling of water exceeds 10 K. This phenomenon is explained by the influence of the oriented electric field, present as a consequence of the ordering of the dipoles of water molecules in the ice block. This field modifies the structure of solvation water (i.e., the layer of water in the immediate vicinity of the ice surface). We showed that the structure of solvation water predetermines the kind of the newly created layer of ice. This effect is temperature-dependent: when the temperature draws nearer to the melting point, the cubic structure becomes the prevailing form. The temperature at which the cubic and the hexagonal structures are formed with the same probabilities is equal to about 260 K. In the case of the prism plane, the new layer that is formed is always the hexagonal one, which is independent of the arrangement of water molecules in the ice block and is in agreement with previous literature data. For the basal plane, as well as for the prism plane, no evident dependence on the ordering of water molecules that constitute the ice block on the rate of crystallization can be observed.

Achieving better-than-3-Å resolution by single-particle cryo-EM at 200 keV

Nearly 98% of single particle cryo-EM structures resolved to better than 4 Å resolution have been determined using 300 keV transmission electron microscopes. We demonstrate that it is possible to obtain reconstructions of macromolecular complexes at a range of sizes to better than 3 Å resolution using a 200 keV transmission electron microscope. These structures are of sufficient quality to unambiguously assign amino acid rotameric conformations and identify ordered water molecules.

Catalytic Mechanism of the Hotdog-Fold Thioesterase PA1618 Revealed by X-ray Structure Determination of a Substrate-Bound Oxygen Ester Analogue Complex.

Thioesterase activity accounts for the majority of the activities in the hotdog-fold superfamily. The structures and mechanisms of catalysis for many hotdog enzymes have been elucidated by X-ray crystallography and kinetics to probe the specific substrate usage and cellular functions. However, structures of hotdog thioesterases in complexes with substrate analogues reported to date utilize ligands that either represent truncations of the substrate or include additional atoms to prevent hydrolysis. Here we present the synthesis of an isosteric and isoelectronic substrate analogue-benzoyl-OdCoA-and the X-ray crystal structure of a complex of the analogue with Pseudomonas aeruginosa hotdog thioesterase PA1618 (at 1.72 Å resolution). The complex is compared with that of the "imperfect" substrate analogue phenacyl-CoA, refined to a resolution of 1.62 Å. Kinetic and structural results are consistent with Glu64 as the catalytic residue and with the involvement of Gln49 in stabilization of the transition state. Structural comparison of the two ligand-bound structures revealed a crucial ordered water molecule coordinated in the active site of the benzoyl-OdCoA structure but not present in the phenacyl-CoA-bound structure. This suggests a general base mechanism of catalysis in which Glu64 activates the coordinated water nucleophile. Together, our findings reveal the importance of a closely similar substrate analogue to determine the true substrate binding and catalytic mechanism.

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.

Testing inhomogeneous solvation theory in structure-based ligand discovery

Binding-site water is often displaced upon ligand recognition, but is commonly neglected in structure-based ligand discovery. Inhomogeneous solvation theory (IST) has become popular for treating this effect, but it has not been tested in controlled experiments at atomic resolution. To do so, we turned to a grid-based version of this method, GIST, readily implemented in molecular docking. Whereas the term only improves docking modestly in retrospective ligand enrichment, it could be added without disrupting performance. We thus turned to prospective docking of large libraries to investigate GIST's impact on ligand discovery, geometry, and water structure in a model cavity site well-suited to exploring these terms. Although top-ranked docked molecules with and without the GIST term often overlapped, many ligands were meaningfully prioritized or deprioritized; some of these were selected for testing. Experimentally, 13/14 molecules prioritized by GIST did bind, whereas none of the molecules that it deprioritized were observed to bind. Nine crystal complexes were determined. In six, the ligand geometry corresponded to that predicted by GIST, for one of these the pose without the GIST term was wrong, and three crystallographic poses differed from both predictions. Notably, in one structure, an ordered water molecule with a high GIST displacement penalty was observed to stay in place. Inclusion of this water-displacement term can substantially improve the hit rates and ligand geometries from docking screens, although the magnitude of its effects can be small and its impact in drug binding sites merits further controlled studies.