(J Am Heart Assoc. 2016;5:e002764 doi: [10.1161/JAHA.115.002764](10.1161/JAHA.115.002764))

An accompanying Video S1 is available at <http://jaha.ahajournals.org/content/5/1/e002764/suppl/DC1>

Introduction {#jah31286-sec-0005}
============

Alterations in coronary wall shear stress (WSS) can interact with the endothelium to affect the distribution, progression, and pathophysiology of atherosclerosis,[1](#jah31286-bib-0001){ref-type="ref"}, [2](#jah31286-bib-0002){ref-type="ref"}, [3](#jah31286-bib-0003){ref-type="ref"}, [4](#jah31286-bib-0004){ref-type="ref"}, [5](#jah31286-bib-0005){ref-type="ref"}, [6](#jah31286-bib-0006){ref-type="ref"}, [7](#jah31286-bib-0007){ref-type="ref"}, [8](#jah31286-bib-0008){ref-type="ref"}, [9](#jah31286-bib-0009){ref-type="ref"} and several studies have demonstrated the association between WSS magnitudes and plaque development.[3](#jah31286-bib-0003){ref-type="ref"}, [4](#jah31286-bib-0004){ref-type="ref"}, [5](#jah31286-bib-0005){ref-type="ref"}, [6](#jah31286-bib-0006){ref-type="ref"}, [7](#jah31286-bib-0007){ref-type="ref"}, [8](#jah31286-bib-0008){ref-type="ref"}, [9](#jah31286-bib-0009){ref-type="ref"}, [10](#jah31286-bib-0010){ref-type="ref"} In particular, low WSS (\<1.0 Pa) has been related to plaque progression,[5](#jah31286-bib-0005){ref-type="ref"}, [6](#jah31286-bib-0006){ref-type="ref"}, [7](#jah31286-bib-0007){ref-type="ref"}, [8](#jah31286-bib-0008){ref-type="ref"}, [9](#jah31286-bib-0009){ref-type="ref"}, [10](#jah31286-bib-0010){ref-type="ref"} whereas high WSS (\>2.5 Pa) has been associated with increased necrotic core area on virtual histology intravascular ultrasound (VH‐IVUS), a marker of plaque vulnerability.[3](#jah31286-bib-0003){ref-type="ref"}, [4](#jah31286-bib-0004){ref-type="ref"}, [9](#jah31286-bib-0009){ref-type="ref"}, [10](#jah31286-bib-0010){ref-type="ref"} The impact of valsartan and simvastatin on the proinflammatory effect of low WSS has been studied in the swine model[11](#jah31286-bib-0011){ref-type="ref"}; however, the effect of an endothelial protective pharmaceutical agent, such as nebivolol, on the interplay between WSS and the progression of human coronary atherosclerosis has not been investigated.

Nebivolol is a third‐generation β1 receptor--selective antagonist with nitric oxide--mediated vasodilatory effect approved by the US Food and Drug Administration for treatment of hypertension and in Europe for left ventricular systolic dysfunction. Previous studies have suggested that nebivolol decreases coronary microvascular resistance through an agonist effect on the endothelial β3 receptor, which mediates nitric oxide-- and endothelium‐dependent vasorelaxation.[12](#jah31286-bib-0012){ref-type="ref"}, [13](#jah31286-bib-0013){ref-type="ref"}, [14](#jah31286-bib-0014){ref-type="ref"}, [15](#jah31286-bib-0015){ref-type="ref"}, [16](#jah31286-bib-0016){ref-type="ref"} Enhanced flow and nitric oxide availability may decrease vascular reactivity and lipid oxidation and prove beneficial in conditions such as atherosclerosis, in which oxidative stress leads to decreased endothelial nitric oxide synthase expression and nitric oxide bioavailability. Potential benefits from nebivolol therapy may be observed in measures of endothelial function, microvascular function, vascular remodeling, or epicardial plaque burden and composition.

Accordingly, we conducted a pilot double‐blinded randomized controlled trial on the effect of 1 year of nebivolol versus atenolol therapy on the changes in the number of thin‐capped fibroatheromas (TCFAs) per vessel and in coronary plaque volumetrics among different WSS categories in patients with nonobstructive coronary artery disease. We hypothesized that nebivolol, by decreasing oxidative stress and improving endothelial function, would favorably affect atherosclerosis progression compared with atenolol, a β‐blocker without nitric oxide--mediated activity.

Methods {#jah31286-sec-0006}
=======

Participants and Study Design {#jah31286-sec-0007}
-----------------------------

This single‐center double‐blinded randomized trial (NCT01230892) enrolled patients with recurrent angina symptoms who presented to the cardiac catheterization laboratory at Emory University Hospital from 2010 to 2012 and were found to have coronary nonobstructive atherosclerotic lesions (\<50% stenosis by angiography or \<70% stenosis with fractional flow reserve \>0.80). Patients provided written informed consent for participation in this study prior to the baseline cardiac catheterization and randomization.

Inclusion criteria included patients being on stable medical therapy, presentation of stable angina or non--ST‐segment elevation acute coronary syndrome, and coronary lesion in the proximal 60 mm of an epicardial vessel deemed significant enough by the operator to warrant further evaluation using physiology and VH‐IVUS. Patients were excluded if they presented with cardiogenic shock, ejection fraction \<30%, or significant hepatic, hematologic, or renal impairment; had a history of coronary artery bypass surgery or severe valvular heart disease; could not provide informed consent; or had any contraindication to β‐blocker therapy. Anatomically, cases were excluded if there were lesions requiring revascularization or significant visual coronary collaterals.

After undergoing baseline coronary angiography with VH‐IVUS and evaluation of coronary microvascular and endothelial function, patients were randomized (1:1) to nebivolol 10 mg daily (first week 5 mg daily) or atenolol 100 mg daily (first week 50 mg daily) for 1 year. Nebivolol 5 mg (Forest Laboratories) and atenolol 50 mg (AstraZeneca Pharmaceuticals) have equivalent efficacy for reducing blood pressure.[17](#jah31286-bib-0017){ref-type="ref"}

The Investigational Drug Service at our institution performed simple randomization using the online program "Research Randomizer" to generate the randomization plan in blocks of 4; dispensed the blinded study medication (same tablet size, shape, and color, with the same bottles and labels); and maintained separate records using WebIDS, a fully Health Insurance Portability and Accountability Act--compliant system that maintains protocol and drug information, patient profiles, randomization assignments, inventory records, and dispensing and patient return information. Enrolled patients were treated for 12±3 months with the study drug plus standard medical therapy (aspirin, statin, and as needed sublingual nitroglycerin) followed by repeat coronary angiography with VH‐IVUS and assessments of microvascular and endothelial function (Figure [1](#jah31286-fig-0001){ref-type="fig"}). Participants, care providers, and study investigators were blinded until after completion of the last follow‐up visit.

![Study flowchart.](JAH3-5-e002764-g001){#jah31286-fig-0001}

The Emory University institutional review board approved the study, and there were no significant changes to the protocol after the trial commenced. The primary end point was the change in the number of TCFAs per vessel between baseline and follow‐up. Secondary end points included inflammatory and oxidative stress biomarkers, endothelial function, microvascular function, static and serial arterial remodeling, and VH‐IVUS plaque area and composition at baseline and 1‐year follow‐up.

Lipids and Biomarkers {#jah31286-sec-0008}
---------------------

Basic metabolic panel, complete blood count, fasting lipid panel, and C‐reactive protein level were collected the day of angiography and measured by the hospital medical laboratory services. Samples for plasma glutathione and cystine measurements were collected immediately after obtaining vascular access in the cardiac catheterization laboratory, were frozen at −80°C, and were measured by the Emory University clinical biomarkers core laboratory. These samples have been shown to be stable under these conditions for 1 year.[18](#jah31286-bib-0018){ref-type="ref"}

Cardiac Catheterization With Assessment of Endothelial and Microvascular Function {#jah31286-sec-0009}
---------------------------------------------------------------------------------

All ß‐blockers were withheld for ≥48 hours and long‐acting nitrates were withheld for ≥24 hours prior to cardiac catheterization (both baseline and follow‐up procedures). Procedures were performed in the morning after patients had fasted for at least 8 hours. Systolic and diastolic blood pressures and heart rate were measured prior to the start of the procedure. Patients underwent angiography in a biplane cardiac catheterization system (Toshiba America Medical Systems) using a standard 6F technique.

Pressure and velocity measurements were obtained using a 0.014‐in pressure and Doppler flow velocity monitoring guidewire (ComboWire: Volcano Corporation). For safety considerations, only the left coronary system was interrogated. The ComboWire was advanced to the guide catheter tip, at which the aortic pressure and guidewire pressures were equalized. Microvascular function was evaluated from pressure and velocity responses to intravenous adenosine infusion (140 μg/kg per minute) for 3 minutes. At maximal hyperemia, fractional flow reserve was measured as the ratio of distal to aortic pressure, and hyperemic microvascular resistance was measured as the ratio of distal pressure to average peak flow velocity. Coronary flow velocity reserve was defined as the ratio of hyperemic to basal average peak flow velocity. Velocity measurements demonstrated good reproducibility, with a concordance correlation coefficient of 0.979 (95% CI 0.966--0.988).[9](#jah31286-bib-0009){ref-type="ref"}

Endothelial function was assessed as the percentage change in coronary diameter and blood flow in response to intracoronary acetylcholine (off‐label use) using Doppler flow velocity measurements and quantitative coronary angiography (QAngio XA 7.3; Medis Medical Imaging Systems). Patients were first evaluated with 10^−8^mol/L acetylcholine for safety prior to proceeding with 10^−6^mol/L acetylcholine. Figure [2](#jah31286-fig-0002){ref-type="fig"} demonstrates Doppler velocity envelopes prior to and during medication infusion in 1 patient at baseline and follow‐up.

![Doppler velocity waveforms obtained from a patient undergoing endothelial and microvascular function assessment at baseline and 1‐year follow‐up. APV indicates average peak velocity; CFR, coronary flow reserve; IC, intracoronary; IV, intravenous.](JAH3-5-e002764-g002){#jah31286-fig-0002}

VH‐IVUS Acquisition and Analysis {#jah31286-sec-0010}
--------------------------------

VH‐IVUS acquisition was performed after administration of 200 μg intracoronary nitroglycerin using phased‐array 20 MHz Eagle Eye catheters and a s5 Imaging System (Volcano Corporation). Automated motorized pullback (0.5 mm/s) was performed, and VH‐IVUS images were continuously acquired up to the guide catheter in the aorta (up to 60 mm of the proximal vessel). Angiography was used to record the catheter start position and its relationship to anatomic branching landmarks to aid with coregistration.

Offline analysis was performed at the Emory cardiovascular imaging and biomechanical core laboratory by experienced investigators who were blinded to the patients\' clinical data according to the criteria of the American College of Cardiology Clinical Consensus document on IVUS,[19](#jah31286-bib-0019){ref-type="ref"}, [20](#jah31286-bib-0020){ref-type="ref"} using VIAS version 3.0 (Volcano Corporation) and echoPlaque 4.0.27 (INDEC Medical Systems). Measurements of the external elastic membrane (EEM), plaque, and lumen cross‐sectional areas were performed for every recorded VH‐IVUS cross‐section (0.5‐mm thickness), defined as a segment in the current analysis. Plaque burden was calculated as plaque area divided by EEM area. Segments involving bifurcating branch points or precluding complete lumen or vessel wall planimetry were excluded from analysis. Intraobserver analysis demonstrated good reproducibility for plaque area (concordance correlation coefficient 0.968 \[95% CI 0.965--0.971\]).[9](#jah31286-bib-0009){ref-type="ref"}

Different nomenclatures have been used to describe arterial remodeling patterns. We defined static remodeling as the ratio of lesion to reference EEM area, also known as the remodeling index. We also used *excessive expansive, compensatory*, and *constrictive* to describe 3 patterns of serial remodeling.[19](#jah31286-bib-0019){ref-type="ref"} For each segment, positive ΔEEM area was defined as positive remodeling, and negative ΔEEM area was defined as constrictive remodeling. Segments with positive remodeling were further subdivided as compensatory if the ratio of ΔEEM area to Δplaque area was between 0.0 and 1.0 or as excessive expansive otherwise.

Finally, lesions were classified based on plaque composition as assessed by VH‐IVUS (fibrous, fibrofatty, necrotic core, and dense calcium) for each segment.[20](#jah31286-bib-0020){ref-type="ref"}, [21](#jah31286-bib-0021){ref-type="ref"} Fibroatheromas were defined as ≥40% plaque burden and ≥10% necrotic core; TCFAs were fibroatheromas with necrotic core abutting the lumen in at least 3 consecutive frames.

Computational Fluid Dynamics and Baseline WSS Analysis {#jah31286-sec-0011}
------------------------------------------------------

WSS analysis methodology using the ANGUS method has been described previously.[1](#jah31286-bib-0001){ref-type="ref"}, [9](#jah31286-bib-0009){ref-type="ref"} Briefly, the 3‐dimensional path of the VH‐IVUS catheter was determined using corresponding biplane angiographic projections acquired at the start of pullback. Each frame was rotated and aligned perpendicular to the catheter core to reconstruct the main artery of interest. Arterial branches were added based on information from angiography and VH‐IVUS. Patient‐specific pulsatile inlet boundary conditions were determined from Doppler wire‐derived waveforms, and the outlets were assumed to be pressure‐free. The reconstructed surface was meshed and imported into Fluent (ANSYS). After computing the pulsatile flow field in the region of interest, WSS was determined as a function of time in the cardiac cycle and spatial location around the lumen (Video S1, Figure [3](#jah31286-fig-0003){ref-type="fig"}) and then averaged over time and circumference at each cross‐section for quantitative analysis. Although WSS is a continuous and dynamic variable, based on previous cell culture and experimental and human data, WSS magnitudes were categorized as low (\<1.0 Pa), intermediate (≥1.0 and \<2.5 Pa), or high (≥2.5 Pa).[3](#jah31286-bib-0003){ref-type="ref"}, [9](#jah31286-bib-0009){ref-type="ref"}, [22](#jah31286-bib-0022){ref-type="ref"}, [23](#jah31286-bib-0023){ref-type="ref"}

![A patient‐specific WSS profile of the left anterior descending coronary artery demonstrates areas of variable WSS magnitudes. WSS indicate wall shear stress.](JAH3-5-e002764-g003){#jah31286-fig-0003}

Statistical Analysis {#jah31286-sec-0012}
--------------------

Sample size was not determined strictly from effect sizes because the purpose of this pilot study was to generate data to power future clinical trials. It was estimated using previous comparative studies between nebivolol and atenolol on endothelial and microvascular function assuming a 0.05 significance level and 0.80 power value.[13](#jah31286-bib-0013){ref-type="ref"}, [15](#jah31286-bib-0015){ref-type="ref"}, [16](#jah31286-bib-0016){ref-type="ref"}

Continuous variables are described as mean±SD or median and interquartile range, as appropriate, and categorical variables are described as counts and proportions. For patient‐level analyses, the chi‐square statistic was used to compare categorical variables (eg, hypertension), and the 2‐sample *t* test was used for continuous variables (eg, age). The differences between baseline and follow‐up patient‐level measures (eg, systolic blood pressure) were computed for each group, and the differences between the 2 groups were compared with the 2‐sample *t* test.

Correlated error is introduced by the clustering of arterial segments within patients. To correct for this, *P* values were adjusted with linear regression for continuous outcomes (eg, Δlumen area) and ordered logistic regression for categorical variables (eg, serial arterial remodeling). In both situations, the Huber White Sandwich Estimator was used to adjust the *P* values for correlated error.[24](#jah31286-bib-0024){ref-type="ref"}, [25](#jah31286-bib-0025){ref-type="ref"} To further examine the association of WSS, drug, and change measures, mixed models with WSS and drug as predictor variables were run. Although there are many observations at the VH‐IVUS segmental level, the effective sample size is 24 participants, and 2 covariates were determined to be a reasonable limit. Analyses were performed using SAS 9.3 (SAS Institute) or Stata 13.1 (StataCorp LP). *P*\<0.05 was established as the level of statistical significance.

Results {#jah31286-sec-0013}
=======

Study Participants, Entire Cohort {#jah31286-sec-0014}
---------------------------------

Fifty patients consented to study participation prior to baseline coronary angiography. Two patients who underwent initial cardiac catheterization developed severe vasospasm during 10^−6^mol/L acetylcholine infusion. They were successfully treated during the procedure and did well clinically but did not participate further in this trial or undergo randomization. Another 19 patients did not meet angiographic inclusion criteria and were considered screen failures. Consequently, 29 patients were enrolled in the study and underwent randomization. A total of 25 patients returned for follow‐up; 1 case was excluded from analysis due to poor data quality (Figure [1](#jah31286-fig-0001){ref-type="fig"}).

For the remaining 24 patients, physiological interrogation was performed in the left circumflex artery for 1 patient and in the left anterior descending artery for 23 patients. Table [1](#jah31286-tbl-0001){ref-type="table-wrap"} displays baseline patient characteristics. At follow‐up, both systolic and diastolic blood pressures decreased numerically (*P* value not significant), but heart rate changed significantly by −5 beats per minute (interquartile range −11 to +1 beats per minute; *P*=0.02). Low‐density lipoprotein cholesterol level also increased by 10 mg/dL (interquartile range −3 to +29 mg/dL; *P*=0.04). There were no significant changes in coronary physiology or endothelial function between baseline and follow‐up (Table [2](#jah31286-tbl-0002){ref-type="table-wrap"}).

###### 

Patient Demographics and Characteristics

                                                    Total (n=24)   Atenolol (n=12)   Nebivolol (n=12)   *P* Value
  ------------------------------------------------- -------------- ----------------- ------------------ -----------
  Age, y                                            52.5±10.3      49.8±10.7         55.2±9.5           0.21
  Male                                              9 (38)         5 (42)            4 (33)             0.67
  White race                                        18 (75)        8 (67)            10 (83)            0.35
  Body mass index                                   29.1±5.8       30.0±6.4          28.2±5.2           0.46
  Follow‐up time period, months                     12.4±1.0       12.5±0.6          12.4±1.3           0.81
  Cardiovascular risk factors                                                                           
  Hypertension                                      16 (67)        9 (75)            7 (58)             0.39
  Diabetes mellitus                                 4 (17)         2 (17)            2 (17)             1.00
  Dyslipidemia                                      18 (75)        8 (67)            10 (83)            0.35
  Prior myocardial infarction                       4 (17)         2 (17)            2 (17)             1.00
  Smoking history                                   12 (50)        7 (58)            5 (42)             0.41
  Family history of coronary disease                11 (46)        4 (33)            7 (58)             0.22
  Presentation                                                                                          
  Non--ST‐segment elevation myocardial infarction   1 (4)          1 (8)             0 (0)              
  Unstable angina                                   5 (21)         2 (17)            3 (25)             0.52
  Stable angina                                     17 (71)        8 (67)            9 (75)             
  Medication use                                                                                        
  Aspirin                                           20 (83)        10 (83)           10 (83)            1.00
  P2Y~12~ inhibitor                                 6 (25)         3 (25)            3 (25)             1.00
  Statin                                            23 (96)        12 (100)          11 (92)            0.31
  Calcium channel blocker                           11 (46)        7 (58)            4 (33)             0.22
  Long‐acting nitrate                               17 (71)        10 (83)           7 (58)             0.18
  ACE inhibitor or ARB                              6 (25)         3 (25)            3 (25)             1.00

Values are expressed as n (%) or mean±SD. ACE indicates angiotensin‐converting enzyme; ARB, angiotensin II receptor blocker.

###### 

Patient and Vessel Characteristics at Baseline and Change After 1 Year

                                                                 Baseline              Change After 1 Year                                                                                                                      
  -------------------------------------------------------------- --------------------- --------------------- ------ --------------------------------------------------- ------------------------------------------------------- -------
  Vital signs                                                                                                                                                                                                                   
  Systolic BP, mm Hg                                             137 \[122, 145\]      143 \[122, 159\]      0.51   −8 \[−25, 10\]                                      −7 \[−20, 1\]                                           0.95
  Diastolic BP, mm Hg                                            73 \[71, 85\]         78 \[69, 92\]         0.63   0 \[−8, 14\]                                        −8 \[−14, 6\]                                           0.08
  Heart rate, bpm                                                71 \[60, 85\]         68 \[62, 78\]         0.91   −5 \[−11, 6\]                                       −6 \[−11, −1\][a](#jah31286-note-0006){ref-type="fn"}   0.66
  Fasting lipid panel                                                                                                                                                                                                           
  Total cholesterol, mg/dL                                       158 \[136, 171\]      149 \[119, 160\]      0.13   15 \[−1, 31\]                                       14 \[−2, 55\]                                           0.80
  HDL cholesterol, mg/dL                                         48 \[38, 54\]         44 \[43, 54\]         0.75   3 \[−5, 16\]                                        −3 \[−8, 2\]                                            0.16
  LDL cholesterol, mg/dL                                         87 \[69, 104\]        81 \[72, 115\]        0.52   6 \[−3, 29\]                                        13 \[−4, 35\]                                           0.89
  Markers of inflammation and oxidative stress                                                                                                                                                                                  
  C‐reactive protein, mg/L                                       1.87 \[0.50, 3.82\]   2.05 \[1.30, 4.06\]   0.48   0.33 \[−1.72, 2.48\]                                0.07 \[−0.78, 0.42\]                                    0.72
  Cystine, μmol/L                                                84 \[72, 112\]        88 \[79, 108\]        0.41   10 \[−4, 20\]                                       17 \[−7, 43\]                                           0.37
  Glutathione, μmol/L                                            1.04 \[0.87, 1.60\]   0.98 \[0.79, 1.26\]   0.74   −0.13 \[−0.56, 0.46\]                               −0.01 \[−0.04, 0.17\]                                   0.52
  Cystine/glutathione ratio                                      84 \[63, 100\]        88 \[69, 139\]        0.53   31 \[−2, 44\]                                       15 \[−24, 48\]                                          0.78
  Physiology, endothelial function, and vessel characteristics                                                                                                                                                                  
  FFR                                                            0.98 \[0.94. 1.00\]   0.91 \[0.90, 0.98\]   0.08   −0.02 \[−0.08, 0.01\]                               −0.01 \[−0.03, 0.05\]                                   0.39
  CFR                                                            2.07 \[1.82, 2.38\]   1.87 \[1.62, 2.33\]   0.74   0.41 \[0.06, 0.48\]                                 0.71 \[−0.58, 1.26\]                                    0.97
  HMR                                                            1.75 \[1.47, 2.66\]   1.91 \[1.46, 2.21\]   0.89   0.07 \[−0.54, 1.07\]                                −0.06 \[−0.24, 0.43\]                                   1.00
  %Diameter change to ACh                                        −5 \[−22, 3\]         −2 \[−8, 15\]         0.47   2 \[−11, 30\]                                       3 \[−12, 25\]                                           0.89
  %CBF change to ACh                                             60 \[−11, 225\]       142 \[37, 181\]       0.81   −5 \[−172, 121\]                                    87 \[30, 385\]                                          0.47
  TCFAs per patient                                              2.2±1.2               1.6±1.0               0.27   −0.4±0.8                                            −0.3±1.0                                                0.94
  Segment characteristics                                                                                                                                                                                                       
  VH‐IVUS run length, mm                                         54.3±10.9             62.5±21.1             0.25                                                                                                               
  Minimum lumen area, mm^2^                                      6.11 \[3.85, 7.15\]   4.15 \[3.87, 7.77\]   0.93                                                                                                               
  Static remodeling                                              0.60±0.19             0.58±0.22             0.50                                                                                                               
  EEM area, mm^2^                                                16.5±7.3              17.5±7.3              0.76   0.20±1.99[a](#jah31286-note-0006){ref-type="fn"}    −0.78±2.46[a](#jah31286-note-0006){ref-type="fn"}       0.02
  Lumen area, mm^2^                                              10.9±4.5              12.7±5.9              0.41   0.33±2.08[a](#jah31286-note-0006){ref-type="fn"}    −1.22±2.68[a](#jah31286-note-0006){ref-type="fn"}       0.004
  Plaque area, mm^2^                                             5.7±4.8               4.9±3.3               0.66   −0.13±2.44[a](#jah31286-note-0006){ref-type="fn"}   0.44±1.95[a](#jah31286-note-0006){ref-type="fn"}        0.27
  FI area, mm^2^                                                 1.67±2.86             1.08±1.66             0.51   −0.16±1.58[a](#jah31286-note-0006){ref-type="fn"}   0.14±1.02[a](#jah31286-note-0006){ref-type="fn"}        0.26
  FF area, mm^2^                                                 0.23±0.46             0.16±0.35             0.57   −0.04±0.38[a](#jah31286-note-0006){ref-type="fn"}   0.02±0.38[a](#jah31286-note-0006){ref-type="fn"}        0.27
  NC area, mm^2^                                                 0.48±1.07             0.36±0.80             0.77   −0.07±0.76[a](#jah31286-note-0006){ref-type="fn"}   0.03±0.49[a](#jah31286-note-0006){ref-type="fn"}        0.35
  DC area, mm^2^                                                 0.09±0.17             0.15±0.38             0.28   0.01±0.17[a](#jah31286-note-0006){ref-type="fn"}    0.04±0.24[a](#jah31286-note-0006){ref-type="fn"}        0.34

Values are expressed as mean±SD or median \[IQR\]. ACh indicates acetylcholine; BP, blood pressure; bpm, beats per minute; CBF, coronary blood flow; CFR, coronary flow reserve; DC indicates dense calcium; EEM, external elastic membrane; FF, fibrofatty; FFR, fractional flow reserve; FI, fibrous; HDL, high‐density lipoprotein; HMR, hyperemic myocardial resistance; LDL, low‐density lipoprotein; NC, necrotic core; RI, remodeling index; TCFAs, thin‐capped fibroatheromas; VH‐IVUS, virtual histology intravascular ultrasound.

Significant change from baseline (*P*\<0.05).

Baseline mean VH‐IVUS run length was 58±17 mm, lumen area was 11.84±5.38 mm^2^, EEM area was 17.08±7.35 mm^2^, and plaque area was 5.24±4.07 mm^2^. At follow‐up, the change in EEM area was −0.33±2.31 mm^2^ (*P*\<0.0001), in lumen area was −0.51±2.54 mm^2^ (*P*\<0.0001), in plaque area was 0.18±2.21 mm^2^ (*P*\<0.0001), in fibrous area was −0.01±1.33 mm^2^ (*P*=0.67), in fibrofatty area was −0.01±0.38 mm^2^ (*P*=0.22), in necrotic core area was −0.02±0.64 mm^2^ (*P*=0.07), and in dense calcium area was 0.03±0.21 mm^2^ (*P*\<0.0001).

Relationship Between Baseline Coronary WSS, Serial Arterial Remodeling, and Plaque Progression {#jah31286-sec-0015}
----------------------------------------------------------------------------------------------

Coronary WSS was calculated in 1843 VH‐IVUS segments using computational fluid dynamics modeling. There were 428 (23%) low‐, 920 (50%) intermediate‐, and 495 (27%) high‐WSS segments. At follow‐up, low‐WSS segments demonstrated an increase in plaque area (*P*\<0.0001) and a decrease in lumen area (*P*\<0.0001) compared with intermediate‐ and high‐WSS segments (Table [3](#jah31286-tbl-0003){ref-type="table-wrap"}) and more constrictive remodeling (*P*=0.04) (Figure [4](#jah31286-fig-0004){ref-type="fig"}).

###### 

Comparison of the Change in Virtual Histology Intravascular Ultrasound Measures After 1 year by Baseline WSS Category

  Characteristic               Low WSS (n=428)   Intermediate WSS (n=920)   High WSS (n=495)   *P* Value
  ---------------------------- ----------------- -------------------------- ------------------ -----------
  ΔLumen area, mm^2^           −1.17±1.24        −0.64±1.61                 0.61±2.15          \<0.0001
  ΔEEM area, mm^2^             −0.70±1.32        −0.37±1.29                 0.05±1.82          0.40
  ΔPlaque area, mm^2^          0.47±1.37         0.27±1.82                  −0.56±2.77         \<0.0001
  ΔFibrous area, mm^2^         0.21±0.76         0.02±1.08                  −0.47±1.64         \<0.0001
  ΔFibrofatty area, mm^2^      0.04±0.30         −0.01±0.24                 −0.05±0.35         \<0.0001
  ΔNecrotic core area, mm^2^   −0.07±0.71        −0.01±0.66                 −0.17±0.88         \<0.0001
  ΔDense calcium area, mm^2^   −0.01±0.22        0.04±0.22                  0.04±0.26          \<0.0001

Values are expressed as mean±SD. EEM indicates external elastic membrane; WSS, wall shear stress.

![Serial arterial remodeling by baseline WSS categories at 1 year follow up. *P* value is for comparison of frequency of 3 remodeling groups across 3 WSS categories. WSS indicate wall shear stress.](JAH3-5-e002764-g004){#jah31286-fig-0004}

Comparisons Between Nebivolol and Atenolol {#jah31286-sec-0016}
------------------------------------------

Between the 2 cohorts, there were no significant differences in clinical characteristics or cardiac medications (Table [1](#jah31286-tbl-0001){ref-type="table-wrap"}) or other baseline and follow‐up variables such as cholesterol levels, inflammatory and oxidative stress biomarkers, coronary flow velocity reserve, hyperemic microvascular resistance, or coronary endothelial function (Table [2](#jah31286-tbl-0002){ref-type="table-wrap"}). There was a trend of more segments in the nebivolol cohort with low WSS (317 versus 111) and fewer segments with high WSS (175 versus 320) compared with the atenolol arm (Figure [5](#jah31286-fig-0005){ref-type="fig"}).

![Distribution of baseline WSS by β‐blocker. *P* value is for comparison of the frequency of 3 WSS categories across 2 β‐blockers. WSS indicates wall shear stress.](JAH3-5-e002764-g005){#jah31286-fig-0005}

On VH‐IVUS, nebivolol segments demonstrated significantly decreased EEM area (*P*=0.02) and lumen area (*P*=0.004) but similar changes in plaque components and number of TCFAs per vessel compared with atenolol (Table [2](#jah31286-tbl-0002){ref-type="table-wrap"}). Nebivolol also showed increased plaque area at follow‐up (*P*\<0.001), whereas atenolol demonstrated plaque regression (*P*=0.047); however, the difference between the 2 arms was not significant (*P*=0.27). With respect to arterial remodeling, static remodeling (remodeling index) was similar between the 2 groups, but serial remodeling was significantly different (*P*=0.002), with nebivolol segments demonstrating more constrictive remodeling (73% versus 46%) and less excessive expansive remodeling (16% versus 39%) than atenolol (Figure [6](#jah31286-fig-0006){ref-type="fig"}).

![Serial arterial remodeling by β‐blocker at 1‐year follow‐up. *P* value is for comparison of the frequency of 3 remodeling groups across 2 β‐blocker categories.](JAH3-5-e002764-g006){#jah31286-fig-0006}

These observations could have been related to differences in baseline WSS patterns between the arms or to pharmacological therapy with the 2 different β‐blockers, so VH‐IVUS results were stratified by categories of WSS (Table [4](#jah31286-tbl-0004){ref-type="table-wrap"}) and β‐blocker (Table [5](#jah31286-tbl-0005){ref-type="table-wrap"}). Among low‐WSS segments (Table [4](#jah31286-tbl-0004){ref-type="table-wrap"}), those treated with nebivolol demonstrated decreased ΔEEM area compared with atenolol (*P*=0.02), consistent with a higher percentage of constrictive serial remodeling (83% versus 63%).

###### 

Comparison of Change in Virtual Histology Intravascular Ultrasound Measures by β‐Blocker, Stratified by Baseline Wall Shear Stress Categories

                        Low WSS      Intermediate WSS   High WSS                                                               
  --------------------- ------------ ------------------ ---------- ------------ ------------ ------- ------------ ------------ ------
  ΔLumen area, mm^2^    −0.80±1.15   −1.29±1.25         0.22       0.21±1.79    −1.24±1.15   0.002   1.15±2.38    −0.39±1.08   0.02
  ΔEEM area, mm^2^      0.01±1.06    −0.95±1.31         0.02       −0.02±1.53   −0.62±1.02   0.15    0.29±1.96    −0.40±1.44   0.43
  ΔPlaque area, mm^2^   0.81±1.05    0.35±1.45          0.30       −0.23±2.38   0.62±1.15    0.28    −0.86±3.23   −0.01±1.45   0.30
  ΔFI area, mm^2^       0.36±0.84    0.16±0.72          0.23       −0.18±1.44   0.16±0.68    0.40    −0.66±1.94   −0.13±0.79   0.31
  ΔFF area, mm^2^       0.04±0.23    0.03±0.32          0.47       −0.03±0.25   0.00±0.24    0.77    −0.09±0.38   0.01±0.27    0.34
  ΔNC area, mm^2^       0.07±0.34    −0.12±0.79         0.45       −0.15±0.92   0.09±0.35    0.14    −0.25±1.06   −0.01±0.30   0.25
  ΔDC area, mm^2^       −0.01±0.12   −0.01±0.24         0.78       0.01±0.16    0.06±0.25    0.19    0.02±0.22    0.09±0.31    0.15

*P*‐values are for comparison between 2 β‐blocker categories. DC indicates dense calcium; EEM, external elastic membrane; FF, fibrofatty; FI, fibrous; NC, necrotic core; WSS, wall shear stress.

Table [5](#jah31286-tbl-0005){ref-type="table-wrap"} compares VH‐IVUS changes by WSS categories within each β‐blocker cohort. Compared with high‐WSS segments, low‐WSS segments demonstrated greater decrease in lumen area (*P*\<0.0001) and increase in plaque area (*P*\<0.01). After adjustment for β‐blocker, WSS category remained an independent predictor of changes in lumen area (*P*\<0.0001) and plaque area (*P*\<0.0001).

###### 

Comparison of Change in VH‐IVUS Measures by Baseline WSS Categories, Stratified by Type of β‐Blocker

                        Atenolol     Nebivolol                                                                   
  --------------------- ------------ ------------ ------------ ---------- ------------ ------------ ------------ ----------
  ΔLumen area, mm^2^    −0.80±1.15   0.21±1.79    1.15±2.38    \<0.0001   −1.29±1.25   −1.24±1.15   −0.39±1.08   \<0.0001
  ΔVessel area, mm^2^   0.01±1.06    −0.02±1.53   0.29±1.96    0.002      −0.95±1.31   −0.62±1.02   −0.40±1.44   \<0.0001
  ΔPlaque area, mm^2^   0.81±1.05    −0.23±2.38   −0.86±3.23   \<0.0001   0.35±1.45    0.62±1.15    −0.01±1.45   0.006
  ΔFI area, mm^2^       0.36±0.84    −0.18±1.44   −0.66±1.94   \<0.0001   0.16±0.72    0.16±0.68    −0.13±0.79   \<0.0001
  ΔFF area, mm^2^       0.04±0.23    −0.03±0.25   −0.09±0.38   \<0.0001   0.03±0.32    0.00±0.24    0.01±0.27    0.02
  ΔNC area, mm^2^       0.07±0.34    −0.15±0.92   −0.25±1.06   \<0.0001   −0.12±0.79   0.09±0.35    −0.01±0.30   0.01
  ΔDC area, mm^2^       −0.01±0.12   0.01±0.16    0.02±0.22    0.57       −0.01±0.24   0.06±0.25    0.09±0.31    \<0.0001

*P*‐values are for comparison among 3 WSS categories. DC indicates dense calcium; FF, fibro‐fatty; FI, fibrous; NC, necrotic core; WSS, wall shear stress.

Discussion {#jah31286-sec-0017}
==========

This investigation is the first randomized controlled clinical trial studying the impact of a β~1~‐blocker with endothelial β~3~‐receptor agonist activity on the proatherosclerotic effect of coronary WSS. The major observations were as follows: First, there was greater plaque progression (*P*\<0.0001) and constrictive remodeling in segments with low WSS (*P*=0.04) and, conversely, greater plaque regression and excessive expansive remodeling in segments with high WSS. Second, there were greater reductions in lumen area (*P*=0.004) and vessel area (*P*=0.02), resulting in more constrictive remodeling, with nebivolol (*P*=0.002), likely driven by a trend toward more low‐WSS segments in the nebivolol cohort (*P*=0.06). Third, after adjusting for the type of β‐blocker therapy, low‐WSS segments remained significantly associated with lumen loss and plaque progression (*P*\<0.0001). Fourth, there were no significant differences in oxidative stress or inflammatory biomarker levels, endothelial function, microvascular function, or number of TCFAs per vessel between nebivolol and atenolol.

Coronary atherosclerosis pathophysiology operates through complex multitude of pathways, including endothelial dysfunction as a conduit for systemic risk factors to affect the vasculature.[18](#jah31286-bib-0018){ref-type="ref"}, [26](#jah31286-bib-0026){ref-type="ref"}, [27](#jah31286-bib-0027){ref-type="ref"} The considerable number of investigations into the effects of WSS magnitudes on plaque development emphasizes this complexity. Low WSS can promote atherogenesis by changing endothelial cell morphology,[2](#jah31286-bib-0002){ref-type="ref"} enhancing production of reactive oxygen species and inflammatory molecules,[22](#jah31286-bib-0022){ref-type="ref"} and stimulating vascular smooth muscle cell migration.[8](#jah31286-bib-0008){ref-type="ref"}, [28](#jah31286-bib-0028){ref-type="ref"} In our study, we observed that low‐WSS regions were associated with more constrictive remodeling (ΔEEM −0.70 mm^2^, *P*=0.04) and an absolute increase in plaque area of 0.47 mm^2^ (*P*\<0.0001). Conversely, high WSS can also contribute to plaque destabilization.[4](#jah31286-bib-0004){ref-type="ref"}, [10](#jah31286-bib-0010){ref-type="ref"}, [29](#jah31286-bib-0029){ref-type="ref"}, [30](#jah31286-bib-0030){ref-type="ref"}, [31](#jah31286-bib-0031){ref-type="ref"} Histological studies have implicated high WSS in smooth muscle cell apoptosis and proteoglycan matrix degradation,[29](#jah31286-bib-0029){ref-type="ref"} whereas animal studies have demonstrated increased mechanical strain, vascular inflammation, macrophage activity, and expansive remodeling in high‐WSS regions.[30](#jah31286-bib-0030){ref-type="ref"}, [32](#jah31286-bib-0032){ref-type="ref"} We previously observed in human coronary arteries that, compared with intermediate‐WSS segments, high‐WSS segments trended toward more excessive expansive arterial remodeling,[9](#jah31286-bib-0009){ref-type="ref"}, [23](#jah31286-bib-0023){ref-type="ref"} an exaggerated remodeling pattern in which both vessel and lumen dimensions increase proportionally more than plaque area. Although this type of remodeling in high‐WSS segments may be partially explained by adaptive arterial enlargement in an attempt to restore WSS to a more physiological level,[31](#jah31286-bib-0031){ref-type="ref"} it can also drive plaque progression through continued lipid accumulation and inflammation. In this study, we demonstrated again that high‐WSS regions were compensated by excessive expansive remodeling (Figure [5](#jah31286-fig-0005){ref-type="fig"}) with fibrous and fibrofatty tissue regression (Table [3](#jah31286-tbl-0003){ref-type="table-wrap"}).

The observation that nebivolol patients had more constrictive remodeling and lumen loss compared with atenolol patients in this study was interesting. Although this was a double‐blinded randomized controlled trial, the nebivolol cohort had a numerically greater number of low coronary WSS segments per patient compared with the atenolol cohort (26.4 versus 9.3, *P*=0.06), despite having similar baseline plaque areas (4.9 mm^2^ versus 5.7 mm^2^, *P*=0.66), static vascular remodeling (remodeling index 0.58 versus 0.60, *P*=0.50), and coronary flow reserve (1.87 versus 2.07, *P*=0.74). After adjusting for type of β‐blocker, low WSS remained significantly associated with plaque progression and lumen loss (both *P*\<0.0001), suggesting that WSS and not type of β‐blocker therapy was driving these findings.

Interestingly, in our study, nebivolol and atenolol had similar changes in lipid profiles, inflammation and oxidative stress biomarkers, and coronary endothelial and microvascular function. In another study of hypertensive patients, we also observed no differences between nebivolol and metoprolol regarding pulse wave velocity, a noninvasive measurement of arterial stiffness and endothelial function, or oxidative stress biomarker levels.[33](#jah31286-bib-0033){ref-type="ref"} In addition, there were no significant differences between nebivolol and atenolol with respect to necrotic core area or the number of TCFAs per vessel. Taken together, the current investigation demonstrates no significant favorable or adverse effects of nebivolol compared with atenolol on comprehensive atherosclerosis phenotyping.

Despite the multidimensional nature of atherosclerosis pathophysiology, previous mechanistic investigations into the impact of pharmacotherapies on atherosclerosis development have generally relied on basic measures of plaque assessment such as the change in carotid intima media thickness[34](#jah31286-bib-0034){ref-type="ref"}, [35](#jah31286-bib-0035){ref-type="ref"} or coronary atheroma volume.[36](#jah31286-bib-0036){ref-type="ref"}, [37](#jah31286-bib-0037){ref-type="ref"}, [38](#jah31286-bib-0038){ref-type="ref"} These measures may underestimate the effect of antiatherosclerotic therapies on the multiple known atherosclerotic pathways and processes. To our knowledge, no prior investigation has comprehensively examined the incremental effect of a potential antiatherosclerotic therapy on epicardial and microvascular endothelial‐dependent and ‐independent function, atherosclerosis burden and composition, static and serial arterial remodeling, and oxidative stress and inflammatory pathways in the context of coronary WSS. This investigation offers a blueprint for how novel cardiovascular therapies can be evaluated, given the numerous possible mechanisms for therapeutic intervention.

Limitations {#jah31286-sec-0018}
-----------

The study was powered for VH‐IVUS end points but may not have been adequately powered to detect changes in endothelial and microvascular function or number of TCFAs per artery. Although the nebivolol dose used in this trial was approved for clinical use, the effective concentration may be too low or the duration of follow‐up insufficient to affect the prespecified end points. Nevertheless, this study is the most comprehensive atherosclerosis phenotyping double‐blinded randomized controlled trial performed thus far and can establish the methodology for future coronary atherosclerosis trials for novel pharmaceutical agents such as protein convertase subtilisin/kexin type 9 inhibitors.

Plaque composition data are derived from VH‐IVUS, which has inherent limitations compared with histology but can be performed in vivo in coronary arteries and has been validated by several studies.[6](#jah31286-bib-0006){ref-type="ref"}, [7](#jah31286-bib-0007){ref-type="ref"}, [8](#jah31286-bib-0008){ref-type="ref"}, [9](#jah31286-bib-0009){ref-type="ref"}, [20](#jah31286-bib-0020){ref-type="ref"}, [21](#jah31286-bib-0021){ref-type="ref"} In addition, although there were many arterial segments, the effective sample size was 24 participants; this limits the number of covariates used in multivariate analysis. Finally, this was a randomized clinical trial between nebivolol and atenolol; however, the constrictive remodeling seen in the nebivolol arm is likely related to the greater number of baseline segments with low WSS. Although randomization can be readily performed with prespecified balance in demographics, it is currently impractical, given the time‐ and resource‐intensive nature of computational fluid dynamic modeling, to randomize patients based on their distribution of coronary WSS segments.

Conclusions {#jah31286-sec-0019}
===========

Compared with atenolol, nebivolol demonstrated greater luminal reduction and constrictive remodeling, but this was likely driven by the higher number of low‐WSS segments in the nebivolol cohort. There were no significant differences in coronary endothelial function or microvascular function between treatment groups. This study suggests that any impact of an endothelial β~3~‐receptor agonist on the effects of WSS magnitudes or on the natural history of epicardial and microvascular atherosclerosis, if present, is minimal. Importantly, it lays down the methodological framework for future mechanistic evaluation of pharmacotherapy on human coronary atherosclerosis.

Disclosures {#jah31286-sec-0022}
===========

McDaniel consults for Medicure, Inc. Samady receives research funding from Forest Laboratories an Allergan Affiliate, Volcano Corporation, St. Jude Medical, Gilead Sciences, Medtronic, Inc and Abbott Vascular.

Supporting information
======================

###### 

**Video S1.** A patient‐specific wall shear stress (WSS) profile of the left anterior descending coronary artery demonstrates areas of variable WSS magnitudes throughout the cardiac cycle.

###### 

Click here for additional data file.

We thank the Emory University Clinical Biomarkers Lab for analyzing the oxidative stress biomarker samples and the Emory University Hospital cardiac catheterization laboratory support staff for their help in the procedures.

Forest Laboratories financed this study. The NIH Ruth L. Kirschstein National Research Service Awards training grant (5T32HL007745) supported Hung and Eshtehardi. The American Heart Association Postdoctoral Fellowships supported Molony and Timmins.
