1def calculate_enzyme_activity(enzyme_conc, substrate_conc, reaction_time, km, vmax):
2    """
3    使用米哈利斯-门腾方程计算酶活性。
4    
5    参数：
6    enzyme_conc (float): 酶浓度（mg/mL）
7    substrate_conc (float): 底物浓度（mM）
8    reaction_time (float): 反应时间（分钟）
9    km (float): 米哈利斯常数（mM）
10    vmax (float): 最大速率（μmol/min）
11    
12    返回：
13    float: 酶活性（U/mg）
14    """
15    reaction_velocity = (vmax * substrate_conc) / (km + substrate_conc)
16    enzyme_activity = reaction_velocity / (enzyme_conc * reaction_time)
17    return enzyme_activity
18
19# 示例用法
20enzyme_conc = 1.0  # mg/mL
21substrate_conc = 10.0  # mM
22reaction_time = 5.0  # min
23km = 5.0  # mM
24vmax = 50.0  # μmol/min
25
26activity = calculate_enzyme_activity(enzyme_conc, substrate_conc, reaction_time, km, vmax)
27print(f"酶活性: {activity:.4f} U/mg")
28

1/**
2 * 使用米哈利斯-门腾方程计算酶活性
3 * @param {number} enzymeConc - 酶浓度（mg/mL）
4 * @param {number} substrateConc - 底物浓度（mM）
5 * @param {number} reactionTime - 反应时间（分钟）
6 * @param {number} km - 米哈利斯常数（mM）
7 * @param {number} vmax - 最大速率（μmol/min）
8 * @returns {number} 酶活性（U/mg）
9 */
10function calculateEnzymeActivity(enzymeConc, substrateConc, reactionTime, km, vmax) {
11  const reactionVelocity = (vmax * substrateConc) / (km + substrateConc);
12  const enzymeActivity = reactionVelocity / (enzymeConc * reactionTime);
13  return enzymeActivity;
14}
15
16// 示例用法
17const enzymeConc = 1.0; // mg/mL
18const substrateConc = 10.0; // mM
19const reactionTime = 5.0; // min
20const km = 5.0; // mM
21const vmax = 50.0; // μmol/min
22
23const activity = calculateEnzymeActivity(enzymeConc, substrateConc, reactionTime, km, vmax);
24console.log(`酶活性: ${activity.toFixed(4)} U/mg`);
25

1public class EnzymeActivityCalculator {
2    /**
3     * 使用米哈利斯-门腾方程计算酶活性
4     * 
5     * @param enzymeConc 酶浓度（mg/mL）
6     * @param substrateConc 底物浓度（mM）
7     * @param reactionTime 反应时间（分钟）
8     * @param km 米哈利斯常数（mM）
9     * @param vmax 最大速率（μmol/min）
10     * @return 酶活性（U/mg）
11     */
12    public static double calculateEnzymeActivity(
13            double enzymeConc, 
14            double substrateConc, 
15            double reactionTime, 
16            double km, 
17            double vmax) {
18        
19        double reactionVelocity = (vmax * substrateConc) / (km + substrateConc);
20        double enzymeActivity = reactionVelocity / (enzymeConc * reactionTime);
21        return enzymeActivity;
22    }
23    
24    public static void main(String[] args) {
25        double enzymeConc = 1.0; // mg/mL
26        double substrateConc = 10.0; // mM
27        double reactionTime = 5.0; // min
28        double km = 5.0; // mM
29        double vmax = 50.0; // μmol/min
30        
31        double activity = calculateEnzymeActivity(
32            enzymeConc, substrateConc, reactionTime, km, vmax);
33        System.out.printf("酶活性: %.4f U/mg%n", activity);
34    }
35}
36

1# R函数用于酶活性计算
2calculate_enzyme_activity <- function(enzyme_conc, substrate_conc, reaction_time, km, vmax) {
3  # 使用米哈利斯-门腾方程计算反应速率
4  reaction_velocity <- (vmax * substrate_conc) / (km + substrate_conc)
5  
6  # 计算酶活性
7  enzyme_activity <- reaction_velocity / (enzyme_conc * reaction_time)
8  
9  return(enzyme_activity)
10}
11
12# 示例用法
13enzyme_conc <- 1.0  # mg/mL
14substrate_conc <- 10.0  # mM
15reaction_time <- 5.0  # min
16km <- 5.0  # mM
17vmax <- 50.0  # μmol/min
18
19activity <- calculate_enzyme_activity(enzyme_conc, substrate_conc, reaction_time, km, vmax)
20cat(sprintf("酶活性: %.4f U/mg", activity))
21