Regulation of PTEN translation by PI3K signaling maintains pathway homeostasis.

The PI3K pathway regulates cell metabolism, proliferation, and migration, and its dysregulation is common in cancer. We now show that both physiologic and oncogenic activation of PI3K signaling increase the expression of its negative regulator PTEN. This limits the duration of the signal and output of the pathway. Physiologic and pharmacologic inhibition of the pathway reduces PTEN and contributes to the rebound in pathway activity in tumors treated with PI3K inhibitors and limits their efficacy. Regulation of PTEN is due to mTOR/4E-BP1-dependent control of its translation and is lost when 4E-BP1 is deleted. Translational regulation of PTEN is therefore a major homeostatic regulator of physiologic PI3K signaling and plays a role in reducing the pathway activation by oncogenic PIK3CA mutants and the antitumor activity of PI3K pathway inhibitors. However, pathway output is hyperactivated in tumor cells with coexistent PI3K mutation and loss of PTEN function.

Computational model for drug research.

This special issue focuses on computational model for drug research regarding drug bioactivity prediction, drug-related interaction prediction, modelling for immunotherapy and modelling for treatment of a specific disease, as conveyed by the following six research and four review articles. Notably, these 10 papers described a wide variety of in-depth drug research from the computational perspective and may represent a snapshot of the wide research landscape.

A Loss of Epigenetic Control Can Promote Cell Death through Reversing the Balance of Pathways in a Signaling Network.

Epigenetic control of regulatory networks is only partially understood. Expression of Insulin-like growth factor-II (IGF2) is controlled by genomic imprinting, mediated by silencing of the maternal allele. Loss of imprinting of IGF2 (LOI) is linked to intestinal and colorectal cancers, causally in murine models and epidemiologically in humans. However, the molecular underpinnings of the LOI phenotype are not clear. Surprisingly, in LOI cells, we find a reversal of the relative activities of two canonical signaling pathways triggered by IGF2, causing further rebalancing between pro- and anti-apoptotic signaling. A predictive mathematical model shows that this network rebalancing quantitatively accounts for the effect of receptor tyrosine kinase inhibition in both WT and LOI cells. This mechanism also quantitatively explains both the stable LOI phenotype and the therapeutic window for selective killing of LOI cells, and thus prevention of epigenetically controlled cancers. These findings suggest a framework for understanding epigenetically modified cell signaling.

Value Estimation versus Effort Mobilization: A General Dissociation between Ventromedial and Dorsomedial Prefrontal Cortex.

Deciding on a course of action requires both an accurate estimation of option values and the right amount of effort invested in deliberation to reach sufficient confidence in the final choice. In a previous study, we have provided evidence, across a series of judgment and choice tasks, for a dissociation between the ventromedial prefrontal cortex (vmPFC), which would represent option values, and the dorsomedial prefrontal cortex (dmPFC), which would represent the duration of deliberation. Here, we first replicate this dissociation and extend it to the case of an instrumental learning task, in which 24 human volunteers (13 women) choose between options associated with probabilistic gains and losses. According to fMRI data recorded during decision-making, vmPFC activity reflects the sum of option values generated by a reinforcement learning model and dmPFC activity the deliberation time. To further generalize the role of the dmPFC in mobilizing effort, we then analyze fMRI data recorded in the same participants while they prepare to perform motor and cognitive tasks (squeezing a handgrip or making numerical comparisons) to maximize gains or minimize losses. In both cases, dmPFC activity is associated with the output of an effort regulation model, and not with response time. Taken together, these results strengthen a general theory of behavioral control that implicates the vmPFC in the estimation of option values and the dmPFC in the energization of relevant motor and cognitive processes.

Antibody inhibition of influenza A virus assembly and release.

