Brain potentials evoked by intraepidermal electrical stimuli reflect the central sensitization of nociceptive pathways

Secondary mechanical punctate hyperalgesia is a cardinal sign of central sensitization (CS), an important mechanism of chronic pain. Our study demonstrates that hyperalgesia from intraepidermal electrical stimulation coexists with mechanical punctate hyperalgesia and elicits electroencephalographic (EEG) potentials that predict the occurrence of punctate hyperalgesia in a human experimental model of CS. These findings inform clinical development of EEG-based biomarkers of CS.

Central sensitization (CS), the increased sensitivity of the central nervous system to somatosensory inputs, accounts for secondary hyperalgesia, a typical sign of several painful clinical conditions. Brain potentials elicited by mechanical punctate stimulation using flat-tip probes can provide neural correlates of CS, but their signal-to-noise ratio is limited by poor synchronization of the afferent nociceptive input. Additionally, mechanical punctate stimulation does not activate nociceptors exclusively. In contrast, low-intensity intraepidermal electrical stimulation (IES) allows selective activation of type II Aδ-mechano-heat nociceptors (II-AMHs) and elicits reproducible brain potentials. However, it is unclear whether hyperalgesia from IES occurs and coexists with secondary mechanical punctate hyperalgesia, and whether the magnitude of the electroencephalographic (EEG) responses evoked by IES within the hyperalgesic area is increased. To address these questions, we explored the modulation of the psychophysical and EEG responses to IES by intraepidermal injection of capsaicin in healthy human subjects. We obtained three main results. First, the intensity of the sensation elicited by IES was significantly increased in participants who developed robust mechanical punctate hyperalgesia after capsaicin injection (i.e., responders), indicating that hyperalgesia from IES coexists with punctate mechanical hyperalgesia. Second, the N2 peak magnitude of the EEG responses elicited by IES was significantly increased after the intraepidermal injection of capsaicin in responders only. Third, a receiver-operator characteristics analysis showed that the N2 peak amplitude is clearly predictive of the presence of CS. These findings suggest that the EEG responses elicited by IES reflect secondary hyperalgesia and therefore represent an objective correlate of CS.

A geometric model of defensive peripersonal space

Potentially harmful stimuli occurring within the defensive peripersonal space (DPPS), a protective area surrounding the body, elicit stronger defensive reactions. The spatial features of the DPPS are poorly defined and limited to descriptive estimates of its extent along a single dimension. Here we postulated a family of geometric models of the DPPS, to address two important questions with respect to its spatial features: What is its fine-grained topography? How does the nervous system represent the body area to be defended? As a measure of the DPPS, we used the strength of the defensive blink reflex elicited by electrical stimulation of the hand (hand-blink reflex, HBR), which is reliably modulated by the position of the stimulated hand in egocentric coordinates. We tested the goodness of fit of the postulated models to HBR data from six experiments in which we systematically explored the HBR modulation by hand position in both head-centered and body-centered coordinates. The best-fitting model indicated that 1) the nervous system's representation of the body area defended by the HBR can be approximated by a half-ellipsoid centered on the face and 2) the DPPS extending from this area has the shape of a bubble elongated along the vertical axis. Finally, the empirical observation that the HBR is modulated by hand position in head-centered coordinates indicates that the DPPS is anchored to the face. The modeling approach described in this article can be generalized to describe the spatial modulation of any defensive response.

Comparison of gSSR and EST-SSR markers for analyzing genetic variability among tomato cultivars (Solanum lycopersicum L.)

In order to study genetic variability and develop better strategies for the utilization of 48 tomato cultivars from America, China, the Netherlands, and Portugal, genomic simple sequence repeat (gSSR) and EST-derived SSR (EST-SSR) markers were applied. In all, 15 of 82 gSSR and 18 of 115 EST-SSR markers showed polymorphic loci. There were 995 and 2072 clear fragments amplified by polymorphic gSSR and EST-SSR markers, respectively. The total and average number of alleles detected by EST-SSRs (75, 4.2) was more than gSSRs (54, 3.6) as a result of some multi-locus EST-SSRs. A lower polymorphism information content value was found in gSSRs (0.529) compared to EST-SSRs (0.620). Similarity coefficient matrixes of the 48 tomato cultivars were established based on the gSSRs and EST-SSRs, and UPGMA dendrograms were constructed from the gSSRs and EST-SSRs similarity coefficient matrixes. A high similarity was observed between the gSSRs and EST-SSRs dendrograms. Genetic variability of four tomato populations from different countries showed that the observed number of alleles and Nei's genetic diversity were highest in the American population, and the effective number of alleles was highest in the Dutch population. The estimated genetic structure showed some tomato cultivars from different countries shared a common genetic background, which might be related to gene flow. It was inferred that both gSSR and EST-SSR markers were effective to assess genetic variability of tomato cultivars, and the combination of both markers could be more effective for genetic diversity analysis in tomato.

