Entropy

MotivationWith the rapid development of AI methods that predict protein structures from sequence, understanding the structure-function relation increasingly depends on quantitative structural descriptors that are both biologically meaningful and scalable to large datasets. Here, we introduce mathematical topology metrics that quantify the entanglement complexity of a tertiary protein structure while respecting uncrossability constraints. ResultsBy employing only three such metrics across all protein structures in the Protein Data Bank, we represent the proteome structural space in a continuous three-dimensional space. Distances within this space capture structural similarity and correlate with functional similarity. We find that the mathematical entanglement based landscape of protein structural space diversifies with the evolutionary expansion of protein function across species. Moreover, this continuous representation reproduces CATH classifications with high accuracy for major structural classes. These results indicate that these metrics efficiently encode structural features linked to protein function and provide a more informative description than conventional metrics. AvailabilityData used in this study are available in the Protein Data Bank. Details of the machine learning model used can be found in https://github.com/roshitac/CATH_Classification-. ContactBanu.Ozkan@asu.edu, Eleni.Panagiotou@asu.edu Supplementary informationSupplementary data are available at Journal Name online.

How morphological complexity evolves is still enigmatic. While there is evidence in algae and plants as well as animals that diversification of the repertoire of transcription factors (TF) is causative for evolution of organismal complexity, there are many examples from lineages that follow their own way of complexity evolution, for example by expansion of particular families. For land plants, correlation of the size of the TF complement with number of cell types (as a proxy for morphological complexity) has been shown, and several families were identified as candidates to drive complexity evolution. Here, we expand a previously available dataset of cell type numbers from 12 to 82 proteomes and introduce a four class body plan scheme. We find that the total TF complement correlates with the number of cell types of Archaeplastida (primary plastid bearing plants and algae). We used TabPFN (Tabular Prior-data Fitted Network) for binary (uni- vs. multicellularity) as well as for four class Bauplan classification. TabPFN is able to predict the morphological complexity with high accuracy. This approach allows to determine organismal complexity based on the gene space of an organism. Based on our results, we can confirm that plant morphological evolution is driven by gain and expansion of TF families.

Network centrality is the workhorse of gene prioritisation, yet what a ranking omits is rarely audited. Scoring each selection against an annotation-count-matched maximum-entropy reference--asking whether a selected gene set covers the genomes functional space or collapses it-reveals that the criterion in standard use has a measurable blind spot in exactly the class it is meant to surface. Degree, the most widely used criterion, returns the cross-module bridges that are also locally dominant--connector hubs--and omits the non-hub connectors: where 26% of the genome occupies these coordinating roles, a degree-ranked list holds 18% and an EDVS-ranked list 55%, and degrees top-1% collapses functional coverage below the reference on all five networks tested. We repurpose EDVS (Entropy of Degree-Vector Sums), an information-theoretic diversity measure, as an annotation-free, partition-free centrality that recovers this omitted class. The coverage it preserves is carried by cross-module participation P, which cannot be computed without a community partition; EDVS matches P-level coverage on all five networks using none, and retains 0.84 of its selection under edge perturbation that leaves partition-based selections at 0.21-0.46. The deficit is general: the collapse holds in the same direction on the two networks built without functional annotation (0.5-1.1 bit; co-expression, physical interaction) as on the three supervised by it (1.6-3.3 bit; RiceNet, AraNet, STRING), so supervision amplifies it rather than creates it. The remedy is bounded: EDVS ceases to preserve coverage on the sparse physical-interaction network. And the class EDVS isolates is organizational, not an importance signal: pre-registered probes--essentiality, transcription-factor identity, tissue-specificity, date/party-hub character, phenotype co-localisation--return null or reversed throughout. The conclusive ones are equivalent to their degree-matched nulls within {+/-}5 percentage points (demonstrated, not merely undetected), and the classical coupling of centrality to importance itself holds only network-dependently. Author SummaryGenes rarely act alone: many diseases and agricultural traits are shaped by genes that coordinate several biological processes rather than specialising in one. The standard way to find such genes in a network of gene interactions is to count each genes connections--its "centrality"--and rank genes by that count. We show this standard approach has a blind spot: it favours genes that dominate one process over genes that quietly bridge several processes without dominating any, and this blind spot appears across rice, thale cress, and yeast gene networks. We repurpose a diversity measure from an unrelated field (originally used to compare citation patterns) as a new way to rank genes that finds these bridging genes from network structure alone, without needing gene-function annotations--which are themselves incomplete and biased toward well-studied genes--or a prior, unstable step of splitting the network into modules. We are careful to show where the new approach also falls short: on sparse, noisy networks it stops working, and the genes it recovers are not shown to be more biologically important than other genes, only differently positioned. What that position is for is a question this work leaves open.

