Spatial distribution functions as a tool in the analysis of ribonucleic acids hydration--molecular dynamics studies

Structural descriptions of ligand interactions to RNA quadruplexes folded from the non-coding region of pseudorabies virus

To rationalise the binding of specific ligands to RNA-quadruplex we investigated several naphthalene diimide ligands that interact with the non-coding region of Pseudorabies virus (PRV). Herein we report on the x-ray structure of the naphthalene diimide ND11 with an RNA G-quadruplex putative forming sequence from rPRV. Consistent with previously observed rPRV sequence it assembles into a bimolecular RNA G-quadruplex consisting of a pair of two tetrads stacked 3' to 5'. We observe that ND11 interacts by binding on both the externally available 5' and 3' quartets. The CUC (loop 1) is structurally altered to enhance the 5' mode of interaction. These loop residues are shifted significantly to generate a new ligand binding pocket whereas the terminal A14 residue is lifted away from the RNA G-quadruplex tetrad plane to be restacked above the bound ND11 ligand NDI core. CD analysis of this family of NDI ligands shows consistency in the spectra between the different ligands in the presence of the rPRV RNA G-quadruplex motif, reflecting a common folded topology and mode of ligand interaction. FRET melt assay confirms the strong stabilising properties of the tetrasubstituted NDI compounds and the contributions length of the substituted groups have on melt temperatures.

Corner Engineering: Tailoring Enzymes for Enhanced Resistance and Thermostability in Deep Eutectic Solvents

Deep eutectic solvents (DESs), heralded for their synthesis simplicity, economic viability, and reduced volatility and flammability, have found increasing application in biocatalysis. However, challenges persist due to a frequent diminution in enzyme activity and stability. Herein, we developed a general protein engineering strategy, termed corner engineering, to acquire DES-resistant and thermostable enzymes via precise tailoring of the transition region in enzyme structure. Employing Bacillus subtilis lipase A (BSLA) as a model, we delineated the engineering process, yielding five multi-DESs resistant variants with highly improved thermostability, such as K88E/N89 K exhibited up to a 10.0-fold catalytic efficiency (kcat /KM ) increase in 30 % (v/v) choline chloride (ChCl): acetamide and 4.1-fold in 95 % (v/v) ChCl: ethylene glycol accompanying 6.7-fold thermal resistance improvement than wild type at ≈50 °C. The generality of the optimized approach was validated by two extra industrial enzymes, endo-β-1,4-glucanase PvCel5A (used for biofuel production) and esterase Bs2Est (used for plastics degradation). The molecular investigations revealed that increased water molecules at substrate binding cleft and finetuned helix formation at the corner region are two dominant determinants governing elevated resistance and thermostability. This study, coupling corner engineering with obtained molecular insights, illuminates enzyme-DES interaction patterns and fosters the rational design of more DES-resistant and thermostable enzymes in biocatalysis and biotransformation.

The influence of the hydrogen-bond network on the structure and dynamics of the RAPRKKG heptapeptide and its mutants

The structural behaviour of the RAPRKKG heptapeptide after individual or multiple mutations was inspected through molecular dynamics simulation. The nature of the mutations provided information on the flexibility of the heptapeptide and on how water molecules establish hydrogen bonds with it. The structural behaviour of the wild-type and the mutated structures were measured through the analysis of protein‒protein and protein‒solvent hydrogen bonds. The conformational behaviours of the different structures were analysed through free energy landscape analysis. The flexibility characteristics of the mutants seem to depend on the reorganization of water molecules and their static or dynamic behaviour around amino acid side chains.

Liquid Structures and Diffusion Dynamics of Alkyl-Pyrene Liquids Studied by Molecular Dynamics Simulations

