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++ Primer Plus习题】13.4

大家好,这里是国中之林! ❥前些天发现了一个巨牛的人工智能学习网站,通俗易懂,风趣幽默,忍不住分享一下给大家。点击跳转到网站。有兴趣的可以点点进去看看← 问题: 解答: main.cpp #include <iostream>#include "port.h"int main() {Port p1;Port p2("Abc", "Bcc", 30);std::cout <<

C++包装器

包装器 在 C++ 中,“包装器”通常指的是一种设计模式或编程技巧,用于封装其他代码或对象,使其更易于使用、管理或扩展。包装器的概念在编程中非常普遍,可以用于函数、类、库等多个方面。下面是几个常见的 “包装器” 类型: 1. 函数包装器 函数包装器用于封装一个或多个函数,使其接口更统一或更便于调用。例如,std::function 是一个通用的函数包装器,它可以存储任意可调用对象(函数、函数

C++11第三弹:lambda表达式 | 新的类功能 | 模板的可变参数

🌈个人主页: 南桥几晴秋 🌈C++专栏: 南桥谈C++ 🌈C语言专栏: C语言学习系列 🌈Linux学习专栏: 南桥谈Linux 🌈数据结构学习专栏: 数据结构杂谈 🌈数据库学习专栏: 南桥谈MySQL 🌈Qt学习专栏: 南桥谈Qt 🌈菜鸡代码练习: 练习随想记录 🌈git学习: 南桥谈Git 🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈🌈�

【C++】_list常用方法解析及模拟实现

相信自己的力量,只要对自己始终保持信心,尽自己最大努力去完成任何事,就算事情最终结果是失败了,努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一:什么是list? •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二:list常用接口 •🌰1.默认成员函数 🔥构造函数(⭐) 🔥析构函数 •🌰2.list对象

06 C++Lambda表达式

lambda表达式的定义 没有显式模版形参的lambda表达式 [捕获] 前属性 (形参列表) 说明符 异常 后属性 尾随类型 约束 {函数体} 有显式模版形参的lambda表达式 [捕获] <模版形参> 模版约束 前属性 (形参列表) 说明符 异常 后属性 尾随类型 约束 {函数体} 含义 捕获:包含零个或者多个捕获符的逗号分隔列表 模板形参:用于泛型lambda提供个模板形参的名

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

6.1.数据结构-c/c++堆详解下篇(堆排序,TopK问题)

上篇:6.1.数据结构-c/c++模拟实现堆上篇(向下,上调整算法,建堆,增删数据)-CSDN博客 本章重点 1.使用堆来完成堆排序 2.使用堆解决TopK问题 目录 一.堆排序 1.1 思路 1.2 代码 1.3 简单测试 二.TopK问题 2.1 思路(求最小): 2.2 C语言代码(手写堆) 2.3 C++代码(使用优先级队列 priority_queue)

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa