Role of a reaction coordinate in free energy calculations

Abstract

The free energy calculation method emerges as a viable technique for ‘in-silico’ calorimetry. Efficient sampling techniques and the good choice of a reaction path connecting the reactant and the product state enable accurate computations of the free energy differences. We argue that in many cases the thermodynamic integration technique has the lowest variance when the transformation between the reactant and the product state proceeds along the natural path of the studied chemical reaction. We provide examples of free energy calculations for the fragmentation of the charged clusters and the swapping reaction of oligomer formation in proteins that follow a tentative reaction mechanism.     

Fluorescence resonance energy transfer as a structural tool for nucleic acids

Fluorescence resonance energy transfer is a spectroscopic method that provides distance information on macromolecules in solution in the range 20-80 A. It is particularly suited to the analysis of the global structure of nucleic acids because the long-range distance information provides constraints when modelling these important structures. The application of fluorescence resonance energy transfer to nucleic acid structure has seen a resurgence of interest in the past decade, which continues to increase. An especially exciting development is the recent extension to single-molecule studies.     

Position-resolved free energy of solvation for amino acids in lipid membranes from molecular dynamics simulations

Studies of insertion and interactions of amino acids in lipid membranes are pivotal to our understanding of membrane protein structure and function. Calculating the insertion cost as a function of transmembrane helix sequence is thus an important step towards improved membrane protein prediction and eventually drug design. Here, we present position-dependent free energies of solvation for all amino acid analogs along the membrane normal. The profiles cover the entire region from bulk water to hydrophobic core, and were produced from all-atom molecular dynamics simulations. Experimental differences corresponding to mutations and costs for entire segments match experimental data well, and in addition the profiles provide the spatial resolution currently not available from experiments. Polar side-chains largely maintain their hydration and assume quite ordered conformations, which indicates the solvation cost is mainly entropic. The cost of solvating charged side-chains is not only significantly lower than for implicit solvation models, but also close to experiments, meaning these could well maintain their protonation states inside the membrane. The single notable exception to the experimental agreement is proline, which is quite expensive to introduce in vivo despite its hydrophobicity--a difference possibly explained by kinks making it harder to insert helices in the translocon.     

Recent developments in multiscale free energy simulations

Physics-based free energy simulations enable the rigorous calculation of properties, such as conformational equilibria, solvation or binding free energies. While historically most applications have occurred at the atomistic level of resolution, a range of advances in the past years make it possible now to reliably cross the temporal, spatial and theory scales for the modeling of complex systems or the efficient prediction of results at the accuracy level of expensive quantum-mechanical calculations. In this mini-review, we discuss recent methodological advances as well as opportunities opened up by the introduction of machine learning approaches, which tackle the diverse challenges across the different scales, improve the accuracy and feasibility, and push the boundaries of multiscale free energy simulations.     

Circulating nucleic acids as a tumor marker

Patients suffering from malignant diseases have been shown to have increased amounts of cell free nucleic acids in their circulation. As genetic and epigenetic alterations are increasingly characterized in different types of tumors, such changes can be used to detect tumor-derived nucleic acids in the circulation. To date, nearly all tumor-associated nucleic acids have been detected in the plasma or serum of cancer patients. Moreover, increased levels of circulating viral nucleic acids have also been demonstrated in patients with certain cancers associated with viral infections. The concentration of these tumor-associated nucleic acid species is generally related to the tumor load and the extent of the disease. Serial monitoring of plasma nucleic acids thus provides a good way to follow disease progress and to predict the outcome of such patients. In this review, different approaches of detecting tumor-related nucleic acids in the circulation and their potential as tumor markers in the screening, monitoring and prognostication of malignant diseases are discussed.     

Circulating nucleic acids as a new diagnostic tool

The discovery of circulating nucleic acids in the 1940s opened up new possibilities for the non-invasive detection, monitoring and screening of various human disorders. Several tumour markers that enable early cancer detection or tumour behaviour prediction have been detected in the plasma of cancer patients. Maternal plasma analysis can be used to detect certain fetal abnormalities, with the quantification of cell-free nucleic acids used to screen for several pregnancy-associated disorders. Some other applications are in transplant monitoring and graft rejection assessment, and in certain medical emergencies such as trauma and burn severity stratification. Many studies have yielded promising results in this field, but the techniques have yet to be applied in routine clinical practice. Large-scale studies using similar technologies and a broad spectrum of patients are still needed to verify the results of the various studies.     

Fluorescence resonance energy transfer in the study of nucleic acids

Recent data are reviewed on the employment of fluorescence resonance energy transfer (FRET) in studying hybridization and higher structures of nucleic acids as well as their enzyme- and ribozyme-catalyzed reactions.     

