Interactions between cyclic nucleotides and common cations: an ab initio molecular dynamics study

Ion-Specific Effects on Ion and Polyelectrolyte Solvation

Ion-specific effects on aqueous solvation of monovalent counter ions, Na+ ${^+ }$ , K+ ${^+ }$ , Cl- ${^- }$ , and Br- ${^- }$ , and two model polyelectrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl- ${^- }$ and Br- ${^- }$ anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.

The development of nucleic acids force fields: From an unchallenged past to a competitive future

Molecular dynamics simulations have strongly matured as a method to study biomolecular processes. Their validity, however, is determined by the accuracy of the underlying force fields that describe the forces between all atoms. In this article, we review the development of nucleic acids force fields. We describe the early attempts in the 1990s and emphasize their strong influence on recent force fields. State-of-the-art force fields still use the same Lennard-Jones parameters derived 25 years ago in spite of the fact that these parameters were in general not fitted for nucleic acids. In addition, electrostatic parameters also are deprecated, which may explain some of the current force field deficiencies. We compare different force fields for various systems and discuss new tests of the recently developed Tumuc1 force field. The OL-force fields and Tumuc1 are arguably the best force fields to describe the DNA double helix. However, no force field is flawless. In particular, the description of sugar-puckering remains a problem for nucleic acids force fields. Future refinements are required, so we review methods for force field refinement and give an outlook to the future of force fields.

Molecular Dynamics Study of Interfacial Micromechanical Behaviors of 6H-SiC/Al Composites under Uniaxial Tensile Deformation

This paper investigated the micromechanical behavior of different 6H-SiC/Al systems during the uniaxial tensile loading by using molecular dynamics simulations. The results showed that the interface models responded diversely to the tensile stress when the four low-index surfaces of the Al were used as the variables of the joint surfaces. In terms of their stress-strain properties, the SiC(0001)/Al(001) models exhibited the highest tensile strength and the smallest elongation, while the other models produced certain deformations to relieve the excessive strain, thus increasing the elongation. The SiC(0001)/Al(110) models exhibited the largest elongations among all the models. From the aspect of their deformation characteristics, the SiC(0001)/Al(001) model performed almost no plastic deformation and dislocations during the tensile process. The deformation of the SiC(0001)/Al(110) model was dominated by the slip of the 1/6 <112> Shockley partial dislocations, which contributed to the intersecting stacking faults in the model. The SiC(0001)/Al(111) model produced a large number of dislocations under the tensile loading. Dislocation entanglement was also found in the model. Meanwhile, a unique defect structure consisting of three 1/6 <110> stair-rod dislocations and three stacking faults were found in the model. The plastic deformation in the SiC(0001)/Al(112) interface model was restricted by the L-C lock and was carried out along the 1/6 <110> stair-rod dislocations' direction. These results reveal the interfacial micromechanical behaviors of the 6H-SiC/Al composites and demonstrate the complexity of the deformation systems of the interfaces under stress.

Recalibrating the calcium trap in amino acid carboxyl groups via classical molecular dynamics simulations

In order to use classical molecular dynamics to complement experiments accurately, it is important to use robust descriptions of the system. The interactions between biomolecules, like aspartic and glutamic acid, and dissolved ions are often studied using standard biomolecular force-fields, where the interactions between biomolecules and cations are often not parameterized explicitly. In this study, we have employed metadynamics simulations to investigate different interactions of Ca with aspartic and glutamic acid and constructed the free energy profiles of Ca²⁺-carboxylate association. Starting from a generally accepted, AMBER-based force field, the association was substantially over and under-estimated, depending on the choice of water model (TIP3P and SPC/fw, respectively). To rectify this discrepancy, we have replaced the default calcium parameters. Additionally, we modified the σ_ij value in the hetero-atomic Lennard-Jones interaction by 0.5% to further improve the interaction between Ca and carboxylate, based on comparison with the experimentally determined association constant for Ca with the carboxylate group of L-aspartic acid. The corrected description retrieved the structural properties of the ion pair in agreement with the original biomolecule - Ca²⁺ interaction in AMBER, whilst also producing an association constant comparable to experimental observations. This refined force field was then used to investigate the interactions between amino acids, calcium and carbonate ions during biogenic and biomimetic calcium carbonate mineralisation.