Antibodies are frontline defenders against influenza virus infection, providing protection through multiple complementary mechanisms. Although a subset of monoclonal antibodies (mAbs) has been shown to restrict replication at the level of virus assembly and release, it remains unclear how potent and pervasive this mechanism of protection is, due in part to the challenge of separating this effect from other aspects of antibody function. To address this question, we developed imaging-based assays to determine how effectively a broad range of mAbs against the IAV surface proteins can specifically restrict viral egress. We find that classically neutralizing antibodies against hemagglutinin are broadly multifunctional, inhibiting virus assembly and release at concentrations 1-20-fold higher than the concentrations at which they inhibit viral entry. These antibodies are also capable of altering the morphological features of shed virions, reducing the proportion of filamentous particles. We find that antibodies against neuraminidase and M2 also restrict viral egress and that inhibition by anti-neuraminidase mAbs is only partly attributable to a loss in enzymatic activity. In all cases, antigen crosslinking-either on the surface of the infected cell, between the viral and cell membrane, or both-plays a critical role in inhibition, and we are able to distinguish between these modes experimentally and through a structure-based computational model. Together, these results provide a framework for dissecting antibody multifunctionality that could help guide the development of improved therapeutic antibodies or vaccines and that can be extended to other viral families and antibody isotypes.IMPORTANCEAntibodies against influenza A virus provide multifaceted protection against infection. Although sensitive and quantitative assays are widely used to measure inhibition of viral attachment and entry, the ability of diverse antibodies to inhibit viral egress is less clear. We address this challenge by developing an imaging-based approach to measure antibody inhibition of virus release across a panel of monoclonal antibodies targeting the influenza A virus surface proteins. Using this approach, we find that inhibition of viral egress is common and can have similar potency to the ability of an antibody to inhibit viral entry. Insights into this understudied aspect of antibody function may help guide the development of improved countermeasures.

Drug-drug interaction prediction: databases, web servers and computational models.

In clinical treatment, two or more drugs (i.e. drug combination) are simultaneously or successively used for therapy with the purpose of primarily enhancing the therapeutic efficacy or reducing drug side effects. However, inappropriate drug combination may not only fail to improve efficacy, but even lead to adverse reactions. Therefore, according to the basic principle of improving the efficacy and/or reducing adverse reactions, we should study drug-drug interactions (DDIs) comprehensively and thoroughly so as to reasonably use drug combination. In this review, we first introduced the basic conception and classification of DDIs. Further, some important publicly available databases and web servers about experimentally verified or predicted DDIs were briefly described. As an effective auxiliary tool, computational models for predicting DDIs can not only save the cost of biological experiments, but also provide relevant guidance for combination therapy to some extent. Therefore, we summarized three types of prediction models (including traditional machine learning-based models, deep learning-based models and score function-based models) proposed during recent years and discussed the advantages as well as limitations of them. Besides, we pointed out the problems that need to be solved in the future research of DDIs prediction and provided corresponding suggestions.

PlDBPred: a novel computational model for discovery of DNA binding proteins in plants.

DNA-binding proteins (DBPs) play crucial roles in numerous cellular processes including nucleotide recognition, transcriptional control and the regulation of gene expression. Majority of the existing computational techniques for identifying DBPs are mainly applicable to human and mouse datasets. Even though some models have been tested on Arabidopsis, they produce poor accuracy when applied to other plant species. Therefore, it is imperative to develop an effective computational model for predicting plant DBPs. In this study, we developed a comprehensive computational model for plant specific DBPs identification. Five shallow learning and six deep learning models were initially used for prediction, where shallow learning methods outperformed deep learning algorithms. In particular, support vector machine achieved highest repeated 5-fold cross-validation accuracy of 94.0% area under receiver operating characteristic curve (AUC-ROC) and 93.5% area under precision recall curve (AUC-PR). With an independent dataset, the developed approach secured 93.8% AUC-ROC and 94.6% AUC-PR. While compared with the state-of-art existing tools by using an independent dataset, the proposed model achieved much higher accuracy. Overall results suggest that the developed computational model is more efficient and reliable as compared to the existing models for the prediction of DBPs in plants. For the convenience of the majority of experimental scientists, the developed prediction server PlDBPred is publicly accessible at https://iasri-sg.icar.gov.in/pldbpred/.The source code is also provided at https://iasri-sg.icar.gov.in/pldbpred/source_code.php for prediction using a large-size dataset.

Enhanced compound-protein binding affinity prediction by representing protein multimodal information via a coevolutionary strategy.