Structural network representations of metastatic dissemination typically focus on static topology without resolving transport dynamics, relaxation timescales, or steady-state behaviour. Here, we formulate a discrete Markovian transport model on a directed higher-order network with transition rates derived from qualitative clinical affinity classes. By constructing a non-Hermitian row-stochastic transfer operator, we characterise the relaxation dynamics through its spectral decomposition. The system exhibits a fast-mixing regime characterised by a spectral gap of {gamma} {approx} 0.67, corresponding to a characteristic relaxation timescale of {tau} {approx} 1.49 discrete steps, with the influence of the primary tumour origin progressively attenuated during dissemination. Convergence towards a non-equilibrium steady state (NESS) is accompanied by a reduction in Shannon entropy, concentrating probability mass within specific topological sinks. This spectral relaxation delineates two distinct dynamical regimes: early transient dissemination (n < {tau}), dominated by local organ-specific transition probabilities (organotropism), and the asymptotic regime (n > {tau}), determined increasingly by the global transport architecture of the network. Comparison with independent clinical and autopsy observations across 21 primary tumours and 23 target organs indicates that the predicted stationary distribution is consistent with the observed hierarchy of metastatic organ involvement.

The folding thermodynamics of proteins are dominated by two opposing forces, the loss in backbone entropy and the packing of hydrophobic groups. The same forces are major contributors to the extension thermodynamics of elastic proteins with the distinction that both processes act in concert, favoring the higher chain and solvent entropy of a relaxed conformation. The relative entropic contributions specify the recoil mechanism; human elastin recoil is primarily driven by hydrophobic forces, whereas fly resilin has a rubber-like mechanism driven by backbone entropy. Despite the importance of elastic proteins to tissue biomechanics, few have been identified, let alone characterized to the same extent as elastin and resilin. We develop a thermodynamic framework that maps proteins by sequence-derived estimates of extension-induced backbone and solvent entropy changes. Putative elastic proteins are proposed and classified by recoil mechanism based on estimated thermodynamic features. Proteins that map to elastic regions are overrepresented by the skin proteome. The set of predicted elastic domains is further extended by incorporating sequence context embedded in protein language models. Protein domains with distinct thermodynamic recoil mechanisms cluster on the latent space manifold. Some of these domains are anticipated to have roles within molecular machines, expanding the scope of elastic protein function beyond mechanical materials like elastin and resilin.

Contextualized protein-protein interaction networks provide crucial insight into diseases and other biological processes, but for a profound understanding of such processes and their distinct effects on individuals, the protein-protein interactions within individual samples must be investigated. A straightforward approach to estimate the PPI network of a sample is to restrict a general network of known PPIs to the proteins that are found in the sample. Although proteomics methods are becoming more accessible and precise, large-scale and single-cell studies still mainly target characterizing the transcriptomics profile of the samples, which is then often used as an approximation of the protein activities. The correlation of gene expression and protein abundance has been addressed in the past, but information about the deviations of the different omics-based estimates of the PPI networks is still lacking. In this study, we performed a comparative analysis of transcriptomic-based and proteomic-based sample-specific PPI network estimates to fill this gap. We created a framework for a comprehensive and transparent comparison of the two omics levels in two independent datasets, with a special focus on time-related network dynamics. We found that the size-adjusted characteristics of the different omics-based networks are very similar; the overall trend of how they change with time is also often the same, but the rate of the changes typically differs. The characteristics of the nodes present in both types of networks also show high similarity and often different time-related rates of change, but this varies among metrics. These results shed light on the properties of PPI network estimations and advise caution in interpreting them appropriately.

Selecting stochastic gene-expression models from single-cell counts requires accurate parameter inference and efficient model selection. Likelihood methods in count space can be costly when full stationary count distributions are unavailable, whereas approximate methods may lose accuracy. Probability generating functions (PGFs) offer a compact analytical alternative, but existing PGF workflows are generally not likelihood based and therefore rely on computationally intensive cross-validation. We develop a likelihood-based PGF framework for both tasks. Correlated empirical PGF values are used to construct a Gaussian quasi-likelihood for parameter inference and PGF-based Bayesian information criterion (BIC) for model selection. We show that the empirical PGF is exactly unbiased and that the parameter estimator is consistent, converges at the inverse-square-root sample-size rate, and is first-order asymptotically unbiased. For large samples and a uniquely preferred model, PGF-BIC selects the same model as leave-one-out cross-validation in PGF space.