A Cr2O3-doped graphene sensor for early diagnosis of liver cirrhosis: a first-principles study

Liver cirrhosis is among the leading causes of death worldwide. Because of its asymptomatic evolution, timely diagnosis of liver cirrhosis via non-invasive techniques is currently under investigation. Among the diagnostic methods employing volatile organic compounds directly detectable from breath, sensing of limonene (C₁₀H₁₆) represents one of the most promising strategies for diagnosing alcohol liver diseases, including cirrhosis. In the present work, by means of state-of-the-art Density Functional Theory calculations including the U correction, we present an investigation on the sensing capabilities of a chromium-oxide-doped graphene (i.e., Cr₂O₃-graphene) structure toward limonene detection. In contrast with other structures such as g-triazobenzol (g-C₆N₆) monolayers and germanane, which revealed their usefulness in detecting limonene via physisorption, the proposed Cr₂O₃-graphene heterostructure is capable of undergoing chemisorption upon molecular approaching of limonene over its surface. In fact, a high adsorption energy is recorded (∼-1.6 eV). Besides, a positive Moss-Burstein effect is observed upon adsorption of limomene on the Cr₂O₃-graphene heterostructure, resulting in a net increase of the bandgap (∼50%), along with a sizeable shift of the Fermi level toward the conduction band. These findings pave the way toward the experimental validation of such predictions and the employment of Cr₂O₃-graphene heterostructures as sensors of key liver cirrhosis biomarkers.

Self-Consistent Parameterization of DNA Residues for the Non-Polarizable AMBER Force Fields

Fixed-charge (non-polarizable) forcefields are accurate and computationally efficient tools for modeling the molecular dynamics of nucleic acid polymers, particularly DNA, well into the µs timescale. The continued utility of these forcefields depends in part on expanding the residue set in step with advancing nucleic acid chemistry and biology. A key step in parameterizing new residues is charge derivation which is self-consistent with the existing residues. As atomic charges are derived by fitting against molecular electrostatic potentials, appropriate structural models are critical. Benchmarking against the existing charge set used in current AMBER nucleic acid forcefields, we report that quantum mechanical models of deoxynucleosides, even at a high level of theory, are not optimal structures for charge derivation. Instead, structures from molecular mechanics minimization yield charges with up to 6-fold lower RMS deviation from the published values, due to the choice of such an approach in the derivation of the original charge set. We present a contemporary protocol for rendering self-consistent charges as well as optimized charges for a panel of nine non-canonical residues that will permit comparison with literature as well as studying the dynamics of novel DNA polymers.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

Structural Analysis of (p)ppGpp Reveals Its Versatile Binding Pattern for Diverse Types of Target Proteins

