Molecular docking and MD simulation approach to... | F1000Research "use strict";function _typeof(t){return(_typeof="function"==typeof Symbol&&"symbol"==typeof Symbol.iterator?function(t){return typeof t}:function(t){return t&&"function"==typeof Symbol&&t.constructor===Symbol&&t!==Symbol.prototype?"symbol":typeof t})(t)}!function(){var t=function(){var t,e,o=[],n=window,r=n;for(;r;){try{if(r.frames.__tcfapiLocator){t=r;break}}catch(t){}if(r===n.top)break;r=r.parent}t||(!function t(){var e=n.document,o=!!n.frames.__tcfapiLocator;if(!o)if(e.body){var r=e.createElement("iframe");r.style.cssText="display:none",r.name="__tcfapiLocator",e.body.appendChild(r)}else setTimeout(t,5);return!o}(),n.__tcfapi=function(){for(var t=arguments.length,n=new Array(t),r=0;r 3&&2===parseInt(n[1],10)&&"boolean"==typeof n[3]&&(e=n[3],"function"==typeof n[2]&&n[2]("set",!0)):"ping"===n[0]?"function"==typeof n[2]&&n[2]({gdprApplies:e,cmpLoaded:!1,cmpStatus:"stub"}):o.push(n)},n.addEventListener("message",(function(t){var e="string"==typeof t.data,o={};if(e)try{o=JSON.parse(t.data)}catch(t){}else o=t.data;var n="object"===_typeof(o)&&null!==o?o.__tcfapiCall:null;n&&window.__tcfapi(n.command,n.version,(function(o,r){var a={__tcfapiReturn:{returnValue:o,success:r,callId:n.callId}};t&&t.source&&t.source.postMessage&&t.source.postMessage(e?JSON.stringify(a):a,"*")}),n.parameter)}),!1))};"undefined"!=typeof module?module.exports=t:t()}(); dataLayer = dataLayer || []; // Standard GTM initialization - Google Consent Mode handles consent automatically (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start': new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0], j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src= 'https://www.googletagmanager.com/gtm.js?id='+i+dl+ '>m_auth=hzk0Vc3qFsQYhCrIoHz68A>m_preview=env-1>m_cookies_win=x';f.parentNode.insertBefore(j,f); })(window,document,'script','dataLayer','GTM-MWFK8L5J'); ;window.NREUM||(NREUM={});NREUM.init={distributed_tracing:{enabled:true},privacy:{cookies_enabled:true},ajax:{deny_list:["bam.nr-data.net"]}}; ;NREUM.loader_config={accountID:"438030",trustKey:"438030",agentID:"772317073",licenseKey:"97f8f67f26",applicationID:"772317073"} ;NREUM.info={beacon:"bam.nr-data.net",errorBeacon:"bam.nr-data.net",licenseKey:"97f8f67f26",applicationID:"772317073",sa:1} ;/*! For license information please see nr-loader-spa-1.236.0.min.js.LICENSE.txt */ (()=>{"use strict";var e,t,r={5763:(e,t,r)=>{r.d(t,{P_:()=>l,Mt:()=>g,C5:()=>s,DL:()=>v,OP:()=>T,lF:()=>D,Yu:()=>y,Dg:()=>h,CX:()=>c,GE:()=>b,sU:()=>_});var n=r(8632),i=r(9567);const o={beacon:n.ce.beacon,errorBeacon:n.ce.errorBeacon,licenseKey:void 0,applicationID:void 0,sa:void 0,queueTime:void 0,applicationTime:void 0,ttGuid:void 0,user:void 0,account:void 0,product:void 0,extra:void 0,jsAttributes:{},userAttributes:void 0,atts:void 0,transactionName:void 0,tNamePlain:void 0},a={};function s(e){if(!e)throw new Error("All info objects require an agent identifier!");if(!a[e])throw new Error("Info for ".concat(e," was never set"));return a[e]}function c(e,t){if(!e)throw new Error("All info objects require an agent identifier!");a[e]=(0,i.D)(t,o),(0,n.Qy)(e,a[e],"info")}var u=r(7056);const d=()=>{const e={blockSelector:"[data-nr-block]",maskInputOptions:{password:!0}};return{allow_bfcache:!0,privacy:{cookies_enabled:!0},ajax:{deny_list:void 0,enabled:!0,harvestTimeSeconds:10},distributed_tracing:{enabled:void 0,exclude_newrelic_header:void 0,cors_use_newrelic_header:void 0,cors_use_tracecontext_headers:void 0,allowed_origins:void 0},session:{domain:void 0,expiresMs:u.oD,inactiveMs:u.Hb},ssl:void 0,obfuscate:void 0,jserrors:{enabled:!0,harvestTimeSeconds:10},metrics:{enabled:!0},page_action:{enabled:!0,harvestTimeSeconds:30},page_view_event:{enabled:!0},page_view_timing:{enabled:!0,harvestTimeSeconds:30,long_task:!1},session_trace:{enabled:!0,harvestTimeSeconds:10},harvest:{tooManyRequestsDelay:60},session_replay:{enabled:!1,harvestTimeSeconds:60,sampleRate:.1,errorSampleRate:.1,maskTextSelector:"*",maskAllInputs:!0,get blockClass(){return"nr-block"},get ignoreClass(){return"nr-ignore"},get maskTextClass(){return"nr-mask"},get blockSelector(){return e.blockSelector},set blockSelector(t){e.blockSelector+=",".concat(t)},get maskInputOptions(){return e.maskInputOptions},set maskInputOptions(t){e.maskInputOptions={...t,password:!0}}},spa:{enabled:!0,harvestTimeSeconds:10}}},f={};function l(e){if(!e)throw new Error("All configuration objects require an agent identifier!");if(!f[e])throw new Error("Configuration for ".concat(e," was never set"));return f[e]}function h(e,t){if(!e)throw new Error("All configuration objects require an agent identifier!");f[e]=(0,i.D)(t,d()),(0,n.Qy)(e,f[e],"config")}function g(e,t){if(!e)throw new Error("All configuration objects require an agent identifier!");var r=l(e);if(r){for(var n=t.split("."),i=0;i {r.d(t,{D:()=>i});var n=r(50);function i(e,t){try{if(!e||"object"!=typeof e)return(0,n.Z)("Setting a Configurable requires an object as input");if(!t||"object"!=typeof t)return(0,n.Z)("Setting a Configurable requires a model to set its initial properties");const r=Object.create(Object.getPrototypeOf(t),Object.getOwnPropertyDescriptors(t)),o=0===Object.keys(r).length?e:r;for(let a in o)if(void 0!==e[a])try{"object"==typeof e[a]&&"object"==typeof t[a]?r[a]=i(e[a],t[a]):r[a]=e[a]}catch(e){(0,n.Z)("An error occurred while setting a property of a Configurable",e)}return r}catch(e){(0,n.Z)("An error occured while setting a Configurable",e)}}},6818:(e,t,r)=>{r.d(t,{Re:()=>i,gF:()=>o,q4:()=>n});const n="1.236.0",i="PROD",o="CDN"},385:(e,t,r)=>{r.d(t,{FN:()=>a,IF:()=>u,Nk:()=>f,Tt:()=>s,_A:()=>o,il:()=>n,pL:()=>c,v6:()=>i,w1:()=>d});const n="undefined"!=typeof window&&!!window.document,i="undefined"!=typeof WorkerGlobalScope&&("undefined"!=typeof self&&self instanceof WorkerGlobalScope&&self.navigator instanceof WorkerNavigator||"undefined"!=typeof globalThis&&globalThis instanceof WorkerGlobalScope&&globalThis.navigator instanceof WorkerNavigator),o=n?window:"undefined"!=typeof WorkerGlobalScope&&("undefined"!=typeof self&&self instanceof WorkerGlobalScope&&self||"undefined"!=typeof globalThis&&globalThis instanceof WorkerGlobalScope&&globalThis),a=""+o?.location,s=/iPad|iPhone|iPod/.test(navigator.userAgent),c=s&&"undefined"==typeof SharedWorker,u=(()=>{const e=navigator.userAgent.match(/Firefox[/\s](\d+\.\d+)/);return Array.isArray(e)&&e.length>=2?+e[1]:0})(),d=Boolean(n&&window.document.documentMode),f=!!navigator.sendBeacon},1117:(e,t,r)=>{r.d(t,{w:()=>o});var n=r(50);const i={agentIdentifier:"",ee:void 0};class o{constructor(e){try{if("object"!=typeof e)return(0,n.Z)("shared context requires an object as input");this.sharedContext={},Object.assign(this.sharedContext,i),Object.entries(e).forEach((e=>{let[t,r]=e;Object.keys(i).includes(t)&&(this.sharedContext[t]=r)}))}catch(e){(0,n.Z)("An error occured while setting SharedContext",e)}}}},8e3:(e,t,r)=>{r.d(t,{L:()=>d,R:()=>c});var n=r(2177),i=r(1284),o=r(4322),a=r(3325);const s={};function c(e,t){const r={staged:!1,priority:a.p[t]||0};u(e),s[e].get(t)||s[e].set(t,r)}function u(e){e&&(s[e]||(s[e]=new Map))}function d(){let e=arguments.length>0&&void 0!==arguments[0]?arguments[0]:"",t=arguments.length>1&&void 0!==arguments[1]?arguments[1]:"feature";if(u(e),!e||!s[e].get(t))return a(t);s[e].get(t).staged=!0;const r=[...s[e]];function a(t){const r=e?n.ee.get(e):n.ee,a=o.X.handlers;if(r.backlog&&a){var s=r.backlog[t],c=a[t];if(c){for(var u=0;s&&u {let[t,r]=e;return r.staged}))&&(r.sort(((e,t)=>e[1].priority-t[1].priority)),r.forEach((e=>{let[t]=e;a(t)})))}function f(e,t){var r=e[1];(0,i.D)(t[r],(function(t,r){var n=e[0];if(r[0]===n){var i=r[1],o=e[3],a=e[2];i.apply(o,a)}}))}},2177:(e,t,r)=>{r.d(t,{c:()=>f,ee:()=>u});var n=r(8632),i=r(2210),o=r(1284),a=r(5763),s="nr@context";let c=(0,n.fP)();var u;function d(){}function f(e){return(0,i.X)(e,s,l)}function l(){return new d}function h(){u.aborted=!0,u.backlog={}}c.ee?u=c.ee:(u=function e(t,r){var n={},c={},f={},g=!1;try{g=16===r.length&&(0,a.OP)(r).isolatedBacklog}catch(e){}var p={on:b,addEventListener:b,removeEventListener:y,emit:v,get:x,listeners:w,context:m,buffer:A,abort:h,aborted:!1,isBuffering:E,debugId:r,backlog:g?{}:t&&"object"==typeof t.backlog?t.backlog:{}};return p;function m(e){return e&&e instanceof d?e:e?(0,i.X)(e,s,l):l()}function v(e,r,n,i,o){if(!1!==o&&(o=!0),!u.aborted||i){t&&o&&t.emit(e,r,n);for(var a=m(n),s=w(e),d=s.length,f=0;fn,p:()=>i});var n=r(2177).ee.get("handle");function i(e,t,r,i,o){o?(o.buffer([e],i),o.emit(e,t,r)):(n.buffer([e],i),n.emit(e,t,r))}},4322:(e,t,r)=>{r.d(t,{X:()=>o});var n=r(5546);o.on=a;var i=o.handlers={};function o(e,t,r,o){a(o||n.E,i,e,t,r)}function a(e,t,r,i,o){o||(o="feature"),e||(e=n.E);var a=t[o]=t[o]||{};(a[r]=a[r]||[]).push([e,i])}},3239:(e,t,r)=>{r.d(t,{bP:()=>s,iz:()=>c,m$:()=>a});var n=r(385);let i=!1,o=!1;try{const e={get passive(){return i=!0,!1},get signal(){return o=!0,!1}};n._A.addEventListener("test",null,e),n._A.removeEventListener("test",null,e)}catch(e){}function a(e,t){return i||o?{capture:!!e,passive:i,signal:t}:!!e}function s(e,t){let r=arguments.length>2&&void 0!==arguments[2]&&arguments[2],n=arguments.length>3?arguments[3]:void 0;window.addEventListener(e,t,a(r,n))}function c(e,t){let r=arguments.length>2&&void 0!==arguments[2]&&arguments[2],n=arguments.length>3?arguments[3]:void 0;document.addEventListener(e,t,a(r,n))}},4402:(e,t,r)=>{r.d(t,{Ht:()=>u,M:()=>c,Rl:()=>a,ky:()=>s});var n=r(385);const i="xxxxxxxx-xxxx-4xxx-yxxx-xxxxxxxxxxxx";function o(e,t){return e?15&e[t]:16*Math.random()|0}function a(){const e=n._A?.crypto||n._A?.msCrypto;let t,r=0;return e&&e.getRandomValues&&(t=e.getRandomValues(new Uint8Array(31))),i.split("").map((e=>"x"===e?o(t,++r).toString(16):"y"===e?(3&o()|8).toString(16):e)).join("")}function s(e){const t=n._A?.crypto||n._A?.msCrypto;let r,i=0;t&&t.getRandomValues&&(r=t.getRandomValues(new Uint8Array(31)));const a=[];for(var s=0;s {r.d(t,{Bq:()=>n,Hb:()=>o,oD:()=>i});const n="NRBA",i=144e5,o=18e5},7894:(e,t,r)=>{function n(){return Math.round(performance.now())}r.d(t,{z:()=>n})},7243:(e,t,r)=>{r.d(t,{e:()=>o});var n=r(385),i={};function o(e){if(e in i)return i[e];if(0===(e||"").indexOf("data:"))return{protocol:"data"};let t;var r=n._A?.location,o={};if(n.il)t=document.createElement("a"),t.href=e;else try{t=new URL(e,r.href)}catch(e){return o}o.port=t.port;var a=t.href.split("://");!o.port&&a[1]&&(o.port=a[1].split("/")[0].split("@").pop().split(":")[1]),o.port&&"0"!==o.port||(o.port="https"===a[0]?"443":"80"),o.hostname=t.hostname||r.hostname,o.pathname=t.pathname,o.protocol=a[0],"/"!==o.pathname.charAt(0)&&(o.pathname="/"+o.pathname);var s=!t.protocol||":"===t.protocol||t.protocol===r.protocol,c=t.hostname===r.hostname&&t.port===r.port;return o.sameOrigin=s&&(!t.hostname||c),"/"===o.pathname&&(i[e]=o),o}},50:(e,t,r)=>{function n(e,t){"function"==typeof console.warn&&(console.warn("New Relic: ".concat(e)),t&&console.warn(t))}r.d(t,{Z:()=>n})},2587:(e,t,r)=>{r.d(t,{N:()=>c,T:()=>u});var n=r(2177),i=r(5546),o=r(8e3),a=r(3325);const s={stn:[a.D.sessionTrace],err:[a.D.jserrors,a.D.metrics],ins:[a.D.pageAction],spa:[a.D.spa],sr:[a.D.sessionReplay,a.D.sessionTrace]};function c(e,t){const r=n.ee.get(t);e&&"object"==typeof e&&(Object.entries(e).forEach((e=>{let[t,n]=e;void 0===u[t]&&(s[t]?s[t].forEach((e=>{n?(0,i.p)("feat-"+t,[],void 0,e,r):(0,i.p)("block-"+t,[],void 0,e,r),(0,i.p)("rumresp-"+t,[Boolean(n)],void 0,e,r)})):n&&(0,i.p)("feat-"+t,[],void 0,void 0,r),u[t]=Boolean(n))})),Object.keys(s).forEach((e=>{void 0===u[e]&&(s[e]?.forEach((t=>(0,i.p)("rumresp-"+e,[!1],void 0,t,r))),u[e]=!1)})),(0,o.L)(t,a.D.pageViewEvent))}const u={}},2210:(e,t,r)=>{r.d(t,{X:()=>i});var n=Object.prototype.hasOwnProperty;function i(e,t,r){if(n.call(e,t))return e[t];var i=r();if(Object.defineProperty&&Object.keys)try{return Object.defineProperty(e,t,{value:i,writable:!0,enumerable:!1}),i}catch(e){}return e[t]=i,i}},1284:(e,t,r)=>{r.d(t,{D:()=>n});const n=(e,t)=>Object.entries(e||{}).map((e=>{let[r,n]=e;return t(r,n)}))},4351:(e,t,r)=>{r.d(t,{P:()=>o});var n=r(2177);const i=()=>{const e=new WeakSet;return(t,r)=>{if("object"==typeof r&&null!==r){if(e.has(r))return;e.add(r)}return r}};function o(e){try{return JSON.stringify(e,i())}catch(e){try{n.ee.emit("internal-error",[e])}catch(e){}}}},3960:(e,t,r)=>{r.d(t,{K:()=>a,b:()=>o});var n=r(3239);function i(){return"undefined"==typeof document||"complete"===document.readyState}function o(e,t){if(i())return e();(0,n.bP)("load",e,t)}function a(e){if(i())return e();(0,n.iz)("DOMContentLoaded",e)}},8632:(e,t,r)=>{r.d(t,{EZ:()=>u,Qy:()=>c,ce:()=>o,fP:()=>a,gG:()=>d,mF:()=>s});var n=r(7894),i=r(385);const o={beacon:"bam.nr-data.net",errorBeacon:"bam.nr-data.net"};function a(){return i._A.NREUM||(i._A.NREUM={}),void 0===i._A.newrelic&&(i._A.newrelic=i._A.NREUM),i._A.NREUM}function s(){let e=a();return e.o||(e.o={ST:i._A.setTimeout,SI:i._A.setImmediate,CT:i._A.clearTimeout,XHR:i._A.XMLHttpRequest,REQ:i._A.Request,EV:i._A.Event,PR:i._A.Promise,MO:i._A.MutationObserver,FETCH:i._A.fetch}),e}function c(e,t,r){let i=a();const o=i.initializedAgents||{},s=o[e]||{};return Object.keys(s).length||(s.initializedAt={ms:(0,n.z)(),date:new Date}),i.initializedAgents={...o,[e]:{...s,[r]:t}},i}function u(e,t){a()[e]=t}function d(){return function(){let e=a();const t=e.info||{};e.info={beacon:o.beacon,errorBeacon:o.errorBeacon,...t}}(),function(){let e=a();const t=e.init||{};e.init={...t}}(),s(),function(){let e=a();const t=e.loader_config||{};e.loader_config={...t}}(),a()}},7956:(e,t,r)=>{r.d(t,{N:()=>i});var n=r(3239);function i(e){let t=arguments.length>1&&void 0!==arguments[1]&&arguments[1],r=arguments.length>2?arguments[2]:void 0,i=arguments.length>3?arguments[3]:void 0;return void(0,n.iz)("visibilitychange",(function(){if(t)return void("hidden"==document.visibilityState&&e());e(document.visibilityState)}),r,i)}},1214:(e,t,r)=>{r.d(t,{em:()=>v,u5:()=>N,QU:()=>S,_L:()=>I,Gm:()=>L,Lg:()=>M,gy:()=>U,BV:()=>Q,Kf:()=>ee});var n=r(2177);const i="nr@original";var o=Object.prototype.hasOwnProperty,a=!1;function s(e,t){return e||(e=n.ee),r.inPlace=function(e,t,n,i,o){n||(n="");var a,s,c,u="-"===n.charAt(0);for(c=0;c 2?n-2:0),o=2;o {r(A[T],e,w),r(E[T],e,w)})),r(l._A,"fetch",y),t.on(y+"end",(function(e,r){var n=this;if(r){var i=r.headers.get("content-length");null!==i&&(n.rxSize=i),t.emit(y+"done",[null,r],n)}else t.emit(y+"done",[e],n)})),t}const O={},j=["pushState","replaceState"];function S(e){const t=function(e){return(e||n.ee).get("history")}(e);return!l.il||O[t.debugId]++||(O[t.debugId]=1,s(t).inPlace(window.history,j,"-")),t}var P=r(3239);const C={},R=["appendChild","insertBefore","replaceChild"];function I(e){const t=function(e){return(e||n.ee).get("jsonp")}(e);if(!l.il||C[t.debugId])return t;C[t.debugId]=!0;var r=s(t),i=/[?&](?:callback|cb)=([^&#]+)/,o=/(.*)\.([^.]+)/,a=/^(\w+)(\.|$)(.*)$/;function c(e,t){var r=e.match(a),n=r[1],i=r[3];return i?c(i,t[n]):t[n]}return r.inPlace(Node.prototype,R,"dom-"),t.on("dom-start",(function(e){!function(e){if(!e||"string"!=typeof e.nodeName||"script"!==e.nodeName.toLowerCase())return;if("function"!=typeof e.addEventListener)return;var n=(a=e.src,s=a.match(i),s?s[1]:null);var a,s;if(!n)return;var u=function(e){var t=e.match(o);if(t&&t.length>=3)return{key:t[2],parent:c(t[1],window)};return{key:e,parent:window}}(n);if("function"!=typeof u.parent[u.key])return;var d={};function f(){t.emit("jsonp-end",[],d),e.removeEventListener("load",f,(0,P.m$)(!1)),e.removeEventListener("error",l,(0,P.m$)(!1))}function l(){t.emit("jsonp-error",[],d),t.emit("jsonp-end",[],d),e.removeEventListener("load",f,(0,P.m$)(!1)),e.removeEventListener("error",l,(0,P.m$)(!1))}r.inPlace(u.parent,[u.key],"cb-",d),e.addEventListener("load",f,(0,P.m$)(!1)),e.addEventListener("error",l,(0,P.m$)(!1)),t.emit("new-jsonp",[e.src],d)}(e[0])})),t}var k=r(5763);const H={};function L(e){const t=function(e){return(e||n.ee).get("mutation")}(e);if(!l.il||H[t.debugId])return t;H[t.debugId]=!0;var r=s(t),i=k.Yu.MO;return i&&(window.MutationObserver=function(e){return this instanceof i?new i(r(e,"fn-")):i.apply(this,arguments)},MutationObserver.prototype=i.prototype),t}const z={};function M(e){const t=function(e){return(e||n.ee).get("promise")}(e);if(z[t.debugId])return t;z[t.debugId]=!0;var r=n.c,o=s(t),a=k.Yu.PR;return a&&function(){function e(r){var n=t.context(),i=o(r,"executor-",n,null,!1);const s=Reflect.construct(a,[i],e);return t.context(s).getCtx=function(){return n},s}l._A.Promise=e,Object.defineProperty(e,"name",{value:"Promise"}),e.toString=function(){return a.toString()},Object.setPrototypeOf(e,a),["all","race"].forEach((function(r){const n=a[r];e[r]=function(e){let i=!1;[...e||[]].forEach((e=>{this.resolve(e).then(a("all"===r),a(!1))}));const o=n.apply(this,arguments);return o;function a(e){return function(){t.emit("propagate",[null,!i],o,!1,!1),i=i||!e}}}})),["resolve","reject"].forEach((function(r){const n=a[r];e[r]=function(e){const r=n.apply(this,arguments);return e!==r&&t.emit("propagate",[e,!0],r,!1,!1),r}})),e.prototype=a.prototype;const n=a.prototype.then;a.prototype.then=function(){var e=this,i=r(e);i.promise=e;for(var a=arguments.length,s=new Array(a),c=0;c e())),t};function m(e,t){i.inPlace(t,["onreadystatechange"],"fn-",E)}function b(){var e=this,t=r.context(e);e.readyState>3&&!t.resolved&&(t.resolved=!0,r.emit("xhr-resolved",[],e)),i.inPlace(e,f,"fn-",E)}if(function(e,t){for(var r in e)t[r]=e[r]}(o,p),p.prototype=o.prototype,i.inPlace(p.prototype,J,"-xhr-",E),r.on("send-xhr-start",(function(e,t){m(e,t),function(e){h.push(e),a&&(y?y.then(A):u?u(A):(w=-w,x.data=w))}(t)})),r.on("open-xhr-start",m),a){var y=c&&c.resolve();if(!u&&!c){var w=1,x=document.createTextNode(w);new a(A).observe(x,{characterData:!0})}}else t.on("fn-end",(function(e){e[0]&&e[0].type===d||A()}));function A(){for(var e=0;e {r.d(t,{t:()=>n});const n=r(3325).D.ajax},6660:(e,t,r)=>{r.d(t,{A:()=>i,t:()=>n});const n=r(3325).D.jserrors,i="nr@seenError"},3081:(e,t,r)=>{r.d(t,{gF:()=>o,mY:()=>i,t9:()=>n,vz:()=>s,xS:()=>a});const n=r(3325).D.metrics,i="sm",o="cm",a="storeSupportabilityMetrics",s="storeEventMetrics"},4649:(e,t,r)=>{r.d(t,{t:()=>n});const n=r(3325).D.pageAction},7633:(e,t,r)=>{r.d(t,{Dz:()=>i,OJ:()=>a,qw:()=>o,t9:()=>n});const n=r(3325).D.pageViewEvent,i="firstbyte",o="domcontent",a="windowload"},9251:(e,t,r)=>{r.d(t,{t:()=>n});const n=r(3325).D.pageViewTiming},3614:(e,t,r)=>{r.d(t,{BST_RESOURCE:()=>i,END:()=>s,FEATURE_NAME:()=>n,FN_END:()=>u,FN_START:()=>c,PUSH_STATE:()=>d,RESOURCE:()=>o,START:()=>a});const n=r(3325).D.sessionTrace,i="bstResource",o="resource",a="-start",s="-end",c="fn"+a,u="fn"+s,d="pushState"},7836:(e,t,r)=>{r.d(t,{BODY:()=>A,CB_END:()=>E,CB_START:()=>u,END:()=>x,FEATURE_NAME:()=>i,FETCH:()=>_,FETCH_BODY:()=>v,FETCH_DONE:()=>m,FETCH_START:()=>p,FN_END:()=>c,FN_START:()=>s,INTERACTION:()=>l,INTERACTION_API:()=>d,INTERACTION_EVENTS:()=>o,JSONP_END:()=>b,JSONP_NODE:()=>g,JS_TIME:()=>T,MAX_TIMER_BUDGET:()=>a,REMAINING:()=>f,SPA_NODE:()=>h,START:()=>w,originalSetTimeout:()=>y});var n=r(5763);const i=r(3325).D.spa,o=["click","submit","keypress","keydown","keyup","change"],a=999,s="fn-start",c="fn-end",u="cb-start",d="api-ixn-",f="remaining",l="interaction",h="spaNode",g="jsonpNode",p="fetch-start",m="fetch-done",v="fetch-body-",b="jsonp-end",y=n.Yu.ST,w="-start",x="-end",A="-body",E="cb"+x,T="jsTime",_="fetch"},5938:(e,t,r)=>{r.d(t,{W:()=>o});var n=r(5763),i=r(2177);class o{constructor(e,t,r){this.agentIdentifier=e,this.aggregator=t,this.ee=i.ee.get(e,(0,n.OP)(this.agentIdentifier).isolatedBacklog),this.featureName=r,this.blocked=!1}}},9144:(e,t,r)=>{r.d(t,{j:()=>m});var n=r(3325),i=r(5763),o=r(5546),a=r(2177),s=r(7894),c=r(8e3),u=r(3960),d=r(385),f=r(50),l=r(3081),h=r(8632);function g(){const e=(0,h.gG)();["setErrorHandler","finished","addToTrace","inlineHit","addRelease","addPageAction","setCurrentRouteName","setPageViewName","setCustomAttribute","interaction","noticeError","setUserId"].forEach((t=>{e[t]=function(){for(var r=arguments.length,n=new Array(r),i=0;i 1?r-1:0),i=1;i {e.exposed&&e.api[t]&&o.push(e.api[t](...n))})),o.length>1?o:o[0]}(t,...n)}}))}var p=r(2587);function m(e){let t=arguments.length>1&&void 0!==arguments[1]?arguments[1]:{},m=arguments.length>2?arguments[2]:void 0,v=arguments.length>3?arguments[3]:void 0,{init:b,info:y,loader_config:w,runtime:x={loaderType:m},exposed:A=!0}=t;const E=(0,h.gG)();y||(b=E.init,y=E.info,w=E.loader_config),(0,i.Dg)(e,b||{}),(0,i.GE)(e,w||{}),(0,i.sU)(e,x),y.jsAttributes??={},d.v6&&(y.jsAttributes.isWorker=!0),(0,i.CX)(e,y),g();const T=function(e,t){t||(0,c.R)(e,"api");const h={};var g=a.ee.get(e),p=g.get("tracer"),m="api-",v=m+"ixn-";function b(t,r,n,o){const a=(0,i.C5)(e);return null===r?delete a.jsAttributes[t]:(0,i.CX)(e,{...a,jsAttributes:{...a.jsAttributes,[t]:r}}),x(m,n,!0,o||null===r?"session":void 0)(t,r)}function y(){}["setErrorHandler","finished","addToTrace","inlineHit","addRelease"].forEach((e=>h[e]=x(m,e,!0,"api"))),h.addPageAction=x(m,"addPageAction",!0,n.D.pageAction),h.setCurrentRouteName=x(m,"routeName",!0,n.D.spa),h.setPageViewName=function(t,r){if("string"==typeof t)return"/"!==t.charAt(0)&&(t="/"+t),(0,i.OP)(e).customTransaction=(r||"http://custom.transaction")+t,x(m,"setPageViewName",!0)()},h.setCustomAttribute=function(e,t){let r=arguments.length>2&&void 0!==arguments[2]&&arguments[2];if("string"==typeof e){if(["string","number"].includes(typeof t)||null===t)return b(e,t,"setCustomAttribute",r);(0,f.Z)("Failed to execute setCustomAttribute.\nNon-null value must be a string or number type, but a type of was provided."))}else(0,f.Z)("Failed to execute setCustomAttribute.\nName must be a string type, but a type of was provided."))},h.setUserId=function(e){if("string"==typeof e||null===e)return b("enduser.id",e,"setUserId",!0);(0,f.Z)("Failed to execute setUserId.\nNon-null value must be a string type, but a type of was provided."))},h.interaction=function(){return(new y).get()};var w=y.prototype={createTracer:function(e,t){var r={},i=this,a="function"==typeof t;return(0,o.p)(v+"tracer",[(0,s.z)(),e,r],i,n.D.spa,g),function(){if(p.emit((a?"":"no-")+"fn-start",[(0,s.z)(),i,a],r),a)try{return t.apply(this,arguments)}catch(e){throw p.emit("fn-err",[arguments,this,"string"==typeof e?new Error(e):e],r),e}finally{p.emit("fn-end",[(0,s.z)()],r)}}}};function x(e,t,r,i){return function(){return(0,o.p)(l.xS,["API/"+t+"/called"],void 0,n.D.metrics,g),i&&(0,o.p)(e+t,[(0,s.z)(),...arguments],r?null:this,i,g),r?void 0:this}}function A(){r.e(439).then(r.bind(r,7438)).then((t=>{let{setAPI:r}=t;r(e),(0,c.L)(e,"api")})).catch((()=>(0,f.Z)("Downloading runtime APIs failed...")))}return["actionText","setName","setAttribute","save","ignore","onEnd","getContext","end","get"].forEach((e=>{w[e]=x(v,e,void 0,n.D.spa)})),h.noticeError=function(e,t){"string"==typeof e&&(e=new Error(e)),(0,o.p)(l.xS,["API/noticeError/called"],void 0,n.D.metrics,g),(0,o.p)("err",[e,(0,s.z)(),!1,t],void 0,n.D.jserrors,g)},d.il?(0,u.b)((()=>A()),!0):A(),h}(e,v);return(0,h.Qy)(e,T,"api"),(0,h.Qy)(e,A,"exposed"),(0,h.EZ)("activatedFeatures",p.T),T}},3325:(e,t,r)=>{r.d(t,{D:()=>n,p:()=>i});const n={ajax:"ajax",jserrors:"jserrors",metrics:"metrics",pageAction:"page_action",pageViewEvent:"page_view_event",pageViewTiming:"page_view_timing",sessionReplay:"session_replay",sessionTrace:"session_trace",spa:"spa"},i={[n.pageViewEvent]:1,[n.pageViewTiming]:2,[n.metrics]:3,[n.jserrors]:4,[n.ajax]:5,[n.sessionTrace]:6,[n.pageAction]:7,[n.spa]:8,[n.sessionReplay]:9}}},n={};function i(e){var t=n[e];if(void 0!==t)return t.exports;var o=n[e]={exports:{}};return r[e](o,o.exports,i),o.exports}i.m=r,i.d=(e,t)=>{for(var r in t)i.o(t,r)&&!i.o(e,r)&&Object.defineProperty(e,r,{enumerable:!0,get:t[r]})},i.f={},i.e=e=>Promise.all(Object.keys(i.f).reduce(((t,r)=>(i.f[r](e,t),t)),[])),i.u=e=>(({78:"page_action-aggregate",147:"metrics-aggregate",242:"session-manager",317:"jserrors-aggregate",348:"page_view_timing-aggregate",412:"lazy-feature-loader",439:"async-api",538:"recorder",590:"session_replay-aggregate",675:"compressor",733:"session_trace-aggregate",786:"page_view_event-aggregate",873:"spa-aggregate",898:"ajax-aggregate"}[e]||e)+"."+{78:"ac76d497",147:"3dc53903",148:"1a20d5fe",242:"2a64278a",317:"49e41428",348:"bd6de33a",412:"2f55ce66",439:"30bd804e",538:"1b18459f",590:"cf0efb30",675:"ae9f91a8",733:"83105561",786:"06482edd",860:"03a8b7a5",873:"e6b09d52",898:"998ef92b"}[e]+"-1.236.0.min.js"),i.o=(e,t)=>Object.prototype.hasOwnProperty.call(e,t),e={},t="NRBA:",i.l=(r,n,o,a)=>{if(e[r])e[r].push(n);else{var s,c;if(void 0!==o)for(var u=document.getElementsByTagName("script"),d=0;d {s.onerror=s.onload=null,clearTimeout(h);var i=e[r];if(delete e[r],s.parentNode&&s.parentNode.removeChild(s),i&&i.forEach((e=>e(n))),t)return t(n)},h=setTimeout(l.bind(null,void 0,{type:"timeout",target:s}),12e4);s.onerror=l.bind(null,s.onerror),s.onload=l.bind(null,s.onload),c&&document.head.appendChild(s)}},i.r=e=>{"undefined"!=typeof Symbol&&Symbol.toStringTag&&Object.defineProperty(e,Symbol.toStringTag,{value:"Module"}),Object.defineProperty(e,"__esModule",{value:!0})},i.j=364,i.p="https://js-agent.newrelic.com/",(()=>{var e={364:0,953:0};i.f.j=(t,r)=>{var n=i.o(e,t)?e[t]:void 0;if(0!==n)if(n)r.push(n[2]);else{var o=new Promise(((r,i)=>n=e[t]=[r,i]));r.push(n[2]=o);var a=i.p+i.u(t),s=new Error;i.l(a,(r=>{if(i.o(e,t)&&(0!==(n=e[t])&&(e[t]=void 0),n)){var o=r&&("load"===r.type?"missing":r.type),a=r&&r.target&&r.target.src;s.message="Loading chunk "+t+" failed.\n("+o+": "+a+")",s.name="ChunkLoadError",s.type=o,s.request=a,n[1](s)}}),"chunk-"+t,t)}};var t=(t,r)=>{var n,o,[a,s,c]=r,u=0;if(a.some((t=>0!==e[t]))){for(n in s)i.o(s,n)&&(i.m[n]=s[n]);if(c)c(i)}for(t&&t(r);u {i.r(o);var e=i(3325),t=i(5763);const r=Object.values(e.D);function n(e){const n={};return r.forEach((r=>{n[r]=function(e,r){return!1!==(0,t.Mt)(r,"".concat(e,".enabled"))}(r,e)})),n}var a=i(9144);var s=i(5546),c=i(385),u=i(8e3),d=i(5938),f=i(3960),l=i(50);class h extends d.W{constructor(e,t,r){let n=!(arguments.length>3&&void 0!==arguments[3])||arguments[3];super(e,t,r),this.auto=n,this.abortHandler,this.featAggregate,this.onAggregateImported,n&&(0,u.R)(e,r)}importAggregator(){let e=arguments.length>0&&void 0!==arguments[0]?arguments[0]:{};if(this.featAggregate||!this.auto)return;const r=c.il&&!0===(0,t.Mt)(this.agentIdentifier,"privacy.cookies_enabled");let n;this.onAggregateImported=new Promise((e=>{n=e}));const o=async()=>{let t;try{if(r){const{setupAgentSession:e}=await Promise.all([i.e(860),i.e(242)]).then(i.bind(i,3228));t=e(this.agentIdentifier)}}catch(e){(0,l.Z)("A problem occurred when starting up session manager. This page will not start or extend any session.",e)}try{if(!this.shouldImportAgg(this.featureName,t))return void(0,u.L)(this.agentIdentifier,this.featureName);const{lazyFeatureLoader:r}=await i.e(412).then(i.bind(i,8582)),{Aggregate:o}=await r(this.featureName,"aggregate");this.featAggregate=new o(this.agentIdentifier,this.aggregator,e),n(!0)}catch(e){(0,l.Z)("Downloading and initializing ".concat(this.featureName," failed..."),e),this.abortHandler?.(),n(!1)}};c.il?(0,f.b)((()=>o()),!0):o()}shouldImportAgg(r,n){return r!==e.D.sessionReplay||!1!==(0,t.Mt)(this.agentIdentifier,"session_trace.enabled")&&(!!n?.isNew||!!n?.state.sessionReplay)}}var g=i(7633),p=i(7894);class m extends h{static featureName=g.t9;constructor(r,n){let i=!(arguments.length>2&&void 0!==arguments[2])||arguments[2];if(super(r,n,g.t9,i),("undefined"==typeof PerformanceNavigationTiming||c.Tt)&&"undefined"!=typeof PerformanceTiming){const n=(0,t.OP)(r);n[g.Dz]=Math.max(Date.now()-n.offset,0),(0,f.K)((()=>n[g.qw]=Math.max((0,p.z)()-n[g.Dz],0))),(0,f.b)((()=>{const t=(0,p.z)();n[g.OJ]=Math.max(t-n[g.Dz],0),(0,s.p)("timing",["load",t],void 0,e.D.pageViewTiming,this.ee)}))}this.importAggregator()}}var v=i(1117),b=i(1284);class y extends v.w{constructor(e){super(e),this.aggregatedData={}}store(e,t,r,n,i){var o=this.getBucket(e,t,r,i);return o.metrics=function(e,t){t||(t={count:0});return t.count+=1,(0,b.D)(e,(function(e,r){t[e]=w(r,t[e])})),t}(n,o.metrics),o}merge(e,t,r,n,i){var o=this.getBucket(e,t,n,i);if(o.metrics){var a=o.metrics;a.count+=r.count,(0,b.D)(r,(function(e,t){if("count"!==e){var n=a[e],i=r[e];i&&!i.c?a[e]=w(i.t,n):a[e]=function(e,t){if(!t)return e;t.c||(t=x(t.t));return t.min=Math.min(e.min,t.min),t.max=Math.max(e.max,t.max),t.t+=e.t,t.sos+=e.sos,t.c+=e.c,t}(i,a[e])}}))}else o.metrics=r}storeMetric(e,t,r,n){var i=this.getBucket(e,t,r);return i.stats=w(n,i.stats),i}getBucket(e,t,r,n){this.aggregatedData[e]||(this.aggregatedData[e]={});var i=this.aggregatedData[e][t];return i||(i=this.aggregatedData[e][t]={params:r||{}},n&&(i.custom=n)),i}get(e,t){return t?this.aggregatedData[e]&&this.aggregatedData[e][t]:this.aggregatedData[e]}take(e){for(var t={},r="",n=!1,i=0;i t.max&&(t.max=e),e 2&&void 0!==arguments[2])||arguments[2];super(e,r,j.t,n),c.il&&((0,t.OP)(e).initHidden=Boolean("hidden"===document.visibilityState),(0,N.N)((()=>(0,s.p)("docHidden",[(0,p.z)()],void 0,j.t,this.ee)),!0),(0,O.bP)("pagehide",(()=>(0,s.p)("winPagehide",[(0,p.z)()],void 0,j.t,this.ee))),this.importAggregator())}}var P=i(3081);class C extends h{static featureName=P.t9;constructor(e,t){let r=!(arguments.length>2&&void 0!==arguments[2])||arguments[2];super(e,t,P.t9,r),this.importAggregator()}}var R,I=i(2210),k=i(1214),H=i(2177),L={};try{R=localStorage.getItem("__nr_flags").split(","),console&&"function"==typeof console.log&&(L.console=!0,-1!==R.indexOf("dev")&&(L.dev=!0),-1!==R.indexOf("nr_dev")&&(L.nrDev=!0))}catch(e){}function z(e){try{L.console&&z(e)}catch(e){}}L.nrDev&&H.ee.on("internal-error",(function(e){z(e.stack)})),L.dev&&H.ee.on("fn-err",(function(e,t,r){z(r.stack)})),L.dev&&(z("NR AGENT IN DEVELOPMENT MODE"),z("flags: "+(0,b.D)(L,(function(e,t){return e})).join(", ")));var M=i(6660);class B extends h{static featureName=M.t;constructor(r,n){let i=!(arguments.length>2&&void 0!==arguments[2])||arguments[2];super(r,n,M.t,i),this.skipNext=0;try{this.removeOnAbort=new AbortController}catch(e){}const o=this;o.ee.on("fn-start",(function(e,t,r){o.abortHandler&&(o.skipNext+=1)})),o.ee.on("fn-err",(function(t,r,n){o.abortHandler&&!n[M.A]&&((0,I.X)(n,M.A,(function(){return!0})),this.thrown=!0,(0,s.p)("err",[n,(0,p.z)()],void 0,e.D.jserrors,o.ee))})),o.ee.on("fn-end",(function(){o.abortHandler&&!this.thrown&&o.skipNext>0&&(o.skipNext-=1)})),o.ee.on("internal-error",(function(t){(0,s.p)("ierr",[t,(0,p.z)(),!0],void 0,e.D.jserrors,o.ee)})),this.origOnerror=c._A.onerror,c._A.onerror=this.onerrorHandler.bind(this),c._A.addEventListener("unhandledrejection",(t=>{const r=function(e){let t="Unhandled Promise Rejection: ";if(e instanceof Error)try{return e.message=t+e.message,e}catch(t){return e}if(void 0===e)return new Error(t);try{return new Error(t+(0,D.P)(e))}catch(e){return new Error(t)}}(t.reason);(0,s.p)("err",[r,(0,p.z)(),!1,{unhandledPromiseRejection:1}],void 0,e.D.jserrors,this.ee)}),(0,O.m$)(!1,this.removeOnAbort?.signal)),(0,k.gy)(this.ee),(0,k.BV)(this.ee),(0,k.em)(this.ee),(0,t.OP)(r).xhrWrappable&&(0,k.Kf)(this.ee),this.abortHandler=this.#e,this.importAggregator()}#e(){this.removeOnAbort?.abort(),this.abortHandler=void 0}onerrorHandler(t,r,n,i,o){"function"==typeof this.origOnerror&&this.origOnerror(...arguments);try{this.skipNext?this.skipNext-=1:(0,s.p)("err",[o||new F(t,r,n),(0,p.z)()],void 0,e.D.jserrors,this.ee)}catch(t){try{(0,s.p)("ierr",[t,(0,p.z)(),!0],void 0,e.D.jserrors,this.ee)}catch(e){}}return!1}}function F(e,t,r){this.message=e||"Uncaught error with no additional information",this.sourceURL=t,this.line=r}let U=1;const q="nr@id";function G(e){const t=typeof e;return!e||"object"!==t&&"function"!==t?-1:e===c._A?0:(0,I.X)(e,q,(function(){return U++}))}function V(e){if("string"==typeof e&&e.length)return e.length;if("object"==typeof e){if("undefined"!=typeof ArrayBuffer&&e instanceof ArrayBuffer&&e.byteLength)return e.byteLength;if("undefined"!=typeof Blob&&e instanceof Blob&&e.size)return e.size;if(!("undefined"!=typeof FormData&&e instanceof FormData))try{return(0,D.P)(e).length}catch(e){return}}}var X=i(7243);class W{constructor(e){this.agentIdentifier=e,this.generateTracePayload=this.generateTracePayload.bind(this),this.shouldGenerateTrace=this.shouldGenerateTrace.bind(this)}generateTracePayload(e){if(!this.shouldGenerateTrace(e))return null;var r=(0,t.DL)(this.agentIdentifier);if(!r)return null;var n=(r.accountID||"").toString()||null,i=(r.agentID||"").toString()||null,o=(r.trustKey||"").toString()||null;if(!n||!i)return null;var a=(0,_.M)(),s=(0,_.Ht)(),c=Date.now(),u={spanId:a,traceId:s,timestamp:c};return(e.sameOrigin||this.isAllowedOrigin(e)&&this.useTraceContextHeadersForCors())&&(u.traceContextParentHeader=this.generateTraceContextParentHeader(a,s),u.traceContextStateHeader=this.generateTraceContextStateHeader(a,c,n,i,o)),(e.sameOrigin&&!this.excludeNewrelicHeader()||!e.sameOrigin&&this.isAllowedOrigin(e)&&this.useNewrelicHeaderForCors())&&(u.newrelicHeader=this.generateTraceHeader(a,s,c,n,i,o)),u}generateTraceContextParentHeader(e,t){return"00-"+t+"-"+e+"-01"}generateTraceContextStateHeader(e,t,r,n,i){return i+"@nr=0-1-"+r+"-"+n+"-"+e+"----"+t}generateTraceHeader(e,t,r,n,i,o){if(!("function"==typeof c._A?.btoa))return null;var a={v:[0,1],d:{ty:"Browser",ac:n,ap:i,id:e,tr:t,ti:r}};return o&&n!==o&&(a.d.tk=o),btoa((0,D.P)(a))}shouldGenerateTrace(e){return this.isDtEnabled()&&this.isAllowedOrigin(e)}isAllowedOrigin(e){var r=!1,n={};if((0,t.Mt)(this.agentIdentifier,"distributed_tracing")&&(n=(0,t.P_)(this.agentIdentifier).distributed_tracing),e.sameOrigin)r=!0;else if(n.allowed_origins instanceof Array)for(var i=0;i 2&&void 0!==arguments[2])||arguments[2];super(r,n,Z.t,i),(0,t.OP)(r).xhrWrappable&&(this.dt=new W(r),this.handler=(e,t,r,n)=>(0,s.p)(e,t,r,n,this.ee),(0,k.u5)(this.ee),(0,k.Kf)(this.ee),function(r,n,i,o){function a(e){var t=this;t.totalCbs=0,t.called=0,t.cbTime=0,t.end=E,t.ended=!1,t.xhrGuids={},t.lastSize=null,t.loadCaptureCalled=!1,t.params=this.params||{},t.metrics=this.metrics||{},e.addEventListener("load",(function(r){_(t,e)}),(0,O.m$)(!1)),c.IF||e.addEventListener("progress",(function(e){t.lastSize=e.loaded}),(0,O.m$)(!1))}function s(e){this.params={method:e[0]},T(this,e[1]),this.metrics={}}function u(e,n){var i=(0,t.DL)(r);i.xpid&&this.sameOrigin&&n.setRequestHeader("X-NewRelic-ID",i.xpid);var a=o.generateTracePayload(this.parsedOrigin);if(a){var s=!1;a.newrelicHeader&&(n.setRequestHeader("newrelic",a.newrelicHeader),s=!0),a.traceContextParentHeader&&(n.setRequestHeader("traceparent",a.traceContextParentHeader),a.traceContextStateHeader&&n.setRequestHeader("tracestate",a.traceContextStateHeader),s=!0),s&&(this.dt=a)}}function d(e,t){var r=this.metrics,i=e[0],o=this;if(r&&i){var a=V(i);a&&(r.txSize=a)}this.startTime=(0,p.z)(),this.listener=function(e){try{"abort"!==e.type||o.loadCaptureCalled||(o.params.aborted=!0),("load"!==e.type||o.called===o.totalCbs&&(o.onloadCalled||"function"!=typeof t.onload)&&"function"==typeof o.end)&&o.end(t)}catch(e){try{n.emit("internal-error",[e])}catch(e){}}};for(var s=0;s 1?e[1]=i:e.push(i)}else e[0]&&e[0].headers&&s(e[0].headers,n)&&(this.dt=n);function s(e,t){var r=!1;return t.newrelicHeader&&(e.set("newrelic",t.newrelicHeader),r=!0),t.traceContextParentHeader&&(e.set("traceparent",t.traceContextParentHeader),t.traceContextStateHeader&&e.set("tracestate",t.traceContextStateHeader),r=!0),r}}function x(e,t){this.params={},this.metrics={},this.startTime=(0,p.z)(),this.dt=t,e.length>=1&&(this.target=e[0]),e.length>=2&&(this.opts=e[1]);var r,n=this.opts||{},i=this.target;"string"==typeof i?r=i:"object"==typeof i&&i instanceof Y?r=i.url:c._A?.URL&&"object"==typeof i&&i instanceof URL&&(r=i.href),T(this,r);var o=(""+(i&&i instanceof Y&&i.method||n.method||"GET")).toUpperCase();this.params.method=o,this.txSize=V(n.body)||0}function A(t,r){var n;this.endTime=(0,p.z)(),this.params||(this.params={}),this.params.status=r?r.status:0,"string"==typeof this.rxSize&&this.rxSize.length>0&&(n=+this.rxSize);var o={txSize:this.txSize,rxSize:n,duration:(0,p.z)()-this.startTime};i("xhr",[this.params,o,this.startTime,this.endTime,"fetch"],this,e.D.ajax)}function E(t){var r=this.params,n=this.metrics;if(!this.ended){this.ended=!0;for(var o=0;o 2&&void 0!==arguments[2])||arguments[2];super(e,t,we.t,r),this.importAggregator()}}new class{constructor(e){let t=arguments.length>1&&void 0!==arguments[1]?arguments[1]:(0,_.ky)(16);c._A?(this.agentIdentifier=t,this.sharedAggregator=new y({agentIdentifier:this.agentIdentifier}),this.features={},this.desiredFeatures=new Set(e.features||[]),this.desiredFeatures.add(m),Object.assign(this,(0,a.j)(this.agentIdentifier,e,e.loaderType||"agent")),this.start()):(0,l.Z)("Failed to initial the agent. Could not determine the runtime environment.")}get config(){return{info:(0,t.C5)(this.agentIdentifier),init:(0,t.P_)(this.agentIdentifier),loader_config:(0,t.DL)(this.agentIdentifier),runtime:(0,t.OP)(this.agentIdentifier)}}start(){const t="features";try{const r=n(this.agentIdentifier),i=[...this.desiredFeatures];i.sort(((t,r)=>e.p[t.featureName]-e.p[r.featureName])),i.forEach((t=>{if(r[t.featureName]||t.featureName===e.D.pageViewEvent){const n=function(t){switch(t){case e.D.ajax:return[e.D.jserrors];case e.D.sessionTrace:return[e.D.ajax,e.D.pageViewEvent];case e.D.sessionReplay:return[e.D.sessionTrace];case e.D.pageViewTiming:return[e.D.pageViewEvent];default:return[]}}(t.featureName);n.every((e=>r[e]))||(0,l.Z)("".concat(t.featureName," is enabled but one or more dependent features has been disabled (").concat((0,D.P)(n),"). This may cause unintended consequences or missing data...")),this.features[t.featureName]=new t(this.agentIdentifier,this.sharedAggregator)}})),(0,T.Qy)(this.agentIdentifier,this.features,t)}catch(e){(0,l.Z)("Failed to initialize all enabled instrument classes (agent aborted) -",e);for(const e in this.features)this.features[e].abortHandler?.();const r=(0,T.fP)();return delete r.initializedAgents[this.agentIdentifier]?.api,delete r.initializedAgents[this.agentIdentifier]?.[t],delete this.sharedAggregator,r.ee?.abort(),delete r.ee?.get(this.agentIdentifier),!1}}}({features:[J,m,S,class extends h{static featureName=oe;constructor(t,r){if(super(t,r,oe,!(arguments.length>2&&void 0!==arguments[2])||arguments[2]),!c.il)return;const n=this.ee;let i;(0,k.QU)(n),this.eventsEE=(0,k.em)(n),this.eventsEE.on(se,(function(e,t){this.bstStart=(0,p.z)()})),this.eventsEE.on(ae,(function(t,r){(0,s.p)("bst",[t[0],r,this.bstStart,(0,p.z)()],void 0,e.D.sessionTrace,n)})),n.on(ce+ne,(function(e){this.time=(0,p.z)(),this.startPath=location.pathname+location.hash})),n.on(ce+ie,(function(t){(0,s.p)("bstHist",[location.pathname+location.hash,this.startPath,this.time],void 0,e.D.sessionTrace,n)}));try{i=new PerformanceObserver((t=>{const r=t.getEntries();(0,s.p)(te,[r],void 0,e.D.sessionTrace,n)})),i.observe({type:re,buffered:!0})}catch(e){}this.importAggregator({resourceObserver:i})}},C,xe,B,class extends h{static featureName=de;constructor(e,r){if(super(e,r,de,!(arguments.length>2&&void 0!==arguments[2])||arguments[2]),!c.il)return;if(!(0,t.OP)(e).xhrWrappable)return;try{this.removeOnAbort=new AbortController}catch(e){}let n,i=0;const o=this.ee.get("tracer"),a=(0,k._L)(this.ee),s=(0,k.Lg)(this.ee),u=(0,k.BV)(this.ee),d=(0,k.Kf)(this.ee),f=this.ee.get("events"),l=(0,k.u5)(this.ee),h=(0,k.QU)(this.ee),g=(0,k.Gm)(this.ee);function m(e,t){h.emit("newURL",[""+window.location,t])}function v(){i++,n=window.location.hash,this[ve]=(0,p.z)()}function b(){i--,window.location.hash!==n&&m(0,!0);var e=(0,p.z)();this[pe]=~~this[pe]+e-this[ve],this[ye]=e}function y(e,t){e.on(t,(function(){this[t]=(0,p.z)()}))}this.ee.on(ve,v),s.on(be,v),a.on(be,v),this.ee.on(ye,b),s.on(ge,b),a.on(ge,b),this.ee.buffer([ve,ye,"xhr-resolved"],this.featureName),f.buffer([ve],this.featureName),u.buffer(["setTimeout"+le,"clearTimeout"+fe,ve],this.featureName),d.buffer([ve,"new-xhr","send-xhr"+fe],this.featureName),l.buffer([me+fe,me+"-done",me+he+fe,me+he+le],this.featureName),h.buffer(["newURL"],this.featureName),g.buffer([ve],this.featureName),s.buffer(["propagate",be,ge,"executor-err","resolve"+fe],this.featureName),o.buffer([ve,"no-"+ve],this.featureName),a.buffer(["new-jsonp","cb-start","jsonp-error","jsonp-end"],this.featureName),y(l,me+fe),y(l,me+"-done"),y(a,"new-jsonp"),y(a,"jsonp-end"),y(a,"cb-start"),h.on("pushState-end",m),h.on("replaceState-end",m),window.addEventListener("hashchange",m,(0,O.m$)(!0,this.removeOnAbort?.signal)),window.addEventListener("load",m,(0,O.m$)(!0,this.removeOnAbort?.signal)),window.addEventListener("popstate",(function(){m(0,i>1)}),(0,O.m$)(!0,this.removeOnAbort?.signal)),this.abortHandler=this.#e,this.importAggregator()}#e(){this.removeOnAbort?.abort(),this.abortHandler=void 0}}],loaderType:"spa"})})(),window.NRBA=o})(); window.jQuery || document.write(' ') CKEDITOR_BASEPATH='https://f1000research.com/js/vendor/ckeditor/' window.reactTheme = 'research'; window.MathJax = { CommonHTML: { linebreaks: { automatic: true } }, 'HTML-CSS': { linebreaks: { automatic: true } }, SVG: { linebreaks: { automatic: true } }, AuthorInit: function() { MathJax.Hub.Register.MessageHook('End Process', function () { let timeout = false; // holder for timeout id const delay = 250; // delay after event is "complete" to run callback const reflowMath = function() { const dispFormulas = document.querySelectorAll('.disp-formula.panel'); if (!dispFormulas) { return; } for (const dispFormula of dispFormulas) { const child = dispFormula.querySelector('.MathJax_Preview').nextSibling.firstChild; const isMultiline = MathJax.Hub.getAllJax(dispFormula)[0].root.isMultiline; if (dispFormula.offsetWidth < child.offsetWidth || isMultiline) { MathJax.Hub.Queue(['Rerender', MathJax.Hub, dispFormula]); } } }; window.addEventListener('resize', function() { clearTimeout(timeout); // clear the timeout timeout = setTimeout(reflowMath, delay); // start timing for event "completion" }); }); }, }; if (window.location.hash == '#_=_'){ window.location = window.location.href.split('#')[0] } !function(f,b,e,v,n,t,s){if(f.fbq)return;n=f.fbq=function() {n.callMethod? n.callMethod.apply(n,arguments):n.queue.push(arguments)} ;if(!f._fbq)f._fbq=n; n.push=n;n.loaded=!0;n.version='2.0';n.queue=[];t=b.createElement(e);t.async=!0; t.src=v;s=b.getElementsByTagName(e)[0];s.parentNode.insertBefore(t,s)}(window, document,'script','https://connect.facebook.net/en_US/fbevents.js'); fbq('init', '1641728616063202'); fbq('track', "PixelInitialized", {}); (function(h,o,t,j,a,r){ h.hj=h.hj||function(){(h.hj.q=h.hj.q||[]).push(arguments)}; h._hjSettings={hjid:2318163,hjsv:6}; a=o.getElementsByTagName('head')[0]; r=o.createElement('script');r.async=1; r.src=t+h._hjSettings.hjid+j+h._hjSettings.hjsv; a.appendChild(r); })(window,document,'https://static.hotjar.com/c/hotjar-','.js?sv='); search file_upload Submit your research search menu close search Browse Gateways & Collections How to Publish Submit your Research My Submissions Article Guidelines Article Guidelines (New Versions) Open Data, Software and Code Guidelines Open Data and Accessible Source Materials Guidelines (HSS) Open Data, Software and Code Guidelines (PSE) Prepublication Checks Production Process Posters and Slides Guidelines Document Guidelines Article Processing Charges Peer Review Finding Article Reviewers About How it Works For Reviewers Our Advisors Policies Glossary FAQs For Developers Newsroom Contact My Research Submissions Content and Tracking Alerts My Details Sign In file_upload Submit your research { "@context": "https://schema.org", "@type": "ScholarlyArticle", "mainEntityOfPage": { "@type": "WebPage", "@id": "https://f1000research.com/articles/13-1271" }, "headline": "Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against...", "datePublished": "2024-10-24T15:09:09", "dateModified": "2025-03-18T10:56:49", "author": [ { "@type": "Person", "name": "Pranaya Sankaranarayanan" }, { "@type": "Person", "name": "Dicky John Davis G" }, { "@type": "Person", "name": "Abhinand PA" }, { "@type": "Person", "name": "M Manikandan" }, { "@type": "Person", "name": "Arabinda Ghosh" } ], "publisher": { "@type": "Organization", "name": "F1000Research", "logo": { "@type": "ImageObject", "url": "https://f1000research.com/img/AMP/F1000Research_image.png", "height": 480, "width": 60 } }, "image": { "@type": "ImageObject", "url": "https://f1000research.com/img/AMP/F1000Research_image.png", "height": 1200, "width": 150 }, "description": " Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, demonstrates high aggressiveness, and has a relatively poor prognosis and clinical outcome. Objective To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified target. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. Result The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo. " } { "@context": "http://schema.org", "@type": "BreadcrumbList", "itemListElement": [ { "@type": "ListItem", "position": "1", "item": { "@id": "https://f1000research.com/", "name": "Home" } }, { "@type": "ListItem", "position": "2", "item": { "@id": "https://f1000research.com/browse/articles", "name": "Browse" } }, { "@type": "ListItem", "position": "3", "item": { "@id": "https://f1000research.com/articles/13-1271", "name": "Molecular docking and MD simulation approach to identify potential..." } } ] } Home Browse Molecular docking and MD simulation approach to identify potential... ALL Metrics - Views Downloads Get PDF Get XML Cite How to cite this article Sankaranarayanan P, G DJD, PA A et al. Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.12688/f1000research.155657.2 ) NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article. Close Copy Citation Details Export Export Citation Sciwheel EndNote Ref. Manager Bibtex ProCite Sente EXPORT Select a format first Track Share ▬ ✚ Research Article Revised Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] Pranaya Sankaranarayanan https://orcid.org/0000-0001-8076-6994 1 , Dicky John Davis G https://orcid.org/0000-0003-2325-7952 1 , Abhinand PA 1 , M Manikandan 2 , Arabinda Ghosh 3 Pranaya Sankaranarayanan https://orcid.org/0000-0001-8076-6994 1 , Dicky John Davis G https://orcid.org/0000-0003-2325-7952 1 , [...] Abhinand PA 1 , M Manikandan 2 , Arabinda Ghosh 3 PUBLISHED 18 Mar 2025 Author details Author details 1 Department of Bioinformatics, Sri Ramachandra Institute of Higher Education and Research (Deemed to be University), Chennai, Tamil Nadu, 600116, India 2 Department of Medical Genetics, Manipal Hospitals, Bengaluru, Karnataka, 560 017, India 3 Department of Botany, Gauhati University, Guwahati, Assam, India Pranaya Sankaranarayanan Roles: Conceptualization, Data Curation, Formal Analysis, Methodology, Writing – Original Draft Preparation Dicky John Davis G Roles: Investigation, Project Administration, Supervision, Validation, Visualization, Writing – Review & Editing Abhinand PA Roles: Conceptualization, Formal Analysis, Software, Validation, Writing – Review & Editing M Manikandan Roles: Data Curation, Formal Analysis, Software, Visualization Arabinda Ghosh Roles: Formal Analysis, Resources, Software, Writing – Original Draft Preparation OPEN PEER REVIEW DETAILS REVIEWER STATUS This article is included in the Bioinformatics gateway. Abstract Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, demonstrates high aggressiveness, and has a relatively poor prognosis and clinical outcome. Objective To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified target. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. Result The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo. READ ALL READ LESS Keywords Triple Negative Breast Cancer, AR target, Phytochemicals, 2–hydroxy naringenin, Virtual screening, Molecular Docking, Molecular dynamics simulation. Corresponding Author(s) Dicky John Davis G ( [email protected] ) Close Corresponding author: Dicky John Davis G Competing interests: No competing interests were disclosed. Grant information: This work was supported by the Founder-Chancellor Shri. N. P. V. Ramasamy Udayar Research Fellowship (U02B160480), Sri Ramachandra Institute of Higher Education and Research. The funders had no role in the study design, data collection and analysis, decision to publish, or manuscript preparation. Copyright: © 2025 Sankaranarayanan P et al . This is an open access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. How to cite: Sankaranarayanan P, G DJD, PA A et al. Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.12688/f1000research.155657.2 ) First published: 24 Oct 2024, 13 :1271 ( https://doi.org/10.12688/f1000research.155657.1 ) Latest published: 18 Mar 2025, 13 :1271 ( https://doi.org/10.12688/f1000research.155657.2 ) Revised Amendments from Version 1 The major differences between the two versions of the article primarily involve refinements to the text, improved clarity, and potentially additional references. However, both versions maintain the same core objectives, methodology, and results. The major differences between the two versions of the article primarily involve refinements to the text, improved clarity, and potentially additional references. However, both versions maintain the same core objectives, methodology, and results. See the authors' detailed response to the review by Shiek S S J Ahmed READ REVIEWER RESPONSES Introduction Breast cancer is the most common type of cancer worldwide, as reported by the World Health Organization (WHO) in 2020 with over 7.8 million women living in the last five years diagnosed with breast cancer. 1 It is responsible for 685,000 deaths worldwide. However, it should be noted that breast cancer is a non-homogenous condition that can be classified into several significant subtypes based on the expression of their genes. Triple-negative breast cancers (TNBC) are characterized by the absence of estrogen, progesterone, and ERBB2 receptors, and are specifically identified as estrogen receptor (ER)-negative, progesterone receptor (PR)-negative, and human epidermal growth factor receptor 2 (HER2). TNBC accounts for 12%–17% of all breast cancers. 2 Sandhu et al. revealed a considerably greater prevalence of TNBC in India than in Western populations. Approximately one in three women diagnosed with breast cancer in India was found to have triple-negative disease. TNBC exhibits unique clinical and pathologic features, is highly aggressive, and has a relatively poor prognosis and clinical outcome. 3 Currently, there is no recognized targeted treatment for TNBC. The primary treatment options for TNBC involve chemotherapy utilizing anthracyclines, taxanes, and/or platinum compounds as the major treatment modalities. A significant proportion of TNBC patients fail to attain Pathological Complete Response (pCR) with standard chemotherapy, prompting concerns about the effectiveness and safety of the chosen chemotherapy. 4 A better understanding of the pathological mechanisms of TNBC onset and progression and the molecular interactions underlying the etiology of the condition can help improve the prophylaxis and design of novel targeted treatment against this cancer type. 5 Gene expression profiling can be invaluable for detecting transcriptional variations between normal and malignant cells and can be extensively used to study gene phenotype associations in breast neoplasms. 6 Protein interaction networks potentially signify patterns in network connectivity between proteins, which can differ between breast cancer subtypes. 7 Phytochemicals are natural, non-toxic compounds found in plants that possess disease-protective or preventive properties. 8 They modulate the molecular pathways associated with cancer growth and progression. 9 The present study aimed to identify a novel therapeutic target protein for TNBC by integrating differential gene expression studies with protein-protein interactions and network topology analysis. Subsequently, phytochemicals with reported anti-breast cancer activities will be subjected to virtual screening by molecular docking against the identified novel target. To validate these findings, Molecular Mechanics with Generalised Born Surface Area solvation and Molecular Dynamics simulations were performed. Based on their binding affinity to the target protein, novel therapeutic phytochemical lead molecules with anti-TNBC activity were identified. Method Gene expression profiling of TNBC microarray datasets A thorough literature mining effort encompassing all eligible studies on gene expression in TNBC was conducted. The search involved querying the Gene Expression Omnibus (GEO) datasets. Gene expression profiling was performed using GEO2R to identify significantly upregulated genes. Figure 1 presents an overview of the methodology. Figure 1. Overview of methodology. During the literature mining process, a microarray dataset was obtained from the NCBI GEO repository using the accession number GSE45498 annotated in the GPL16299 platform. This dataset encompasses 40 samples from healthy normal tissues, 160 from individuals with cancer, and 54 from metastatic cases. NGS datasets were obtained from the NCBI GEO repository using accession number GSE214101 annotated in the GPL20301 platform. This dataset included 24 samples derived from the MDA-MB-231 and MDA-MB-436 cell lines. Gene expression profiling values underwent log (base2) transformation and percentage shift normalization was applied. To assess the differences in gene expression between normal and diseased samples, the fold change for each gene was individually calculated. A threshold of 1.25-fold change was used to categorize genes as being upregulated. Gene expression profiling followed the protocol reported previously. 10 Study of protein–protein interactions The selected genes were subjected to the Bisogenet plug-in of Cytoscape to identify protein-protein interactions of all genes differentially regulated in TNBC. STRING is an open-source bioinformatics platform integrated in Cytoscape, designed for the study of both predicted and known protein-protein interactions. This database gathers, evaluates, and integrates information on protein-protein interactions from all publicly available sources. Additionally, it augments these data with computational predictions. 11 These interactions encompass both indirect (functional) and direct (physical) associations. 12 The genes were uploaded and a string network was built. Molecular Complex Detection (MCODE) detects Protein-Protein Interactions subnetworks and highly interconnected clusters within the PPI network. 13 PPI networks were broken down into top-ranked dense cliques (sub-clusters) using the MCODE plugin. The top-ranked dense clique was selected for further analysis. Building a library of phytochemicals with anti-breast cancer activity Phytochemicals are naturally occurring biologically active chemical compounds found in plants that serve as medicinal ingredients and nutrients, offering health benefits to humans. 14 Many natural products and their analogs have been identified as potent anticancer agents and the anticancer properties of various plants and phytochemicals. 15 Phytochemicals were identified through a systematic literature search indicating anti-breast cancer activity were selected, and their 3D structures in SDF format were retrieved from PubChem database. Subsequently, phytochemicals that did not conform to Lipinski’s rule of five were excluded, and the remaining compounds were subjected to further analyses. Virtual screening Understanding the fundamental principles governing how protein receptors recognize, interact, and form associations with molecular substrates and inhibitors is crucial for drug discovery. PyRx v0.8 software 16 with an inbuilt AutoDock Vina 1.2.5 17 for molecular docking was used to scan phytochemicals conforming to Lipinski’s rule of 5. AutoDock Vina uses a semi-empirical free-energy force field to predict the binding free energies of small molecules to macromolecular targets. The human Androgen Receptor (PDB ID: 1E3G) was sourced from the RCSB Protein Data Bank. Initially, the protein structure underwent a curation process to remove any crystallographic water molecules and heteroatoms that might interfere with docking simulations. Subsequently, energy minimization was performed using UCSF Chimera vs 1.54 ( https://www.cgl.ucsf.edu/chimera/ ) to optimize the geometry of the protein. The steepest descent algorithm was applied for 100 steps, which is a common approach to relieve steric clashes and achieve a more stable conformation. Partial charges were then assigned to the protein using the AMBER ff14SB force field, which is well known for accurately modeling protein dynamics and interactions. The co-crystallized ligand metribolone (R18) was used as the control, and the ligands were docked at its active site. ADMET - ProTox II The development of high-quality in silico ADMET models will enable compound efficacy and druggability features to be optimized concurrently, thereby improving the overall quality of drug candidates. 18 ProTox-II was used to experimentally validate the chemical toxicity and their combination. It uses machine learning models, the most common features, pharmacophore-based, fragment propensities, and chemical similarity to forecast different toxicity endpoints. 19 Based on the virtual screening results, the top ten phytochemical compounds were chosen for ADMET analysis. Induced fit docking Induced fit docking was carried out using Schrodinger vs. 2020.3, which takes into account the flexibility of both the protein receptor and ligand, allowing for conformational changes to occur upon binding. The energy-minimized ligands were saved in PDB format for compatibility with the Schrodinger software, and the partial charges of the ligands were assigned, such as Gasteiger charges, which estimate the distribution of charges on the molecule based on its structure. Similarly, the protein charges may also be assigned using OPLS_2005 force fields to accurately capture its electrostatic properties. The grid box is a crucial parameter in docking simulations, as it defines the search space where the ligand can orient itself around the protein receptor. The dimensions of the grid box are typically specified in terms of the number of grid points along each axis (nx, ny, nz) and the grid spacing (Å) around the binding cavity residues LEU701, LEU707, MET742, MET745, ARG752, MET780, MET787, ALA748, LEU880, LEU873, PHE876, MET895, ILE899, THR877, GLN774, PHE764, LEU746, GLY708, GLN711, TRP741, ASN705. The dimensions were set to (58, 64, and 52 Å), providing a sufficient volume to explore potential binding modes of the ligand within the protein’s active site with a charge cutoff polarity set for a charge cutoff of 0.25 Å. Molecular dynamics simulation Molecular dynamics (MD) simulations were conducted for the docked complex of the human Androgen Receptor with the best-docked molecule, employing Schrodinger Desmond 2020.1. 20 The OPLS-2005 force field, 21 along with an explicit solvent model using SPC water molecules, 22 were employed in this system. The simulation was performed in a periodic boundary solvation box with dimensions of 10 × 10 × 10 Å. To neutralize the charge, Na+ ions were added, and a 0.15 M NaCl solution was added to mimic the physiological environment. The initial equilibration was carried out using an NVT ensemble for 10 ns to allow the system to relax over the protein-ligand complexes. Subsequently, a short run of equilibration and minimization was conducted using an NPT ensemble for 12 ns. The NPT ensemble utilized the Nose-Hoover chain coupling scheme 23 with a temperature set at 37 °C, relaxation time of 1.0 ps, and pressure maintained at 1 bar in all simulations. A time step of 2 fs was used. Pressure control was achieved using the Martyna-Tuckerman-Klein chain coupling scheme 24 with a relaxation time of 2 ps. The long-range electrostatic interactions were calculated using the particle mesh Ewald method, 25 and the Coulomb interaction radius was fixed at 9 Å. A RESPA integrator with a time step of 2 fs was used for each trajectory to calculate the bonded forces. The final production run was extended for 100 ns for the Human Androgen Receptor with the best-docked molecule complex. To track the stability of the MD simulations, a variety of parameters were computed, including the number of hydrogen bonds, radius of gyration (Rg), root-mean-square fluctuation (RMSF), and root-mean-square deviation (RMSD). Binding free energy analysis Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approaches are less computationally intensive than biochemical free energy methods and more precise than most molecular docking scoring systems. This method is useful for predicting the binding free energy in molecular systems. MM-GBSA is a useful technique for comprehending the impact of mutations on large biomolecular systems. 26 Biomolecular research has been utilized in investigations of protein folding, protein-ligand binding, protein-protein interactions etc. 27 The MM-GBSA approach was used to determine the binding free energies of the ligand-protein complexes. The MM-GBSA binding free energy was computed using the Python script thermal mmgbsa.py in the simulation trajectory with the VSGB solvation model and OPLS5 force field over the last 50 frames with a 1 step sampling size. The binding free energy of MM-GBSA (kcal/mol) was calculated using the additivity principle, wherein the differences in free energies, GBSA solvation energies, and surface area energies of ligand-protein complexes compared to their respective total energies of them individually were calculated. Results Differentially expressed genes (DEGs) analysis Gene expression in TNBC and normal microarray datasets was compared to assess the underlying molecular pathways driving TNBC, and further network analysis was performed. Boolean operators and relevant filters were used to filter the microarray datasets using the GEO2R. The Benjamini-Hochberg-Yekutieli approach was used to adjust the P-value for the DEGs, and only the top 10% of the upregulated genes (P-value 1.25 and P-value 1.25 and P-value <0.05. Gene ID Description log 2 FC p-Value ESR1 Estrogen Receptor 1 3.45098 8.51E-14 IGFBP6 Insulin-like growth factors binding protein-6 3.115311 1.71E-14 NGFR Nerve growth factor receptor 3.069617 3.26E-10 DLC1 Deleted in liver cancer 1 2.833933 1.03E-12 TGFBR3 Transforming Growth Factor Beta Receptor 3 2.631049 2.84E-10 EGR1 Early growth response factor 1 2.31673 5.84E-11 NTRK2 Neurotrophic Tyrosine Receptor Kinase 2.19261 1.77E-06 PPARG Peroxisome proliferator-activated receptor gamma 2.151492 3.32E-10 CD34 CD34 1.887035 5.93E-09 IGF1 Insulin-Like Growth Factor-1 1.870246 1.53E-10 FOS FOS 1.734574 5.27E-08 CAV1 Caveolin 1 1.694425 6.72E-07 FGF2 Fibroblast Growth Factor 2 1.61343 4.41E-04 KIT KIT 1.547563 2.93E-05 AR Androgen Receptor 1.381295 2.51E-04 Table 2. The list of upregulated genes in dataset GSE214101 with LogFC > 1.25 and P-value <0.05. Gene ID Description log 2 FC p-value CDH4 cadherin 4 2.805 2.26E-06 MAP 2K6 mitogen-activated protein kinase kinase 6 2.659 2.16E-16 SHANK2 SH3 and multiple ankyrin repeat domains 2 2.62 7.80E-08 NEGR1 neuronal growth regulator 1 2.388 2.80E-03 AKAP6 A-kinase anchoring protein 6 2.26 7.91E-04 AR androgen receptor 2.15 7.05E-02 MAP 2 microtubule associated protein 2 2.116 3.90E-08 NCAM2 neural cell adhesion molecule 2 2.091 2.00E-03 NLGN1 neuroligin 1 2.074 1.92E-04 ADGRL3 adhesion G protein-coupled receptor L3 2.049 1.37E-03 PRKG1 protein kinase cGMP-dependent 1 1.976 7.03E-05 PDE11A phosphodiesterase 11A 1.895 1.30E-04 FAM78B family with sequence similarity 78 member B 1.705 1.30E-04 PLXDC2 plexin domain containing 2 1.685 3.61E-11 SEMA3D semaphorin 3D 1.657 2.45E-06 ID1 inhibitor of DNA binding 1 1.637 3.42E-03 The STRING tool was used to identify potential connections between DEGs in different tissues. 12 To build PPI networks, active interaction sources such as databases, co-expression, text mining, experiments, and species restricted to “Homo sapiens” were used, along with an interaction score greater than 0.4. The PPI network was displayed using Cytoscape v3.6.1 software as depicted in Figure 2 . Figure 2. Protein–protein interaction network where Androgen receptor (AR) is the central hub gene. Table 3. List of compounds used for Tox prediction. Compound name Docking score Predicted LD 50 Hepatotoxicity Carcinogenicity Immunotoxicity Mutagenicity Cytotoxicity Chrysin 7-O-beta-D-glucopyranuronoside -7.2 5000 mg/kg Inactive Inactive Inactive Inactive Inactive 0.73 0.51 0.96 0.74 0.81 Atalantoflavone -7.1 2570 mg/kg Inactive Inactive Inactive Inactive Inactive 0.77 0.5 0.51 0.62 0.83 8-Prenyldaidzein -6.8 2500 mg/kg Inactive Inactive Inactive Inactive - 0.7 0.66 0.8 0.65 6-Prenylnaringenin -6.7 2000 mg/kg Inactive Inactive Inactive Inactive Inactive 0.69 0.69 0.5 0.64 0.79 alpha-Isowighteone -6.6 2875 mg/kg Inactive Inactive Inactive Inactive Inactive 0.71 0.61 0.85 0.55 0.81 2-Hydroxynaringenin -6.5 2000 mg/kg Inactive Inactive Inactive Inactive Inactive 0.71 0.57 0.8 0.77 0.55 Carpachromene -6.5 4000 mg/kg Inactive Inactive Inactive Inactive Inactive 0.77 0.5 0.61 0.62 0.83 8-Demethyleucalyptin -6.3 3919 mg/kg Inactive Inactive Inactive Inactive Inactive 0.71 0.54 0.83 0.73 0.93 5-Hydroxy-7-acetoxy-8-methoxyflavone -6.3 5000 mg/kg Inactive Inactive Inactive Inactive Inactive 0.76 0.54 0.87 0.7 0.83 Apigenin -6.3 2500 mg/kg Inactive Inactive Inactive Inactive Inactive 0.68 0.62 0.99 0.57 0.87 The MCode plugin was employed to identify the highly linked regions inside the PPI network, while the CentiScape plugin was utilized to calculate the network topology parameters. Using degree and betweenness as the primary parameters, hub genes were identified. A complete set of algorithms, called CentiScape, was used to analyze the centrality of the network nodes. It can calculate multiple centralities for weighted, directed, and undirected networks. 28 The human Androgen Receptor was determined to be an appropriate hub gene in the protein-protein interaction network consisting of DEG genes. Virtual screening of phytochemical library The human Androgen Receptor (hAR), covering the C-terminal amino acid residues (1E3G) with the co-crystallized ligand metribolone (R18), consists of 250 amino acid residues arranged in a three-layered α-helical sandwich structure. The ligand-binding pocket is located within the hydrophobic cavity formed by helices. A total of 1358 compounds were initially identified through systematic literature search, and their structures were retrieved from the PubChem database. Of these, only 543 compounds met the criteria outlined by Lipinski’s rule of five. These 543 compounds were then selected for the initial virtual screening against human Androgen Receptor using PyRx, and their binding affinities were tabulated 29 (refer to extended data Table S1). The top 50 ranked compounds were subjected to ADMET analysis on the ProTox II server. Only the top 10 ranked compounds that showed favorable binding affinity towards hAR based on their docking interaction and ideal ADMET properties were chosen for further analysis. The initial docking results and ADMET properties are shown in extended data. Induced fit docking and the molecular interactions Molecular interaction studies of the binding cavity of the human Androgen Receptor and molecules are listed in extended data This was compared with the co-crystallized ligand associated with hAR protein R18 and analyzed by Schrodinger-induced fit docking. The ligand 2-hydroxynaringenin demonstrated high affinity for flexible residues within the binding pocket of the Human Androgen receptor protein. The calculated free energy of binding (ΔG) was determined to be -8.59 kcal/mol, indicating a strong binding interaction. While couple of other molecules 8-Prenyldaidzein and 5-Hydroxy-7-acetoxy-8-methoxyflavone also exhibited significant binding with hAR having ΔG = -8.54 kcal/mol and -8.26 kcal/mol, respectively. The highest affinity with a low negative binding energy was observed for 2-hydroxynaringenin, where the ligand formed conventional hydrogen bonds with Leu704, Asn705, Gln711, Met745, Arg752, and Thr877. Leu707, Met780, Leu873, and Phe876 were found to be involved in pi-alkyl and alkyl interactions with the 2-Hydroxynaringenin ligand. The binding energies of 2-Hydroxynaringenin and protein-ligand interactions are displayed in Figure 4 and the binding energies of other molecules are depicted in extended data. Figure 3. Role of Androgen receptor (Source modified from Ref. 30 ). Figure 4. Induced fit docking pose of the ligand (A) 2-Hydroxynaringenin and co-crystallized (B) R18 molecules with hAR (PDB ID: 1E3G) displaying the ribbon shaped 3D protein and ligand interaction, 3D image of binding cavity residues and 2D interaction profile of bidning cavity residues with the respective ligands. Molecular dynamics simulation studies Molecular dynamics simulation (MD) investigations were performed to ascertain the convergence and stability of 1E3G-Apo (no ligand hAR protein), 1E3G+R18 (R18 co-crystallized ligand) and 1E3G+2-Hydroxynaringenin complexes. When comparing the root mean square deviation (RMSD) measurements, the 100 ns simulation showed a stable conformation. The Apo protein’s Cα-backbone’s RMSD showed a 3.0 Å divergence ( Figure 5A ). While 1E3G+R18 and 1E3G+2-Hydroxynaringenin both showed 2.9 Å, the overall RMSD is shown to be 2.9 Å ( Figure 5A ). The root mean square fluctuations (RMSF) plot of the 1E3G+2-Hydroxynaringenin complex protein revealed notable variations at residues 60–70, 110–120, and 180–185, which may have been caused by the residues’ increased flexibility. The rest of the residues fluctuated less during the course of the 100 ns simulation ( Figure 5B ). Radius of gyration (Rg) in this study, 1E3G Cα-backbone bound to Apo protein displayed increment of Rg values indicating lesser compactness while stable Rg was observed from 20.2 to 20.3 Å in 1E3G+R18 ( Figure 5C ). The number of hydrogen bonds was significantly different between 1E3G+2-Hydroxynaringenin, throughout the simulation time of 100 ns ( Figure 5D ). The average number of hydrogen bonds observed in 1E3G+2-Hydroxynaringenin was two on average in MD simulation studies ( Figure 5D , red color). Figure 5. MD simulation analysis of 100 ns trajectories of (A) Cα backbone RMSD of 1E3G+2-Hydroxynaringenin (red), RMSD of 1E3GApo (black), and 1E3G+R18 (blue) (B) RMSF of Cα backbone RMSD of 1E3G+2-Hydroxynaringenin (red), RMSD of 1E3GApo (black), and 1E3G+R18 (blue) (C) Radius of gyration (Rg) of Cα backbone of Cα backbone RMSD of 1E3G+2-Hydroxynaringenin (red), RMSD of 1E3GApo (black), and 1E3G+R18 (blue) (D) Formation of hydrogen bonds in 1E3G+2-hydroxynaringenin (red) and R18 (black). Mechanics generalized born surface area (MM-GBSA) calculations Utilizing the MD simulation trajectory, the binding free energy along with other contributing energy in form of MM-GBSA were determined for hAR+2-hydroxynaringenin. The results ( Table 4 ) suggested that the maximum contribution to ΔG bind in the stability of the simulated complexes were due to ΔG bind Coulomb, ΔG bind vdW and ΔG bind Lipo, while, ΔG bind Covalent and ΔG bind SolvGB contributed to the instability of the corresponding complexes. The hAR+2-hydroxynaringenin complex showed significantly higher binding free energies. These results supported the potential of 2-hydroxynaringenin, showed the efficiency in binding to the selected protein and the ability to form stable protein-ligand complexes. Table 4. Binding free energy components for the 1E3G+2-hydroxynaringenin and 1E3G+R18 calculated by MM-GBSA. Energies (kcal/mol) 1E3G+2-hydroxynaringenin 1E3G+R18 ΔG bind -31.53±5.3 -29.95±4.1 ΔG bind Lipo -29.83±3.2 -23.51±3.2 ΔG bind vdW -22.68±3.22 -16.27±1.21 ΔG bind Coulomb -5.22±2.11 -7.45±2.8 ΔG bind H bond -0.9±0.1 -0.6±0.2 ΔG bind SolvGB 33.91±1.27 41.27±1.76 ΔG bind Covalent 0.79±0.3 1.24±0.23 Discussion The integrated analysis of gene expression and protein-protein interactions (PPI) would help to identify candidates that could serve as therapeutic targets. In this study, we compared TNBC datasets to normal datasets to assess the underlying molecular pathways that drive TNBC. Differential gene expression profiling of the selected datasets using the Benjamini-Hochberg-Yekutieli approach was used to adjust the P-value, which controls the rate of false discovery under positive dependence assumptions. Then, using STRING, which incorporates both known and anticipated PPIs, the protein-protein interactions between the previously mentioned genes were investigated using Cytoscape. CentiScape was used to analyze the centrality of network nodes, and the Human Androgen Receptor was determined to be an appropriate hub gene in the protein-protein interaction network consisting of DEG genes. The Androgen Receptor (AR) pathway is becoming a viable therapeutic target in breast cancer. 31 12-55% of TNBC cases, which provides a chance for targeted treatment. The “Luminal AR (LAR)” molecular subtype of TNBC is where AR is most prevalent. 32 The LAR subtype exhibits the highest amount of AR expression amongst the many molecular subtypes of TNBC in which it is present. All AR+ TNBC primary tumors that were evaluated showed nuclear localization of AR, a sign of transcriptionally active receptors. Many investigations have shown that AR expression in breast cancer, particularly in the TNBC subtype, has been linked to an overall better outcome. Considering that > 70% of AR expression is consistent between primary and metastatic breast cancers, AR may be a novel diagnostic and therapeutic target for patients with AR-positive breast cancer. 33 In luminal mammary carcinomas, a high percentage of cases express androgen receptors (AR), and the ratio of AR to estrogen receptors (ER) or progesterone receptors (PR) is considered a potential prognostic factor. However, in estrogen receptor-negative (ER-) tumors, AR expression is associated with a poorer prognosis. Androgen receptor (AR) expression has demonstrated predictive value for potential response to adjuvant hormonal therapy in estrogen receptor-positive (ER+) breast cancers. Additionally, AR expression has been associated with predicting responses to neoadjuvant chemotherapy in triple-negative breast cancer (TNBC). The role of the AR is shown in Figure 3 . The human Androgen Receptor, which has 920 amino acid residues, was identified as the primary therapeutic target for TNBC. The 3D structure of Androgen Receptor (PDB ID: 1E3G) does not cover the entire protein and contains only the amino acids 671-920. This region encompasses the nuclear receptor ligand-binding domain (NR LBD) of human AR. It consists of 250 amino acid residues, arranged in a three-layered α-helical sandwich structure. The ligand-binding pocket is located within the hydrophobic cavity formed by helices. Virtual screening of 543 ligands against human AR was performed using PyRx at the co-crystallized ligand-binding site. The top 10 ranked compounds that showed favorable binding affinity towards human AR and ideal ADMET properties were chosen for induced fit docking. Unlike rigid docking, induced fit docking treats the ligand and protein as typically flexible entities allowing for conformational changes to occur upon binding. The ligand 2-hydroxynaringenin demonstrated a high affinity for the flexible residues within the binding pocket of human AR, with an interaction binding energy of-8.59 kcal/mol with six conventional hydrogen bonds, indicating a strong binding interaction. Interestingly, the interaction binding energy of the human AR protein with R18 was observed to be -7.8 kcal/mol and only one conventional hydrogen bond formed between R18 and Arg752 ( Figure 4 ). No other potential interactions were observed, except for van der Waal’s instructions. For both 2-Hydroxynaringenin and R18, it was observed that Arg752 is the key residue for ligand binding and could play an active role in protein function. Molecular dynamics (MD) simulation studies of 100 ns showed stable conformations with 1E3G+2-Hydroxynaringenin complexes. The RMSD of the Cα-backbone of the Apo protein exhibited a deviation of 3.0 Å. While 1E3G+R18 exhibited 2.9 Å and simlarly 1E3G+2-Hydroxynaringenin also exhibited the total RMSD is depicted to be 2.9 Å ( Figure 5A ). All RMSD values were below the acceptable range of 3 Å. 34 Stable RMSD plots of apo-1E3G, 1E3G+R18 and 1E3G+2-Hydroxynaringenin were observed to be less than 3 Å. Therefore, it can be suggested that apo-1E3G, 1E3G+R18 and 1E3G+2-Hydroxynaringenin complexes are well converged and equilibrated. The RMSF of the 1E3G+2-Hydroxynaringenin complex protein exhibited notable fluctuation spikes at residues 60–70, 110–120, and 180–185, which may have been brought on by the residues’ increased flexibility. During the course of the 100 ns simulation, the remaining residues fluctuated less. A more rigid conformation with fewer fluctuations was observed in the Apo-protein and 1E3G+R18 complex. Therefore, from the RMSF plots, it can be suggested that the structures of 1E3G+2-Hydroxynaringenin are more flexible during simulation in ligand-bound conformations. The radius of gyration (Rg) is a measure of protein compactness. Lowering and stable of radius of gyration (Rg) from 20.0 to 20.02 Å in 1E3G+2-Hydroxynaringenin was observed. The quantity of hydrogen bonds forming between the ligand and protein indicates a strong connection and stability of the complex. Over the course of the 100 ns simulation, there was a considerable difference in the amount of hydrogen bonds between 1E3G+2-Hydroxynaringenin ( Figure 5D ). The average number of hydrogen bonds observed in 1E3G+2-Hydroxynaringenin was two on average in MD simulation studies ( Figure 5D , red). Using the MD simulation trajectory, the binding free energy and additional contributing energies in the form of MM-GBSA were found for hAR+2-hydroxynaringenin. The findings ( Table 4 ) show that ΔGbindCoulomb, ΔGbindvdW, and ΔGbindLipo were the main contributors to ΔGbind in the simulated complexes’ stability, whereas ΔGbindCovalent and ΔGbindSolvGB were responsible for the corresponding complexes’ instability. hAR+2-hydroxynaringenin complex showed significantly higher binding free energies. The capacity of 2-hydroxynaringenin to bind to the chosen protein efficiently and form stable protein-ligand complexes was demonstrated by these data, which further validated the compound’s potential. Conclusion In recent years, bioinformatic analysis has become essential for studying the pathogenesis of human diseases. Differential gene expression studies, protein–protein interactions, and network topology analyses were performed. The current study identified the human Androgen Receptor (AR) as a potential drug target to combat TNBC. This was concluded based on gene expression profiling, protein-protein interaction, and network topology analysis. The specific role of the Androgen Receptor in breast cancer growth and progression remains uncertain, although the AR is expressed in approximately 77% of all breast cancers, even higher than Estrogen Receptors (ERs). 31 A more luminal, well-differentiated, and less aggressive tumor may be indicated by high expression of Androgen Receptor in breast cancer, which could improve prognosis. 32 AR inhibition tends to be well-tolerated, and patients with TNBC may benefit from it when paired with other medications, as its toxicity is much lower than that of chemotherapy. Combinations involving mTOR inhibitors, EGFR and other ErbB inhibitors, PIK3 inhibitors, anti-PDL1 antibodies, paclitaxel, and other chemotherapeutic drugs are supported by preclinical results. Randomized clinical trials would be required to ascertain the clinical utility of AR inhibitors. 32 , 36 , 37 Flavonoids are a class of natural compounds found in various fruits, vegetables, and plants and have been extensively studied for their potential therapeutic effects, including their ability to combat cancer. Naringenin, specifically categorized as a flavanone, is a flavonoid present in grapefruit and tomatoes, among other dietary sources. 38 The antioxidant and anti-inflammatory properties of naringenin have led to its exploration for various potential use in the pharmaceutical industry. 39 Limitations of the study The current study identified the human Androgen Receptor as a potential candidate drug target to combat TNBC and recognized 2-hydroxynaringenin as a potential lead molecule. The in vitro and in vivo efficacies of 2-hydroxynaringenin require further investigation. Safety, pharmacokinetics, and pharmacodynamics tests need to be performed to further develop hydroxynaringenin for clinical use. Data availability statement Underlying data 1. GEO DATASET 1 - Accession number- GSE45498 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE45498 Platform–GPL16299 2. GEO DATASET 2 - Accession number- GSE214101 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE214101 Extended data Supplementary data 1. Figshare: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer - Supplementary Figures.docx - DOI: 10.6084/m9.figshare.26967880.v1 38 Data are available under the terms of the Creative Commons Zero “No rights reserved” data waiver (CC0 1.0 Public domain dedication). 2. Figshare: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer - Supplementary Table.docx - DOI: 10.6084/m9.figshare.26967733.v1 39 Data are available under the terms of the Creative Commons Zero “No rights reserved” data waiver (CC0 1.0 Public domain dedication). References 1. Arnold M, Morgan E, Rumgay H, et al. : Current and future burden of breast cancer: Global statistics for 2020 and 2040. The Breast. 2022 Dec 1; 66 : 15–23. PubMed Abstract | Publisher Full Text | Free Full Text 2. Bergin AR, Loi S: Triple-negative breast cancer: recent treatment advances. F1000Research. 2019; 8 : 1342. PubMed Abstract | Publisher Full Text | Free Full Text 3. Kalimutho M, Parsons K, Mittal D, et al. : Targeted therapies for triple-negative breast cancer: combating a stubborn disease. Trends in Pharmacological Sciences. 2015 Dec 1; 36 (12): 822–846. PubMed Abstract | Publisher Full Text 4. Bianchini G, Balko JM, Mayer IA, et al. : Triple-negative breast cancer: challenges and opportunities of a heterogeneous disease. Nature Reviews Clinical Oncology. 2016 Nov; 13 (11): 674–690. PubMed Abstract | Publisher Full Text | Free Full Text 5. Ramadan E, Alinsaif S, Hassan MR: Network topology measures for identifying disease-gene association in breast cancer. BMC Bioinformatics. 2016 Jul; 17 : 473–480. Publisher Full Text 6. Price PD, Palmer Droguett DH, Taylor JA, et al. : Detecting signatures of selection on gene expression. Nature Ecology & Evolution. 2022 Jul; 6 (7): 1035–1045. Publisher Full Text 7. Hozhabri H, Ghasemi Dehkohneh RS, Razavi SM, et al. : Comparative analysis of protein-protein interaction networks in metastatic breast cancer. PLoS One. 2022 Jan 19; 17 (1): e0260584. PubMed Abstract | Publisher Full Text | Free Full Text 8. Israel BE, Tilghman SL, Parker-Lemieux K, et al. : Phytochemicals: Current strategies for treating breast cancer. Oncology Letters. 2018 May 1; 15 (5): 7471–7478. PubMed Abstract | Publisher Full Text 9. Choudhari AS, Mandave PC, Deshpande M, et al. : Phytochemicals in cancer treatment: From preclinical studies to clinical practice. Frontiers in Pharmacology. 2020 Jan 28; 10 : 497776. Publisher Full Text 10. Pranaya S, Ragunath PK, Venkatesan P: Diagnosis of triple negative breast cancer using expression data with several machine learning tools. Bioinformation. 2022; 18 (4): 325–330. PubMed Abstract | Publisher Full Text | Free Full Text 11. Szklarczyk D, Gable AL, Lyon D, et al. : STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Research. 2019 Jan 8; 47 (D1): D607–D613. PubMed Abstract | Publisher Full Text | Free Full Text 12. Szklarczyk D, Gable AL, Nastou KC, et al. : The STRING database in 2021: customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Research. 2021 Jan 8; 49 (D1): D605–D612. PubMed Abstract | Publisher Full Text | Free Full Text 13. Dashti S, Taheri M, Ghafouri-Fard S: An in-silico method leads to recognition of hub genes and crucial pathways in survival of patients with breast cancer. Scientific Reports. 2020 Oct 30; 10 (1): 18770. PubMed Abstract | Publisher Full Text | Free Full Text 14. Kawale MA: AN OVERERVIEW OF MAJOR CLASSES OF PHYTOCHEMICALS: THEIR TYPES AND ROLE IN DISEASE PREVENTION.2018 [cited 2023 Dec 7]. Reference Source 15. Shukla S, Mehta A: Anticancer potential of medicinal plants and their phytochemicals: a review. Brazilian Journal of Botany. 2015 Jun; 38 : 199–210. Publisher Full Text 16. Dallakyan S, Olson AJ: Small-molecule library screening by docking with PyRx. Chemical Biology: Methods and Protocols. 2015; 1263 : 243–250. PubMed Abstract | Publisher Full Text 17. Trott O, Olson AJ: AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. Journal of Computational Chemistry. 2010 Jan 30; 31 (2): 455–461. PubMed Abstract | Publisher Full Text | Free Full Text 18. Wang Y, Xing J, Xu Y, et al. : In silico ADME/T modelling for rational drug design. Quarterly Reviews of Biophysics. 2015 Nov; 48 (4): 488–515. PubMed Abstract | Publisher Full Text 19. Banerjee P, Eckert AO, Schrey AK, et al. : ProTox-II: a webserver for the prediction of toxicity of chemicals. Nucleic Acids Ressearch. 2018 Jul 2 [cited 2023 Jan 9]; 46 (W1): W257–W263. PubMed Abstract | Publisher Full Text | Free Full Text Reference Source 20. Shaw DE, Maragakis P, Lindorff-Larsen K, et al. : Atomic-level characterization of the structural dynamics of proteins. Science. 2010 Oct 15; 330 (6002): 341–346. Publisher Full Text 21. Shivakumar D, Williams J, Wu Y, et al. : Prediction of absolute solvation free energies using molecular dynamics free energy perturbation and the OPLS force field. Journal of Chemical Theory and Computation. 2010 May 11; 6 (5): 1509–1519. PubMed Abstract | Publisher Full Text 22. Jorgensen WL, Chandrasekhar J, Madura JD, et al. : Comparison of simple potential functions for simulating liquid water. The Journal of Chemical Physics. 1983 Jul 15; 79 (2): 926–935. 23. Martyna GJ, Tobias DJ, Klein ML: Constant pressure molecular dynamics algorithms. The Journal of Chemical Physics. 1994 Sep 1; 101 (5): 4177–4189. Publisher Full Text 24. Martyna GJ, Klein ML, Tuckerman M: Nosé–Hoover chains: The canonical ensemble via continuous dynamics. The Journal of Chemical Physics. 1992 Aug 15; 97 (4): 2635–2643. Publisher Full Text 25. Toukmaji AY, Board JA Jr: Ewald summation techniques in perspective: a survey. Computer Physics Communications. 1996 Jun 1; 95 (2-3): 73–92. Publisher Full Text 26. Kollman PA, Massova I, Reyes C, et al. : Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. Accounts of Chemical Research. 2000 Dec 19; 33 (12): 889–897. PubMed Abstract | Publisher Full Text 27. Wang E, Sun H, Wang J, et al. : End-point binding free energy calculation with MM/PBSA and MM/GBSA: strategies and applications in drug design. Chemical Reviews. 2019 Jun 24; 119 (16): 9478–9508. PubMed Abstract | Publisher Full Text 28. Scardoni G, Tosadori G, Faizan M, et al. : Biological network analysis with CentiScaPe: centralities and experimental dataset integration. F1000Research. 2014; 3 : 3. Publisher Full Text 29. Prabhakar L: Meta-analysis of lean and obese RNA-seq datasets to identify genes targeting obesity. Bioinformation. 2023; 19 (3): 331–335. PubMed Abstract | Publisher Full Text | Free Full Text 30. Shorning BY, Dass MS, Smalley MJ, et al. : The PI3K-AKT-mTOR Pathway and Prostate Cancer: At the Crossroads of AR, MAPK, and WNT Signaling. International Journal of Molecular Sciences. 2020; 21 : 4507. Publisher Full Text 31. Anestis A, Zoi I, Papavassiliou AG, et al. : Androgen receptor in breast cancer—clinical and preclinical research insights. Molecules. 2020 Jan 15; 25 (2): 358. PubMed Abstract | Publisher Full Text | Free Full Text 32. Barton VN, D’Amato NC, Gordon MA, et al. : Multiple molecular subtypes of triple-negative breast cancer critically rely on androgen receptor and respond to enzalutamide in vivo. Molecular Cancer Therapeutics. 2015 Mar 1; 14 (3): 769–778. PubMed Abstract | Publisher Full Text | Free Full Text 33. Vasiliou SK, Diamandis EP: Androgen receptor: A promising therapeutic target in breast cancer. Critical Reviews in Clinical Laboratory Sciences. 2019 Apr 3 [cited 2025 Feb 24]; 56 (3): 200–223. PubMed Abstract | Publisher Full Text 34. Ramírez D, Caballero J: Is It Reliable to Take the Molecular Docking Top Scoring Position as the Best Solution without Considering Available Structural Data? Molecules: A Journal of Synthetic Chemistry and Natural Product Chemistry [Internet]. 2018 [cited 2025 Feb 24]; 23 (5): 1038. PubMed Abstract | Publisher Full Text | Free Full Text 35. Barton VN, Gordon MA, Richer JK, et al. : Anti-androgen therapy in triple-negative breast cancer. Therapeutic Advances in Medical Oncology. 2016 Jul; 8 (4): 305–308. PubMed Abstract | Publisher Full Text | Free Full Text 36. Kriegsmann M, Endris V, Wolf T, et al. : Mutational profiles in triple-negative breast cancer defined by ultradeep multigene sequencing show high rates of PI3K pathway alterations and clinically relevant entity subgroup specific differences. Oncotarget. 2014 Oct; 5 (20): 9952–9965. PubMed Abstract | Publisher Full Text | Free Full Text 37. Tung NM, Garber JE, Torous V, et al. : Prevalence and predictors of androgen receptor (AR) and programmed death-ligand 1 (PD-L1) expression in BRCA1-associated and sporadic triple negative breast cancer (TNBC). 38. Muralidharan S, Velanganni AA, Shanmugam K: Inhibition of Breast Cancer Proteins by the Flavonoid Naringenin and its Derivative: A Molecular Docking Study. Journal of Natural Remedies. 2022 Jan 22; 51–64. Publisher Full Text 39. Zhao Z, Jin G, Ge Y, et al. : Naringenin inhibits migration of breast cancer cells via inflammatory and apoptosis cell signaling pathways. Inflammopharmacology. 2019 Oct; 27 : 1021–1036. PubMed Abstract | Publisher Full Text 40. Sankaranarayanan P, John DGD, Abhinand PA, et al. : Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer - Supplementary Figures.docx. figshare. 2024. Reference Source 41. Sankaranarayanan P, John DGD, Abhinand PA, et al. : Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer - Supplementary Table.docx. figshare. 2024. Reference Source Comments on this article Comments (0) Version 2 VERSION 2 PUBLISHED 24 Oct 2024 ADD YOUR COMMENT Comment Author details Author details 1 Department of Bioinformatics, Sri Ramachandra Institute of Higher Education and Research (Deemed to be University), Chennai, Tamil Nadu, 600116, India 2 Department of Medical Genetics, Manipal Hospitals, Bengaluru, Karnataka, 560 017, India 3 Department of Botany, Gauhati University, Guwahati, Assam, India Pranaya Sankaranarayanan Roles: Conceptualization, Data Curation, Formal Analysis, Methodology, Writing – Original Draft Preparation Dicky John Davis G Roles: Investigation, Project Administration, Supervision, Validation, Visualization, Writing – Review & Editing Abhinand PA Roles: Conceptualization, Formal Analysis, Software, Validation, Writing – Review & Editing M Manikandan Roles: Data Curation, Formal Analysis, Software, Visualization Arabinda Ghosh Roles: Formal Analysis, Resources, Software, Writing – Original Draft Preparation Competing interests No competing interests were disclosed. Grant information This work was supported by the Founder-Chancellor Shri. N. P. V. Ramasamy Udayar Research Fellowship (U02B160480), Sri Ramachandra Institute of Higher Education and Research. The funders had no role in the study design, data collection and analysis, decision to publish, or manuscript preparation. Article Versions (2) version 2 Revised Published: 18 Mar 2025, 13:1271 https://doi.org/10.12688/f1000research.155657.2 version 1 Published: 24 Oct 2024, 13:1271 https://doi.org/10.12688/f1000research.155657.1 Copyright © 2025 Sankaranarayanan P et al . This is an open access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Download Export To Sciwheel Bibtex EndNote ProCite Ref. Manager (RIS) Sente metrics Views Downloads F1000Research - - PubMed Central info_outline Data from PMC are received and updated monthly. - - Citations open_in_new 0 open_in_new 0 open_in_new SEE MORE DETAILS CITE how to cite this article Sankaranarayanan P, G DJD, PA A et al. Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.12688/f1000research.155657.2 ) NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS track receive updates on this article Track an article to receive email alerts on any updates to this article. TRACK THIS ARTICLE Share Open Peer Review Current Reviewer Status: ? Key to Reviewer Statuses VIEW HIDE Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions Version 2 VERSION 2 PUBLISHED 18 Mar 2025 Revised Views 0 Cite How to cite this report: Ahmed SSSJ. Reviewer Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.178643.r371506 ) The direct URL for this report is: https://f1000research.com/articles/13-1271/v2#referee-response-371506 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 24 Aug 2025 Shiek S S J Ahmed , Chettinad Academy of Research and Education, Kelambakkam 603103, Tamil Nadu,, India Approved VIEWS 0 https://doi.org/10.5256/f1000research.178643.r371506 The author responded to all comments. If there is no negative ... Continue reading READ ALL The author responded to all comments. If there is no negative report from other reviewers, it can be accepted Competing Interests: No competing interests were disclosed. Reviewer Expertise: Neuroscience, Cancer research, immunoinformatic, systems biology, NGS, artificial intelligence, omics research, Big data analytics I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Ahmed SSSJ. Reviewer Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.178643.r371506 ) The direct URL for this report is: https://f1000research.com/articles/13-1271/v2#referee-response-371506 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Respond or Comment COMMENT ON THIS REPORT Views 0 Cite How to cite this report: Muthiah I. Reviewer Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.178643.r393020 ) The direct URL for this report is: https://f1000research.com/articles/13-1271/v2#referee-response-393020 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 25 Jul 2025 Indiraleka Muthiah , Mepco Schlenk Engineering College, Sivakasi, Tamil Nadu, India Approved with Reservations VIEWS 0 https://doi.org/10.5256/f1000research.178643.r393020 The research titled “Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer” effectively presents an in silico strategy to identify potential phytochemical ligands targeting TNBC molecular markers. The study concludes that the ... Continue reading READ ALL The research titled “Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer” effectively presents an in silico strategy to identify potential phytochemical ligands targeting TNBC molecular markers. The study concludes that the androgen receptor (AR) is a promising target and highlights 2-Hydroxynaringenin as a potential lead molecule for AR. The manuscript is well-written, with a clearly defined methodology that effectively addresses the research objective. However, the following suggestions could enhance the manuscript further: The significance of the current research should be emphasised in the background section of the abstract to convey its impact and relevance. Numerous molecular targets and phytochemical ligands for TNBC have been explored through in silico studies. In the introduction, it would strengthen the manuscript to differentiate this study from existing research by outlining its novelty and contribution, supported by up-to-date literature. Additionally, comparing the current results with previously reported studies would provide valuable context. The ligand selection process requires more detailed explanation, including the specific criteria and rationale for choosing the phytochemicals investigated. Since many phytochemicals have been studied against breast cancer, clarifying why only a limited subset was selected here will add clarity. The docking results section could benefit from identifying and discussing the binding sites, particularly in comparison to standard AR inhibitors. The explanation of molecular docking outcomes is currently insufficient, making it challenging to interpret and compare the ligands presented in Table 3. Clarification on whether the ligands are ranked by binding affinity and the reasoning behind selecting 2-Hydroxynaringenin for further study should be included. The discussion section should be expanded to provide a more comprehensive interpretation of the findings, leading to clearer conclusions. It is important to note that while this in silico approach can propose potential lead molecules for AR, confirmation of their efficacy and potential requires further validation through in vitro and in vivo assays. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? I cannot comment. A qualified statistician is required. Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly Competing Interests: No competing interests were disclosed. Reviewer Expertise: Research on the design and delivery of drugs for TNBC I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Muthiah I. Reviewer Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.178643.r393020 ) The direct URL for this report is: https://f1000research.com/articles/13-1271/v2#referee-response-393020 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Respond or Comment COMMENT ON THIS REPORT Views 0 Cite How to cite this report: Zou Y. Reviewer Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.178643.r372224 ) The direct URL for this report is: https://f1000research.com/articles/13-1271/v2#referee-response-372224 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 25 Apr 2025 Yutian Zou , Sun Yat-sen University Cancer Center, Guangzhou, China Approved with Reservations VIEWS 0 https://doi.org/10.5256/f1000research.178643.r372224 The manuscript titled “Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer” presents an in silico approach to identify novel therapeutic candidates for triple-negative breast cancer (TNBC). Using differential gene expression analysis ... Continue reading READ ALL The manuscript titled “Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer” presents an in silico approach to identify novel therapeutic candidates for triple-negative breast cancer (TNBC). Using differential gene expression analysis and protein-protein interaction networks, the androgen receptor (AR) was identified as a potential target. Subsequent virtual screening and molecular dynamics simulations suggested that 2-hydroxynaringenin may serve as a promising phytochemical inhibitor of AR. The study demonstrates the structural stability and favorable binding affinity of the 2-hydroxynaringenin–AR complex through MM-GBSA calculations. However, the following issues are required for explaining: The computational identification of 2-hydroxynaringenin as a potential AR inhibitor is interesting, but it is essential to experimentally validate the binding affinity and functional relevance of this interaction in vitro. Cell-based assays using AR-expressing TNBC cell lines should be conducted to confirm the direct binding and biological effect of 2-hydroxynaringenin. The inhibitory effect of 2-hydroxynaringenin on AR activity should be quantified using standard biochemical assays to determine its IC₅₀ value. This quantitative measurement is crucial to evaluate its potency compared to known AR inhibitors. The authors should consider validating their docking and simulation workflow by applying it to known AR inhibitors such as enzalutamide, bicalutamide, and apalutamide. Comparing the binding affinity and stability of these compounds with 2-hydroxynaringenin using the same computational pipeline is important to demonstrate the predictive power of the proposed approach. The manuscript would benefit from a more detailed discussion on how the computational findings could inform future drug development. Specifically, how can the identified phytochemical be optimized or modified for enhanced efficacy and bioavailability? What are the steps for translating this lead compound into preclinical or clinical testing? Essential details of data and analysis in this study remain unclear. For example, there is a lack of details on how docking data were processed and normalized? A reproducible detail should be provided. Some similar studies regarding the docking and cancer should be cited and discussed. For example, PMID: 39912365, 39081282. The current study lacks a discussion on the limitations of in silico-only studies. The authors should acknowledge the absence of experimental validation and outline clear next steps, including in vitro assays, pharmacokinetic evaluation, and possible structural modification of the lead compound to improve drug-like properties. The authors are recommended to consider engaging a professional language editing service to ensure the clarity and coherence of the manuscript. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Partly If applicable, is the statistical analysis and its interpretation appropriate? Yes Are all the source data underlying the results available to ensure full reproducibility? Partly Are the conclusions drawn adequately supported by the results? Partly References 1. Sultan S, Tawfik SS, Selim KB, Nasr MNA: Design, Synthesis and Molecular Docking of New Thieno[2,3‑d]Pyrimidin-4-One Derivatives as Dual EGFR and FGFR Inhibitors. Drug Dev Res . 2025; 86 (1): e70061 PubMed Abstract | Publisher Full Text 2. Dong XD, Lu Q, Li YD, Cai CY, et al.: RN486, a Bruton's Tyrosine Kinase inhibitor, antagonizes multidrug resistance in ABCG2-overexpressing cancer cells. J Transl Int Med . 2024; 12 (3): 288-298 PubMed Abstract | Publisher Full Text Competing Interests: No competing interests were disclosed. Reviewer Expertise: breast cancer I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Zou Y. Reviewer Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.178643.r372224 ) The direct URL for this report is: https://f1000research.com/articles/13-1271/v2#referee-response-372224 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Respond or Comment COMMENT ON THIS REPORT Version 1 VERSION 1 PUBLISHED 24 Oct 2024 Views 0 Cite How to cite this report: Ahmed SSSJ. Reviewer Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.170851.r335168 ) The direct URL for this report is: https://f1000research.com/articles/13-1271/v1#referee-response-335168 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 10 Dec 2024 Shiek S S J Ahmed , Chettinad Academy of Research and Education, Kelambakkam 603103, Tamil Nadu,, India Approved with Reservations VIEWS 0 https://doi.org/10.5256/f1000research.170851.r335168 This manuscript provides valuable insights into drug discovery in the field of oncology. The authors have done excellent work in exploring breast cancer treatment by utilizing efficient computational methods such as differential gene expression, protein network, molecular docking, molecular dynamics ... Continue reading READ ALL This manuscript provides valuable insights into drug discovery in the field of oncology. The authors have done excellent work in exploring breast cancer treatment by utilizing efficient computational methods such as differential gene expression, protein network, molecular docking, molecular dynamics simulations, and MM-GBSA. The study aims to identify potential lead compounds against a candidate target for triple-negative breast cancer. However, the following issues need to be addressed before indexing. 1. Abstract needs to reframed { ( #Note : Can be re-constructed based to authors preference) Example: Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, high aggressiveness, and has a relatively poor prognosis and clinical outcome. Objective To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified targets. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. Results The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo.} please refer to following link for more details https://f1000research.s3.amazonaws.com/linked/695061.155657-F1000_review_1_.docx 2. Please do not use scientific term multiple times. rather use the abbreviation. 3. In the methodology section, the calculation of MD trajectories is unclear. Was the Desmond plug-in employed, or manually? Additionally, there appears to be a contrast between the described MD simulation methodology and the reported results. The MD simulation results might have been obtained using GROMACS rather than Desmond. This assumption arises because Desmond typically calculates parameters such as the Rg and intramolecular hydrogen bonds (intraHB) for the ligand's extendness and internal hydrogen bonding, but not for the entire complex. 4. The authors should clearly define the criteria used for selecting those two datasets (gene expression) in the methodology section. The datasets seems heterogeneous one consists of human sample (expression array) while the other consists of cell line samples (RNA-sequencing). How the authors address and account for batch effects arising from these differences? 4a. The basis for the adopting the LogFC > 1.25 as a threshold while in gene expression analysis should be explained 5. How the 2D interaction of protein and ligand was visualized after docking. 6. In Fig. 5B, the RMSF does not cover the entire protein (263 amino acids). 7. Some important references are missing in the discussion section, particularly in AR association with Breast cancer (second paragraph). 8. In Fig 2, highlighting the hub target protein or constructing the sub-network by centering AR would enhance visualization and helpful for the new readers to understand. 9. The following sentence needs referencing “All RMSD values were below the acceptable range of 3Å”. 10. The interpretation of the MM-GBSA results needs to be clarified. Specifically, the outcomes of the MM-GBSA analysis comparing the test and control compounds should be explained in detail. Is the work clearly and accurately presented and does it cite the current literature? Yes Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? I cannot comment. A qualified statistician is required. Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly Competing Interests: No competing interests were disclosed. Reviewer Expertise: Neuroscience, Cancer research, immunoinformatic, systems biology, NGS, artificial intelligence, omics research, Big data analytics I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Ahmed SSSJ. Reviewer Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.170851.r335168 ) The direct URL for this report is: https://f1000research.com/articles/13-1271/v1#referee-response-335168 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Author Response 18 Mar 2025 Pranaya Sankaranarayanan , Department of Bioinformatics, Sri Ramachandra Institute of Higher Education and Research (Deemed to be University), Chennai, 600116, India 18 Mar 2025 Author Response 1. Abstract needs to reframed Response: Thank you for your suggestion. The abstract has been restructured. Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of ... Continue reading 1. Abstract needs to reframed Response: Thank you for your suggestion. The abstract has been restructured. Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, demonstrates high aggressiveness, and has a relatively poor prognosis and clinical outcome. Objective To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified targets. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. Results The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo. 2. Please do not use scientific term multiple times. rather use the abbreviation. Response: Thank you for the input. The changes have been made. 3. In the methodology section, the calculation of MD trajectories is unclear. Was the Desmond plug-in employed, or manually? Additionally, there appears to be a contrast between the described MD simulation methodology and the reported results. The MD simulation results might have been obtained using GROMACS rather than Desmond. This assumption arises because Desmond typically calculates parameters such as the Rg and intramolecular hydrogen bonds (intraHB) for the ligand's extendness and internal hydrogen bonding, but not for the entire complex. Response: In the methodology section it was quite clearly mentioned that Desmond 2020.1 engine was used manually for the study. Moreover, the concern was raised the superiority of GROMACS over Desmond is somehow confusing and as per the suggestion in order to obtain the exact parameters there are options such as simulation event analysis where the hydrogen bonding calculates between the protein and ligand complex. Desmond is also well documented and parameterized simulation engine like GROMACS and there are many reports and previous literatures successfully attributed the specificity of the tool. Fig. 5B, RMSF extracted from the tool outcome and regions having spikes with high and moderate significance exhibited in the figure. The outcome of MMGBSA is compared with the control R18 compound which was co-crytallized inhibitor in the pdb structure of the receptor. 4. The authors should clearly define the criteria used for selecting those two datasets (gene expression) in the methodology section. The datasets seem heterogeneous one consists of human sample (expression array) while the other consists of cell line samples (RNA-sequencing). How the authors address and account for batch effects arising from these differences? Response: We appreciate your insightful comments regarding the dataset selection and batch effect considerations. In our study, we selected the two datasets—GSE45498 (microarray from human tissue samples) and GSE214101 (RNA-seq from cell lines)—based on their relevance to triple-negative breast cancer (TNBC) and their comprehensive gene expression profiles. These datasets were chosen to enable a broader investigation of TNBC-related gene expression variations across different biological contexts. To address potential batch effects arising from differences in platform technology and sample types, we implemented standard normalization and transformation techniques. 4a. The basis for the adopting the LogFC > 1.25 as a threshold while in gene expression analysis should be explained Response: A log₂FC of 1.25 corresponds to a 2.4-fold change in expression, meaning the gene expression is increased by ~140% (if positive) or decreased to ~42% of its original level (if negative). Many biological processes, such as regulatory networks and signalling pathways, involve moderate gene expression changes rather than extreme shifts. Traditional thresholds like log₂FC ≥ 2 (4-fold change) are quite conservative and may exclude relevant genes, especially for genes with small but functional changes. A log₂FC threshold of 1.25 provides a balance—capturing biologically significant changes while avoiding detection of genes with minor fluctuations. Maintaining a balance between sensitivity and specificity : Setting too high a threshold (e.g., log₂FC ≥ 2) might miss key regulatory genes with moderate but meaningful expression changes. A threshold of log₂FC ≥ 1.25, combined with an adjusted p-value (e.g., FDR ≤ 0.05), ensures that identified genes are statistically reliable and biologically relevant. Reducing false positives while capturing meaningful expression shifts : Log₂FC values <1 (e.g., 0.5, which is ~1.4-fold change) might capture noise or small fluctuations due to experimental variability. Using 1.25 as a threshold ensures changes are more likely due to biological regulation rather than technical variation. 5. How the 2D interaction of protein and ligand was visualized after docking Response: By using discovery studio, the 2D interaction of protein and ligand was visualized. 6. In Fig. 5B, the RMSF does not cover the entire protein (263 amino acids). Response: The 3D structure of Androgen Receptor (PDB ID: 1E3G) does not cover the entire protein and contains only the amino acids 671-920 (comprising of 250 amino acids). The manuscript also has been modified to reflect this. 7. Some important references are missing in the discussion section, particularly in AR association with Breast cancer (second paragraph). Response: References have been added 8. In Fig 2, highlighting the hub target protein or constructing the sub-network by centering AR would enhance visualization and helpful for the new readers to understand. Response: Thank you for your valuable suggestion. We have incorporated the changes. https://f1000research.s3.amazonaws.com/linked/707287.155657-Comment_to_reviewer_Sheik_by_Author_Pranaya.pdf 9. The following sentence needs referencing “All RMSD values were below the acceptable range of 3Å”. Response: Reference has been added 10. The interpretation of the MM-GBSA results needs to be clarified. Specifically, the outcomes of the MM-GBSA analysis comparing the test and control compounds should be explained in detail. Response: Utilizing the MD simulation trajectory, the binding free energy along with other contributing energy in form of MM-GBSA were determined for HAR+2-hydroxynaringenin. The results (Table 4) suggested that the maximum contribution to ΔG bind in the stability of the simulated complexes were due to ΔG bind Coulomb, ΔG bind vdW and ΔG bind Lipo, while, ΔG bind Covalent and ΔG bind SolvGB contributed to the instability of the corresponding complexes. The HAR+2-hydroxynaringenin complex showed significantly higher binding free energies. These results supported the potential of 2-hydroxynaringenin, showed the efficiency in binding to the selected protein and the ability to form stable protein-ligand complexes. Binding free energy components for the 1E3G+2-hydroxynaringenin and 1E3G+R18 calculated by MM-GBSA. Table is attached in the manuscript. 1. Abstract needs to reframed Response: Thank you for your suggestion. The abstract has been restructured. Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, demonstrates high aggressiveness, and has a relatively poor prognosis and clinical outcome. Objective To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified targets. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. Results The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo. 2. Please do not use scientific term multiple times. rather use the abbreviation. Response: Thank you for the input. The changes have been made. 3. In the methodology section, the calculation of MD trajectories is unclear. Was the Desmond plug-in employed, or manually? Additionally, there appears to be a contrast between the described MD simulation methodology and the reported results. The MD simulation results might have been obtained using GROMACS rather than Desmond. This assumption arises because Desmond typically calculates parameters such as the Rg and intramolecular hydrogen bonds (intraHB) for the ligand's extendness and internal hydrogen bonding, but not for the entire complex. Response: In the methodology section it was quite clearly mentioned that Desmond 2020.1 engine was used manually for the study. Moreover, the concern was raised the superiority of GROMACS over Desmond is somehow confusing and as per the suggestion in order to obtain the exact parameters there are options such as simulation event analysis where the hydrogen bonding calculates between the protein and ligand complex. Desmond is also well documented and parameterized simulation engine like GROMACS and there are many reports and previous literatures successfully attributed the specificity of the tool. Fig. 5B, RMSF extracted from the tool outcome and regions having spikes with high and moderate significance exhibited in the figure. The outcome of MMGBSA is compared with the control R18 compound which was co-crytallized inhibitor in the pdb structure of the receptor. 4. The authors should clearly define the criteria used for selecting those two datasets (gene expression) in the methodology section. The datasets seem heterogeneous one consists of human sample (expression array) while the other consists of cell line samples (RNA-sequencing). How the authors address and account for batch effects arising from these differences? Response: We appreciate your insightful comments regarding the dataset selection and batch effect considerations. In our study, we selected the two datasets—GSE45498 (microarray from human tissue samples) and GSE214101 (RNA-seq from cell lines)—based on their relevance to triple-negative breast cancer (TNBC) and their comprehensive gene expression profiles. These datasets were chosen to enable a broader investigation of TNBC-related gene expression variations across different biological contexts. To address potential batch effects arising from differences in platform technology and sample types, we implemented standard normalization and transformation techniques. 4a. The basis for the adopting the LogFC > 1.25 as a threshold while in gene expression analysis should be explained Response: A log₂FC of 1.25 corresponds to a 2.4-fold change in expression, meaning the gene expression is increased by ~140% (if positive) or decreased to ~42% of its original level (if negative). Many biological processes, such as regulatory networks and signalling pathways, involve moderate gene expression changes rather than extreme shifts. Traditional thresholds like log₂FC ≥ 2 (4-fold change) are quite conservative and may exclude relevant genes, especially for genes with small but functional changes. A log₂FC threshold of 1.25 provides a balance—capturing biologically significant changes while avoiding detection of genes with minor fluctuations. Maintaining a balance between sensitivity and specificity : Setting too high a threshold (e.g., log₂FC ≥ 2) might miss key regulatory genes with moderate but meaningful expression changes. A threshold of log₂FC ≥ 1.25, combined with an adjusted p-value (e.g., FDR ≤ 0.05), ensures that identified genes are statistically reliable and biologically relevant. Reducing false positives while capturing meaningful expression shifts : Log₂FC values <1 (e.g., 0.5, which is ~1.4-fold change) might capture noise or small fluctuations due to experimental variability. Using 1.25 as a threshold ensures changes are more likely due to biological regulation rather than technical variation. 5. How the 2D interaction of protein and ligand was visualized after docking Response: By using discovery studio, the 2D interaction of protein and ligand was visualized. 6. In Fig. 5B, the RMSF does not cover the entire protein (263 amino acids). Response: The 3D structure of Androgen Receptor (PDB ID: 1E3G) does not cover the entire protein and contains only the amino acids 671-920 (comprising of 250 amino acids). The manuscript also has been modified to reflect this. 7. Some important references are missing in the discussion section, particularly in AR association with Breast cancer (second paragraph). Response: References have been added 8. In Fig 2, highlighting the hub target protein or constructing the sub-network by centering AR would enhance visualization and helpful for the new readers to understand. Response: Thank you for your valuable suggestion. We have incorporated the changes. https://f1000research.s3.amazonaws.com/linked/707287.155657-Comment_to_reviewer_Sheik_by_Author_Pranaya.pdf 9. The following sentence needs referencing “All RMSD values were below the acceptable range of 3Å”. Response: Reference has been added 10. The interpretation of the MM-GBSA results needs to be clarified. Specifically, the outcomes of the MM-GBSA analysis comparing the test and control compounds should be explained in detail. Response: Utilizing the MD simulation trajectory, the binding free energy along with other contributing energy in form of MM-GBSA were determined for HAR+2-hydroxynaringenin. The results (Table 4) suggested that the maximum contribution to ΔG bind in the stability of the simulated complexes were due to ΔG bind Coulomb, ΔG bind vdW and ΔG bind Lipo, while, ΔG bind Covalent and ΔG bind SolvGB contributed to the instability of the corresponding complexes. The HAR+2-hydroxynaringenin complex showed significantly higher binding free energies. These results supported the potential of 2-hydroxynaringenin, showed the efficiency in binding to the selected protein and the ability to form stable protein-ligand complexes. Binding free energy components for the 1E3G+2-hydroxynaringenin and 1E3G+R18 calculated by MM-GBSA. Table is attached in the manuscript. Competing Interests: No competing interests were disclosed. Close Report a concern Respond or Comment COMMENTS ON THIS REPORT Author Response 18 Mar 2025 Pranaya Sankaranarayanan , Department of Bioinformatics, Sri Ramachandra Institute of Higher Education and Research (Deemed to be University), Chennai, 600116, India 18 Mar 2025 Author Response 1. Abstract needs to reframed Response: Thank you for your suggestion. The abstract has been restructured. Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of ... Continue reading 1. Abstract needs to reframed Response: Thank you for your suggestion. The abstract has been restructured. Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, demonstrates high aggressiveness, and has a relatively poor prognosis and clinical outcome. Objective To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified targets. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. Results The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo. 2. Please do not use scientific term multiple times. rather use the abbreviation. Response: Thank you for the input. The changes have been made. 3. In the methodology section, the calculation of MD trajectories is unclear. Was the Desmond plug-in employed, or manually? Additionally, there appears to be a contrast between the described MD simulation methodology and the reported results. The MD simulation results might have been obtained using GROMACS rather than Desmond. This assumption arises because Desmond typically calculates parameters such as the Rg and intramolecular hydrogen bonds (intraHB) for the ligand's extendness and internal hydrogen bonding, but not for the entire complex. Response: In the methodology section it was quite clearly mentioned that Desmond 2020.1 engine was used manually for the study. Moreover, the concern was raised the superiority of GROMACS over Desmond is somehow confusing and as per the suggestion in order to obtain the exact parameters there are options such as simulation event analysis where the hydrogen bonding calculates between the protein and ligand complex. Desmond is also well documented and parameterized simulation engine like GROMACS and there are many reports and previous literatures successfully attributed the specificity of the tool. Fig. 5B, RMSF extracted from the tool outcome and regions having spikes with high and moderate significance exhibited in the figure. The outcome of MMGBSA is compared with the control R18 compound which was co-crytallized inhibitor in the pdb structure of the receptor. 4. The authors should clearly define the criteria used for selecting those two datasets (gene expression) in the methodology section. The datasets seem heterogeneous one consists of human sample (expression array) while the other consists of cell line samples (RNA-sequencing). How the authors address and account for batch effects arising from these differences? Response: We appreciate your insightful comments regarding the dataset selection and batch effect considerations. In our study, we selected the two datasets—GSE45498 (microarray from human tissue samples) and GSE214101 (RNA-seq from cell lines)—based on their relevance to triple-negative breast cancer (TNBC) and their comprehensive gene expression profiles. These datasets were chosen to enable a broader investigation of TNBC-related gene expression variations across different biological contexts. To address potential batch effects arising from differences in platform technology and sample types, we implemented standard normalization and transformation techniques. 4a. The basis for the adopting the LogFC > 1.25 as a threshold while in gene expression analysis should be explained Response: A log₂FC of 1.25 corresponds to a 2.4-fold change in expression, meaning the gene expression is increased by ~140% (if positive) or decreased to ~42% of its original level (if negative). Many biological processes, such as regulatory networks and signalling pathways, involve moderate gene expression changes rather than extreme shifts. Traditional thresholds like log₂FC ≥ 2 (4-fold change) are quite conservative and may exclude relevant genes, especially for genes with small but functional changes. A log₂FC threshold of 1.25 provides a balance—capturing biologically significant changes while avoiding detection of genes with minor fluctuations. Maintaining a balance between sensitivity and specificity : Setting too high a threshold (e.g., log₂FC ≥ 2) might miss key regulatory genes with moderate but meaningful expression changes. A threshold of log₂FC ≥ 1.25, combined with an adjusted p-value (e.g., FDR ≤ 0.05), ensures that identified genes are statistically reliable and biologically relevant. Reducing false positives while capturing meaningful expression shifts : Log₂FC values <1 (e.g., 0.5, which is ~1.4-fold change) might capture noise or small fluctuations due to experimental variability. Using 1.25 as a threshold ensures changes are more likely due to biological regulation rather than technical variation. 5. How the 2D interaction of protein and ligand was visualized after docking Response: By using discovery studio, the 2D interaction of protein and ligand was visualized. 6. In Fig. 5B, the RMSF does not cover the entire protein (263 amino acids). Response: The 3D structure of Androgen Receptor (PDB ID: 1E3G) does not cover the entire protein and contains only the amino acids 671-920 (comprising of 250 amino acids). The manuscript also has been modified to reflect this. 7. Some important references are missing in the discussion section, particularly in AR association with Breast cancer (second paragraph). Response: References have been added 8. In Fig 2, highlighting the hub target protein or constructing the sub-network by centering AR would enhance visualization and helpful for the new readers to understand. Response: Thank you for your valuable suggestion. We have incorporated the changes. https://f1000research.s3.amazonaws.com/linked/707287.155657-Comment_to_reviewer_Sheik_by_Author_Pranaya.pdf 9. The following sentence needs referencing “All RMSD values were below the acceptable range of 3Å”. Response: Reference has been added 10. The interpretation of the MM-GBSA results needs to be clarified. Specifically, the outcomes of the MM-GBSA analysis comparing the test and control compounds should be explained in detail. Response: Utilizing the MD simulation trajectory, the binding free energy along with other contributing energy in form of MM-GBSA were determined for HAR+2-hydroxynaringenin. The results (Table 4) suggested that the maximum contribution to ΔG bind in the stability of the simulated complexes were due to ΔG bind Coulomb, ΔG bind vdW and ΔG bind Lipo, while, ΔG bind Covalent and ΔG bind SolvGB contributed to the instability of the corresponding complexes. The HAR+2-hydroxynaringenin complex showed significantly higher binding free energies. These results supported the potential of 2-hydroxynaringenin, showed the efficiency in binding to the selected protein and the ability to form stable protein-ligand complexes. Binding free energy components for the 1E3G+2-hydroxynaringenin and 1E3G+R18 calculated by MM-GBSA. Table is attached in the manuscript. 1. Abstract needs to reframed Response: Thank you for your suggestion. The abstract has been restructured. Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, demonstrates high aggressiveness, and has a relatively poor prognosis and clinical outcome. Objective To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified targets. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. Results The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo. 2. Please do not use scientific term multiple times. rather use the abbreviation. Response: Thank you for the input. The changes have been made. 3. In the methodology section, the calculation of MD trajectories is unclear. Was the Desmond plug-in employed, or manually? Additionally, there appears to be a contrast between the described MD simulation methodology and the reported results. The MD simulation results might have been obtained using GROMACS rather than Desmond. This assumption arises because Desmond typically calculates parameters such as the Rg and intramolecular hydrogen bonds (intraHB) for the ligand's extendness and internal hydrogen bonding, but not for the entire complex. Response: In the methodology section it was quite clearly mentioned that Desmond 2020.1 engine was used manually for the study. Moreover, the concern was raised the superiority of GROMACS over Desmond is somehow confusing and as per the suggestion in order to obtain the exact parameters there are options such as simulation event analysis where the hydrogen bonding calculates between the protein and ligand complex. Desmond is also well documented and parameterized simulation engine like GROMACS and there are many reports and previous literatures successfully attributed the specificity of the tool. Fig. 5B, RMSF extracted from the tool outcome and regions having spikes with high and moderate significance exhibited in the figure. The outcome of MMGBSA is compared with the control R18 compound which was co-crytallized inhibitor in the pdb structure of the receptor. 4. The authors should clearly define the criteria used for selecting those two datasets (gene expression) in the methodology section. The datasets seem heterogeneous one consists of human sample (expression array) while the other consists of cell line samples (RNA-sequencing). How the authors address and account for batch effects arising from these differences? Response: We appreciate your insightful comments regarding the dataset selection and batch effect considerations. In our study, we selected the two datasets—GSE45498 (microarray from human tissue samples) and GSE214101 (RNA-seq from cell lines)—based on their relevance to triple-negative breast cancer (TNBC) and their comprehensive gene expression profiles. These datasets were chosen to enable a broader investigation of TNBC-related gene expression variations across different biological contexts. To address potential batch effects arising from differences in platform technology and sample types, we implemented standard normalization and transformation techniques. 4a. The basis for the adopting the LogFC > 1.25 as a threshold while in gene expression analysis should be explained Response: A log₂FC of 1.25 corresponds to a 2.4-fold change in expression, meaning the gene expression is increased by ~140% (if positive) or decreased to ~42% of its original level (if negative). Many biological processes, such as regulatory networks and signalling pathways, involve moderate gene expression changes rather than extreme shifts. Traditional thresholds like log₂FC ≥ 2 (4-fold change) are quite conservative and may exclude relevant genes, especially for genes with small but functional changes. A log₂FC threshold of 1.25 provides a balance—capturing biologically significant changes while avoiding detection of genes with minor fluctuations. Maintaining a balance between sensitivity and specificity : Setting too high a threshold (e.g., log₂FC ≥ 2) might miss key regulatory genes with moderate but meaningful expression changes. A threshold of log₂FC ≥ 1.25, combined with an adjusted p-value (e.g., FDR ≤ 0.05), ensures that identified genes are statistically reliable and biologically relevant. Reducing false positives while capturing meaningful expression shifts : Log₂FC values <1 (e.g., 0.5, which is ~1.4-fold change) might capture noise or small fluctuations due to experimental variability. Using 1.25 as a threshold ensures changes are more likely due to biological regulation rather than technical variation. 5. How the 2D interaction of protein and ligand was visualized after docking Response: By using discovery studio, the 2D interaction of protein and ligand was visualized. 6. In Fig. 5B, the RMSF does not cover the entire protein (263 amino acids). Response: The 3D structure of Androgen Receptor (PDB ID: 1E3G) does not cover the entire protein and contains only the amino acids 671-920 (comprising of 250 amino acids). The manuscript also has been modified to reflect this. 7. Some important references are missing in the discussion section, particularly in AR association with Breast cancer (second paragraph). Response: References have been added 8. In Fig 2, highlighting the hub target protein or constructing the sub-network by centering AR would enhance visualization and helpful for the new readers to understand. Response: Thank you for your valuable suggestion. We have incorporated the changes. https://f1000research.s3.amazonaws.com/linked/707287.155657-Comment_to_reviewer_Sheik_by_Author_Pranaya.pdf 9. The following sentence needs referencing “All RMSD values were below the acceptable range of 3Å”. Response: Reference has been added 10. The interpretation of the MM-GBSA results needs to be clarified. Specifically, the outcomes of the MM-GBSA analysis comparing the test and control compounds should be explained in detail. Response: Utilizing the MD simulation trajectory, the binding free energy along with other contributing energy in form of MM-GBSA were determined for HAR+2-hydroxynaringenin. The results (Table 4) suggested that the maximum contribution to ΔG bind in the stability of the simulated complexes were due to ΔG bind Coulomb, ΔG bind vdW and ΔG bind Lipo, while, ΔG bind Covalent and ΔG bind SolvGB contributed to the instability of the corresponding complexes. The HAR+2-hydroxynaringenin complex showed significantly higher binding free energies. These results supported the potential of 2-hydroxynaringenin, showed the efficiency in binding to the selected protein and the ability to form stable protein-ligand complexes. Binding free energy components for the 1E3G+2-hydroxynaringenin and 1E3G+R18 calculated by MM-GBSA. Table is attached in the manuscript. Competing Interests: No competing interests were disclosed. Close Report a concern COMMENT ON THIS REPORT Comments on this article Comments (0) Version 2 VERSION 2 PUBLISHED 24 Oct 2024 ADD YOUR COMMENT Comment keyboard_arrow_left keyboard_arrow_right Open Peer Review Reviewer Status info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions Reviewer Reports Invited Reviewers 1 2 3 Version 2 (revision) 18 Mar 25 read read read Version 1 24 Oct 24 read Shiek S S J Ahmed , Chettinad Academy of Research and Education, Kelambakkam 603103, Tamil Nadu,, India Yutian Zou , Sun Yat-sen University Cancer Center, Guangzhou, China Indiraleka Muthiah , Mepco Schlenk Engineering College, Sivakasi, India Comments on this article All Comments (0) Add a comment Sign up for content alerts Sign Up You are now signed up to receive this alert Browse by related subjects keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2025 Ahmed S. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 24 Aug 2025 | for Version 2 Shiek S S J Ahmed , Chettinad Academy of Research and Education, Kelambakkam 603103, Tamil Nadu,, India 0 Views copyright © 2025 Ahmed S. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (0) Approved info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions The author responded to all comments. If there is no negative report from other reviewers, it can be accepted Competing Interests No competing interests were disclosed. Reviewer Expertise Neuroscience, Cancer research, immunoinformatic, systems biology, NGS, artificial intelligence, omics research, Big data analytics I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. reply Respond to this report Responses (0) Ahmed SSSJ. Peer Review Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.178643.r371506) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/13-1271/v2#referee-response-371506 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2025 Muthiah I. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 25 Jul 2025 | for Version 2 Indiraleka Muthiah , Mepco Schlenk Engineering College, Sivakasi, Tamil Nadu, India 0 Views copyright © 2025 Muthiah I. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (0) Approved With Reservations info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions The research titled “Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer” effectively presents an in silico strategy to identify potential phytochemical ligands targeting TNBC molecular markers. The study concludes that the androgen receptor (AR) is a promising target and highlights 2-Hydroxynaringenin as a potential lead molecule for AR. The manuscript is well-written, with a clearly defined methodology that effectively addresses the research objective. However, the following suggestions could enhance the manuscript further: The significance of the current research should be emphasised in the background section of the abstract to convey its impact and relevance. Numerous molecular targets and phytochemical ligands for TNBC have been explored through in silico studies. In the introduction, it would strengthen the manuscript to differentiate this study from existing research by outlining its novelty and contribution, supported by up-to-date literature. Additionally, comparing the current results with previously reported studies would provide valuable context. The ligand selection process requires more detailed explanation, including the specific criteria and rationale for choosing the phytochemicals investigated. Since many phytochemicals have been studied against breast cancer, clarifying why only a limited subset was selected here will add clarity. The docking results section could benefit from identifying and discussing the binding sites, particularly in comparison to standard AR inhibitors. The explanation of molecular docking outcomes is currently insufficient, making it challenging to interpret and compare the ligands presented in Table 3. Clarification on whether the ligands are ranked by binding affinity and the reasoning behind selecting 2-Hydroxynaringenin for further study should be included. The discussion section should be expanded to provide a more comprehensive interpretation of the findings, leading to clearer conclusions. It is important to note that while this in silico approach can propose potential lead molecules for AR, confirmation of their efficacy and potential requires further validation through in vitro and in vivo assays. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? I cannot comment. A qualified statistician is required. Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly Competing Interests No competing interests were disclosed. Reviewer Expertise Research on the design and delivery of drugs for TNBC I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above. reply Respond to this report Responses (0) Muthiah I. Peer Review Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.178643.r393020) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/13-1271/v2#referee-response-393020 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2025 Zou Y. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 25 Apr 2025 | for Version 2 Yutian Zou , Sun Yat-sen University Cancer Center, Guangzhou, China 0 Views copyright © 2025 Zou Y. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (0) Approved With Reservations info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions The manuscript titled “Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer” presents an in silico approach to identify novel therapeutic candidates for triple-negative breast cancer (TNBC). Using differential gene expression analysis and protein-protein interaction networks, the androgen receptor (AR) was identified as a potential target. Subsequent virtual screening and molecular dynamics simulations suggested that 2-hydroxynaringenin may serve as a promising phytochemical inhibitor of AR. The study demonstrates the structural stability and favorable binding affinity of the 2-hydroxynaringenin–AR complex through MM-GBSA calculations. However, the following issues are required for explaining: The computational identification of 2-hydroxynaringenin as a potential AR inhibitor is interesting, but it is essential to experimentally validate the binding affinity and functional relevance of this interaction in vitro. Cell-based assays using AR-expressing TNBC cell lines should be conducted to confirm the direct binding and biological effect of 2-hydroxynaringenin. The inhibitory effect of 2-hydroxynaringenin on AR activity should be quantified using standard biochemical assays to determine its IC₅₀ value. This quantitative measurement is crucial to evaluate its potency compared to known AR inhibitors. The authors should consider validating their docking and simulation workflow by applying it to known AR inhibitors such as enzalutamide, bicalutamide, and apalutamide. Comparing the binding affinity and stability of these compounds with 2-hydroxynaringenin using the same computational pipeline is important to demonstrate the predictive power of the proposed approach. The manuscript would benefit from a more detailed discussion on how the computational findings could inform future drug development. Specifically, how can the identified phytochemical be optimized or modified for enhanced efficacy and bioavailability? What are the steps for translating this lead compound into preclinical or clinical testing? Essential details of data and analysis in this study remain unclear. For example, there is a lack of details on how docking data were processed and normalized? A reproducible detail should be provided. Some similar studies regarding the docking and cancer should be cited and discussed. For example, PMID: 39912365, 39081282. The current study lacks a discussion on the limitations of in silico-only studies. The authors should acknowledge the absence of experimental validation and outline clear next steps, including in vitro assays, pharmacokinetic evaluation, and possible structural modification of the lead compound to improve drug-like properties. The authors are recommended to consider engaging a professional language editing service to ensure the clarity and coherence of the manuscript. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Partly If applicable, is the statistical analysis and its interpretation appropriate? Yes Are all the source data underlying the results available to ensure full reproducibility? Partly Are the conclusions drawn adequately supported by the results? Partly References 1. Sultan S, Tawfik SS, Selim KB, Nasr MNA: Design, Synthesis and Molecular Docking of New Thieno[2,3‑d]Pyrimidin-4-One Derivatives as Dual EGFR and FGFR Inhibitors. Drug Dev Res . 2025; 86 (1): e70061 PubMed Abstract | Publisher Full Text 2. Dong XD, Lu Q, Li YD, Cai CY, et al.: RN486, a Bruton's Tyrosine Kinase inhibitor, antagonizes multidrug resistance in ABCG2-overexpressing cancer cells. J Transl Int Med . 2024; 12 (3): 288-298 PubMed Abstract | Publisher Full Text Competing Interests No competing interests were disclosed. Reviewer Expertise breast cancer I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above. reply Respond to this report Responses (0) Zou Y. Peer Review Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.178643.r372224) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/13-1271/v2#referee-response-372224 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2024 Ahmed S. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 10 Dec 2024 | for Version 1 Shiek S S J Ahmed , Chettinad Academy of Research and Education, Kelambakkam 603103, Tamil Nadu,, India 0 Views copyright © 2024 Ahmed S. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (1) Approved With Reservations info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions This manuscript provides valuable insights into drug discovery in the field of oncology. The authors have done excellent work in exploring breast cancer treatment by utilizing efficient computational methods such as differential gene expression, protein network, molecular docking, molecular dynamics simulations, and MM-GBSA. The study aims to identify potential lead compounds against a candidate target for triple-negative breast cancer. However, the following issues need to be addressed before indexing. 1. Abstract needs to reframed { ( #Note : Can be re-constructed based to authors preference) Example: Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, high aggressiveness, and has a relatively poor prognosis and clinical outcome. Objective To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified targets. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. Results The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo.} please refer to following link for more details https://f1000research.s3.amazonaws.com/linked/695061.155657-F1000_review_1_.docx 2. Please do not use scientific term multiple times. rather use the abbreviation. 3. In the methodology section, the calculation of MD trajectories is unclear. Was the Desmond plug-in employed, or manually? Additionally, there appears to be a contrast between the described MD simulation methodology and the reported results. The MD simulation results might have been obtained using GROMACS rather than Desmond. This assumption arises because Desmond typically calculates parameters such as the Rg and intramolecular hydrogen bonds (intraHB) for the ligand's extendness and internal hydrogen bonding, but not for the entire complex. 4. The authors should clearly define the criteria used for selecting those two datasets (gene expression) in the methodology section. The datasets seems heterogeneous one consists of human sample (expression array) while the other consists of cell line samples (RNA-sequencing). How the authors address and account for batch effects arising from these differences? 4a. The basis for the adopting the LogFC > 1.25 as a threshold while in gene expression analysis should be explained 5. How the 2D interaction of protein and ligand was visualized after docking. 6. In Fig. 5B, the RMSF does not cover the entire protein (263 amino acids). 7. Some important references are missing in the discussion section, particularly in AR association with Breast cancer (second paragraph). 8. In Fig 2, highlighting the hub target protein or constructing the sub-network by centering AR would enhance visualization and helpful for the new readers to understand. 9. The following sentence needs referencing “All RMSD values were below the acceptable range of 3Å”. 10. The interpretation of the MM-GBSA results needs to be clarified. Specifically, the outcomes of the MM-GBSA analysis comparing the test and control compounds should be explained in detail. Is the work clearly and accurately presented and does it cite the current literature? Yes Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? I cannot comment. A qualified statistician is required. Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly Competing Interests No competing interests were disclosed. Reviewer Expertise Neuroscience, Cancer research, immunoinformatic, systems biology, NGS, artificial intelligence, omics research, Big data analytics I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above. reply Respond to this report Responses (1) Author Response 18 Mar 2025 Pranaya Sankaranarayanan , Department of Bioinformatics, Sri Ramachandra Institute of Higher Education and Research (Deemed to be University), Chennai, 600116, India 1. Abstract needs to reframed Response: Thank you for your suggestion. The abstract has been restructured. Background Triple-negative breast cancers (TNBC) are defined as tumors that lack the expression of the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). It exhibits unique clinical and pathological features, demonstrates high aggressiveness, and has a relatively poor prognosis and clinical outcome. Objective To identify a novel drug target protein against TNBC and potential phytochemical lead molecules against the identified targets. Methods In this study, we retrieved TNBC samples from NGS and microarray datasets in the Gene Expression Omnibus database. We employed a combination of differential gene expression studies, protein-protein interaction analysis, and network topology investigation to identify the target protein. Additionally, the molecular docking and molecular dynamics (MD) simulation studies followed by Molecular Mechanics with Generalised Born Surface Area salvation was used to identify potential lead molecule. Results The upregulated genes with LogFC > 1.25 and P-value < 0.05 from the TNBC gene expression dataset were identified. Androgen receptor (AR) was found to be an appropriate hub target in the protein-protein interaction network. Phytochemicals that inhibit breast cancer target were retrieved from the PubChem database and virtual screening was performed using PyRx against the AR protein. Thereby, the AR was found to be the target protein and 2-hydroxynaringenin was discovered to be a possible phytochemical lead molecule for combating TNBC. Moreover, the AR and the 2-hydroxynaringenin complex showed structural stability and higher binding affinity through molecular dynamics and MM-GBSA studies. Conclusion AR was identified as a hub protein that is highly expressed in breast cancer and 2-hydroxynaringenin efficacy of counter TNBC requires further investigation both in vitro and in vivo. 2. Please do not use scientific term multiple times. rather use the abbreviation. Response: Thank you for the input. The changes have been made. 3. In the methodology section, the calculation of MD trajectories is unclear. Was the Desmond plug-in employed, or manually? Additionally, there appears to be a contrast between the described MD simulation methodology and the reported results. The MD simulation results might have been obtained using GROMACS rather than Desmond. This assumption arises because Desmond typically calculates parameters such as the Rg and intramolecular hydrogen bonds (intraHB) for the ligand's extendness and internal hydrogen bonding, but not for the entire complex. Response: In the methodology section it was quite clearly mentioned that Desmond 2020.1 engine was used manually for the study. Moreover, the concern was raised the superiority of GROMACS over Desmond is somehow confusing and as per the suggestion in order to obtain the exact parameters there are options such as simulation event analysis where the hydrogen bonding calculates between the protein and ligand complex. Desmond is also well documented and parameterized simulation engine like GROMACS and there are many reports and previous literatures successfully attributed the specificity of the tool. Fig. 5B, RMSF extracted from the tool outcome and regions having spikes with high and moderate significance exhibited in the figure. The outcome of MMGBSA is compared with the control R18 compound which was co-crytallized inhibitor in the pdb structure of the receptor. 4. The authors should clearly define the criteria used for selecting those two datasets (gene expression) in the methodology section. The datasets seem heterogeneous one consists of human sample (expression array) while the other consists of cell line samples (RNA-sequencing). How the authors address and account for batch effects arising from these differences? Response: We appreciate your insightful comments regarding the dataset selection and batch effect considerations. In our study, we selected the two datasets—GSE45498 (microarray from human tissue samples) and GSE214101 (RNA-seq from cell lines)—based on their relevance to triple-negative breast cancer (TNBC) and their comprehensive gene expression profiles. These datasets were chosen to enable a broader investigation of TNBC-related gene expression variations across different biological contexts. To address potential batch effects arising from differences in platform technology and sample types, we implemented standard normalization and transformation techniques. 4a. The basis for the adopting the LogFC > 1.25 as a threshold while in gene expression analysis should be explained Response: A log₂FC of 1.25 corresponds to a 2.4-fold change in expression, meaning the gene expression is increased by ~140% (if positive) or decreased to ~42% of its original level (if negative). Many biological processes, such as regulatory networks and signalling pathways, involve moderate gene expression changes rather than extreme shifts. Traditional thresholds like log₂FC ≥ 2 (4-fold change) are quite conservative and may exclude relevant genes, especially for genes with small but functional changes. A log₂FC threshold of 1.25 provides a balance—capturing biologically significant changes while avoiding detection of genes with minor fluctuations. Maintaining a balance between sensitivity and specificity : Setting too high a threshold (e.g., log₂FC ≥ 2) might miss key regulatory genes with moderate but meaningful expression changes. A threshold of log₂FC ≥ 1.25, combined with an adjusted p-value (e.g., FDR ≤ 0.05), ensures that identified genes are statistically reliable and biologically relevant. Reducing false positives while capturing meaningful expression shifts : Log₂FC values <1 (e.g., 0.5, which is ~1.4-fold change) might capture noise or small fluctuations due to experimental variability. Using 1.25 as a threshold ensures changes are more likely due to biological regulation rather than technical variation. 5. How the 2D interaction of protein and ligand was visualized after docking Response: By using discovery studio, the 2D interaction of protein and ligand was visualized. 6. In Fig. 5B, the RMSF does not cover the entire protein (263 amino acids). Response: The 3D structure of Androgen Receptor (PDB ID: 1E3G) does not cover the entire protein and contains only the amino acids 671-920 (comprising of 250 amino acids). The manuscript also has been modified to reflect this. 7. Some important references are missing in the discussion section, particularly in AR association with Breast cancer (second paragraph). Response: References have been added 8. In Fig 2, highlighting the hub target protein or constructing the sub-network by centering AR would enhance visualization and helpful for the new readers to understand. Response: Thank you for your valuable suggestion. We have incorporated the changes. https://f1000research.s3.amazonaws.com/linked/707287.155657-Comment_to_reviewer_Sheik_by_Author_Pranaya.pdf 9. The following sentence needs referencing “All RMSD values were below the acceptable range of 3Å”. Response: Reference has been added 10. The interpretation of the MM-GBSA results needs to be clarified. Specifically, the outcomes of the MM-GBSA analysis comparing the test and control compounds should be explained in detail. Response: Utilizing the MD simulation trajectory, the binding free energy along with other contributing energy in form of MM-GBSA were determined for HAR+2-hydroxynaringenin. The results (Table 4) suggested that the maximum contribution to ΔG bind in the stability of the simulated complexes were due to ΔG bind Coulomb, ΔG bind vdW and ΔG bind Lipo, while, ΔG bind Covalent and ΔG bind SolvGB contributed to the instability of the corresponding complexes. The HAR+2-hydroxynaringenin complex showed significantly higher binding free energies. These results supported the potential of 2-hydroxynaringenin, showed the efficiency in binding to the selected protein and the ability to form stable protein-ligand complexes. Binding free energy components for the 1E3G+2-hydroxynaringenin and 1E3G+R18 calculated by MM-GBSA. Table is attached in the manuscript. View more View less Competing Interests No competing interests were disclosed. reply Respond Report a concern Ahmed SSSJ. Peer Review Report For: Molecular docking and MD simulation approach to identify potential phytochemical lead molecule against triple negative breast cancer [version 2; peer review: 1 approved, 2 approved with reservations] . F1000Research 2025, 13 :1271 ( https://doi.org/10.5256/f1000research.170851.r335168) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/13-1271/v1#referee-response-335168 Alongside their report, reviewers assign a status to the article: Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions Adjust parameters to alter display View on desktop for interactive features Includes Interactive Elements View on desktop for interactive features Competing Interests Policy Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list: Examples of 'Non-Financial Competing Interests' Within the past 4 years, you have held joint grants, published or collaborated with any of the authors of the selected paper. You have a close personal relationship (e.g. parent, spouse, sibling, or domestic partner) with any of the authors. You are a close professional associate of any of the authors (e.g. scientific mentor, recent student). You work at the same institute as any of the authors. You hope/expect to benefit (e.g. favour or employment) as a result of your submission. You are an Editor for the journal in which the article is published. Examples of 'Financial Competing Interests' You expect to receive, or in the past 4 years have received, any of the following from any commercial organisation that may gain financially from your submission: a salary, fees, funding, reimbursements. You expect to receive, or in the past 4 years have received, shared grant support or other funding with any of the authors. You hold, or are currently applying for, any patents or significant stocks/shares relating to the subject matter of the paper you are commenting on. Stay Updated Sign up for content alerts and receive a weekly or monthly email with all newly published articles Register with F1000Research Already registered? Sign in Not now, thanks close PLEASE NOTE If you are an AUTHOR of this article, please check that you signed in with the account associated with this article otherwise we cannot automatically identify your role as an author and your comment will be labelled as a “User Comment”. If you are a REVIEWER of this article, please check that you have signed in with the account associated with this article and then go to your account to submit your report, please do not post your review here. If you do not have access to your original account, please contact us . All commenters must hold a formal affiliation as per our Policies . The information that you give us will be displayed next to your comment. User comments must be in English, comprehensible and relevant to the article under discussion. We reserve the right to remove any comments that we consider to be inappropriate, offensive or otherwise in breach of the User Comment Terms and Conditions . Commenters must not use a comment for personal attacks. When criticisms of the article are based on unpublished data, the data should be made available. I accept the User Comment Terms and Conditions Please confirm that you accept the User Comment Terms and Conditions. Affiliation ✕ refresh Please enter your institution. Note: To add your institution or organisation, start typing the name and then select the correct name from the list. Where applicable, the name will appear in both the original language and in English. Do not paste in the name. If the name does not appear in the drop-down list, we will display the information you have entered. ✕ refresh Country/Region * USA UK Canada China France Germany Afghanistan Aland Islands Albania Algeria American Samoa Andorra Angola Anguilla Antarctica Antigua and Barbuda Argentina Armenia Aruba Australia Austria Azerbaijan Bahamas Bahrain Bangladesh Barbados Belarus Belgium Belize Benin Bermuda Bhutan Bolivia Bosnia and Herzegovina Botswana Bouvet Island Brazil British Indian Ocean Territory British Virgin Islands Brunei Bulgaria Burkina Faso Burundi Cambodia Cameroon Canada Cape Verde Cayman Islands Central African Republic Chad Chile China Christmas Island Cocos (Keeling) Islands Colombia Comoros Congo Cook Islands Costa Rica Cote d'Ivoire Croatia Cuba Cyprus Czech Republic Democratic Republic of the Congo Denmark Djibouti Dominica Dominican Republic Ecuador Egypt El Salvador Equatorial Guinea Eritrea Estonia Ethiopia Falkland Islands Faroe Islands Federated States of Micronesia Fiji Finland France French Guiana French Polynesia French Southern Territories Gabon Georgia Germany Ghana Gibraltar Greece Greenland Grenada Guadeloupe Guam Guatemala Guernsey Guinea Guinea-Bissau Guyana Haiti Heard Island and Mcdonald Islands Holy See (Vatican City State) Honduras Hong Kong Hungary Iceland India Indonesia Iran Iraq Ireland Israel Italy Jamaica Japan Jersey Jordan Kazakhstan Kenya Kiribati Kosovo (Serbia and Montenegro) Kuwait Kyrgyzstan Lao People's Democratic Republic Latvia Lebanon Lesotho Liberia Libya Liechtenstein Lithuania Luxembourg Macao Madagascar Malawi Malaysia Maldives Mali Malta Marshall Islands Martinique Mauritania Mauritius Mayotte Mexico Minor Outlying Islands of the United States Moldova Monaco Mongolia Montenegro Montserrat Morocco Mozambique Myanmar Namibia Nauru Nepal Netherlands Antilles New Caledonia New Zealand Nicaragua Niger Nigeria Niue Norfolk Island North Korea North Macedonia Northern Mariana Islands Norway Oman Pakistan Palau Palestinian Territory Panama Papua New Guinea Paraguay Peru Philippines Pitcairn Poland Portugal Puerto Rico Qatar Reunion Romania Russian Federation Rwanda Saint Helena Saint Kitts and Nevis Saint Lucia Saint Pierre and Miquelon Saint Vincent and the Grenadines Samoa San Marino Sao Tome and Principe Saudi Arabia Senegal Serbia Seychelles Sierra Leone Singapore Slovakia Slovenia Solomon Islands Somalia South Africa South Georgia and the South Sandwich Is South Korea South Sudan Spain Sri Lanka Sudan Suriname Svalbard and Jan Mayen Swaziland Sweden Switzerland Syria Taiwan Tajikistan Tanzania Thailand The Gambia The Netherlands Timor-Leste Togo Tokelau Tonga Trinidad and Tobago Tunisia Turkey Turkmenistan Turks and Caicos Islands Tuvalu UK USA Uganda Ukraine United Arab Emirates United States Virgin Islands Uruguay Uzbekistan Vanuatu Venezuela Vietnam Wallis and Futuna West Bank and Gaza Strip Western Sahara Yemen Zambia Zimbabwe Please select your country/region. You must enter a comment. Competing Interests Please disclose any competing interests that might be construed to influence your judgment of the article's or peer review report's validity or importance. Competing Interests Policy Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list: Examples of 'Non-Financial Competing Interests' Within the past 4 years, you have held joint grants, published or collaborated with any of the authors of the selected paper. You have a close personal relationship (e.g. parent, spouse, sibling, or domestic partner) with any of the authors. You are a close professional associate of any of the authors (e.g. scientific mentor, recent student). You work at the same institute as any of the authors. You hope/expect to benefit (e.g. favour or employment) as a result of your submission. You are an Editor for the journal in which the article is published. Examples of 'Financial Competing Interests' You expect to receive, or in the past 4 years have received, any of the following from any commercial organisation that may gain financially from your submission: a salary, fees, funding, reimbursements. You expect to receive, or in the past 4 years have received, shared grant support or other funding with any of the authors. You hold, or are currently applying for, any patents or significant stocks/shares relating to the subject matter of the paper you are commenting on. Please state your competing interests The comment has been saved. An error has occurred. Please try again. Cancel Post var lTitle = "Molecular docking and MD simulation approach...".replace("'", ''); var linkedInUrl = "http://www.linkedin.com/shareArticle?url=https://f1000research.com/articles/13-1271/v2" + "&title=" + encodeURIComponent(lTitle) + "&summary=" + encodeURIComponent('Read the article by '); var deliciousUrl = "https://del.icio.us/post?url=https://f1000research.com/articles/13-1271/v2&title=" + encodeURIComponent(lTitle); var redditUrl = "http://reddit.com/submit?url=https://f1000research.com/articles/13-1271/v2" + "&title=" + encodeURIComponent(lTitle); linkedInUrl += encodeURIComponent('Sankaranarayanan P et al.'); var offsetTop = /chrome/i.test( navigator.userAgent ) ? 4 : -10; var addthis_config = { ui_offset_top: offsetTop, services_compact : "facebook,twitter,www.linkedin.com,www.mendeley.com,reddit.com", services_expanded : "facebook,twitter,www.linkedin.com,www.mendeley.com,reddit.com", services_custom : [ { name: "LinkedIn", url: linkedInUrl, icon:"/img/icon/at_linkedin.svg" }, { name: "Mendeley", url: "http://www.mendeley.com/import/?url=https://f1000research.com/articles/13-1271/v2/mendeley", icon:"/img/icon/at_mendeley.svg" }, { name: "Reddit", url: redditUrl, icon:"/img/icon/at_reddit.svg" }, ] }; var addthis_share = { url: "https://f1000research.com/articles/13-1271", templates : { twitter : "Molecular docking and MD simulation approach to identify potential.... Sankaranarayanan P et al., published by " + "@F1000Research" + ", https://f1000research.com/articles/13-1271/v2" } }; if (typeof(addthis) != "undefined"){ addthis.addEventListener('addthis.ready', checkCount); addthis.addEventListener('addthis.menu.share', checkCount); } $(".f1r-shares-twitter").attr("href", "https://twitter.com/intent/tweet?text=" + addthis_share.templates.twitter); $(".f1r-shares-facebook").attr("href", "https://www.facebook.com/sharer/sharer.php?u=" + addthis_share.url); $(".f1r-shares-linkedin").attr("href", addthis_config.services_custom[0].url); $(".f1r-shares-reddit").attr("href", addthis_config.services_custom[2].url); $(".f1r-shares-mendelay").attr("href", addthis_config.services_custom[1].url); function checkCount(){ setTimeout(function(){ $(".addthis_button_expanded").each(function(){ var count = $(this).text(); if (count !== "" && count != "0") $(this).removeClass("is-hidden"); else $(this).addClass("is-hidden"); }); }, 1000); } close How to cite this report {{reportCitation}} Cancel Copy Citation Details $(function(){R.ui.buttonDropdowns('.dropdown-for-downloads');}); $(function(){R.ui.toolbarDropdowns('.toolbar-dropdown-for-downloads');}); $.get("/articles/acj/155657/178643") new F1000.Clipboard(); new F1000.ThesaurusTermsDisplay("articles", "article", "178643"); $(document).ready(function() { $( "#frame1" ).on('load', function() { var mydiv = $(this).contents().find("div"); var h = mydiv.height(); console.log(h) }); var tooltipLivingFigure = jQuery(".interactive-living-figure-label .icon-more-info"), titleLivingFigure = tooltipLivingFigure.attr("title"); tooltipLivingFigure.simpletip({ fixed: true, position: ["-115", "30"], baseClass: 'small-tooltip', content:titleLivingFigure + " " }); tooltipLivingFigure.removeAttr("title"); $("body").on("click", ".cite-living-figure", function(e) { e.preventDefault(); var ref = $(this).attr("data-ref"); $(this).closest(".living-figure-list-container").find("#" + ref).fadeIn(200); }); $("body").on("click", ".close-cite-living-figure", function(e) { e.preventDefault(); $(this).closest(".popup-window-wrapper").fadeOut(200); }); $(document).on("mouseup", function(e) { var metricsContainer = $(".article-metrics-popover-wrapper"); if (!metricsContainer.is(e.target) && metricsContainer.has(e.target).length === 0) { $(".article-metrics-close-button").click(); } }); var articleId = $('#articleId').val(); if($("#main-article-count-box").attachArticleMetrics) { $("#main-article-count-box").attachArticleMetrics(articleId, { articleMetricsView: true }); } }); var figshareWidget = $(".new_figshare_widget"); if (figshareWidget.length > 0) { window.figshare.load("f1000", function(Widget) { // Select a tag/tags defined in your page. In this tag we will place the widget. _.map(figshareWidget, function(el){ var widget = new Widget({ articleId: $(el).attr("figshare_articleId") //height:300 // this is the height of the viewer part. [Default: 550] }); widget.initialize(); // initialize the widget widget.mount(el); // mount it in a tag that's on your page // this will save the widget on the global scope for later use from // your JS scripts. This line is optional. //window.widget = widget; }); }); } close Error Close Add Reset F1000.MICROSERVICES.AFFILIATION = ''; $(document).ready(function () { $('.js-affiliations-form').each((index, form) => { new AffiliationForm({ formId: form.id, institutionErrorSelector: '.comment-enter-institution', departmentErrorSelector: '.comment-enter-department', placeSelector: '.js-add-comment-place', stateSelector: '.js-add-comment-state', zipCodeSelector: '.js-add-comment-zipcode', countrySelector: '.js-add-comment-country', countryErrorSelector: '.comment-enter-country', }); }); }); $(document).ready(function () { var reportIds = { "341253": 0, "372229": 0, "341252": 0, "372228": 0, "341255": 0, "372231": 0, "341254": 0, "372230": 0, "372225": 0, "372224": 16, "341251": 0, "372227": 0, "341250": 0, "372226": 0, "341257": 0, "341256": 0, "341259": 0, "341258": 0, "381213": 0, "381212": 0, "381215": 0, "381214": 0, "381209": 0, "381211": 0, "381210": 0, "381217": 0, "381216": 0, "381218": 0, "393013": 0, "393012": 0, "393015": 0, "393014": 0, "355377": 0, "371506": 2, "335165": 0, "355389": 0, "393021": 0, "355388": 0, "393020": 10, "335167": 0, "355391": 0, "353087": 0, "335166": 0, "355390": 0, "393017": 0, "393016": 0, "393019": 0, "393018": 0, "335173": 0, "353093": 0, "365893": 0, "335172": 0, "355396": 0, "353092": 0, "335174": 0, "353094": 0, "335169": 0, "355393": 0, "353089": 0, "355392": 0, "335168": 38, "353088": 0, "335171": 0, "355395": 0, "353091": 0, "335170": 0, "355394": 0, "353090": 0, "353096": 0, "338317": 0, "338316": 0, "338319": 0, "383375": 0, "338318": 0, "383374": 0, "338315": 0, "383381": 0, "338324": 0, "383380": 0, "383383": 0, "383382": 0, "338321": 0, "383377": 0, "338320": 0, "383376": 0, "338323": 0, "383379": 0, "338322": 0, "383378": 0, "374445": 0, "374447": 0, "374446": 0, "378805": 0, "374453": 0, "378804": 0, "374452": 0, "378807": 0, "378806": 0, "374449": 0, "374448": 0, "374451": 0, "374450": 0, "378813": 0, "378812": 0, "378809": 0, "378808": 0, "378811": 0, "378810": 0, "402414": 0, "407534": 0, "402415": 0, "407535": 0, "402412": 0, "407532": 0, "402413": 0, "407533": 0, "407530": 0, "402411": 0, "407531": 0, "407529": 0, "402418": 0, "407538": 0, "402416": 0, "407536": 0, "402417": 0, "407537": 0, "372223": 0, "372222": 0, }; $(".referee-response-container,.js-referee-report").each(function(index, el) { var reportId = $(el).attr("data-reportid"), reportCount = reportIds[reportId] || 0; $(el).find(".comments-count-container,.js-referee-report-views").html(reportCount); }); var uuidInput = $("#article_uuid"), oldUUId = uuidInput.val(), newUUId = "06961b46-4133-476e-804f-b58deced2674"; uuidInput.val(newUUId); $("a[href*='article_uuid=']").each(function(index, el) { var newHref = $(el).attr("href").replace(oldUUId, newUUId); $(el).attr("href", newHref); }); }); An innovative open access publishing platform offering rapid publication and open peer review, whilst supporting data deposition and sharing. Browse Gateways Collections How it Works Contact For Developers Cookie Notice Privacy Notice RSS Submit Your Research Follow us © 2012-2026 F1000 Research Ltd. ISSN 2046-1402 | Legal | Partner of Research4Life • CrossRef • ORCID • FAIRSharing R.templateTests.simpleTemplate = R.template(' $text $text $text $text $text '); R.templateTests.runTests(); var F1000platform = new F1000.Platform({ name: "f1000research", displayName: "F1000Research", hostName: "f1000research.com", id: "1", editorialEmail: "
[email protected]", infoEmail: "
[email protected]", usePmcStats: true }); $(function(){R.ui.dropdowns('.dropdown-for-authors, .dropdown-for-about, .dropdown-for-myresearch');}); // $(function(){R.ui.dropdowns('.dropdown-for-referees');}); $(document).ready(function () { if ($(".cookie-warning").is(":visible")) { $(".sticky").css("margin-bottom", "35px"); $(".devices").addClass("devices-and-cookie-warning"); } $(".cookie-warning .close-button").click(function (e) { $(".devices").removeClass("devices-and-cookie-warning"); $(".sticky").css("margin-bottom", "0"); }); $("#tweeter-feed .tweet-message").each(function (i, message) { var self = $(message); self.html(linkify(self.html())); }); $(".partner").on("mouseenter mouseleave", function() { $(this).find(".gray-scale, .colour").toggleClass("is-hidden"); }); }); Sign In Remember me Forgotten your password? Sign In Cancel Email or password not correct. Please try again Please wait... $(function(){ // Note: All the setup needs to run against a name attribute and *not* the id due the clonish // nature of facebox... $("a[id=googleSignInButton]").click(function(event){ event.preventDefault(); $("input[id=oAuthSystem]").val("GOOGLE"); $("form[id=oAuthForm]").submit(); }); $("a[id=facebookSignInButton]").click(function(event){ event.preventDefault(); $("input[id=oAuthSystem]").val("FACEBOOK"); $("form[id=oAuthForm]").submit(); }); $("a[id=orcidSignInButton]").click(function(event){ event.preventDefault(); $("input[id=oAuthSystem]").val("ORCID"); $("form[id=oAuthForm]").submit(); }); }); If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password. The email address should be the one you originally registered with F1000. Email address not valid, please try again You registered with F1000 via Google, so we cannot reset your password. To sign in, please click here . If you still need help with your Google account password, please click here . You registered with F1000 via Facebook, so we cannot reset your password. To sign in, please click here . If you still need help with your Facebook account password, please click here . Code not correct, please try again Reset password Cancel Email us for further assistance. Server error, please try again. If your email address is registered with us, we will email you instructions to reset your password. If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance. Please wait... Register $(document).ready(function () { signIn.createSignInAsRow($("#sign-in-form-gfb-popup")); $(".target-field").each(function () { var uris = $(this).val().split("/"); if (uris.pop() === "login") { $(this).val(uris.toString().replace(",","/")); } }); });
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.