Within the studies of stochastic gene expression, apart from the variability of copy number of gene products, the problems of threshold crossing of those products are biologically important as they often lead to terminal cellular events. Here, we study the threshold crossing problem of the messenger ribonucleic acid (mRNA) and present an exact probability distribution of first passage times in Laplace space. The function furnishes moments of any order and also predicts the characteristic time of the exponential tail of the distribution, which we match against Gillespie simulations. We find that all the measures of relative fluctuations of the threshold crossing times show U-shapes within this simple model of mRNA, as was found earlier in more mathematically involved models of threshold crossing time statistics of proteins. Furthermore, we extend the exact formula to include the phenomenon of DNA duplication and the corresponding doubling of transcription rate. As expected, the distribution varies considerably depending on the onset of the duplication stage within the cell cycle.

Based on our experimental observation of activation state oscillations of different frequencies used to communicate information by the master human stress response regulator protein p38 MAPK 1, we develop a simple graphical approach for understanding and predicting the behavior of complex biological networks acting upon signals carrying information as different frequency waves of chemicals. This approach uses the same techniques used for analyzing and understanding information transmission using waves of electrical currents and fields used in electrical alternating current (AC) circuits. We show how biological components can be organized to behave as standard components found in electronic telecommunications circuits. Finally, we demonstrate how such components can be organized into complex biological signaling cascades whose behavior can be qualitatively and quantitatively understood, and whose output resembles that experimentally observed in p38.

A fundamental biological process of transcription occurs in the cell nucleus, which is a complex medium that also contains multiple heterogeneous structures known as biomolecular condensates. Interestingly, some of these condensates contain mRNA molecules in addition to proteins, suggesting an important cellular role in transcription that is not yet well understood. In this work, we develop a minimal theoretical framework for quantitative investigation of the role of reversible mRNA condensation in transcription. Our discrete-state stochastic approach accounts for the most relevant processes, allowing us to explicitly evaluate the properties of the system and clarify the effects of condensation. Analytical calculations supported by computer simulations suggest that reversible mRNA condensation influences the transcription processes by maintaining a constant level of free mRNA in the nucleoplasm while lowering the degree of stochastic noise and increasing the robustness against external perturbations. Physicochemical arguments are presented to explain these observations. The proposed theoretical framework elucidates important microscopic aspects of transcription, providing a convenient quantitative tool for investigating complex biological phenomena.

Time compaction theory is a general framework explaining how a brain can efficiently deal with dynamic situations occurring in, e.g., sports games. It involves a geometric representation of the time dimension, which enables effective learning and strategic action planning. The theory has recently received experimental support in humans. However, its current computational model has an important limitation: it does not account for deliberate waiting and speed modulation, behaviors ubiquitous in natural environments. This work substantially extends the original model formulation by a dimensional lifting of an n-D workspace into (n + 1)-D mental space, where time remains geometrically embedded. The proposed biologically inspired computational model can generate adaptive behavior across increasingly complex situations, from navigation in everyday social environments to competitive sports. Furthermore, by actively conditioning the expected responses of other agents and stabilizing future predictions, we introduce the concept of uncertainty points in sequences of generalized cognitive maps to support the generation of adaptive strategies in interactive environments, where future prediction has a limited time horizon. Thus, we provide a mechanism for chaining short-term solutions into long-term strategies, which is illustrated by simulating the behavior of a player in a real football game. Author summaryHumans often anticipate future interactions in dynamic environments. Many behaviors, such as avoiding other pedestrians, letting someone pass through a narrow corridor, or reproducing the kind of dribbling maneuvers performed by elite football players, require deciding not only where to move but also when to move. Existing theories suggest that the brain simplifies such situations by representing future interactions as static spatial maps, making them easier to learn and recall. However, current computational models cannot naturally account for common behaviors such as waiting, slowing down, or modulating speed. Here we show that these behaviors readily emerge if the model space is extended by an additional virtual coordinate that encodes accumulated waiting rather than physical time. The proposed model simultaneously admits a wide variety of behaviors, including speed modulation, multigoal decisions, and compound actions, while preserving the principles of time compaction. We illustrate the model in everyday situations and by reproducing two real football plays, comparing the observed behaviors with model simulations. Our results suggest computational principles through which the human brain may efficiently represent, memorize, and exploit dynamic situations.