(p)ppGpp, highly phosphorylated guanosine, are global regulatory nucleotides that modulate several biochemical events in bacterial physiology ranging from core central dogma to various metabolic pathways. Conventionally, (p)ppGpp collectively refers to two nucleotides, ppGpp, and pppGpp in the literature. Initially, (p)ppGpp has been discovered as a transcription regulatory molecule as it binds to RNA polymerase and regulates transcriptional gene regulation. During the past decade, several other target proteins of (p)ppGpp have been discovered and as of now, more than 30 proteins have been reported to be regulated by the binding of these two signaling nucleotides. The regulation of diverse biochemical activities by (p)ppGpp requires fine-tuned molecular interactions with various classes of proteins so that it can moderate varied functions. Here we report a structural dynamics of (p)ppGpp in the unbound state using well-defined computational tools and its interactions with target proteins to understand the differential regulation by (p)ppGpp at the molecular level. We carried out replica exchange molecular dynamics simulation studies to enhance sampling of conformations during (p)ppGpp simulation. The detailed comparative analysis of torsion angle conformation of ribose sugar of unbound (p)ppGpp and bound states of (p)ppGpp was carried out. The structural dynamics shows that two linear phosphate chains provide plasticity to (p)ppGpp nucleotides for the binding to diverse proteins. Moreover, the intermolecular interactions between (p)ppGpp and target proteins were characterized through various physicochemical parameters including, hydrogen bonds, van der Waal’s interactions, aromatic stacking, and side chains of interacting residues of proteins. Surprisingly, we observed that interactions of (p)ppGpp to target protein have a consensus binding pattern for a particular functional class of enzymes. For example, the binding of (p)ppGpp to RNA polymerase is significantly different from the binding of (p)ppGpp to the proteins involved in the ribosome biogenesis pathway. Whereas, (p)ppGpp binding to enzymes involved in nucleotide metabolism facilitates the functional regulation through oligomerization. Analysis of these datasets revealed that guanine base-specific contacts are key determinants to discriminate functional class of protein. Altogether, our studies provide significant information to understand the differential interaction pattern of (p)ppGpp to its target and this information may be useful to design antibacterial compounds based on (p)ppGpp analogs.

Arsenic-nucleotides interactions: an experimental and computational investigation

Albeit arsenic As(iii) is a well-known carcinogenic contaminant, the modalities by which it interacts with living organisms are still elusive. Details pertaining to the binding properties of As(iii) by common nucleotides such as AMP, ADP and ATP are indeed mostly unknown. Here we present an investigation, conducted via experimental and quantum-based computational approaches, on the stability of the complexes formed by arsenic with those nucleotides. By means of potentiometric and calorimetric measurements, the relative stability of AMP, ADP and ATP has been evaluated as a function of the pH. It turns out that ATP forms more stable structures with As(iii) than ADP which, in turn, better chelates arsenic than AMP. Such a stability sequestration capability of arsenic (ATP > ADP > AMP) has been interpreted on a twofold basis via state-of-the-art ab initio molecular dynamics (AIMD) and metadynamics (MetD) simulations performed on aqueous solutions of As(iii) chelated by AMP and ATP. In fact, we demonstrate that ATP offers a larger number of effective binding sites than AMP, thus indicating a higher statistical probability for chelating arsenic. Moreover, an evaluation of the free energy associated with the interactions that As(iii) establishes with the nucleotide atoms responsible for the binding quantitatively proves the greater effectiveness of ATP as a chelating agent.

Removal of As(III) from Biological Fluids: Mono- versus Dithiolic Ligands

Arsenic is one of the inorganic pollutants typically found in natural waters, and its toxic effects on the human body are currently of great concern. For this reason, the search for detoxifying agents that can be used in a so-called “chelation therapy” is of primary importance. However, to the aim of finding the thermodynamic behavior of efficient chelating agents, extensive speciation studies, capable of reproducing physiological conditions in terms of pH, temperature, and ionic strength, are in order. Here, we report on the acid–base properties of meso-2,3-dimercaptosuccinic acid (DMSA) at different temperatures (i.e., T = 288.15, 298.15, 310.15, and 318.15 K). In particular, its capability to interact with As(III) has been investigated by experimentally evaluating some crucial thermodynamic parameters (ΔH and TΔS), stability constants, and its speciation model. Additionally, in order to gather information on the microscopic coordination modalities of As(III) with the functional groups of DMSA and, at the same time, to better interpret the experimental results, a series of state-of-the-art ab initio molecular dynamics simulations have been performed. For the sake of completeness, the sequestering capabilities of DMSA—a simple dithiol ligand—toward As(III) are directly compared with those recently emerged from similar analyses reported on monothiol ligands.