Laser-Evoked Vertex Potentials Predict Defensive Motor Actions

The vertex potential is the largest response that can be recorded in the electroencephalogram of an awake, healthy human. It is elicited by sudden and intense stimuli, and is composed by a negative–positive deflection. The stimulus properties that determine the vertex potential amplitude have been well characterized. Nonetheless, its functional significance remains elusive. The dominant interpretation is that it reflects neural activities related to the detection of salient stimuli. However, given that threatening stimuli elicit both vertex potentials and defensive movements, we hypothesized that the vertex potential is related to the execution of defensive actions. Here, we directly compared the salience and motoric interpretations by investigating the relationship between the amplitude of laser-evoked potentials (LEPs) and the response time of movements with different defensive values. First, we show that a larger LEP negative wave (N2 wave) predicts faster motor response times. Second, this prediction is significantly stronger when the motor response is defensive in nature. Third, the relation between the N2 wave and motor response time depends not only on the kinematic form of the movement, but also on whether that kinematic form serves as a functional defense of the body. Therefore, the N2 wave of the LEP encodes key defensive reactions to threats.

Electronic damage in S atoms in a native protein crystal induced by an intense X-ray free-electron laser pulse

Current hard X-ray free-electron laser (XFEL) sources can deliver doses to biological macromolecules well exceeding 1 GGy, in timescales of a few tens of femtoseconds. During the pulse, photoionization can reach the point of saturation in which certain atomic species in the sample lose most of their electrons. This electronic radiation damage causes the atomic scattering factors to change, affecting, in particular, the heavy atoms, due to their higher photoabsorption cross sections. Here, it is shown that experimental serial femtosecond crystallography data collected with an extremely bright XFEL source exhibit a reduction of the effective scattering power of the sulfur atoms in a native protein. Quantitative methods are developed to retrieve information on the effective ionization of the damaged atomic species from experimental data, and the implications of utilizing new phasing methods which can take advantage of this localized radiation damage are discussed.

Separated response function ratios in exclusive, forward π(±) electroproduction

The study of exclusive π(±) electroproduction on the nucleon, including separation of the various structure functions, is of interest for a number of reasons. The ratio RL=σL(π-)/σL(π+) is sensitive to isoscalar contamination to the dominant isovector pion exchange amplitude, which is the basis for the determination of the charged pion form factor from electroproduction data. A change in the value of RT=σT(π-)/σT(π+) from unity at small -t, to 1/4 at large -t, would suggest a transition from coupling to a (virtual) pion to coupling to individual quarks. Furthermore, the mentioned ratios may show an earlier approach to perturbative QCD than the individual cross sections. We have performed the first complete separation of the four unpolarized electromagnetic structure functions above the dominant resonances in forward, exclusive π(±) electroproduction on the deuteron at central Q(2) values of 0.6, 1.0, 1.6  GeV(2) at W=1.95  GeV, and Q(2)=2.45  GeV(2) at W=2.22  GeV. Here, we present the L and T cross sections, with emphasis on RL and RT, and compare them with theoretical calculations. Results for the separated ratio RL indicate dominance of the pion-pole diagram at low -t, while results for RT are consistent with a transition between pion knockout and quark knockout mechanisms.

Akt and mTOR mediate programmed necrosis in neurons