A novel mathematical framework to define the threshold of action potentials in excitable cells is presented. Unlike previously applied methods that rely on approximations or bifurcations, the approach focuses on the geometry of membrane potential trajectories. The changes in concavity during the upstroke of an action potential can be directly obtained from a time series of voltages. The concavity criterion is then extended to models based on autonomous dynamical systems where the changes in concavity can be obtained analytically from a curve of inflection points in phase space. The inflection point manifold defines a region required for excitability: all the orbits that cross it contain action potentials, and all the trajectories that contain action potentials are in it. This analytical principle can then be used to define excitability in a dynamical system, and also a measure of excitability that enables quantification and comparisons of excitability across dynamical system. The measure provides a way to compare the excitabilities of systems that model neurons with different electrophysiological phenotypes and consider different stimulus conditions. The traditionally vague physiological concept of electrical excitability is transformed into a rigorous analytical description by considering the time-dependent curvature of the membrane potential. The criterion is robust across smooth, single compartment models of electrical excitability and can be can be extended to single compartment models in higher dimensions, and multicompartment models as well.

This paper presents and critically reviews the results of original and some literature based experimental studies conducted by the authors last years on the structural organization of DNA in dormant (starvation stress), anabiotic dormant (4 HR treatment) E. coli cells, as well as the K12 {Delta}dps strain, which lacks the Dps protein (Dps null E. coli). The experimental data includes small-angle synchrotron radiation diffraction (SAXS) and transmission electron microscopy (TEM) data. Synchrotron radiation diffraction experiments on K12{Delta}dps cells allowed us to conclude that peaks at 44.3, 22.1, and 14.8 angstrom resolutions are associated exclusively with ordered DNA organization. Peaks at 44.3, 22.1, and 14.8 angstrom resolutions are also observed for samples of dormant (starvation stress) cells and anabiotically dormant cells. Therefore, this ordered DNA organization also applies to samples of dormant and anabiotically dormant cells. A model is proposed that considers the ordered DNA organization in the cell as a cholesteric liquid crystal. The powder diffraction pattern calculated based on this model is compared with experimental small angle X ray scattering (SAXS) data obtained on Dps-null cell samples. The model completely reproduces the key features of the experimental diffraction pattern from Dps-null cell samples. Accordingly, the cholesteric liquid crystal model corresponds to DNA packaging in dormant and anabiotically dormant cells. Cholesteric liquid crystal ordering should be further considered in all models of cellular DNA packaging. To address the question of which structural organization of DNA predominates in the cell: the cholesteric liquid crystal or nanocrystalline or whether they coexist and fully manifest themselves under different external conditions, it is necessary to utilize the latest methodological advances in structural analysis.

Neutron scattering is a powerful technique to investigate atomic structures and molecular dynamics of proteins at the nano-scale. When it comes to dynamics, incoherent and coherent scattering respectively provide information on the single and collective dynamics of nuclei. In proteins, hydrogen has the highest incoherent cross-section, and it is common practice to overlook the contribution of coherent terms stemming from all nuclei. However, the fast collective dynamics of heavier nuclei could also be studied if coherent scattering and incoherent scattering were experimentally separated. The recent advent of polarized neutron spectroscopy with sufficient flux and energy resolution has made it possible, and opens new perspectives to investigate the relative importance of coherent scattering and the information it provides on biological samples. The present study reports on the use of polarized quasi-elastic neutron scattering (QENS) and the application of a minimalistic model adapted to both individual and collective dynamics. Using a perdeuterated green fluorescent protein as a model globular protein, the study provides an interpretation of the dynamical parameters obtained with QENS, and a comparative study of the Elastic Coherent and Incoherent Scattering Factor. Based on both experiments and calculations, we discuss the relative importance of distinct and self components of coherent scattering, which is often wrongly assumed to be representative of collective dynamics only. The results highlight the current impediments rendering complicated a straightforward analysis of fast collective dynamics in hydrated protein samples.