Functional molecular liquids (FMLs) based on alkylated π-conjugated molecules have attracted attention as solvent-free and nonvolatile liquid materials with prominent optoelectronic features. Recently, novel FML compounds containing pyrene as the functional core were synthesized, and their rheological and photochemical properties were investigated. Although the molecules differ only in the number of alkyl chain substituents and their substitution positions, their viscosity coefficients are largely different beyond the Stokes-Einstein relation on the assumption of identical microscopic friction, indicating that local microscopic molecular interactions are crucial for the macroscopic rheological properties. Here, we report a theoretical study on the rheological properties of the alkyl-pyrene liquids by means of atomistic molecular dynamics (MD) simulations. We performed long-time MD simulations for tens of microseconds to obtain ample statistical samples of the alkyl-pyrene liquids and analyzed their liquid structures and diffusion dynamics based on spatiotemporal correlation functions. We found the formation of characteristic local liquid structures of π-π stacking of the pyrene moieties and locally anisotropic and anomalous diffusion dynamics, which remarkably vary depending on the alkyl substituent patterns. The present results provide an atomistic insight into the macroscopic rheological properties of alkyl-π FMLs and molecular design strategy for them.

Investigations of CO2 Capture from Gas Mixtures Using Porous Liquids

Porous liquids, a new porous material with fluidity, can be applied in numerous fields, such as gas storage and/or separation. In this work, the separation of binary gas mixtures CO₂/N₂ and CO₂/CH₄ with porous liquids was examined by molecular dynamics (MD) simulations. The pure gas adsorption capacity was analyzed with different concentrations of porous liquids. The dependence of the separation effect of a gas mixture on the total pressure and temperature was investigated. Meanwhile, for both CO₂/N₂ and CO₂/CH₄ systems, the adsorption and separation effects of porous liquids with a cage:solvent ratio of 1:12 are better than those of 1:91 and 1:170. The results of the spatial distribution function and/or trajectories indicated that porous liquids prefer CO₂, leading to the location of CO₂ in the channels formed in porous liquids. However, N₂ and CH₄ are hardly adsorbed into the bulk. The diffusion of gas molecules follows the order of CO₂ > N₂ (for CO₂/N₂) and CH₄ > CO₂ (for CO₂/CH₄) in the bulk and N₂ > CO₂ (for CO₂/N₂) and CH₄ > CO₂ (for CO₂/CH₄) at the interface of porous liquids. Upon increasing the concentrations of porous liquids, the working capacities of CO₂ show small decreases in CO₂/N₂ and CO₂/CH₄ systems, but the sorbent selection parameters are higher in pressure- and temperature-swing adsorption processes. The porous liquid with a cage:solvent ratio of 1:12 is more suitable for the separation of CO₂/N₂ and CO₂/CH₄ systems than ratios of 1:91 and 1:170.

Investigation of bound and unbound phosphoserine phosphatase conformations through elastic network models and molecular dynamics simulations

The human phosphoserine phosphatase (hPSP) catalyses the last step in the biosynthesis of L-serine. It involves conformational changes of the enzyme lid once the substrate, phosphoserine (PSer), is bound in the active site. Here, Elastic Network Model (ENM) is applied to the crystal structure of hPSP to probe the transition between open and closed conformations of hPSP. Molecular Dynamics (MD) simulations are carried out on several PSer-hPSP systems to characterise the intermolecular interactions and their effect on the dynamics of the enzyme lid. Systems involving either Ca⁺⁺ or Mg⁺⁺ are considered. The first ENM normal mode shows that an open-closed transition can be explained from a simple description of the enzyme in terms of harmonic potentials. Principal Component Analyses applied to the MD trajectories also highlight a trend for a closing/opening motion. Different PSer orientations inside the enzyme cavity are identified, i.e. either the carboxylate, the phosphate group of PSer, or both, are oriented towards the cation. The interaction patterns are analysed in terms of hydrogen bonds, electrostatics, and bond critical points of the electron density distributions. The latter approach yields a global description of the bonding intermolecular interactions. The PSer orientation determines the content of the cation coordination shell and the mobility of the substrate, while Lys158 and Thr182, involved in the reaction mechanism, are always in interaction with the substrate. Closed enzyme conformations involve Met52-Gln204, Arg49-Glu29, and Arg50-Glu29 interactions. Met52, as well as Arg49 and Arg50, also stabilize PSer inside the cavity. Communicated by Ramaswamy H. Sarma.