Due to the lack of a method to efficiently represent the multimodal information of a protein, including its structure and sequence information, predicting compound-protein binding affinity (CPA) still suffers from low accuracy when applying machine-learning methods. To overcome this limitation, in a novel end-to-end architecture (named FeatNN), we develop a coevolutionary strategy to jointly represent the structure and sequence features of proteins and ultimately optimize the mathematical models for predicting CPA. Furthermore, from the perspective of data-driven approach, we proposed a rational method that can utilize both high- and low-quality databases to optimize the accuracy and generalization ability of FeatNN in CPA prediction tasks. Notably, we visually interpret the feature interaction process between sequence and structure in the rationally designed architecture. As a result, FeatNN considerably outperforms the state-of-the-art (SOTA) baseline in virtual drug evaluation tasks, indicating the feasibility of this approach for practical use. FeatNN provides an outstanding method for higher CPA prediction accuracy and better generalization ability by efficiently representing multimodal information of proteins via a coevolutionary strategy.

Hinge point emergence in mammalian spinal neurulation.

Neurulation is the process in early vertebrate embryonic development during which the neural plate folds to form the neural tube. Spinal neural tube folding in the posterior neuropore changes over time, first showing a median hinge point, then both the median hinge point and dorsolateral hinge points, followed by dorsolateral hinge points only. The biomechanical mechanism of hinge point formation in the mammalian neural tube is poorly understood. Here we employ a mechanical finite element model to study neural tube formation. The computational model mimics the mammalian neural tube using microscopy data from mouse and human embryos. While intrinsic curvature at the neural plate midline has been hypothesized to drive neural tube folding, intrinsic curvature was not sufficient for tube closure in our simulations. We achieved neural tube closure with an alternative model combining mesoderm expansion, nonneural ectoderm expansion, and neural plate adhesion to the notochord. Dorsolateral hinge points emerged in simulations with low mesoderm expansion and zippering. We propose that zippering provides the biomechanical force for dorsolateral hinge point formation in settings where the neural plate lateral sides extend above the mesoderm. Together, these results provide a perspective on the biomechanical and molecular mechanism of mammalian spinal neurulation.

Synchronous inhibitory pathways create both efficiency and diversity in the retina.

Sensory receptive fields combine features that originate in different neural pathways. Retinal ganglion cell receptive fields compute intensity changes across space and time using a peripheral region known as the surround, a property that improves information transmission about natural scenes. The visual features that construct this fundamental property have not been quantitatively assigned to specific interneurons. Here, we describe a generalizable approach using simultaneous intracellular and multielectrode recording to directly measure and manipulate the sensory feature conveyed by a neural pathway to a downstream neuron. By directly controlling the gain of individual interneurons in the circuit, we show that rather than transmitting different temporal features, inhibitory horizontal cells and linear amacrine cells synchronously create the linear surround at different spatial scales and that these two components fully account for the surround. By analyzing a large population of ganglion cells, we observe substantial diversity in the relative contribution of amacrine and horizontal cell visual features while still allowing individual cells to increase information transmission under the statistics of natural scenes. Established theories of efficient coding have shown that optimal information transmission under natural scenes allows a diverse set of receptive fields. Our results give a mechanism for this theory, showing how distinct neural pathways synthesize a sensory computation and how this architecture both generates computational diversity and achieves the objective of high information transmission.

Network modeling predicts personalized gene expression and drug responses in valve myofibroblasts cultured with patient sera.

Aortic valve stenosis (AVS) patients experience pathogenic valve leaflet stiffening due to excessive extracellular matrix (ECM) remodeling. Numerous microenvironmental cues influence pathogenic expression of ECM remodeling genes in tissue-resident valvular myofibroblasts, and the regulation of complex myofibroblast signaling networks depends on patient-specific extracellular factors. Here, we combined a manually curated myofibroblast signaling network with a data-driven transcription factor network to predict patient-specific myofibroblast gene expression signatures and drug responses. Using transcriptomic data from myofibroblasts cultured with AVS patient sera, we produced a large-scale, logic-gated differential equation model in which 11 biochemical and biomechanical signals were transduced via a network of 334 signaling and transcription reactions to accurately predict the expression of 27 fibrosis-related genes. Correlations were found between personalized model-predicted gene expression and AVS patient echocardiography data, suggesting links between fibrosis-related signaling and patient-specific AVS severity. Further, global network perturbation analyses revealed signaling molecules with the most influence over network-wide activity, including endothelin 1 (ET1), interleukin 6 (IL6), and transforming growth factor ß (TGFß), along with downstream mediators c-Jun N-terminal kinase (JNK), signal transducer and activator of transcription (STAT), and reactive oxygen species (ROS). Lastly, we performed virtual drug screening to identify patient-specific drug responses, which were experimentally validated via fibrotic gene expression measurements in valvular interstitial cells cultured with AVS patient sera and treated with or without bosentan-a clinically approved ET1 receptor inhibitor. In sum, our work advances the ability of computational approaches to provide a mechanistic basis for clinical decisions including patient stratification and personalized drug screening.

