c++计算DNA探针的熔解温度

2024-04-23 05:44
文章标签 dna c++ 计算 温度 探针 熔解

本文主要是介绍c++计算DNA探针的熔解温度,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

DNA探针的熔解温度(Tm)是指DNA双链在解离过程中的温度,可以用来估计DNA探针与靶序列的结合强度。

DNA探针富集实验中使用的盐浓度通常是在高盐条件下进行的,以帮助DNA与探针结合并提高富集效率。一般来说,盐浓度在0.5 M到1 M之间是常见的范围,但具体的盐浓度会根据实验的具体条件和目的而有所不同。例如,常用的盐包括NaCl和NaI。

DNA探针富集实验中常用的二价阳离子是Mg^2+。Mg^2+离子在DNA结合和酶活性中起着重要作用,因此通常需要在实验中添加适量的Mg^2+。一般来说,Mg^2+的浓度可以在0.5 mM到10 mM范围内变化,具体的浓度取决于实验条件、所用的DNA探针和靶序列的特性以及所需的结合强度等因素。

//
//  main.cpp
//  test3
//
//  Created by  zhengxueming on 2024/4/22.
//#include <iostream>
#include <math.h>
#include <string.h>#define BAIT_ERROR -1using namespace std;double divalentToMonovalent(double divalent, double dntp) {if(divalent==0) dntp=0;if(divalent<0 || dntp<0) return BAIT_ERROR;if(divalent<dntp)/* According to theory, melting temperature does not depend ondivalent cations */divalent=dntp;return 120*(sqrt(divalent-dntp));
}/*
const string seq: The bait DNA sequencedouble salt_conc: Concentration of divalent cations (millimolar).
一般来说,盐浓度在0.5 M到1 M之间是常见的范围常用的盐包括NaCl和NaIdouble dna_conc: DNA concentration (nanomolar).
double divalent_conc: Concentration of divalent cations (millimolar)
double dntp_conc: Concentration of dNTPs (millimolar)*/double
calBaitTm(const string seq, double salt_conc=50,double divalent_conc=0.5, double dntp_conc=0)
{int GC_count = 0;//const char *p, *end;if (divalentToMonovalent(divalent_conc, dntp_conc) == BAIT_ERROR)return BAIT_ERROR;salt_conc = salt_conc + divalentToMonovalent(divalent_conc, dntp_conc);//std::cout << salt_conc << std::endl;/* Length <= 0 is nonsensical. */for (auto base: seq) {if (base == 'G' || base == 'C')GC_count++;}int64_t len = seq.size();//std::cout << len << std::endl;return81.5+ (16.6 * log10(salt_conc / 1000.0))+ (41.0 * (((double) GC_count) / len))- (600.0 / len);
}int main(int argc, const char * argv[]) {// insert code here...const string normal_bait = "ACAAAGAAATCTCATATTGCTAATGAAGTTGAAGAAAATGACAGCATCTTTGTAAAGCTTCTTAAGATATCAGGAATTATTCTTAAAACGGGAGAGAGTCAGAATCAACTAGGTAATATT";const string high_gc_bait = "GCCCGGGGCGCGCGCGCGCGCGCGCAGTTGAAGAAAATGCCCGGGGCGCGCGCGCGCGCGCGCATGCCCGGGGCGCGCGCGCGCGCGCGCATC";const string high_at_bait = "TTGAATAAAATATTGAATAAAATATTGAATAAAATATTGAATAAAATATTGAATAAAATATTGAATAAAATATTGAATAAAATTTGAATAAAATA";std::cout << "normal_bait tm: " << calBaitTm(normal_bait) << std::endl;std::cout << "high_gc_bait tm: " << calBaitTm(high_gc_bait)<< std::endl;std::cout << "high_at_bait tm: " << calBaitTm(high_at_bait)<< std::endl;return 0;
}

这篇关于c++计算DNA探针的熔解温度的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/927935