Necroptosis is a newly described form of regulated necrosis that contributes to neuronal death in experimental models of stroke and brain trauma. Although much work has been done elucidating initiating mechanisms, signaling events governing necroptosis remain largely unexplored. Akt is known to inhibit apoptotic neuronal cell death. Mechanistic target of rapamycin (mTOR) is a downstream effector of Akt that controls protein synthesis. We previously reported that dual inhibition of Akt and mTOR reduced acute cell death and improved long term cognitive deficits after controlled-cortical impact in mice. These findings raised the possibility that Akt/mTOR might regulate necroptosis. To test this hypothesis, we induced necroptosis in the hippocampal neuronal cell line HT22 using concomitant treatment with tumor necrosis factor α (TNFα) and the pan-caspase inhibitor N-benzyloxycarbonyl-Val-Ala-Asp-fluoromethylketone. TNFα/zVAD treatment induced cell death within 4 h. Cell death was preceded by RIPK1–RIPK3–pAkt assembly, and phosphorylation of Thr-308 and Thr473 of AKT and its direct substrate glycogen synthase kinase-3β, as well as mTOR and its direct substrate S6 ribosomal protein (S6), suggesting activation of Akt/mTOR pathways. Pretreatment with Akt inhibitor viii and rapamycin inhibited Akt and S6 phosphorylation events, mitochondrial reactive oxygen species production, and necroptosis by over 50% without affecting RIPK1–RIPK3 complex assembly. These data were confirmed using small inhibitory ribonucleic acid-mediated knockdown of AKT1/2 and mTOR. All of the aforementioned biochemical events were inhibited by necrostatin-1, including Akt and mTOR phosphorylation, generation of oxidative stress, and RIPK1–RIPK3–pAkt complex assembly. The data suggest a novel, heretofore unexpected role for Akt and mTOR downstream of RIPK1 activation in neuronal cell death.

Protein engineering and preclinical development of a GM-CSF receptor antibody for the treatment of rheumatoid arthritis

BACKGROUND AND PURPOSE

For antibody therapies against receptor targets, in vivo outcomes can be difficult to predict because of target-mediated clearance or antigen 'sink' effects. The purpose of this work was to engineer an antibody to the GM-CSF receptor α (GM-CSFRα) with pharmacological properties optimized for chronic, s.c. treatment of rheumatoid arthritis (RA) patients.

EXPERIMENTAL APPROACH

We used an in silico model of receptor occupancy to guide the target affinity and a combinatorial phage display approach for affinity maturation. Mechanism of action and internalization assays were performed on the optimized antibody in vitro before refining the modelling predictions of the eventual dosing in man. Finally, in vivo pharmacology studies in cynomolgus monkeys were carried out to inform the predictions and support future clinical development.

KEY RESULTS

Antibody potency was improved 8600-fold, and the target affinity was reached. The refined model predicted pharmacodynamic effects at doses as low as 1 mg kg(-1) and a study in cynomolgus monkeys confirmed in vivo efficacy at 1 mg kg(-1) dosing.

CONCLUSIONS AND IMPLICATIONS

This rational approach to antibody drug discovery enabled the isolation of a potent molecule compatible with chronic, s.c. self-administration by RA patients. We believe this general approach enables the development of optimal biopharmaceuticals.

Population Structures of 'Candidatus Liberibacter asiaticus' in Southern China

Huanglongbing (HLB) is a highly detrimental citrus disease associated with 'Candidatus Liberibacter asiaticus', a nonculturable alpha-proteobacterium. Characterization of the bacterial populations is important for development of disease management strategies. In this study, the 'Ca. L. asiaticus' populations in eight provinces in southern China where HLB is endemic were analyzed based on tandem repeat number (TRN) variations in a previously characterized genomic locus CLIBASIA_01645. Of the 224 HLB samples collected, 175 (78.3%) samples yielded single polymerase chain reaction (PCR) amplicons (the single amplicon group, SAG) and 49 (21.7%) samples produced multiple PCR amplicons (the multiple amplicon group, MAG). Variations in SAG are summarized by Nei's diversity index (H) and ratio of TRN ≤ 10/TRN > 10 genotypes (R10). Variations in the MAG are described by the percentage of occurrence (PMAG). At an orchard-level comparison, the 'Ca. L. asiaticus' population from a Guangdong orchard (n = 24) showed H = 0.50, R10 = 23, and PMAG = 0, significantly different from that of the non-Guangdong orchards in Yunnan (n = 23), H = 0.83, R10 = 2.3, and PMAG = 11.5, and in Hainan (n = 35), H = 0.88, R10 = 1.5, and PMAG = 16.7. In a region-level consideration, the Guangdong 'Ca. L. asiaticus' population (n = 78) was H = 0.77, R10 = 25, and PMAG = 1.3, whereas the non-Guangdong population (n = 84) was H = 0.91, R10 = 1.6, and PMAG = 26.9. Overall, significant differences were observed between the 'Ca. L. asiaticus' population from Guangdong Province and those from the other provinces. A strong aggregation of TRN = 6, 7, and 8 genotypes is characteristic to the 'Ca. L. asiaticus' population in Guangdong. Referenced to genome annotation, we propose that rearrangement of tandem repeats at locus CLIBASIA_01645 could be associated with bacterial environmental adaptation.