Closing the loop between microstructure and charge transport in conjugated polymers by combining microscopy and simulation.

A grand challenge in materials science is to identify the impact of molecular composition and structure across a range of length scales on macroscopic properties. We demonstrate a unified experimental-theoretical framework that coordinates experimental measurements of mesoscale structure with molecular-level physical modeling to bridge multiple scales of physical behavior. Here we apply this framework to understand charge transport in a semiconducting polymer. Spatially-resolved nanodiffraction in a transmission electron microscope is combined with a self-consistent framework of the polymer chain statistics to yield a detailed picture of the polymer microstructure ranging from the molecular to device relevant scale. Using these data as inputs for charge transport calculations, the combined multiscale approach highlights the underrepresented role of defects in existing transport models. Short-range transport is shown to be more chaotic than is often pictured, with the drift velocity accounting for a small portion of overall charge motion. Local transport is sensitive to the alignment and geometry of polymer chains. At longer length scales, large domains and gradual grain boundaries funnel charges preferentially to certain regions, creating inhomogeneous charge distributions. While alignment generally improves mobility, these funneling effects negatively impact mobility. The microstructure is modified in silico to explore possible design rules, showing chain stiffness and alignment to be beneficial while local homogeneity has no positive effect. This combined approach creates a flexible and extensible pipeline for analyzing multiscale functional properties and a general strategy for extending the accesible length scales of experimental and theoretical probes by harnessing their combined strengths.

Geometry of neural computation unifies working memory and planning.

Real-world tasks require coordination of working memory, decision-making, and planning, yet these cognitive functions have disproportionately been studied as independent modular processes in the brain. Here, we propose that contingency representations, defined as mappings for how future behaviors depend on upcoming events, can unify working memory and planning computations. We designed a task capable of disambiguating distinct types of representations. In task-optimized recurrent neural networks, we investigated possible circuit mechanisms for contingency representations and found that these representations can explain neurophysiological observations from the prefrontal cortex during working memory tasks. Our experiments revealed that human behavior is consistent with contingency representations and not with traditional sensory models of working memory. Finally, we generated falsifiable predictions for neural data to identify contingency representations in neural data and to dissociate different models of working memory. Our findings characterize a neural representational strategy that can unify working memory, planning, and context-dependent decision-making.

Neural signature of flexible coding in prefrontal cortex.

The ability of prefrontal cortex to quickly encode novel associations is crucial for adaptive behavior and central to working memory. Fast Hebbian changes in synaptic strength permit forming new associations, but neuronal signatures of this have been elusive. We devised a trialwise index of pattern similarity to look for rapid changes in population codes. Based on a computational model of working memory, we hypothesized that synaptic strength-and consequently, the tuning of neurons-could change if features of a subsequent stimulus need to be "reassociated," i.e., if bindings between features need to be broken to encode the new item. As a result, identical stimuli might elicit different neural responses. As predicted, neural response similarity dropped following rebinding, but only in prefrontal cortex. The history-dependent changes were expressed on top of traditional, fixed selectivity and were not explainable by carryover of previous firing into the current trial or by neural adaptation.

Emotional Memory Processing during REM Sleep with Implications for Post-Traumatic Stress Disorder.