ShortThe SwissProt database contains a stable 20,418 human protein-coding genes and 42,541 human protein sequences. Ribo-Seq suggests about 7,000 additional, non-canonical Open Reading Frames (ORFs) are present in humans, though only a few of them are confirmed by Mass Spectrometry (MS). Detecting these proteins requires extensive database searches, increasing computational load and inflating False Discovery Rates (FDR). Using the ionbot search engine with the OpenProt database allows for reliable detection of non-canonical proteins while controlling FDR. Ionbot surpasses the Trans-Proteomics Pipeline (TPP) in reproducibility, identifying more peptides and proteins supported by multiple spectra. In addition, open modification searches yield better PSMs compared to closed searches. This work highlights the importance of employing cutting-edge search engines in non-canonical protein research, as well as the value of open modification search in correcting errors in non-canonical protein detection. LongO_ST_ABSBackgroundC_ST_ABSThe SwissProt database reports a quite stable 20,418 human protein-coding genes and 42,541 human protein sequences, figures that have remained stable. New techniques like Ribo-Seq indicate that approximately 7,000 additional, non-canonical Open Reading Frames (ORFs) are translated in humans, few of which have been confirmed by Mass Spectrometry (MS). Detecting these non-canonical proteins requires comprehensive database searches, which increase computational load and False Discovery Rate (FDR). Here, we use the open search engine ionbot in combination with the OpenProt proteogenomics database to reproducibly detect non-canonical proteins while maintaining a well-controlled FDR. ResultsCompared to the current gold standard, the Trans-Proteomics Pipeline (TPP), ionbot shows higher reproducibility, with a higher number of peptides and proteins supported by multiple spectra, and across multiple samples. We observe that PSMs from the open modification search against OpenProt have higher fragment ion intensity correlation compared to PSMs obtained from the closed search, or by only searching canonical proteins. ConclusionsIn this work, we show the potential for open modification searching to correct potential mistakes in non-canonical proteins detection by preventing modified canonical peptides or variants from being incorrectly identified as non-canonical peptides. We also highlight the importance of assessing the FDR of non-canonical identifications separately from canonical ones, as global FDR calculations are biased by the scarcity of non-canonical identifications in each dataset.

Standard methods for ancestral sequence reconstruction (ASR) rely on substitution models for the residues in a biological sequence and assume independent evolution across these sites, ignoring the epistatic interactions that shape molecular evolution. In contrast, deep learning models like variational autoencoders (VAEs) can learn low-dimensional representations ("embeddings") of sequences in a protein family that may implicitly handle these dependencies, raising the possibility of performing more accurate ASR by interpolating between extant sequence embeddings within the VAE's latent space. In this study, we test this hypothesis by developing and evaluating a VAE-based ASR pipeline. Benchmarking this approach against established likelihood-based and parsimony methods using various simulations of protein evolution, including scenarios with and without epistasis, we find that the VAE-based approach is consistently and significantly outperformed by standard methods, even in epistatic regimes where it was hypothesized to have an advantage. We further show that this failure is not due to a lack of phylogenetic structure in the latent space, which does contain evolutionary signal. Rather, the primary limitation is the information loss inherent to the autoencoding process: the VAE's decoder cannot generate sequences with sufficient fidelity for the precise demands of ASR.

Determining the amount of information a neuron receives per unit of time is key to understanding brain connectivity and how neural networks encode and transmit information. In particular, estimating this information flow by distinguishing between excitatory and inhibitory synaptic contributions is critical to understanding neural network function, as maintaining the excitation-inhibition (E/I) balance regulates neuronal excitability and circuit stability, whereas its disruption can lead to a plethora of brain disorders, including neurodegenerative and psychiatric conditions. However, because synaptic conductances cannot be measured directly, inverse methods are required to infer them from the membrane potential --a readily measurable quantity. Although partial solutions have been proposed, accurately estimating these conductances remains a significant challenge due to the complexity and diversity of the inputs. This is particularly true in the spiking regime, where neurons actively fire. In this work, we introduce a novel computational strategy that combines two critical metrics extracted from the time course of the membrane potential recording: the amplitude of the spike and the interspike interval. By using these quantities, the proposed method enables the accurate separation of excitatory and inhibitory contributions, yielding highly favorable results in the spiking regime. Author summaryQuantifying the continuous stream of inputs a neuron receives is key to understanding brain connectivity. Inside the brain, individual cells must maintain a tight balance between excitation and inhibition (E/I) to process information correctly, as any disruption in this equilibrium can impair its functionality. However, directly measuring the underlying excitatory and inhibitory synaptic conductances is technically challenging, and existing mathematical tools often fail when neurons enter their active firing regime. In this work, we introduce a novel computational strategy designed to extract and separate these time-varying conductances directly from the neurons spiking activity. By dynamically tracking just two accessible metrics - the amplitude of the spikes and the time intervals between them - our algorithm estimates both conductance profiles with high precision. Furthermore, we demonstrate that this procedure is highly robust against realistic experimental noise and data variability, providing an accessible framework that does not require complex hardware or an unfeasible number of repetitive experimental trials. By tracking changes in the E/I ratio of the synaptic input, this method provides an efficient approach to detecting pathological imbalances and understanding how local connectivity shapes cellular functionality.