Primary sensory cortices contain distinguishable spatial patterns of activity for each sense

Whether primary sensory cortices are essentially multisensory or whether they respond to only one sense is an emerging debate in neuroscience. Here we use a multivariate pattern analysis of functional magnetic resonance imaging data in humans to demonstrate that simple and isolated stimuli of one sense elicit distinguishable spatial patterns of neuronal responses, not only in their corresponding primary sensory cortex, but in other primary sensory cortices. These results indicate that primary sensory cortices, traditionally regarded as unisensory, contain unique signatures of other senses and, thereby, prompt a reconsideration of how sensory information is coded in the human brain.

Human primary sensory cortices are traditionally regarded as being able to process only one sensory modality. Liang and colleagues use brain imaging to show that, as well as being processed in typically corresponding cortical areas, different sensory modalities are also processed in atypical cortical areas.

The primary somatosensory cortex contributes to the latest part of the cortical response elicited by nociceptive somatosensory stimuli in humans

Nociceptive laser pulses elicit temporally-distinct cortical responses (the N1, N2 and P2 waves of laser-evoked potentials, LEPs) mainly reflecting the activity of the primary somatosensory cortex (S1) contralateral to the stimulated side, and of the bilateral operculoinsular and cingulate cortices. Here, by performing two different EEG experiments and applying a range of analysis approaches (microstate analysis, scalp topography, single-trial estimation), we describe a distinct component in the last part of the human LEP response (P4 wave). We obtained three main results. First, the LEP is reliably decomposed in four main and distinct functional microstates, corresponding to the N1, N2, P2, and P4 waves, regardless of stimulus territory. Second, the scalp and source configurations of the P4 wave follow a clear somatotopical organization, indicating that this response is likely to be partly generated in contralateral S1. Third, single-trial latencies and amplitudes of the P4 are tightly coupled with those of the N1, and are similarly sensitive to experimental manipulations (e.g., to crossing the hands over the body midline), suggesting that the P4 and N1 may have common neural sources. These results indicate that the P4 wave is a clear and distinct LEP component, which should be considered in LEP studies to achieve a comprehensive understanding of the brain response to nociceptive stimulation.

Single-particle structure determination by correlations of snapshot X-ray diffraction patterns

Diffractive imaging with free-electron lasers allows structure determination from ensembles of weakly scattering identical nanoparticles. The ultra-short, ultra-bright X-ray pulses provide snapshots of the randomly oriented particles frozen in time, and terminate before the onset of structural damage. As signal strength diminishes for small particles, the synthesis of a three-dimensional diffraction volume requires simultaneous involvement of all data. Here we report the first application of a three-dimensional spatial frequency correlation analysis to carry out this synthesis from noisy single-particle femtosecond X-ray diffraction patterns of nearly identical samples in random and unknown orientations, collected at the Linac Coherent Light Source. Our demonstration uses unsupported test particles created via aerosol self-assembly, and composed of two polystyrene spheres of equal diameter. The correlation analysis avoids the need for orientation determination entirely. This method may be applied to the structural determination of biological macromolecules in solution.

Chimerism but not neonatal antigen exposure induces transplant tolerance

Acquired transplant tolerance could be readily induced during foetal or neonatal period through donor cell infusion, but it is not the case in adults. This phenomenon has been attributed to the variation of immune system development in neonatal and adult periods. To investigate the role of immature immune system and chimerism in neonatal transplant tolerance, irradiated spleen cells or cell fraction from F1 (BALB/c × C57BL/6) or GFP-F1 mice were injected intravenously into neonatal C57BL/6 mice to induce tolerance. Irradiated cells or cell fraction could not induce chimerism and transplant tolerance in neonatal mice, even increasing the dose of donor cells to 5 × 10(7). Living donor cells induced tolerance in neonatal mice, and the quantity of living cells was correlated with the degree of chimerism and tolerance. At the amount of 3 × 10(7) F1 spleen cells, skin grafts were survived permanently in more than 80% of treated mice. However, the amount of 0.7 × 10(7) F1 spleen cells could only slightly prolong allografts' survival. The more donor cells were infused, the higher level of chimerism was achieved, and the higher frequency of alloreactive T cells was deleted. Chimerism is prerequisite for the induction and maintenance of tolerance. Chimerism in long-term tolerant mice was significantly higher than that in chronic graft rejected mice, with 6.48 ± 4.02% versus 1.41 ± 0.77%. It implies that transplant tolerance depends on the establishment of chimerism, but not on antigen exposure to immature immune system in foetus or neonates.