Methyl green-pyronin; stoichiometry of reaction with nucleic acids

Study of the stoichiometry of the DNA-methyl green reaction by dialysis, precipitation of stain-nucleic acid mixtures, and the staining of nuclei of known DNA content, indicate that the compound consists of one dye molecule per 10 P. The significance of this result was discussed in the preceding paper (1). Histone and lanthanum (and probably other multivalent cations (3)) compete with the dye for the nucleic acid molecule, indicating a common site of attachment, presumably the phosphoric acid groups. With care in the avoidance of procedures which might depolymerize DNA, and the use of a buffer at about pH 4.1, a quantitative histochemical method for DNA by the use of methyl green is possible. Pyronin staining appears to be of qualitative significance only. Slight differences in degree of polymerization, as between the shad and mammalian DNA appear to have no effect on methyl green staining. It may be that a critical level of polymerization for DNA staining exists. This level must exceed 20 nucleotides to account for the 10 P to 1 dye molecule and the effect on the methyl green absorption spectrum; but it may be considerably greater. Beyond this critical level, whatever it may be, further polymerization probably has no influence on staining.     

Using estimative reaction free energy to predict splice sites and their flanking competitors

A crucial part in the gene structure prediction is to identify the accurate splice sites, not only constitutive but also alternative ones. Here, we use the maximum information principle (MIP) to analyze the conservative segments around splice sites. According to the MIP, a reaction free energy (RFE) expression is deduced, which can be employed to estimate the free energy change during splicing reaction involving a donor or acceptor site. The expression contains not only the background probability factors, but also all kinds of dependencies among both adjacent and non-adjacent bases. We apply the RFE expression to recognize splice sites and their flanking competitors in human genes, the results show high sensitivity and specificity, so the RFE expression accords well with the splicing reaction process. Moreover, the RFE expression is better than previous methods for predicting competitors of splice sites, and it outperforms the reaction free energy subtraction (RFES), that implies RFE competition between a given splice site and its flanking competitor may not be an only primary factor for alternative splice site selection. The work is helpful to not only the understanding of splicing reaction from its relation to MIP, but also the research on computational recognition of splicing sites and alternative splice events.     

Cellular nucleic acid binding protein binds G-rich single-stranded nucleic acids and may function as a nucleic acid chaperone

Cellular nucleic acid binding protein (CNBP) is a small single-stranded nucleic acid binding protein made of seven Zn knuckles and an Arg-Gly rich box. CNBP is strikingly conserved among vertebrates and was reported to play broad-spectrum functions in eukaryotic cells biology. Neither its biological function nor its mechanisms of action were elucidated yet. The main goal of this work was to gain further insights into the CNBP biochemical and molecular features. We studied Bufo arenarum CNBP (bCNBP) binding to single-stranded nucleic acid probes representing the main reported CNBP putative targets. We report that, although bCNBP is able to bind RNA and single-stranded DNA (ssDNA) probes in vitro, it binds RNA as a preformed dimer whereas both monomer and dimer are able to bind to ssDNA. A systematic analysis of variant probes shows that the preferred bCNBP targets contain unpaired guanosine-rich stretches. These data expand the knowledge about CNBP binding stoichiometry and begins to dissect the main features of CNBP nucleic acid targets. Besides, we show that bCNBP presents a highly disordered predicted structure and promotes the annealing and melting of nucleic acids in vitro. These features are typical of proteins that function as nucleic acid chaperones. Based on these data, we propose that CNBP may function as a nucleic acid chaperone through binding, remodeling, and stabilizing nucleic acids secondary structures. This novel CNBP biochemical activity broadens the field of study about its biological function and may be the basis to understand the diverse ways in which CNBP controls gene expression.     

A catalytic mechanism for cysteine N-terminal nucleophile hydrolases, as revealed by free energy simulations

The N-terminal nucleophile (Ntn) hydrolases are a superfamily of enzymes specialized in the hydrolytic cleavage of amide bonds. Even though several members of this family are emerging as innovative drug targets for cancer, inflammation, and pain, the processes through which they catalyze amide hydrolysis remains poorly understood. In particular, the catalytic reactions of cysteine Ntn-hydrolases have never been investigated from a mechanistic point of view. In the present study, we used free energy simulations in the quantum mechanics/molecular mechanics framework to determine the reaction mechanism of amide hydrolysis catalyzed by the prototypical cysteine Ntn-hydrolase, conjugated bile acid hydrolase (CBAH). The computational analyses, which were confirmed in water and using different CBAH mutants, revealed the existence of a chair-like transition state, which might be one of the specific features of the catalytic cycle of Ntn-hydrolases. Our results offer new insights on Ntn-mediated hydrolysis and suggest possible strategies for the creation of therapeutically useful inhibitors.