REM sleep is important for the processing of emotional memories, including fear memories. Rhythmic interactions, especially in the theta band, between the medial prefrontal cortex (mPFC) and limbic structures are thought to play an important role, but the ways in which memory processing occurs at a mechanistic and circuits level are largely unknown. To investigate how rhythmic interactions lead to fear extinction during REM sleep, we used a biophysically based model that included the infralimbic cortex (IL), a part of the mPFC with a critical role in suppressing fear memories. Theta frequency (4-12 Hz) inputs to a given cell assembly in IL, representing an emotional memory, resulted in the strengthening of connections from the IL to the amygdala and the weakening of connections from the amygdala to the IL, resulting in the suppression of the activity of fear expression cells for the associated memory. Lower frequency (4 Hz) theta inputs effected these changes over a wider range of input strengths. In contrast, inputs at other frequencies were ineffective at causing these synaptic changes and did not suppress fear memories. Under post-traumatic stress disorder (PTSD) REM sleep conditions, rhythmic activity dissipated, and 4 Hz theta inputs to IL were ineffective, but higher-frequency (10 Hz) theta inputs to IL induced changes similar to those seen with 4 Hz inputs under normal REM sleep conditions, resulting in the suppression of fear expression cells. These results suggest why PTSD patients may repeatedly experience the same emotionally charged dreams and suggest potential neuromodulatory therapies for the amelioration of PTSD symptoms.SIGNIFICANCE STATEMENT Rhythmic interactions in the theta band between the mPFC and limbic structures are thought to play an important role in processing emotional memories, including fear memories, during REM sleep. The infralimbic cortex (IL) in the mPFC is thought to play a critical role in suppressing fear memories. We show that theta inputs to the IL, unlike other frequency inputs, are effective in producing synaptic changes that suppress the activity of fear expression cells associated with a given memory. Under PTSD REM sleep conditions, lower-frequency (4 Hz) theta inputs to the IL do not suppress the activity of fear expression cells associated with the given memory but, surprisingly, 10 Hz inputs do. These results suggest potential neuromodulatory therapies for PTSD.

Creation of Neuronal Ensembles and Cell-Specific Homeostatic Plasticity through Chronic Sparse Optogenetic Stimulation.

Cortical computations emerge from the dynamics of neurons embedded in complex cortical circuits. Within these circuits, neuronal ensembles, which represent subnetworks with shared functional connectivity, emerge in an experience-dependent manner. Here we induced ensembles in ex vivo cortical circuits from mice of either sex by differentially activating subpopulations through chronic optogenetic stimulation. We observed a decrease in voltage correlation, and importantly a synaptic decoupling between the stimulated and nonstimulated populations. We also observed a decrease in firing rate during Up-states in the stimulated population. These ensemble-specific changes were accompanied by decreases in intrinsic excitability in the stimulated population, and a decrease in connectivity between stimulated and nonstimulated pyramidal neurons. By incorporating the empirically observed changes in intrinsic excitability and connectivity into a spiking neural network model, we were able to demonstrate that changes in both intrinsic excitability and connectivity accounted for the decreased firing rate, but only changes in connectivity accounted for the observed decorrelation. Our findings help ascertain the mechanisms underlying the ability of chronic patterned stimulation to create ensembles within cortical circuits and, importantly, show that while Up-states are a global network-wide phenomenon, functionally distinct ensembles can preserve their identity during Up-states through differential firing rates and correlations.SIGNIFICANCE STATEMENT The connectivity and activity patterns of local cortical circuits are shaped by experience. This experience-dependent reorganization of cortical circuits is driven by complex interactions between different local learning rules, external input, and reciprocal feedback between many distinct brain areas. Here we used an ex vivo approach to demonstrate how simple forms of chronic external stimulation can shape local cortical circuits in terms of their correlated activity and functional connectivity. The absence of feedback between different brain areas and full control of external input allowed for a tractable system to study the underlying mechanisms and development of a computational model. Results show that differential stimulation of subpopulations of neurons significantly reshapes cortical circuits and forms subnetworks referred to as neuronal ensembles.

The acquired dyad inclination and decreased interpersonal brain communication in the pursuit of collective benefit.

People perform better collectively than individually, a phenomenon known as the collective benefit. To pursue the benefit, they may learn from previous behaviors, come to know whose initial opinion should be valued, and develop the inclination to take it as the collective one. Such learning may affect interpersonal brain communication. To test these hypotheses, this study recruited participant dyads to conduct a perceptual task on which they made individual decisions first and then the collective one. The enhanced interpersonal brain synchronization (IBS) between participants was explored when individual decisions were in disagreement vs. agreement. Computational modeling revealed that participant dyads developed the dyad inclination of taking the higher-able participants', not the lower-able ones' decisions as their collective ones. Brain analyses unveiled the enhanced IBS at frontopolar areas, premotor areas, supramarginal gyri, and right temporal-parietal junctions. The premotor IBS correlated negatively with dyad inclination and collective benefit in the absence of correction. The Granger causality analyses further supported the negative relation of dyad inclination with inter-brain communication. This study highlights that dyads learn to weigh individuals' decisions, resulting in dyad inclinations, and explores associated inter-brain communication, offering insights into the dynamics of collective decision-making.