Fractal morphology, imaging and mass spectrometry of single aerosol particles in flight

The morphology of micrometre-size particulate matter is of critical importance in fields ranging from toxicology to climate science, yet these properties are surprisingly difficult to measure in the particles' native environment. Electron microscopy requires collection of particles on a substrate; visible light scattering provides insufficient resolution; and X-ray synchrotron studies have been limited to ensembles of particles. Here we demonstrate an in situ method for imaging individual sub-micrometre particles to nanometre resolution in their native environment, using intense, coherent X-ray pulses from the Linac Coherent Light Source free-electron laser. We introduced individual aerosol particles into the pulsed X-ray beam, which is sufficiently intense that diffraction from individual particles can be measured for morphological analysis. At the same time, ion fragments ejected from the beam were analysed using mass spectrometry, to determine the composition of single aerosol particles. Our results show the extent of internal dilation symmetry of individual soot particles subject to non-equilibrium aggregation, and the surprisingly large variability in their fractal dimensions. More broadly, our methods can be extended to resolve both static and dynamic morphology of general ensembles of disordered particles. Such general morphology has implications in topics such as solvent accessibilities in proteins, vibrational energy transfer by the hydrodynamic interaction of amino acids, and large-scale production of nanoscale structures by flame synthesis.

Femtosecond dark-field imaging with an X-ray free electron laser

The emergence of femtosecond diffractive imaging with X-ray lasers has enabled pioneering structural studies of isolated particles, such as viruses, at nanometer length scales. However, the issue of missing low frequency data significantly limits the potential of X-ray lasers to reveal sub-nanometer details of micrometer-sized samples. We have developed a new technique of dark-field coherent diffractive imaging to simultaneously overcome the missing data issue and enable us to harness the unique contrast mechanisms available in dark-field microscopy. Images of airborne particulate matter (soot) up to two microns in length were obtained using single-shot diffraction patterns obtained at the Linac Coherent Light Source, four times the size of objects previously imaged in similar experiments. This technique opens the door to femtosecond diffractive imaging of a wide range of micrometer-sized materials that exhibit irreproducible complexity down to the nanoscale, including airborne particulate matter, small cells, bacteria and gold-labeled biological samples.

Defensive peripersonal space: the blink reflex evoked by hand stimulation is increased when the hand is near the face

Electrical stimulation of the median nerve at the wrist may elicit a blink reflex [hand blink reflex (HBR)] mediated by a neural circuit at brain stem level. As, in a Sherringtonian sense, the blink reflex is a defensive response, in a series of experiments we tested, in healthy volunteers, whether and how the HBR is modulated by the proximity of the stimulated hand to the face. Electromyographic activity was recorded from the orbicularis oculi, bilaterally. We observed that the HBR is enhanced when the stimulated hand is inside the peripersonal space of the face, compared with when it is outside, irrespective of whether the proximity of the hand to the face is manipulated by changing the position of the arm ( experiment 1) or by rotating the head while keeping the arm position constant ( experiment 3). Experiment 2 showed that such HBR enhancement has similar magnitude when the participants have their eyes closed. Experiments 4 and 5 showed, respectively, that the blink reflex elicited by the electrical stimulation of the supraorbital nerve, as well as the N20 wave of the somatosensory evoked potentials elicited by the median nerve stimulation, are entirely unaffected by hand position. Taken together, our results provide compelling evidence that the brain stem circuits mediating the HBR in humans undergo tonic and selective top-down modulation from higher order cortical areas responsible for encoding the location of somatosensory stimuli in external space coordinates. These findings support the existence of a “defensive” peripersonal space, representing a safety margin advantageous for survival.

Bypassing primary sensory cortices--a direct thalamocortical pathway for transmitting salient sensory information

Detection and appropriate reaction to sudden and intense events happening in the sensory environment is crucial for survival. By combining Bayesian model selection with dynamic causal modeling of functional magnetic resonance imaging data, a novel analysis approach that allows inferring the causality between neural activities in different brain areas, we demonstrate that salient sensory information reaches the multimodal cortical areas responsible for its detection directly from the thalamus, without being first processed in primary and secondary sensory-specific areas. This direct thalamocortical transmission of multimodal salient information is parallel to the processing of finer stimulus attributes, which are transmitted in a modality-specific fashion from the thalamus to the relevant primary sensory areas. Such direct thalamocortical connections bypassing primary sensory cortices provide a fast and efficient way for transmitting information from subcortical structures to multimodal cortical areas, to allow the early detection of salient events and, thereby, trigger immediate and appropriate behavior.