Role of Mono- and Divalent Ions in Peptide Glu-Asp-Arg-DNA Interaction

The interaction of the regulatory biologically active peptide Glu-Asp-Arg (EDR) with DNA is considered by spectral, NMR, viscosimetry, and molecular dynamics methods. It was shown that EDR can partly penetrate into the major groove of DNA and affect the base atoms, mainly the N7 and O6 of guanine. It was observed that Mg²⁺ ions can promote DNA-EDR interaction due to their effective screening of the negatively charged phosphate groups of DNA. This action of Mg²⁺ remains in salted solution as well.

Representation of DNA environment: Spiral staircase distribution function

In the present study, we investigated the local structure of DNA and its environment using a new visualization technique. The spiral staircase distribution function (SSDF) is determined as two-dimensional density distribution of atoms of water and ligands in local reference frames linked with each base pair of poly-DNA molecule, either GC or AT. This property of SSDF provides opportunity to study sequence-specific binding of ions, peptides, and other agents derived from a molecular dynamics computer simulation. The spatial structure of double-stranded DNA environment in water solution containing either Mg²⁺ or Na⁺ ions was investigated using of SSDF. The distributions of ions around GC and AT base pairs are shown separately. It is observed that Mg²⁺ ions interact with DNA atoms by means of the layer of water molecules and penetrate into the major groove only. Na⁺ ions have a direct contact with DNA atoms and penetrate both into the major and minor grooves of DNA. © 2018 Wiley Periodicals, Inc.

Critical hydrogen bond formation for activation of the angiotensin II type 1 receptor

G protein-coupled receptors contain selectively important residues that play central roles in the conformational changes that occur during receptor activation. Asparagine 111 (N111(3.35)) is such a residue within the angiotensin II type 1 (AT(1)) receptor. Substitution of N111(3.35) for glycine leads to a constitutively active receptor, whereas substitution for tryptophan leads to an inactivable receptor. Here, we analyzed the AT(1) receptor and two mutants (N111G and N111W) by molecular dynamics simulations, which revealed a novel molecular switch involving the strictly conserved residue D74(2.50). Indeed, D74(2.50) forms a stable hydrogen bond (H-bond) with the residue in position 111(3.35) in the wild-type and the inactivable receptor. However, in the constitutively active mutant N111G-AT(1) receptor, residue D74 is reoriented to form a new H-bond with another strictly conserved residue, N46(1.50). When expressed in HEK293 cells, the mutant N46G-AT(1) receptor was poorly activable, although it retained a high binding affinity. Interestingly, the mutant N46G/N111G-AT(1) receptor was also inactivable. Molecular dynamics simulations also revealed the presence of a cluster of hydrophobic residues from transmembrane domains 2, 3, and 7 that appears to stabilize the inactive form of the receptor. Whereas this hydrophobic cluster and the H-bond between D74(2.50) and W111(3.35) are more stable in the inactivable N111W-AT(1) receptor, the mutant N111W/F77A-AT(1) receptor, designed to weaken the hydrophobic core, showed significant agonist-induced signaling. These results support the potential for the formation of an H-bond between residues D74(2.50) and N46(1.50) in the activation of the AT(1) receptor.

Conformational study of bovine lactoferricin in membrane-micking conditions by molecular dynamics simulation and circular dichroism

Lactoferricins are potent antimicrobial peptides released by pepsin cleavage of Lactoferrins. Bovine Lactoferricin (LfcinB) has higher activity than the intact bovine Lactoferrin, and is the most active among the other Lactoferricins of human, murine and caprine origin. In the intact protein the fragment corresponding to LfcinB is in an helical conformation, while in water LfcinB adopts an amphipathic β-hairpin structure. However, whether any of these structural motifs is the antibacterial active conformation, i.e., the one interacting with bacterial membrane components, remains to be seen. Here we present Circular Dichroism (CD) spectra and Molecular Dynamics (MD) simulations indicating that in membrane-mimicking solvents the LfcinB adopts an amphipathic β-hairpin structure similar to that observed in water, but differing in the dynamic behavior of the side-chains of the two tryptophan residues. In the membrane-mimicking solvent these side-chains show a high propensity to point towards the hydrophobic environment, rather than being in the hydrophobic core as seen in water, while the backbone preserves the hairpin conformation as found in water. These results suggest that the tryptophans might act as anchors pulling the stable, solvent-invariant hairpin structure into the membrane.