Multiformity of extracellular microelectrode recordings from Aδ neurons in the dorsal root ganglia: a computational modeling study.

Microelectrodes serve as a fundamental tool in electrophysiology research throughout the nervous system, providing a means of exploring neural function with a high resolution of neural firing information. We constructed a hybrid computational model using the finite element method and multicompartment cable models to explore factors that contribute to extracellular voltage waveforms that are produced by sensory pseudounipolar neurons, specifically smaller A-type neurons, and that are recorded by microelectrodes in dorsal root ganglia. The finite element method model included a dorsal root ganglion, surrounding tissues, and a planar microelectrode array. We built a multicompartment neuron model with multiple trajectories of the glomerular initial segment found in many A-type sensory neurons. Our model replicated both the somatic intracellular voltage profile of Aδ low-threshold mechanoreceptor neurons and the unique extracellular voltage waveform shapes that are observed in experimental settings. Results from this model indicated that tortuous glomerular initial segment geometries can introduce distinct multiphasic properties into a neuron's recorded waveform. Our model also demonstrated how recording location relative to specific microanatomical components of these neurons, and recording distance from these components, can contribute to additional changes in the multiphasic characteristics and peak-to-peak voltage amplitude of the waveform. This knowledge may provide context for research employing microelectrode recordings of pseudounipolar neurons in sensory ganglia, including functional mapping and closed-loop neuromodulation. Furthermore, our simulations gave insight into the neurophysiology of pseudounipolar neurons by demonstrating how the glomerular initial segment aids in increasing the resistance of the stem axon and mitigating rebounding somatic action potentials.NEW & NOTEWORTHY We built a computational model of sensory neurons in the dorsal root ganglia to investigate factors that influence the extracellular waveforms recorded by microelectrodes. Our model demonstrates how the unique structure of these neurons can lead to diverse and often multiphasic waveform profiles depending on the location of the recording contact relative to microanatomical neural components. Our model also provides insight into the neurophysiological function of axon glomeruli that are often present in these neurons.

Estimation of the transmission delays in the basal ganglia of the macaque monkey and subsequent predictions about oscillatory activity under dopamine depletion.

The timescales of the dynamics of a system depend on the combination of the timescales of its components and of its transmission delays between components. Here, we combine experimental stimulation data from 10 studies in macaque monkeys that reveal the timing of excitatory and inhibitory events in the basal ganglia circuit, to estimate its set of transmission delays. In doing so, we reveal possible inconsistencies in the existing data, calling for replications, and we propose two possible sets of transmission delays. We then integrate these delays in a model of the primate basal ganglia that does not rely on direct and indirect pathways' segregation and show that extrastriatal dopaminergic depletion in the external part of the globus pallidus and in the subthalamic nucleus is sufficient to generate ß-band oscillations (in the high part, 20-35 Hz, of the band). More specifically, we show that D2 and D5 dopamine receptors in these nuclei play opposing roles in the emergence of these oscillations, thereby explaining how completely deactivating D5 receptors in the subthalamic nucleus can, paradoxically, cancel oscillations.

A large-scale neuronal network modelling study: Stimulus size modulates gamma oscillations in the primary visual cortex by long-range connections.

Stimulus size modulation of neuronal firing activity is a fundamental property of the primary visual cortex. Numerous biological experiments have shown that stimulus size modulation is affected by multiple factors at different spatiotemporal scales, but the exact pathways and mechanisms remain incompletely understood. In this paper, we establish a large-scale neuronal network model of primary visual cortex with layer 2/3 to study how gamma oscillation properties are modulated by stimulus size and especially how long-range connections affect the modulation as realistic neuronal properties and spatial distributions of synaptic connections are considered. It is shown that long-range horizontal synaptic connections are sufficient to produce dimensional modulation of firing rates and gamma oscillations. In particular, with increasing grating stimulus size, the firing rate increases and then decreases, the peak frequency of gamma oscillations decreases and the spectral power increases. These are consistent with biological experimental observations. Furthermore, we explain in detail how the number and spatial distribution of long-range connections affect the size modulation of gamma oscillations by using the analysis of neuronal firing activity and synaptic current fluctuations. Our results provide a mechanism explanation for size modulation of gamma oscillations in the primary visual cortex and reveal the important and unique role played by long-range connections, which contributes to a deeper understanding of the cognitive function of gamma oscillations in visual cortex.