MEDI-573 as a novel approach to IGF-1R and IR-A signaling inhibition by blocking IGF ligands: Phase I PK/PD, safety data, and disease linkage studies in breast cancer

271

Background: MEDI-573 is a dual-targeting human antibody that neutralizes IGF-I/-II ligands and inhibits IGF-1R and insulin receptor-A (IR-A) signaling pathways, which play a role in breast and other epithelial cancers. By sparing IR-B and its hybrid receptors, MEDI-573 is expected to achieve antitumor activity without perturbing glucose homeostasis.

Methods: In a 3 + 3 dose-escalation trial ( NCT00816361 ), patients with advanced solid tumors, KPS ≥ 60, received MEDI-573 as a weekly 1-hr IV infusion at 0.5, 1.5, 5, 10, or 15 mg/kg every 21 days. Objectives included determination of MTD, safety profile, pharmacokinetics (PK), pharmacodynamics, and tumor response. In a separate biomarker study, mRNA expression of IR-A and IR-B was compared in HR+/HER- breast cancers.

Results: Patients included 13M/12F (17 evaluable for efficacy at data cutoff); median age 59 yrs (range 37-78). No DLTs, drug-related SAEs, or significant toxicity patterns were reported. Forty-five AEs were considered treatment-related (% of patients): decreased appetite (28%), fatigue (24%), nausea (20%), diarrhea (16%), and anemia (12%). Hyperglycemia occurred in 2 patients, 1 with evidence of prestudy insulin resistance. No other significant changes in plasma glucose levels occurred. MEDI-573 had dose-proportionality, with AUC of 415 ± 165, 597 ± 298, and 1940 ± 904 d*μg/mL at 5, 10, and 15 mg/kg, respectively. Suppression of free IGF-1 and -2 was fully achieved at doses ≥ 5 mg/kg. Stable disease spanning ≥ 12 wks (range 12-36+) was seen in 6/17 patients. No antidrug antibodies were reported. In a study to identify breast cancer subsets that are dependent on IGF signaling, the IR-A:IR-B ratio was significantly higher in HR+/HER2- breast cancer than matched adjacent normal tissues and in proliferative luminal B than in luminal A breast cancer.

Conclusions: MEDI-573 showed acceptable safety and favorable PK profiles without significant changes in glucose levels. The IR-A:IR-B ratio supports a biomarker approach for MEDI-573 patient stratification. A biomarker-rich phase Ib/II study of MEDI-573 with an aromatase inhibitor in patients with advanced HR+/HER2- breast cancer opens in 2011.

Taking into account latency, amplitude, and morphology: improved estimation of single-trial ERPs by wavelet filtering and multiple linear regression

Across-trial averaging is a widely used approach to enhance the signal-to-noise ratio (SNR) of event-related potentials (ERPs). However, across-trial variability of ERP latency and amplitude may contain physiologically relevant information that is lost by across-trial averaging. Hence, we aimed to develop a novel method that uses 1) wavelet filtering (WF) to enhance the SNR of ERPs and 2) a multiple linear regression with a dispersion term (MLR(d)) that takes into account shape distortions to estimate the single-trial latency and amplitude of ERP peaks. Using simulated ERP data sets containing different levels of noise, we provide evidence that, compared with other approaches, the proposed WF+MLR(d) method yields the most accurate estimate of single-trial ERP features. When applied to a real laser-evoked potential data set, the WF+MLR(d) approach provides reliable estimation of single-trial latency, amplitude, and morphology of ERPs and thereby allows performing meaningful correlations at single-trial level. We obtained three main findings. First, WF significantly enhances the SNR of single-trial ERPs. Second, MLR(d) effectively captures and measures the variability in the morphology of single-trial ERPs, thus providing an accurate and unbiased estimate of their peak latency and amplitude. Third, intensity of pain perception significantly correlates with the single-trial estimates of N2 and P2 amplitude. These results indicate that WF+MLR(d) can be used to explore the dynamics between different ERP features, behavioral variables, and other neuroimaging measures of brain activity, thus providing new insights into the functional significance of the different brain processes underlying the brain responses to sensory stimuli.