Structural and dynamic effects of cholesterol at preferred sites of interaction with rhodopsin identified from microsecond length molecular dynamics simulations

An unresolved question about GPCR function is the role of membrane components in receptor stability and activation. In particular, cholesterol is known to affect the function of membrane proteins, but the details of its effect on GPCRs are still elusive. Here, we describe how cholesterol modulates the behavior of the TM1-TM2-TM7-helix 8(H8) functional network that comprises the highly conserved NPxxY(x)(5,6)F motif, through specific interactions with the receptor. The inferences are based on the analysis of microsecond length molecular dynamics (MD) simulations of rhodopsin in an explicit membrane environment. Three regions on the rhodopsin exhibit the highest cholesterol density throughout the trajectory: the extracellular end of TM7, a location resembling the high-density sterol area from the electron microscopy data; the intracellular parts of TM1, TM2, and TM4, a region suggested as the cholesterol binding site in the recent X-ray crystallography data on beta(2)-adrenergic GPCR; and the intracellular ends of TM2-TM3, a location that was categorized as the high cholesterol density area in multiple independent 100 ns MD simulations of the same system. We found that cholesterol primarily affects specific local perturbations of the helical TM domains such as the kinks in TM1, TM2, and TM7. These local distortions, in turn, relate to rigid-body motions of the TMs in the TM1-TM2-TM7-H8 bundle. The specificity of the effects stems from the nonuniform distribution of cholesterol around the protein. Through correlation analysis we connect local effects of cholesterol on structural perturbations with a regulatory role of cholesterol in the structural rearrangements involved in GPCR function.

RNA solvation: a molecular dynamics simulation perspective

With the availability of accurate methods to treat the electrostatic long-range interactions, molecular dynamics simulations have resulted in refined dynamical models of the structure of the hydration shell around RNA motifs. The models reviewed here range from basic Watson-Crick to more specific noncanonical base pairs, from "simple" double helices to RNA molecules displaying more complex tertiary folds, and from DNA/RNA hybrid double helices to RNA hybrids formed with a chemically modified strand.

The 1.19 A X-ray structure of 2'-O-Me(CGCGCG)(2) duplex shows dehydrated RNA with 2-methyl-2,4-pentanediol in the minor groove

The crystal and molecular structure of 2'-O-Me(CGCGCG)(2) has been determined at 1.19 A resolution, at 100 K, using synchrotron radiation. The structure in space group P3(2)12 is a half-turn right-handed helix that includes two 2-methyl-2,4-pentanediol (MPD) molecules bound in the minor groove. The structure deviates from A-form RNA. The duplex is overwound with an average value of 9.7 bp per turn, characterised as having a C3'-endo sugar pucker, very low base pair rise and high helical twist and inclination angles. The structure includes 65 ordered water molecules. Only a single row of water molecules is observed in the minor groove due to the presence of hydrophobic 2'-O-methyl groups. As many as five magnesium ions are located in the structure. Two are in the major groove and interact with O(6) and N(7) of guanosine and N(4) of cytidine residues through their hydration spheres. This work provides the first example of molecular interactions of nucleic acids with MPD, which was used as a precipitant, cryo-solvent and resolution enhancing agent. The two MPD molecules intrude into the hydration network in the minor groove, each forming hydrogen bonds between their secondary hydroxyl group and exo-amino functions of guanosine residues. Comparison of the 2'-O-Me(CGCGCG)(2) structure in the P3(2)12 and P6(1)22 crystals delineates stability of the water network within the minor groove to dehydration by MPD and is of interest for evaluating factors governing small molecule binding to RNA. Intrusion of MPD into the minor groove of 2'-O-Me(CGCGCG)(2) is discussed with respect to RNA dehydration, a prerequisite of Z-RNA formation.