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++如何通过Qt反射机制实现数据类序列化

《C++如何通过Qt反射机制实现数据类序列化》在C++工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作,所以本文就来聊聊C++如何通过Qt反射机制实现数据类序列化吧... 目录设计预期设计思路代码实现使用方法在 C++ 工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作。由于数据类

Linux下如何使用C++获取硬件信息

《Linux下如何使用C++获取硬件信息》这篇文章主要为大家详细介绍了如何使用C++实现获取CPU,主板,磁盘,BIOS信息等硬件信息,文中的示例代码讲解详细,感兴趣的小伙伴可以了解下... 目录方法获取CPU信息:读取"/proc/cpuinfo"文件获取磁盘信息:读取"/proc/diskstats"文

C++使用printf语句实现进制转换的示例代码

《C++使用printf语句实现进制转换的示例代码》在C语言中,printf函数可以直接实现部分进制转换功能,通过格式说明符(formatspecifier)快速输出不同进制的数值,下面给大家分享C+... 目录一、printf 原生支持的进制转换1. 十进制、八进制、十六进制转换2. 显示进制前缀3. 指

C++中初始化二维数组的几种常见方法

《C++中初始化二维数组的几种常见方法》本文详细介绍了在C++中初始化二维数组的不同方式,包括静态初始化、循环、全部为零、部分初始化、std::array和std::vector,以及std::vec... 目录1. 静态初始化2. 使用循环初始化3. 全部初始化为零4. 部分初始化5. 使用 std::a

C++ vector的常见用法超详细讲解

《C++vector的常见用法超详细讲解》:本文主要介绍C++vector的常见用法,包括C++中vector容器的定义、初始化方法、访问元素、常用函数及其时间复杂度,通过代码介绍的非常详细,... 目录1、vector的定义2、vector常用初始化方法1、使编程用花括号直接赋值2、使用圆括号赋值3、ve

如何高效移除C++关联容器中的元素

《如何高效移除C++关联容器中的元素》关联容器和顺序容器有着很大不同,关联容器中的元素是按照关键字来保存和访问的,而顺序容器中的元素是按它们在容器中的位置来顺序保存和访问的,本文介绍了如何高效移除C+... 目录一、简介二、移除给定位置的元素三、移除与特定键值等价的元素四、移除满足特android定条件的元

Python获取C++中返回的char*字段的两种思路

《Python获取C++中返回的char*字段的两种思路》有时候需要获取C++函数中返回来的不定长的char*字符串,本文小编为大家找到了两种解决问题的思路,感兴趣的小伙伴可以跟随小编一起学习一下... 有时候需要获取C++函数中返回来的不定长的char*字符串,目前我找到两种解决问题的思路,具体实现如下:

C++ Sort函数使用场景分析

《C++Sort函数使用场景分析》sort函数是algorithm库下的一个函数,sort函数是不稳定的,即大小相同的元素在排序后相对顺序可能发生改变,如果某些场景需要保持相同元素间的相对顺序,可使... 目录C++ Sort函数详解一、sort函数调用的两种方式二、sort函数使用场景三、sort函数排序

Java调用C++动态库超详细步骤讲解(附源码)

《Java调用C++动态库超详细步骤讲解(附源码)》C语言因其高效和接近硬件的特性,时常会被用在性能要求较高或者需要直接操作硬件的场合,:本文主要介绍Java调用C++动态库的相关资料,文中通过代... 目录一、直接调用C++库第一步:动态库生成(vs2017+qt5.12.10)第二步:Java调用C++

C/C++错误信息处理的常见方法及函数

《C/C++错误信息处理的常见方法及函数》C/C++是两种广泛使用的编程语言,特别是在系统编程、嵌入式开发以及高性能计算领域,:本文主要介绍C/C++错误信息处理的常见方法及函数,文中通过代码介绍... 目录前言1. errno 和 perror()示例:2. strerror()示例:3. perror(