With 302 neurons and a rigorously characterized connectome, the nematode Caenorhabditis elegans represents a powerful model organism to study the fundamental roles of neuronal circuits in behavior. However, despite the breadth of research, many questions remain unanswered regarding how these organisms are able to successfully navigate their environment. Here, we present a biologically grounded dynamical circuit model for the investigation of sensory-guided behavior during C. elegans chemotaxis. Our mathematical model consists of the chemosensory neuron AWA, interneurons RIM and RIA, motor neurons, including SMDs and RMDs, and body wall muscles that provide proprioceptive feedback through stretch receptors. After optimization with an evolutionary algorithm, the model locomotes effectively toward a chemical attractant, realistically capturing nematode chemotactic behavior. Chemotaxis is ensured by sharp turns, which resemble the omega turns of living nematodes, as a key emergent property of the model. The sharp turning behavior is triggered by decreases in the concentration of the attractant. These result in reduced AWA activity, which in turn triggers disinhibition of RIM and subsequent changes in RIA oscillations. The ensuing coordinated changes in downstream motor neurons activity patterns produce sharp turns, which correct the nematodes path, so that the model worm heads toward the attractant, and remains at its proximity, after it reaches the gradient peak. The proposed framework, along with its emergent dynamics, provides new insights into the minimum requirements for C. elegans circuitry to display major features of its chemotactic behavior, including omega turns. In parallel, it generates experimentally testable hypotheses with respect to the participating neuronal elements.

The study of tissue dynamics has been stimulated during the last decades thanks to the use of quantitative descriptions, with the development of several theoretical and computational frameworks, many of them revolving around the notion of reaction-diffusion systems, eventually with additional structure variables beyond time and space. The use of structure variables can accommodate phenotypic traits. In this work, we study a family of competition models, where a given population depends on a resource (e.g. oxygen) and several populations are competing for it. Our quantitative description incorporates phenotypic traits and heterogeneity at the level of cell cycle variations, which influence replication rates via oxygen consumption. This enables us to replicate the fitness of specific subpopulations to environmental conditions (e.g. oxygen shortage or external influences). Using numerical simulations, we show that such models display dynamical pattern formation in the form of coupled travelling wave profiles that expand or retreat at the same wave speed. The full theoretical analysis of such dynamics is quite involved; to circumvent this difficulty, we introduce a quasi-stationary approximation for the resource dynamics. We find that this approximation can reproduce the overall behaviour very accurately, with the additional benefit of allowing theoretical treatment of the reduced model. In this way, we provide estimates on the wave speed which are numerically shown to be robust across a wide range of macroscopic parameters of the full model. The wave speeds are thus found to depend strongly on the proliferation rate of the fittest population, resembling a winner-takes-all dynamics.

Classical neurogeometric models describe the primary visual cortex as a fibered structure in which retinal position and local orientation are coupled through the geometry of the roto-translation group. We extend this approach to the visuomotor cortex by modeling it as an assemblage of visual and motor cortical geometries. The model combines orientation-selective representations, analogous to those of the primary visual cortex, with movement-direction-selective representations, analogous to those of the primary motor cortex, in order to describe the mixed visual and motor selectivity observed in the visuomotor cortex. We introduce a coupled visuomotor structure in which visual orientation and motor direction coexist over a common spatial plane and interact through a relative-orientation constraint. Neural responses are modeled by orientation- and direction-dependent profile functions, and preference maps are obtained from vectorized population responses. Numerical simulations generate visual, motor, and mixed visuomotor response maps. A competition rule between visual and motor responses produces incidence ratios close to experimental observations in macaque visuomotor cortex. This framework provides a first neurogeometric approximation of visuomotor functional architecture and a mathematical setting for studying visually guided action.