相关文章

C++ 中的 if-constexpr语法和作用

《C++中的if-constexpr语法和作用》if-constexpr语法是C++17引入的新语法特性,也被称为常量if表达式或静态if(staticif),:本文主要介绍C++中的if-c... 目录1 if-constexpr 语法1.1 基本语法1.2 扩展说明1.2.1 条件表达式1.2.2 fa

C++中::SHCreateDirectoryEx函数使用方法

《C++中::SHCreateDirectoryEx函数使用方法》::SHCreateDirectoryEx用于创建多级目录,类似于mkdir-p命令,本文主要介绍了C++中::SHCreateDir... 目录1. 函数原型与依赖项2. 基本使用示例示例 1:创建单层目录示例 2:创建多级目录3. 关键注

C++从序列容器中删除元素的四种方法

《C++从序列容器中删除元素的四种方法》删除元素的方法在序列容器和关联容器之间是非常不同的,在序列容器中,vector和string是最常用的,但这里也会介绍deque和list以供全面了解,尽管在一... 目录一、简介二、移除给定位置的元素三、移除与某个值相等的元素3.1、序列容器vector、deque

C++常见容器获取头元素的方法大全

《C++常见容器获取头元素的方法大全》在C++编程中,容器是存储和管理数据集合的重要工具,不同的容器提供了不同的接口来访问和操作其中的元素,获取容器的头元素(即第一个元素)是常见的操作之一,本文将详细... 目录一、std::vector二、std::list三、std::deque四、std::forwa

C++字符串提取和分割的多种方法

《C++字符串提取和分割的多种方法》在C++编程中,字符串处理是一个常见的任务,尤其是在需要从字符串中提取特定数据时,本文将详细探讨如何使用C++标准库中的工具来提取和分割字符串,并分析不同方法的适用... 目录1. 字符串提取的基本方法1.1 使用 std::istringstream 和 >> 操作符示

C++原地删除有序数组重复项的N种方法

《C++原地删除有序数组重复项的N种方法》给定一个排序数组,你需要在原地删除重复出现的元素,使得每个元素只出现一次,返回移除后数组的新长度,不要使用额外的数组空间,你必须在原地修改输入数组并在使用O(... 目录一、问题二、问题分析三、算法实现四、问题变体:最多保留两次五、分析和代码实现5.1、问题分析5.

C++ 各种map特点对比分析

《C++各种map特点对比分析》文章比较了C++中不同类型的map(如std::map,std::unordered_map,std::multimap,std::unordered_multima... 目录特点比较C++ 示例代码 ​​​​​​代码解释特点比较1. std::map底层实现:基于红黑

C++中函数模板与类模板的简单使用及区别介绍

《C++中函数模板与类模板的简单使用及区别介绍》这篇文章介绍了C++中的模板机制,包括函数模板和类模板的概念、语法和实际应用,函数模板通过类型参数实现泛型操作,而类模板允许创建可处理多种数据类型的类,... 目录一、函数模板定义语法真实示例二、类模板三、关键区别四、注意事项 ‌在C++中,模板是实现泛型编程

利用Python和C++解析gltf文件的示例详解

《利用Python和C++解析gltf文件的示例详解》gltf,全称是GLTransmissionFormat,是一种开放的3D文件格式,Python和C++是两个非常强大的工具,下面我们就来看看如何... 目录什么是gltf文件选择语言的原因安装必要的库解析gltf文件的步骤1. 读取gltf文件2. 提

C++快速排序超详细讲解

《C++快速排序超详细讲解》快速排序是一种高效的排序算法,通过分治法将数组划分为两部分,递归排序,直到整个数组有序,通过代码解析和示例,详细解释了快速排序的工作原理和实现过程,需要的朋友可以参考下... 目录一、快速排序原理二、快速排序标准代码三、代码解析四、使用while循环的快速排序1.代码代码1.由快